谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
4 k# ^+ c: p, H7 J. `
>> clear
- |3 X) c9 ^8 X) C% t* Z>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')

y =

. B, L1 m; G/ g2 Y) 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))
& J: k0 S' |% ?  X( I6 b/ L& H
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
- Y: v3 \8 R7 z( \+ O8 F
>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
% z* J3 X4 I" N% U>> y=eval(y)

: E4 p( H- \7 v) G( B4 }y =
% S( G( B: d9 \9 p
   1.1932 + 0.0000i
这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
syms a b c d e;
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
$ @5 L) }. n& X, ~( Q3 s: m+ }x=200:0.1:400;%取值范围，步长
a=1;b=1;c=1;d=1;e=1;%带入参数
y=eval(y);
plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
5 M9 P. V$ s$ Q8 @4 v: v- tvb做界面了，matlab做比较轻松愉快~！~