谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
5 ]) g3 _/ j4 b0 L$ |
>> clear
  p' \/ r4 ^4 J5 c- b>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
4 |: u. l; h: y1 R- j
( _) x5 k- g2 e; j1 `& P3 c3 vy =

(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;%若常数已知

>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
>> y=eval(y)
9 V* U' b1 A4 y: Z
* l; x- Z2 L/ j! N9 g# ~' t( X2 Z7 sy =

   1.1932 + 0.0000i
' o) u) E5 Y( \这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
syms a b c d e;
% _, K5 v* A% C/ u; `7 }% s% \  u) Ty=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
x=200:0.1:400;%取值范围，步长
' j% t. [) B" @0 ?, r' [$ Za=1;b=1;c=1;d=1;e=1;%带入参数
+ L: U4 W' E0 A1 k1 m6 gy=eval(y);
: Q/ J4 ?5 P# E0 S7 R& ?2 tplot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
" K9 T, d' K1 i5 Z& e9 Svb做界面了，matlab做比较轻松愉快~！~