谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
  m1 o( u+ C  `
>> clear
>> syms a b c d e;
/ ]* Q9 v/ J2 Y>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
5 W& ~5 M9 b6 T+ m  @  q) O& I: v
! }5 P+ P0 Y& s) x# o% \4 _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))
& A& L7 O" m" j# C0 r
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
5 X/ w) o' n9 {* Z( I0 a( E
>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
8 ^; M( L, F0 P8 u* C# T2 N" ]+ d>> y=eval(y)

6 p- N; h! q* q. D  s. J' z9 qy =

   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');
x=200:0.1:400;%取值范围，步长
a=1;b=1;c=1;d=1;e=1;%带入参数
. Y, s" ?) t* q0 g' d) S) J" N: ~* uy=eval(y);
5 v4 d, r! ~" h  [+ {* ~plot(x,y)

youhuizhong

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