找回密码
 加入华同
搜索
黄金广告位联系EMAIL:[email protected] 黄金广告[email protected]
查看: 1814|回复: 3

谁会编程解二阶微分方程

[复制链接]
发表于 2013-10-11 23:06:08 | 显示全部楼层 |阅读模式
vb能做吗
发表于 2013-10-12 06:56:57 | 显示全部楼层
matlab做
回复

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程
2 ~  f8 s  \. Y( ?. G& V* \9 l# b" l1 p" [3 ], n
>> clear
) _. g# F. H" P% r9 X- B>> syms a b c d e;2 S- G  T- d2 J7 i; h2 X
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')1 `! K" y" R% {7 I1 E

# J! f! w0 c% K. D5 Ty =
  V& T! n9 q/ Y  \
5 l9 _; B  _3 a. X5 c) 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))
5 G! ?6 N/ C6 U7 y7 b # a- w' w7 W( u; q6 e: R
>> a=1;b=1;c=1;d=1;e=1;%若常数已知6 Z) Q: T7 ~+ O0 e" K

) ?5 D9 d# h9 A* A' N% b>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值& v: j* g& q, }' Z/ i* K/ S
>> y=eval(y) # _+ e; m; ]" l
' R0 q+ z& W8 z+ }; L
y =
* S8 o5 L- z: T6 `% o  W1 _7 ^( }* t1 o0 X" E+ u1 B& u
   1.1932 + 0.0000i9 |- s1 G3 F$ a7 N6 Y& ^/ A( K: q- J
这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
+ ?. k" r  K& ]1 S+ P7 ^5 X; S9 n. ysyms a b c d e;
; o6 z8 a, D; x% b& H3 Qy=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
& h% ^& R$ d+ {1 T# t5 q) Qx=200:0.1:400;%取值范围,步长
* n  T& V* _2 K' W2 o3 u1 L/ [a=1;b=1;c=1;d=1;e=1;%带入参数) S8 l: _/ F& f9 Q" T8 i' z
y=eval(y);" |$ o! A8 l  e
plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~, K2 f- n/ [6 V: g5 D
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 加入华同

本版积分规则

Archiver|手机版|小黑屋|华人同志

GMT+8, 2025-4-4 09:37 , Processed in 0.069287 second(s), 6 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表