微分方程式の解法 | |
時間領域解法 |
s領域解法 |
y"+y'+y=0 , y(0)=-1,
y'(0)=0 特性方程式: u^2+u+1=0, u=(-1±i √3)/2 y(t)=e^(-t/2) { Ca e^(i √3)t/2 +Cb e^(-i √3)t/2 } =e^(-t/2) { A cos((√3)t/2) +B sin(-(√3)t/2) } y(0)=A=-1 y=e^(-t/2) { -cos((√3)t/2) +B sin((√3)t/2) } y'=-(1/2)e^(-t/2) { -cos((√3)t/2) +B sin((√3)t/2) }+ +(√3/2)e^(-t/2) { sin((√3)t/2) +B cos((√3)t/2) } y'(0)=(1/2)+(√3/2)B=0 B=-1/√3 ∴y=-e^(-t/2) { cos((√3)t/2) + (1/√3)sin(-(√3)t/2) } |
一般の微分方程式 |
y=xuと置換して変数分離
|
x*tan(y/x)-y+ x
y'=0の一般解を求めよ。 ただし,x≠0, |y/x|≦π/2とする。 |
[Mapale]による解法 > g:=x*tan(y(x)/x)-y(x)+x*diff(y(x),x); dsolve(g=0); y(x)=arcsin(1/x_C1)*x [解答] y=x*arcsin(C/x) [y=ux置換による解答] dy/dx=u+du/dx を代入 (x^2)du/dx+x*tan(u)=0 x≠0ゆえ xdu/dx+tan(u)=0 du/tan(u)=-(1/x)dx ln|sin(u)|=ln|(C/x)|, Cは積分定数 sin(u)=±C/x ∴sin(y/x)=±C/x または y=±x*arcsin(C/x) |
非線形微分方程式:2次 |
|
y''(x)+Ay'(x)+B=(y'(x))2
|
y(x)=(1/2)(Ax-2ln(cos(c1√(-4B-A2)+(1/2)x√(-4B-A2))))+c2 [解法] y'(x)=v(x)とおくと y''(x)=v'(x) v'(x)+Av(x)+B=(v(x))2 v'(x)=-B+(v(x)) 2 -Av(x) v'(x)/{-B+(v(x)) 2 -Av(x)}=1 ∫ v'(x)/{-B+(v(x)) 2 -Av(x)}dx=∫1dx 2tan-1((2v(x)-A)/√(-4B-A2))/√(-4B-A2)=x+c1 v(x)=(1/2)(A+√(-4B-A2)tan((1/2)√(-4B-A2)(x+c1))) y'(x)=(1/2)(A+√(-4B-A2)tan((1/2)√(-4B-A2)(x+c1))) y(x)=∫(1/2)(A+√(-4B-A2)tan((1/2)√(-4B-A2)(x+c1)))dx =(1/2)(Ax-2ln(cos((1/2)√(-2B-A2) x+(1/2)√(-4B-A2) c1)))+c2 任意定数を簡単化すると y(x)=(1/2)(Ax-2ln(cos((1/2)√(-2B-A2) x+ c1)))+c2 |
初版:2007.7.3
Update:2008.01.02
Update:2012.03.29