$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) 時刻が経つにつれて、質点の位置はどのようにふるまうかを答える。

応用数学微分方程式運動方程式速度加速度積分減衰振動終端速度
2025/5/23
## 【問1】

1. 問題の内容

xx軸上を運動する質点の時刻 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. 解き方の手順

(i) グラフの概形:
v(t)=et2sin(2t)v(t) = e^{-\frac{t}{2}} \sin(2t) について、0t2π0 \le t \le 2\pi で考える。
* t=0t=0 のとき、v(0)=e0sin(0)=0v(0) = e^0 \sin(0) = 0
* v(t)=0v(t)=0 となる tt の値を求める。et2>0e^{-\frac{t}{2}} > 0 であるから、sin(2t)=0\sin(2t) = 0 となる tt を探す。
2t=0,π,2π,3π,4π2t = 0, \pi, 2\pi, 3\pi, 4\pi
t=0,π2,π,3π2,2πt = 0, \frac{\pi}{2}, \pi, \frac{3\pi}{2}, 2\pi
* 0<t<π20 < t < \frac{\pi}{2} のとき sin(2t)>0\sin(2t) > 0, et2>0e^{-\frac{t}{2}}>0 なので v(t)>0v(t) > 0
* π2<t<π\frac{\pi}{2} < t < \pi のとき sin(2t)<0\sin(2t) < 0, et2>0e^{-\frac{t}{2}}>0 なので v(t)<0v(t) < 0
* π<t<3π2\pi < t < \frac{3\pi}{2} のとき sin(2t)>0\sin(2t) > 0, et2>0e^{-\frac{t}{2}}>0 なので v(t)>0v(t) > 0
* 3π2<t<2π\frac{3\pi}{2} < t < 2\pi のとき sin(2t)<0\sin(2t) < 0, et2>0e^{-\frac{t}{2}}>0 なので v(t)<0v(t) < 0
* et2e^{-\frac{t}{2}}tt が増加するとともに減少する。
以上を考慮すると、tt が大きくなるにつれて振幅が小さくなる振動となる。
(ii) 加速度 a(t)a(t) の計算:
加速度は速度の時間微分で求められる。
a(t)=dv(t)dt=ddt(et2sin(2t))a(t) = \frac{dv(t)}{dt} = \frac{d}{dt} \left( e^{-\frac{t}{2}} \sin(2t) \right)
積の微分公式より、
a(t)=12et2sin(2t)+et2(2cos(2t))a(t) = -\frac{1}{2}e^{-\frac{t}{2}} \sin(2t) + e^{-\frac{t}{2}} (2 \cos(2t))
a(t)=et2(2cos(2t)12sin(2t))a(t) = e^{-\frac{t}{2}} \left( 2\cos(2t) - \frac{1}{2}\sin(2t) \right)
(iii) 位置 x(t)x(t) の計算:
位置は速度の積分で求められる。x(0)=0x(0) = 0 であることに注意する。
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=(2et2)sin(2t)dt=2et2sin(2t)2et22cos(2t)dtI = \int (-2e^{-\frac{t}{2}})' \sin(2t) dt = -2e^{-\frac{t}{2}} \sin(2t) - \int -2e^{-\frac{t}{2}} 2\cos(2t) dt
I=2et2sin(2t)+4et2cos(2t)dtI = -2e^{-\frac{t}{2}} \sin(2t) + 4 \int e^{-\frac{t}{2}} \cos(2t) dt
I=2et2sin(2t)+4(2et2)cos(2t)dt=2et2sin(2t)+4(2et2cos(2t)2et2(sin(2t))2dt)I = -2e^{-\frac{t}{2}} \sin(2t) + 4 \int (-2e^{-\frac{t}{2}})' \cos(2t) dt = -2e^{-\frac{t}{2}} \sin(2t) + 4(-2e^{-\frac{t}{2}} \cos(2t) - \int -2e^{-\frac{t}{2}} (-\sin(2t))2dt)
I=2et2sin(2t)8et2cos(2t)16et2sin(2t)dtI = -2e^{-\frac{t}{2}} \sin(2t) - 8 e^{-\frac{t}{2}} \cos(2t) - 16\int e^{-\frac{t}{2}} \sin(2t) dt
I=2et2sin(2t)8et2cos(2t)16II = -2e^{-\frac{t}{2}} \sin(2t) - 8 e^{-\frac{t}{2}} \cos(2t) - 16I
17I=2et2sin(2t)8et2cos(2t)17I = -2e^{-\frac{t}{2}} \sin(2t) - 8 e^{-\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
x(0)=0x(0) = 0 より、
0=217e0sin(0)817e0cos(0)+C0 = -\frac{2}{17}e^{0} \sin(0) - \frac{8}{17} e^{0} \cos(0) + C
0=0817+C0 = 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 のときの質点の位置:
tt \to \infty のとき、et20e^{-\frac{t}{2}} \to 0 なので、
limtx(t)=limt(217et2sin(2t)817et2cos(2t)+817)=817\lim_{t \to \infty} x(t) = \lim_{t \to \infty} \left( -\frac{2}{17}e^{-\frac{t}{2}} \sin(2t) - \frac{8}{17} e^{-\frac{t}{2}} \cos(2t) + \frac{8}{17} \right) = \frac{8}{17}
質点の位置は 817\frac{8}{17} に近づく。

3. 最終的な答え

(i) グラフの概形:振幅が減衰する振動
(ii) 加速度:a(t)=et2(2cos(2t)12sin(2t))a(t) = e^{-\frac{t}{2}} \left( 2\cos(2t) - \frac{1}{2}\sin(2t) \right)
(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) 質点の位置:817\frac{8}{17} に近づく
## 【問2】

1. 問題の内容

質量 mm の物体が、重力と空気抵抗を受けて落下する運動について、以下の問いに答える。重力加速度の大きさを gg とし、鉛直上向きを yy 軸の正の向きとする。空気抵抗は速度 vv に比例する粘性抵抗 (大きさ 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. 解き方の手順

(i) 運動方程式:
運動方程式は、物体の質量 mm と加速度 aa の積が、物体に作用する力の和に等しいというものである。
この場合、物体に作用する力は、鉛直下向きの重力 mgmg と、鉛直上向きの空気抵抗 bvbv である。
鉛直上向きを正としているので、重力は負の方向に働く。
よって、運動方程式は
ma=mg+bvma = -mg + bv
mdvdt=mg+bvm \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}
dvdt=bm(v0mgb)ebmt\frac{dv}{dt} = -\frac{b}{m} (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t}
運動方程式に代入する。
mdvdt=mg+bvm \frac{dv}{dt} = -mg + bv
m(bm(v0mgb)ebmt)=mg+b((v0mgb)ebmt+mgb)m \left( -\frac{b}{m} (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t} \right) = -mg + b \left( (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t} + \frac{mg}{b} \right)
b(v0mgb)ebmt=mg+b(v0mgb)ebmt+mg-b (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t} = -mg + b (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t} + mg
b(v0mgb)ebmt=b(v0mgb)ebmt-b (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t} = b (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t}
t=0t=0の時v=v0v=v_0を確認する。
v(0)=(v0mgb)ebm0+mgb=(v0mgb)+mgb=v0v(0) = (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}0} + \frac{mg}{b} = (v_0 - \frac{mg}{b}) + \frac{mg}{b} = v_0
(iii) 終端速度:
十分時間が経過したとき、tt \to \infty とすると、ebmt0e^{-\frac{b}{m}t} \to 0 となる。
したがって、
limtv(t)=limt((v0mgb)ebmt+mgb)=mgb\lim_{t \to \infty} v(t) = \lim_{t \to \infty} \left( (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t} + \frac{mg}{b} \right) = \frac{mg}{b}
よって、速度は mgb\frac{mg}{b} に漸近する。これが終端速度である。

3. 最終的な答え

(i) 運動方程式:mdvdt=mg+bvm \frac{dv}{dt} = -mg + bv
(ii) 解の確認:省略(上記参照)
(iii) 終端速度:mgb\frac{mg}{b}

「応用数学」の関連問題

図1のように、初速度 $v_0$ で角度 $\theta$ で原点から発射された物体が、座標 $(X, Y)$ に到達するような投射角 $\tan \theta$ を求める問題です。

力学物理運動方程式放物運動二次方程式微分三角関数
2025/6/18

初速度 $v_0$ で角度 $\theta$ で投射された物体が、座標 $(X, Y)$ に命中するような $\tan \theta$ の値を求める問題です。重力加速度を $g$ とします。

力学放物運動軌道微分積分最適化
2025/6/18

(1) 水平方向に速さ40 m/sで飛んできた質量0.15 kgのボールをバットで打ち返したところ、打ち返したボールが水平逆向きに同じ速さで飛んでいった。このとき、バットとの衝突の間にボールが受ける力...

力学運動量力積物理
2025/6/18

球座標系 $(r, \theta, \phi)$ におけるベクトル場 $\vec{A}$ の発散を、微小体積における流出量から求める。

ベクトル解析発散球座標系微分
2025/6/18

この問題は、与えられたマクロ経済モデルにおいて、IS-LM分析を用いて均衡実質GDP($Y^*$)と均衡実質利子率($r^*$)を求める問題です。物価水準は$P=3$で一定であり、期待物価上昇率は$0...

マクロ経済学IS-LM分析連立方程式経済モデル
2025/6/18

長さ $3l$ [m] の糸に質量 $m = 5$ [kg] の重りを取り付け、糸が鉛直と $60^\circ$ の角度をなして同一水平面を等速円運動している。円周率は $\pi$ 、重力加速度の大き...

力学円運動物理ベクトル運動方程式
2025/6/18

ヘリウム(He)をファンデルワールスの実在気体として扱ったとき、ファンデルワールス定数 $b$ が教科書表8-1より $0.0238 \text{ L/mol}$ である。Heを球と仮定したときの半径...

物理化学ファンデルワールス力気体体積半径アボガドロ定数計算
2025/6/18

問題1:周期 $T = 0.2$ [s] の単振動を行っている質量 $m = 0.5$ [kg] の物体の最大の速さが $4\pi$ [m/s] であるとき、(1)この単振動の振幅と(2)物体が受ける...

単振動物理角振動数振幅力の大きさ
2025/6/18

ある分子の基底状態と第一励起状態の縮退度がそれぞれ1と3である。これらのエネルギー準位間の差は $365 \text{ cm}^{-1}$ に相当する。$25^\circ\text{C}$ では、第一...

ボルツマン分布熱力学指数関数物理化学
2025/6/18

ヘリウム(He)をファンデルワールスの実在気体として扱ったとき、ファンデルワールス定数 $b$ は教科書の表8-1より $0.0238 \ L \ mol^{-1}$ である。この$b$の値から、ヘリ...

物理化学ファンデルワールス力原子半径体積
2025/6/18