谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程

>> clear
>> syms a b c d e;
, `" W6 x! J3 `% r>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
3 N& X  K- T' y2 P
y =
. m: x  y. u# P( T6 x* @
7 O% |( D, k- ?& v(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))
6 K; e) |/ h: o# y# C# n# s
% B0 }5 B. c+ Z/ O; v>> a=1;b=1;c=1;d=1;e=1;%若常数已知
! v; r1 _" H- q* V- M6 f# |; |
>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
>> y=eval(y)
: n5 D# H& t  s' `* H
y =
$ t6 L: d  w1 z( k6 T1 Y- \9 O5 L
   1.1932 + 0.0000i
5 W' G( K4 c9 u' \& j这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
" R* \- q  f- B2 }2 g4 j9 _syms a b c d e;
7 Z7 E, U  x7 ^1 g/ T# Z9 P) ey=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
+ j, N, F  a" X& @3 zx=200:0.1:400;%取值范围，步长
a=1;b=1;c=1;d=1;e=1;%带入参数
8 ^, w3 c" [# }6 n3 k$ By=eval(y);
plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
+ k9 @2 |* f- K* w; r3 Y0 q! nvb做界面了，matlab做比较轻松愉快~！~