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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程
$ S( n: C% j$ b
) j/ t- e( ~9 I) \, O# d% |' c3 j>> clear! j! ^3 V! K  O" T. ~2 M& p+ S
>> syms a b c d e;
9 N5 l# r( o/ p7 |7 E- F>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')0 y! L/ |( O6 f0 |  e) Q8 Q, o: y
8 ]' E* B! l, }0 ?7 k
y =
* V+ o# g. P0 |  B$ j* v( K: q4 D
. `# D5 D! R" U$ M+ W9 V3 [(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 i) Y, \6 q$ E: O+ i( d& Y. Y; Z
/ W$ k$ t# g7 n0 e) m- F>> a=1;b=1;c=1;d=1;e=1;%若常数已知; W% O8 f& W8 S% w, i
( ~2 a" v; _9 t5 N  O. \
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值2 l, n- ]; M0 W
>> y=eval(y) * E7 P5 V, S; j& w1 g: t
1 N0 y* F/ {) k6 `4 q; N
y =
; w) o2 Z' j6 q, O$ U& {6 D! l
1 x) a- E3 B, Y, H: C   1.1932 + 0.0000i
4 A: z, F4 q: w这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
& Y9 l# d7 q$ o9 ~syms a b c d e;
- ?5 L" [3 E1 [0 |2 D+ cy=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');% }1 w. ?: e6 {
x=200:0.1:400;%取值范围,步长
+ k: L1 A. {; l/ _4 W  j& ^a=1;b=1;c=1;d=1;e=1;%带入参数$ \9 b) {! ]. [. B$ ?/ G4 q. w
y=eval(y);
) l2 j# A7 }) l/ J0 g2 F. tplot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~
& z: m% Z* ?9 y! Z- o5 Tvb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2025-2-19 06:16 , Processed in 0.066527 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

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