谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程

>> clear
>> syms a b c d e;
5 V. `& X- p% l' ~9 ~/ d/ [>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')

y =

1 t$ j5 Y; ^: O: p(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))

>> a=1;b=1;c=1;d=1;e=1;%若常数已知
$ a" D1 N; K, _: ?5 t. E: M
3 s4 C, G, \; Z/ D>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
0 P4 j6 g5 n6 a8 M3 g$ i>> y=eval(y)

y =

! e- s: n# U/ s. Q' ?   1.1932 + 0.0000i
这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
& Y+ A. y. _& v/ M+ M; Isyms a b c d e;
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
5 M" a  T, J3 W, Nx=200:0.1:400;%取值范围，步长
8 A9 o$ f* i; O8 {1 {a=1;b=1;c=1;d=1;e=1;%带入参数
: r* }& t; @2 @y=eval(y);
plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
7 x( J- ^9 d3 p! ]6 ~8 L3 zvb做界面了，matlab做比较轻松愉快~！~