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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程( w8 e. g$ R+ T, i' u+ E: u2 R$ v
9 @! D. m& |# V: b
>> clear: d4 {  j$ x3 U$ o
>> syms a b c d e;1 L: o' P# w& l7 T2 u
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')% X% y, G; F8 m

+ V; z$ p6 U5 a0 D- q) }y =1 o, S% ?& F0 @; x; X
* x6 {! a: [0 l4 U) F+ \
(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' L7 c: \, H4 R: {
" r5 n8 u  G1 ~
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
9 z, g+ ^6 k2 u# U, F/ N. }* _  Z2 T, U. q9 N8 e( m" z
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值9 [- c0 T7 p8 z) {0 }: O4 ~8 }1 n5 ~
>> y=eval(y) % F1 h: L4 u( Q# s0 N* n7 O

4 b7 w; m+ e  my =# p9 i8 ^5 a' c5 V! T. H1 J! G
* ~& T" P4 L8 j% }& ?
   1.1932 + 0.0000i
! |" Y+ |, R4 u7 L4 r这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。7 q& q7 K! d, N& _
syms a b c d e;8 Q5 n! ^" N6 @5 f$ C# x. |9 q
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');9 n# w+ G" p! ~- h
x=200:0.1:400;%取值范围,步长) r4 h2 J' |5 O' x+ ~! u
a=1;b=1;c=1;d=1;e=1;%带入参数
: [) v  b6 }9 {% W- Q$ @y=eval(y);. S& {$ i3 `. m7 q: ^" F. v
plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~1 z% c8 b& U. r- d. K7 p% n# x" r
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-10-18 12:05 , Processed in 0.056808 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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