|
发表于 2014-7-14 00:28:29
|
显示全部楼层
matlab 二阶微分方程6 t) e. }# O5 Z8 z& h/ x" I6 V
. O. \$ ]" R3 T( F>> clear8 H1 |+ P; `9 A
>> syms a b c d e;
/ ?# `" F. W5 R- L: H>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')( o( g+ L: O5 k4 J1 s# Q) V
: c6 x( U- f& w2 h+ ]9 v- h/ }y =
/ |, d4 i! r! {1 R
2 [- l% `7 o# S* v% c0 {, B1 ]6 @: o(2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b + (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2))0 w& b" k7 Q) \0 J
. G: G6 n d8 q {>> a=1;b=1;c=1;d=1;e=1;%若常数已知, X2 w2 @) a0 F
. E: z: ~4 J# N' [* y7 M0 d
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值; d2 a- d b" C- C" n
>> y=eval(y) 7 G. Q5 Q$ O/ m! O
1 @* V1 d4 t1 Ry =
: n7 D7 i7 q9 J5 e' Y$ { g
5 s1 p3 Q) U2 }: ?1 A1 R9 ]$ M 1.1932 + 0.0000i
6 a6 J4 X+ E( ^. S0 T1 x这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
4 a: h K0 V! }! gsyms a b c d e;
! ?; j5 ^7 K1 ?, }y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
; p! c- Q# | qx=200:0.1:400;%取值范围,步长' y' u; x. `- h2 q' `
a=1;b=1;c=1;d=1;e=1;%带入参数
: K' M4 @9 V5 c5 L' uy=eval(y);
3 n" Z9 f6 K# c8 {plot(x,y) |
|