谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程
- W5 E7 r  W" p! ]" J0 H5 J
: z: o7 P3 I* ~) w>> clear
>> syms a b c d e;
3 Y; T- y8 Z  h( U" V+ R+ K, T& m>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')

+ K6 e3 q9 `2 [4 F1 Fy =

(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))
7 `1 D1 [, @' K* O
4 w% S' L0 U, S>> a=1;b=1;c=1;d=1;e=1;%若常数已知
% u- e9 Z* ?+ P& i6 \/ C& P. b. L
>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
8 u" y9 m# w0 k( A>> y=eval(y)
% Z8 V# ]- v) K3 D, k/ i
! Q0 X6 {8 U5 {3 Gy =

3 t3 I4 k' S9 F# u; u! V+ K* x   1.1932 + 0.0000i
这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
syms a b c d e;
! T$ ~# k/ b* [! Iy=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
" K. e' L0 y4 v; g: p  K: |x=200:0.1:400;%取值范围，步长
a=1;b=1;c=1;d=1;e=1;%带入参数
y=eval(y);
# o' S* q2 g' e4 k: B5 S1 Nplot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
/ B6 ~& h: m+ t' N/ evb做界面了，matlab做比较轻松愉快~！~