与えられた微分方程式 $y'' + 4y = \sin{t}$ を、初期条件 $y(0) = 0$ および $y'(0) = 0$ の下で、ラプラス変換を用いて解く問題です。

応用数学微分方程式ラプラス変換初期条件逆ラプラス変換
2025/5/31

1. 問題の内容

与えられた微分方程式 y+4y=sinty'' + 4y = \sin{t} を、初期条件 y(0)=0y(0) = 0 および y(0)=0y'(0) = 0 の下で、ラプラス変換を用いて解く問題です。

2. 解き方の手順

まず、与えられた微分方程式の両辺をラプラス変換します。Y(s)Y(s)y(t)y(t) のラプラス変換とすると、
L{y(t)}=s2Y(s)sy(0)y(0)\mathcal{L}\{y''(t)\} = s^2Y(s) - sy(0) - y'(0)
L{y(t)}=Y(s)\mathcal{L}\{y(t)\} = Y(s)
L{sint}=1s2+1\mathcal{L}\{\sin{t}\} = \frac{1}{s^2 + 1}
初期条件 y(0)=0y(0) = 0y(0)=0y'(0) = 0 を代入すると、
L{y(t)}=s2Y(s)\mathcal{L}\{y''(t)\} = s^2Y(s)
したがって、ラプラス変換された微分方程式は、
s2Y(s)+4Y(s)=1s2+1s^2Y(s) + 4Y(s) = \frac{1}{s^2 + 1}
Y(s)Y(s) について解くと、
Y(s)(s2+4)=1s2+1Y(s)(s^2 + 4) = \frac{1}{s^2 + 1}
Y(s)=1(s2+4)(s2+1)Y(s) = \frac{1}{(s^2 + 4)(s^2 + 1)}
次に、Y(s)Y(s) を部分分数分解します。
1(s2+4)(s2+1)=As2+4+Bs2+1\frac{1}{(s^2 + 4)(s^2 + 1)} = \frac{A}{s^2 + 4} + \frac{B}{s^2 + 1}
1=A(s2+1)+B(s2+4)1 = A(s^2 + 1) + B(s^2 + 4)
1=(A+B)s2+(A+4B)1 = (A+B)s^2 + (A+4B)
両辺の係数を比較すると、
A+B=0A+B = 0
A+4B=1A+4B = 1
これらの式を解くと、A=13A = -\frac{1}{3} および B=13B = \frac{1}{3} となります。したがって、
Y(s)=131s2+4+131s2+1Y(s) = -\frac{1}{3} \cdot \frac{1}{s^2 + 4} + \frac{1}{3} \cdot \frac{1}{s^2 + 1}
Y(s)=162s2+22+131s2+12Y(s) = -\frac{1}{6} \cdot \frac{2}{s^2 + 2^2} + \frac{1}{3} \cdot \frac{1}{s^2 + 1^2}
最後に、Y(s)Y(s) を逆ラプラス変換して、y(t)y(t) を求めます。
L1{2s2+22}=sin2t\mathcal{L}^{-1}\left\{\frac{2}{s^2 + 2^2}\right\} = \sin{2t}
L1{1s2+12}=sint\mathcal{L}^{-1}\left\{\frac{1}{s^2 + 1^2}\right\} = \sin{t}
したがって、
y(t)=16sin2t+13sinty(t) = -\frac{1}{6} \sin{2t} + \frac{1}{3} \sin{t}

3. 最終的な答え

y(t)=13sint16sin2ty(t) = \frac{1}{3}\sin{t} - \frac{1}{6}\sin{2t}

「応用数学」の関連問題

長さ $l$ の弦が張力 $S$ で張られている。弦の一端から $x$ の位置に質量 $m$ のおもりをつけ、水平面内で糸に垂直な方向に微小振動させた。このときの周期を $x$ の関数として求める。 ...

力学振動微分方程式物理
2025/6/2

長さ $l$ の弦が張力 $S$ で張られています。弦の一端から $x$ の位置に質量 $m$ のおもりをつけて、水平面内で糸に垂直な方向に微小振動させたときの周期 $T$ を $x$ の関数として求...

力学振動微分方程式単振動物理
2025/6/2

与えられた例を参考に、(1)から(4)の酸の水溶液のpHを求める問題です。pHは水素イオン濃度[H+]を用いて、$pH = -log_{10}[H^+]$と表されます。水素イオン濃度が $1 \tim...

pH対数化学酸・塩基水素イオン濃度電離度
2025/6/2

断面積が変化する配管があり、断面①から水が $u_1 = 7 \text{ m/s}$ の速度で噴出している。配管途中の断面②の半径は $30 \text{ mm}$ で、流速は $u_2 = 4 \...

流体力学連続の式物理計算
2025/6/2

断面が滑らかに拡大している配管があり、流量 $Q = 9.2 \, \text{l/min}$ である。断面1, 2, 3 の直径はそれぞれ $50 \, \text{mm}$, $70 \, \te...

流体力学流量流速断面積円の面積単位変換
2025/6/2

図3に示すマノメータにおいて、タンクAとB内の圧力 $p_A$ と $p_B$ の差 $p_A - p_B$ を、水の密度 $\rho$、水銀の密度 $\rho_g$、重力加速度 $g$、高さ $H$...

流体力学圧力物理マノメータ計算
2025/6/2

図1のようなマノメータにおいて、配管の図に示す位置での圧力 $p$ を求めます。水の密度 $\rho_w = 1000 \, \text{kg/m}^3$、水銀の密度 $\rho_{Hg} = 136...

圧力流体力学マノメータ物理
2025/6/2

図に示すようなマノメータにおいて、配管内の圧力 $P_A$ を求める問題です。与えられている値は、水の密度 $\rho_w = 1000 \, \text{kg/m}^3$、水銀の密度 $\rho_{...

流体圧力マノメータ物理
2025/6/2

1. ベクトル場 $\mathbf{A} = (x^2, 2, z)$ に対して、原点を中心とする半径 $a$ の球面 $S$ 上での面積分 $\oint_S \mathbf{A} \cdot d\...

ベクトル解析面積分ガウスの発散定理球面座標
2025/6/2

問題1は、ベクトル場 $\mathbf{A} = (x^2, 2, z)$ の閉曲面 $S$ 上の面積分 $\oint_S \mathbf{A} \cdot d\mathbf{S}$ を求める問題です...

ベクトル解析面積分ガウスの発散定理球面座標多重積分
2025/6/2