谁会编程解二阶微分方程

andren

vb能做吗

freddy

matlab做

sunhaiyan

matlab 二阶微分方程

>> clear
>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
! U; @  A  R3 }2 Y1 X  z* b
y =

9 _; @/ D1 S: }  y; T(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))
/ Q# U9 u3 {0 V' p" q  M8 _# m: |2 n
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
  V9 w7 K) |2 q' q! b
% k8 y) s, m4 V+ p/ B>> t=1;%desolve中没指定自变量x，这里默认为t，带入任意一个x值
- J9 `) b' _7 x>> y=eval(y)
# b  I( f# c! ^
y =
  `7 U; v- i7 B  t/ O' x
  v6 u5 Z$ o! N0 Y, K8 b' W   1.1932 + 0.0000i
% a5 a7 V1 l" {1 ~% f9 j8 G4 m" W这是解方程的方法，画图也差不多，我在命令窗写的改麻烦，下面是我在editor写的。
. E2 N, \9 m. H# [6 ~( t5 osyms a b c d e;
; i0 T. Q& b9 g0 }, ky=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
! H6 P, H+ F, j0 f4 s+ I$ `4 nx=200:0.1:400;%取值范围，步长
5 p  r9 e/ M% x2 V: l$ ^a=1;b=1;c=1;d=1;e=1;%带入参数
y=eval(y);
7 r: S* n  `. A$ p% ?plot(x,y)

youhuizhong

大神跑这来问问题。。厉害~！~
$ T1 U& H2 V2 ^& ]2 Hvb做界面了，matlab做比较轻松愉快~！~