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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程- `" J4 d0 w3 x8 f4 E  B4 {

8 R9 v3 G2 s: T9 h0 t2 p: U>> clear- r- A) m( Z; H4 [3 K
>> syms a b c d e;0 }' f  ?( h9 k: F
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')  w+ ~1 j5 J' `$ E3 D! S' W

! r0 a8 e5 j7 ty =1 S& z% J; t9 Q5 G

1 g7 D8 [1 d  r: D4 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))7 I3 m8 c/ g3 u) b6 J+ _

9 f3 D4 w5 q4 }  J. f>> a=1;b=1;c=1;d=1;e=1;%若常数已知
" l" x) K, e8 `: k/ @$ k: C
# w+ O8 z1 p, v8 f- O' T( _>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值! e! ^- v) W3 g& n
>> y=eval(y) / h9 H& N* w0 b
* p: v1 i% f# |- U5 i
y =
( c" i! e" n4 ]# O. d; i( X# i. Z, Q! n8 T& f, H- i4 N
   1.1932 + 0.0000i7 k6 F; x! L1 d7 H  _. c# ?
这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。3 I# U/ Q7 H& m& _* A& [1 R; \
syms a b c d e;) q2 C! v/ o1 X8 w* `2 C- ~" j
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');- U' F: V) P' C2 L6 O/ B
x=200:0.1:400;%取值范围,步长8 X( s: r" d6 n4 H, J
a=1;b=1;c=1;d=1;e=1;%带入参数
+ G9 z( T% Y  O, ?y=eval(y);0 o% n' x- o$ @$ F: f
plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~
# e/ S( ~& }+ o: ^2 Ivb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-7-3 09:31 , Processed in 0.071168 second(s), 6 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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