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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程8 S1 A* L) c$ H+ P! }
: Y+ E; V: K: y% f2 r2 m( G
>> clear; t' a  n. Q. E' Z# t: @- G
>> syms a b c d e;
' |1 B3 ~7 q6 p, O2 I% r* T>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
7 r5 U$ S" M( p* m5 _& C1 z
) q2 U( G5 D/ Q/ u# T9 `2 X9 U( dy =
; j; b, w% r0 E9 I, o% L4 i0 b  Z* X! | * S2 z) B( Y6 E" Y$ X
(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))- u. f5 H# ]4 e' b; W
2 ~$ A0 B! t* P( `# h$ k& n
>> a=1;b=1;c=1;d=1;e=1;%若常数已知1 R; \2 \) T5 h: b0 G

' `  `" y! f! V/ z$ J) i/ ^) s>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
; r$ H5 c9 o. ~' ~3 P. E>> y=eval(y) & M3 f' u  R+ H

9 h# U  k9 S0 H4 @( `! |! g+ B3 b) s9 my =2 ?# m* {* b. o3 F& s+ M7 j1 L
* H3 d1 q8 T# @" Q' p
   1.1932 + 0.0000i
/ }% M' O# i0 i* @这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
) o* g5 d8 t3 |/ \syms a b c d e;" m, ]9 m+ O; ^8 z
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');$ h9 o. P$ }5 s
x=200:0.1:400;%取值范围,步长
" k- j, H& S8 K8 ba=1;b=1;c=1;d=1;e=1;%带入参数
2 T4 Z. e5 i: [1 o# ~3 n+ J! my=eval(y);
: j8 y( ~. W4 ~  \plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~
& |! Q  d, x( {6 ^: `; {2 Cvb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2025-4-28 15:46 , Processed in 0.066565 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

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