コンテンツにスキップ

定数変化法

必ず標準形 y′′+p(x)y′+q(x)y=r(x)y''+p(x)y'+q(x)y=r(x) へ直し、同次方程式の基本解系を y1,y2y_1,y_2 とします。yp=u1y1+u2y2y_p=u_1y_1+u_2y_2 と置きます。未知関数が2つある自由度を整理するため u1′y1+u2′y2=0u_1'y_1+u_2'y_2=0 を補助条件にします。微分して方程式へ代入すると u1′y1′+u2′y2′=r.u_1'y_1'+u_2'y_2'=r. したがって (y1y2y1′y2′)(u1′u2′)=(0r).\begin{pmatrix}y_1&y_2\\y_1'&y_2'\end{pmatrix}\begin{pmatrix}u_1'\\u_2'\end{pmatrix}=\begin{pmatrix}0\\r\end{pmatrix}.

基本解の順序を交換すると WW の符号も変わり、最終的な特殊解は同値になります。初期値問題では、x0x_0 から xx までの定積分を使うと特殊解とその導関数が x0x_0 で0になり、定数決定が整理しやすくなります。これは後にグリーン関数へつながります。

Q基本:公式の確認★☆☆
y′′+y=r(x)y''+y=r(x)、y1=cos⁡x,y2=sin⁡xy_1=\cos x,y_2=\sin x の u1′,u2′u_1',u_2' を求めてください。
HINTヒントを見る
W=1W=1。
ANS解答を見る
u1′=−r(x)sin⁡x,u2′=r(x)cos⁡xu_1'=-r(x)\sin x, u_2'=r(x)\cos x。
Q標準:指数基本解★★☆
y′′−y=1xy''-y=\frac{1}{x}(x>0x>0)の特殊解を積分表示してください。
HINTヒントを見る
y1=ex,y2=e−x,W=−2y_1=e^x,y_2=e^{-x},W=-2。
ANS解答を見る
u1′=e−x2x,u2′=−ex2xu_1'=\frac{e^{-x}}{2x}, u_2'=-\frac{e^x}{2x}。したがって yp=ex2∫e−xx dx−e−x2∫exx dxy_p=\frac{e^x}{2}\int \frac{e^{-x}}{x}\,dx-\frac{e^{-x}}2\int \frac{e^x}{x}\,dx。初等関数に積分できなくても有効です。
Q発展:完全計算★★★
y′′+y=csc⁡xy''+y=\csc x(0<x<π0<x<\pi)の特殊解を求めてください。
HINTヒントを見る
u1′=−1,u2′=cot⁡xu_1'=-1, u_2'=\cot x。
ANS解答を見る
u1=−x,u2=log⁡(sin⁡x)u_1=-x,u_2=\log(\sin x) と選び、yp=−xcos⁡x+sin⁡xlog⁡(sin⁡x)y_p=-x\cos x+\sin x\log(\sin x)。公式から構成したので代入すると yp′′+yp=csc⁡xy_p''+y_p=\csc x。