与えられた微分方程式は、以下の通りです。 $y'' + y = t$ 初期条件は、$t=0$のとき、$y(0) = 1$、$y'(0) = -2$です。

応用数学微分方程式ラプラス変換初期条件逆ラプラス変換
2025/5/31
はい、承知いたしました。微分方程式の問題ですね。今回は、問題(1) y+y=ty'' + y = t のラプラス変換を利用した解法を説明します。

1. 問題の内容

与えられた微分方程式は、以下の通りです。
y+y=ty'' + y = t
初期条件は、t=0t=0のとき、y(0)=1y(0) = 1y(0)=2y'(0) = -2です。

2. 解き方の手順

(1) ラプラス変換の適用:
微分方程式の両辺をラプラス変換します。ラプラス変換の線形性より、
L[y+y]=L[t]L[y'' + y] = L[t]
ここで、Y(s)=L[y(t)]Y(s) = L[y(t)]とおくと、ラプラス変換の性質から、
L[y]=s2Y(s)sy(0)y(0)L[y''] = s^2Y(s) - sy(0) - y'(0)
L[y]=Y(s)L[y] = Y(s)
L[t]=1s2L[t] = \frac{1}{s^2}
したがって、
s2Y(s)sy(0)y(0)+Y(s)=1s2s^2Y(s) - sy(0) - y'(0) + Y(s) = \frac{1}{s^2}
(2) 初期条件の代入:
初期条件y(0)=1y(0) = 1y(0)=2y'(0) = -2を代入します。
s2Y(s)s(1)(2)+Y(s)=1s2s^2Y(s) - s(1) - (-2) + Y(s) = \frac{1}{s^2}
s2Y(s)s+2+Y(s)=1s2s^2Y(s) - s + 2 + Y(s) = \frac{1}{s^2}
(3) Y(s)Y(s)について解く:
Y(s)Y(s)について式を整理します。
(s2+1)Y(s)=1s2+s2(s^2 + 1)Y(s) = \frac{1}{s^2} + s - 2
Y(s)=1s2(s2+1)+ss2+12s2+1Y(s) = \frac{1}{s^2(s^2 + 1)} + \frac{s}{s^2 + 1} - \frac{2}{s^2 + 1}
(4) 部分分数分解:
1s2(s2+1)\frac{1}{s^2(s^2 + 1)}を部分分数分解します。
1s2(s2+1)=As+Bs2+Cs+Ds2+1\frac{1}{s^2(s^2 + 1)} = \frac{A}{s} + \frac{B}{s^2} + \frac{Cs + D}{s^2 + 1}
1=As(s2+1)+B(s2+1)+(Cs+D)s21 = As(s^2 + 1) + B(s^2 + 1) + (Cs + D)s^2
1=As3+As+Bs2+B+Cs3+Ds21 = As^3 + As + Bs^2 + B + Cs^3 + Ds^2
係数を比較して、
A+C=0A + C = 0
B+D=0B + D = 0
A=0A = 0
B=1B = 1
したがって、A=0,B=1,C=0,D=1A=0, B=1, C=0, D=-1
1s2(s2+1)=1s21s2+1\frac{1}{s^2(s^2 + 1)} = \frac{1}{s^2} - \frac{1}{s^2 + 1}
(5) 逆ラプラス変換:
Y(s)Y(s)を逆ラプラス変換します。
y(t)=L1[Y(s)]=L1[1s21s2+1+ss2+12s2+1]y(t) = L^{-1}[Y(s)] = L^{-1}[\frac{1}{s^2} - \frac{1}{s^2 + 1} + \frac{s}{s^2 + 1} - \frac{2}{s^2 + 1}]
y(t)=L1[1s2]L1[1s2+1]+L1[ss2+1]2L1[1s2+1]y(t) = L^{-1}[\frac{1}{s^2}] - L^{-1}[\frac{1}{s^2 + 1}] + L^{-1}[\frac{s}{s^2 + 1}] - 2L^{-1}[\frac{1}{s^2 + 1}]
y(t)=tsint+cost2sinty(t) = t - \sin t + \cos t - 2\sin t
y(t)=t3sint+costy(t) = t - 3\sin t + \cos t

3. 最終的な答え

y(t)=t3sint+costy(t) = t - 3\sin t + \cos t

「応用数学」の関連問題

質量 $m$ の質点が、重力と速度の二乗に比例する抵抗力(慣性抵抗)を受けながら落下運動をする。鉛直下向きを正として、質点の位置を $y(t)$、速度を $v(t)$ とする。重力加速度の大きさを $...

微分方程式力学運動終端速度積分
2025/6/2

図に示すような形状の容器に海水を入れたときの、点A, B, C, Dにおける圧力(Pa)をそれぞれ求める問題です。 海水の密度 $\rho = 1025 \text{ kg/m}^3$、重力加速度 $...

圧力流体物理学密度
2025/6/2

質量 $m_e = 9.11 \times 10^{-31} \text{ kg}$、速度 $v = 2.20 \times 10^6 \text{ m/s}$ で飛ぶ電子のド・ブロイ波長を計算します...

物理波長ド・ブロイ波長計算
2025/6/2

幅 $b = 20$ mm を一定とした矩形断面の片持ち梁(長さ $L = 0.6$ m)に先端荷重 $P = 1.0 \times 10^3$ N を作用させる。梁全体の曲げ応力 $\sigma =...

力学構造力学曲げ応力片持ち梁
2025/6/2

問題1:長さ $L = 1.0 \, \text{m}$ の片持ち梁に先端荷重 $P = 1.0 \times 10^3 \, \text{N}$ が作用する。梁材料のヤング率を $E = 200 \...

材料力学たわみヤング率断面2次モーメント
2025/6/2

電子の質量 $m_e = 9.11 \times 10^{-31} \text{ kg}$、速度 $v = 2.20 \times 10^6 \text{ m/s}$ で飛ぶ電子のド・ブロイ波長 $\...

物理学ド・ブロイ波長計算科学
2025/6/2

電子の質量 $m_e = 9.11 \times 10^{-31} \text{ kg}$、速度 $v = 2.20 \times 10^6 \text{ m/s}$ で飛んでいる電子のド・ブロイ波長...

物理量子力学ド・ブロイ波長計算
2025/6/2

質量 $50.0 kg$ の人が速度 $10.0 km/h$ で動いているときのド・ブロイ波長を計算する問題です。 ド・ブロイ波長の公式 $\lambda = h/p$ が与えられています。 ここで、...

物理波長運動量ド・ブロイ波長物理定数
2025/6/2

1993年の半導体素子と集積回路の合計生産額の前年比増加率が5%であるとき、1992年の半導体素子と集積回路の合計生産額を求める。

割合計算経済
2025/6/2

GaAs中の電子の移動度 $\mu$ が $8500 \text{ cm}^2/\text{Vs}$、電子密度 $n$ が $1.0 \times 10^{15} \text{ cm}^{-3}$、電...

物理電子工学移動度抵抗率半導体
2025/6/2