谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
7 j9 ^9 \8 h; Q2 K
>> clear
>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
( _: m7 `; p0 e/ ?" ]
y =
. W9 `+ y) p8 ^
(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))

3 i7 h+ V/ r) {3 E5 Z  b>> a=1;b=1;c=1;d=1;e=1;%若常数已知
( G7 v$ T2 P; L$ j% y4 n# b) f6 _
>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
3 m' e6 i) M8 v0 \' y* Y>> y=eval(y)
% c5 j" @3 }" N
y =
( q1 W! v" g, g7 C  @# V* n+ U* I
   1.1932 + 0.0000i
这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
syms a b c d e;
, K+ d. z" t- Y# ], |y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
x=200:0.1:400;%取值范围，步长
/ b3 S9 v% M# D; {1 J/ p; @a=1;b=1;c=1;d=1;e=1;%带入参数
y=eval(y);
7 B9 I$ }9 T# h  g8 wplot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
" j6 C( |5 u& ^; rvb做界面了，matlab做比较轻松愉快~！~