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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程. u. S  F& c3 N( ?% q

# q' f- ^0 e1 b4 M) [>> clear
+ _8 v1 W! I: G* ?>> syms a b c d e;
. V1 a9 y0 \3 I8 L$ T2 V>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')4 j5 |1 c4 K$ Q- N: [, N! [+ X' K
/ I* h1 C  D+ J7 t* J3 x. Z: s" t
y =$ n/ O# O0 Y( V1 f- R

+ K3 Z$ e6 d  m. _(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))
! b6 s+ {1 g- v; G$ o3 w0 F 6 @  `3 L' X: ~2 }' S, S9 r
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
* P  ]% }, f* c/ [/ j. R+ `# F# y7 r/ ]9 @3 t7 X! m
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
; F7 t$ e: C. o) p>> y=eval(y)
. U0 Y# n( p( J  o0 b( b3 G7 F- n4 K6 ]
y =
7 `' H5 b, S( s$ g2 r
( n1 K1 f$ S& S   1.1932 + 0.0000i3 J2 Q2 c' N( {$ n% z* Z  i
这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
6 Z" b+ L! F& hsyms a b c d e;
: u# y; a2 _9 M" N" _# ey=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');1 {' U* Y( B4 `/ w/ B/ O, R6 `: z
x=200:0.1:400;%取值范围,步长
/ q2 I. A& j, L7 |1 ya=1;b=1;c=1;d=1;e=1;%带入参数
; c2 ]7 `4 X' d6 N+ ?- [8 ]y=eval(y);
/ V2 _+ c, d& j3 A* nplot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~: ?/ Q" Q3 b$ n  D# ~% L
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-6-27 13:06 , Processed in 0.068854 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

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