谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
# D2 L% {0 Q) s* l+ K: R
>> clear
>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
* }% ^" k, |8 U/ P
y =
( N% M" R! r$ K) ?" d! L
(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;%若常数已知

& b* n& Y* @4 g>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
>> y=eval(y)
% @2 H: K$ d- A1 k2 S" e" h1 K
9 X2 o) k/ [- T( n9 e* @- ny =

   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');
' K5 o4 t- X0 [2 Kx=200:0.1:400;%取值范围，步长
, @5 Q: o9 u. ]+ V: Z! ~1 {a=1;b=1;c=1;d=1;e=1;%带入参数
/ F5 J3 q4 B; S5 z/ n* S2 U  N& ty=eval(y);
- _( H0 K) M' X8 \plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
* |1 A6 s8 A) d% R& Hvb做界面了，matlab做比较轻松愉快~！~