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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程
6 V! f7 i+ K5 E* N0 D( l  }9 S  f- H6 Q" J
>> clear
9 [/ ~/ A5 H( C# P>> syms a b c d e;
, t7 o4 ?: s7 f* K# ~+ L4 |>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')+ p+ \# l: m( X
2 G& G  P3 h% T4 s9 k- K
y =- z" K- {' D; `* v& \/ Z$ R0 `

( C1 J2 M: V4 z' y0 e. [8 O(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))
" @. y: _8 \- h$ @
4 |& }/ z, Z# v6 t/ |  o5 B>> a=1;b=1;c=1;d=1;e=1;%若常数已知7 S% w& s: \. F6 h& U, ^

; o6 L- b$ F5 B* b$ e>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
, w% r7 }. Q7 R>> y=eval(y) ' A0 d, K- p; ?: {: ]
. z8 H1 E2 T0 d" E4 I# V
y =
' l: z) _1 S5 d  w
  w! W) T; P3 S- T- g! v8 Z   1.1932 + 0.0000i
, R, C+ d' u$ q0 }这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。$ W; ^$ |& q% ?' k
syms a b c d e;
2 f( m+ L0 r" y" W! A6 w. Ay=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
% b; S( J4 `/ A+ t3 _/ x& H' v6 xx=200:0.1:400;%取值范围,步长
* I4 Q( ~' y- w, a) Sa=1;b=1;c=1;d=1;e=1;%带入参数" a3 \4 Q9 `7 _( |( F& H# k  z
y=eval(y);1 ~) j  \9 c/ y) G
plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~0 N: D; Y2 X( r
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2025-4-2 16:08 , Processed in 0.068203 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

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