谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
1 }8 r& t0 X& F/ J
2 V! r/ B5 o. X# f>> clear
>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
. Y$ }$ ^, q3 \5 C
/ B$ g! Y+ k6 J( @y =

(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))
1 b5 O, o3 }( h2 {2 {' g1 u/ L
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
( B/ N# d9 |( @5 e/ O
9 j# W0 H8 k  b. Q* B>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
>> y=eval(y)
' g6 M" P' |! [8 y
y =
: ?4 ~$ q7 G# e4 x. T
   1.1932 + 0.0000i
这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
: t1 g8 X8 ~* m( F: \5 `6 dsyms a b c d e;
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
x=200:0.1:400;%取值范围，步长
a=1;b=1;c=1;d=1;e=1;%带入参数
3 b. P" @9 o3 K! }: Hy=eval(y);
plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
vb做界面了，matlab做比较轻松愉快~！~