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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程! x# k. j+ m  ~) ]

- b/ t5 ?$ Y# P) b5 O. ^>> clear
+ u0 G+ y. R# j& u5 t; Z4 n>> syms a b c d e;* o$ J8 h) @1 H; v  \% F0 s
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')/ T7 E& M5 a8 Z

' C- U5 {& {  iy =
/ g- a9 A$ J+ W  d$ j$ r! i
  K0 q1 c) `! ~3 {/ F: U, q3 A(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))( f$ r# l, f  ?4 S) [1 E
$ q# Q: H0 H$ M# w5 m1 T% a9 C8 L
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
& @! `9 E2 s+ p7 e
- c. E1 G8 X6 i% k) x>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
0 W% ~. E- p; G# x+ V4 }5 v>> y=eval(y) , U- _8 L! p2 n3 S5 a1 T
9 }( N8 B) [$ H: h( _( m% u# x
y =# O2 v. F4 i% d6 R0 V: b( x7 T& E
/ a/ y- k$ S. n; F
   1.1932 + 0.0000i
* J+ x/ K' P4 B8 y. a3 \这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
- @" b1 m& p3 w6 o5 gsyms a b c d e;
/ Y# V$ z5 z/ d6 b- y& t! ~y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');7 i5 C% t3 b7 Z1 n- a
x=200:0.1:400;%取值范围,步长
2 W+ ^1 y7 u; z, k: Q2 v( Ma=1;b=1;c=1;d=1;e=1;%带入参数
9 Y& }! t) a5 z* z7 p* [' my=eval(y);
3 T5 s9 K1 a  E: t+ k9 |! splot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~
# {8 S1 J$ h) S; {. C; pvb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2025-4-2 01:35 , Processed in 0.066329 second(s), 6 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

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