谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
+ E' y* \- X9 g  S- n
>> clear
' ~, K  {* R. }% x8 Z: K, L( y% B>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')

1 H0 y( K$ {9 R3 p: Y8 d" f1 Vy =
7 a* J& F% F# k! ~5 J: p. G
(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 ?( H  m) b! D' ^9 V( q>> a=1;b=1;c=1;d=1;e=1;%若常数已知
4 ~5 O$ h7 z( t% k/ J# V( j
' J0 D! ~2 R- V+ k% [>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
: i0 b* I, `" r) A0 i+ C; Y>> y=eval(y)
3 z- m1 G6 y8 x7 l- h! u
" U' a2 r- Z6 D# Fy =

# \$ V9 i5 {% C/ D   1.1932 + 0.0000i
这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
- \+ U4 X+ T5 P. Q+ L2 `2 usyms 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=eval(y);
plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
. E5 r3 [- g1 C4 J+ Ivb做界面了，matlab做比较轻松愉快~！~