谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
, g# x6 x  c+ {( s' N% Q. s( t  Q( }
>> clear
7 g1 S+ m9 |5 p) c0 Z>> syms a b c d e;
+ U5 f: c2 [3 P+ @( h5 M>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')

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))

7 l9 @5 o  y7 j( t0 d: J>> a=1;b=1;c=1;d=1;e=1;%若常数已知

>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
>> y=eval(y)

y =

   1.1932 + 0.0000i
4 \- g4 {/ k( s3 Y5 N0 E0 Q这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
3 \  t9 |6 {' Y. gsyms a b c d e;
4 h7 t$ {8 E4 r3 t2 v- i. c6 l( _- iy=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
x=200:0.1:400;%取值范围，步长
1 n2 s9 _* M5 o* y. ba=1;b=1;c=1;d=1;e=1;%带入参数
; [4 o" U* Y& O! l# \y=eval(y);
plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
0 {, X' i7 T) t, u! V( L- @% m. Rvb做界面了，matlab做比较轻松愉快~！~