問1: x軸上を運動する質点の時刻 $t$ における速度が $v(t) = e^{-\frac{t}{2}} \sin(2t)$ で与えられている。 (i) $0 \le t \le 2\pi$ の範囲で $v(t)$ の概形を描け。どのように考えて描いたかも説明する。 (ii) 時刻 $t$ における加速度 $a(t)$ を求めよ。 (iii) 時刻 $t$ における位置 $x(t)$ を求めよ。ただし、$t=0$ のとき $x=0$ である。 (iv) 時刻が経つにつれて、質点の位置はどのように振る舞うか答えよ。 問2: 地表面付近で空気抵抗を受ける質量 $m$ の物体の落下運動を考える。重力加速度の大きさを $g$ とし、鉛直上向きを $y$ 軸とする。物体に作用する力は、鉛直下向きの重力 (大きさ $mg$) と、速度に比例する粘性抵抗 (大きさ $bv$, $b>0$) である。 (i) 物体が満たす運動方程式を立てよ。 (ii) $t=0$ で $v=v_0$ を満たす運動方程式の解が $v(t) = (v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b}$ となることを確かめよ。 (iii) 十分時間が経過したとき、速度が一定の速度に漸近することを示し、終端速度を求めよ。

解析学微分積分運動微分方程式減衰振動
2025/5/24

1. 問題の内容

問1:
x軸上を運動する質点の時刻 tt における速度が v(t)=et2sin(2t)v(t) = e^{-\frac{t}{2}} \sin(2t) で与えられている。
(i) 0t2π0 \le t \le 2\pi の範囲で v(t)v(t) の概形を描け。どのように考えて描いたかも説明する。
(ii) 時刻 tt における加速度 a(t)a(t) を求めよ。
(iii) 時刻 tt における位置 x(t)x(t) を求めよ。ただし、t=0t=0 のとき x=0x=0 である。
(iv) 時刻が経つにつれて、質点の位置はどのように振る舞うか答えよ。
問2:
地表面付近で空気抵抗を受ける質量 mm の物体の落下運動を考える。重力加速度の大きさを gg とし、鉛直上向きを yy 軸とする。物体に作用する力は、鉛直下向きの重力 (大きさ mgmg) と、速度に比例する粘性抵抗 (大きさ bvbv, b>0b>0) である。
(i) 物体が満たす運動方程式を立てよ。
(ii) t=0t=0v=v0v=v_0 を満たす運動方程式の解が v(t)=(v0mgb)ebmt+mgbv(t) = (v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b} となることを確かめよ。
(iii) 十分時間が経過したとき、速度が一定の速度に漸近することを示し、終端速度を求めよ。

2. 解き方の手順

問1:
(i) 速度の概形:
まず、v(t)=et2sin(2t)v(t) = e^{-\frac{t}{2}} \sin(2t) のグラフを考える。
- et2e^{-\frac{t}{2}} は減衰関数であり、tt が大きくなるにつれて0に近づく。
- sin(2t)\sin(2t) は周期 π\pi のサイン関数。
- v(t)v(t) は、サイン関数が減衰関数で減衰していくようなグラフになる。
- t=0t=0 のとき、v(0)=0v(0) = 0
- v(t)=0v(t) = 0 となるのは、sin(2t)=0\sin(2t) = 0 のとき、すなわち 2t=nπ2t = n\pi (nn は整数) のとき。t=nπ2t = \frac{n\pi}{2}
- 0t2π0 \le t \le 2\pi の範囲では、t=0,π2,π,3π2,2πt = 0, \frac{\pi}{2}, \pi, \frac{3\pi}{2}, 2\piv(t)=0v(t) = 0
- グラフは、t=π2t = \frac{\pi}{2} ごとにx軸と交わり、振幅は減衰していく。
(ii) 加速度:
加速度 a(t)a(t) は速度 v(t)v(t) の時間微分である。
a(t)=dv(t)dt=ddt(et2sin(2t))a(t) = \frac{dv(t)}{dt} = \frac{d}{dt} (e^{-\frac{t}{2}} \sin(2t))
積の微分公式を用いて、
a(t)=12et2sin(2t)+et2(2cos(2t))=et2(2cos(2t)12sin(2t))a(t) = -\frac{1}{2} e^{-\frac{t}{2}} \sin(2t) + e^{-\frac{t}{2}} (2\cos(2t)) = e^{-\frac{t}{2}}(2\cos(2t) - \frac{1}{2}\sin(2t))
(iii) 位置:
位置 x(t)x(t) は速度 v(t)v(t) の時間積分である。
x(t)=v(t)dt=et2sin(2t)dtx(t) = \int v(t) dt = \int e^{-\frac{t}{2}} \sin(2t) dt
この積分は部分積分を2回行うことで計算できる。
まず、I=et2sin(2t)dtI = \int e^{-\frac{t}{2}} \sin(2t) dt とおく。
I=et2sin(2t)dt=2et2sin(2t)+2et2(2cos(2t))dt=2et2sin(2t)+4et2cos(2t)dtI = \int e^{-\frac{t}{2}} \sin(2t) dt = -2e^{-\frac{t}{2}} \sin(2t) + \int 2e^{-\frac{t}{2}} (2\cos(2t)) dt = -2e^{-\frac{t}{2}} \sin(2t) + 4 \int e^{-\frac{t}{2}} \cos(2t) dt
さらに部分積分を行う。
et2cos(2t)dt=2et2cos(2t)2et2(2sin(2t))dt=2et2cos(2t)4et2sin(2t)dt=2et2cos(2t)4I\int e^{-\frac{t}{2}} \cos(2t) dt = -2e^{-\frac{t}{2}} \cos(2t) - \int -2e^{-\frac{t}{2}} (-2\sin(2t)) dt = -2e^{-\frac{t}{2}} \cos(2t) - 4 \int e^{-\frac{t}{2}} \sin(2t) dt = -2e^{-\frac{t}{2}} \cos(2t) - 4I
よって、I=2et2sin(2t)+4(2et2cos(2t)4I)=2et2sin(2t)8et2cos(2t)16II = -2e^{-\frac{t}{2}} \sin(2t) + 4(-2e^{-\frac{t}{2}} \cos(2t) - 4I) = -2e^{-\frac{t}{2}} \sin(2t) - 8e^{-\frac{t}{2}} \cos(2t) - 16I
17I=2et2sin(2t)8et2cos(2t)17I = -2e^{-\frac{t}{2}} \sin(2t) - 8e^{-\frac{t}{2}} \cos(2t)
I=217et2sin(2t)817et2cos(2t)+CI = -\frac{2}{17} e^{-\frac{t}{2}} \sin(2t) - \frac{8}{17} e^{-\frac{t}{2}} \cos(2t) + C
x(t)=217et2sin(2t)817et2cos(2t)+Cx(t) = -\frac{2}{17} e^{-\frac{t}{2}} \sin(2t) - \frac{8}{17} e^{-\frac{t}{2}} \cos(2t) + C
t=0t=0 のとき x(0)=0x(0) = 0 より、
0=217e0sin(0)817e0cos(0)+C=0817+C0 = -\frac{2}{17} e^{0} \sin(0) - \frac{8}{17} e^{0} \cos(0) + C = 0 - \frac{8}{17} + C
C=817C = \frac{8}{17}
したがって、x(t)=217et2sin(2t)817et2cos(2t)+817x(t) = -\frac{2}{17} e^{-\frac{t}{2}} \sin(2t) - \frac{8}{17} e^{-\frac{t}{2}} \cos(2t) + \frac{8}{17}
(iv) 時刻が経つにつれて:
tt \to \infty のとき、et20e^{-\frac{t}{2}} \to 0 なので、
limtx(t)=817\lim_{t \to \infty} x(t) = \frac{8}{17}
質点の位置は 817\frac{8}{17} に近づく。
問2:
(i) 運動方程式:
yy 軸上向きを正とする。
質量 mm の物体に働く力は、重力 mg-mg と粘性抵抗 bv-bv である。ニュートンの運動方程式は、
mdvdt=mgbvm\frac{dv}{dt} = -mg - bv
(ii) 解の確認:
v(t)=(v0mgb)ebmt+mgbv(t) = (v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b} を運動方程式に代入する。
dvdt=(v0mgb)(bm)ebmt\frac{dv}{dt} = (v_0 - \frac{mg}{b}) (-\frac{b}{m}) e^{-\frac{b}{m}t}
mdvdt=m(v0mgb)(bm)ebmt=b(v0mgb)ebmtm \frac{dv}{dt} = m(v_0 - \frac{mg}{b}) (-\frac{b}{m}) e^{-\frac{b}{m}t} = -b(v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t}
mgbv=mgb((v0mgb)ebmt+mgb)=mgb(v0mgb)ebmtmg=2mgb(v0mgb)ebmt-mg - bv = -mg - b((v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b}) = -mg - b(v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} - mg = -2mg - b(v_0 - \frac{mg}{b})e^{-\frac{b}{m}t}
しかし、このままでは等しくならない。
正しくは、mdvdt=mgbvm\frac{dv}{dt} = -mg - bv なので、mdvdt+bv=mgm\frac{dv}{dt} + bv = -mg である。
mdvdt+bv=b(v0mgb)ebmt+b((v0mgb)ebmt+mgb)=bmgb=mgm\frac{dv}{dt} + bv = -b(v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t} + b((v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b}) = b\frac{mg}{b} = mg
(iii) 終端速度:
tt \to \infty のとき、ebmt0e^{-\frac{b}{m}t} \to 0 なので、
limtv(t)=mgb\lim_{t \to \infty} v(t) = \frac{mg}{b}
終端速度は mgb\frac{mg}{b} である。

3. 最終的な答え

問1:
(i) グラフの説明は上記の通り。グラフは省略。
(ii) a(t)=et2(2cos(2t)12sin(2t))a(t) = e^{-\frac{t}{2}}(2\cos(2t) - \frac{1}{2}\sin(2t))
(iii) x(t)=217et2sin(2t)817et2cos(2t)+817x(t) = -\frac{2}{17} e^{-\frac{t}{2}} \sin(2t) - \frac{8}{17} e^{-\frac{t}{2}} \cos(2t) + \frac{8}{17}
(iv) limtx(t)=817\lim_{t \to \infty} x(t) = \frac{8}{17}
問2:
(i) mdvdt=mgbvm\frac{dv}{dt} = -mg - bv
(ii) 問題文に誤りがあり、mdvdt+bv=mgm\frac{dv}{dt} + bv = -mg と変形すれば解となる。
(iii) limtv(t)=mgb\lim_{t \to \infty} v(t) = \frac{mg}{b}

「解析学」の関連問題

数列 $\{a_n\}$ が与えられたとき、$n \ge N$ ならば $|a_n - \alpha| < 10^{-4}$ が成り立つような最小の自然数 $N$ を求める問題です。ここで、$\alp...

数列極限収束不等式
2025/5/24

与えられた三角関数 $y = 2\sin(\frac{\theta}{2} - \frac{\pi}{3}) + 1$ の周期を求め、さらに、関数 $y = 2\sin\frac{\theta}{2}...

三角関数周期グラフの平行移動振幅
2025/5/24

領域$D$上で、関数$e^{x+y}$の重積分を計算します。ここで領域$D$は、$y \ge 0$, $y \le x$, $x+y \le 2$によって定義されます。

重積分多重積分積分領域
2025/5/24

領域 $D$ が $y \ge 1-x$, $x \le 1$, $y \le 1$ で定義されるとき、二重積分 $\iint_D x^2 y \, dx dy$ の値を計算します。

二重積分領域積分計算
2025/5/24

領域 $D$ 上の重積分 $\iint_D x^2 y \,dx\,dy$ を計算します。領域 $D$ は $y \ge 1-x$, $x \le 1$, $y \le 1$ で定義されています。

重積分積分領域二重積分
2025/5/24

半径2の円周上を運動する質点AとBについて、それぞれの時刻$t$における位置が与えられています。 $r^A(t) = 2(\cos(\frac{\pi t}{3} - \frac{\pi}{6})i ...

円運動軌跡角速度加速度速度微分
2025/5/24

半径2の円周上を運動する2つの質点AとBについて、与えられた位置ベクトル $\vec{r}^A(t)$ と $\vec{r}^B(t)$ をもとに、以下の問いに答える。 * (i) $0 \le ...

ベクトル円運動軌跡角速度加速度微分
2025/5/24

半径2の円周上を運動する質点A, Bの位置がそれぞれ $r^A(t) = 2 (\cos(\frac{\pi t}{3} - \frac{\pi}{6})i + \sin(\frac{\pi t}{3...

ベクトル解析円運動角速度加速度速度
2025/5/24

2変数関数 $z = f(x, y) = x^2 - 6xy + 2y^3$ について、極値があればその値と極大値または極小値を求め、極値がなければ「なし」と答える。

多変数関数偏微分極値ヘッセ行列
2025/5/24

半径2の円周上を運動する質点A, Bについて、それぞれの時刻 $t$ における位置が与えられています。 - 質点Aの位置: $\vec{r}^A(t) = 2(\cos(\frac{\pi t}{3}...

円運動ベクトル速度加速度微分軌跡
2025/5/24