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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程
; K( G, D* ~9 {4 k
1 [; H) Z6 W$ u0 j# `3 J>> clear3 H- y2 q& I% {% a0 f
>> syms a b c d e;. S  l; g* h+ ]/ z+ I2 t
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
; X- L( ~$ o. E$ F. w; N( ?, p
' }- l  I8 {8 U7 s6 Ry =1 Q3 ?$ ~( Q% f1 g3 G
& j9 y% Q8 w" r+ u3 R/ `( P, 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))
+ Z6 P1 c/ I/ D$ T5 v* W9 `/ ]
; p0 z; U: T9 Q8 V* a>> a=1;b=1;c=1;d=1;e=1;%若常数已知
5 B4 l( J, x/ H. P% z! F- W
/ X7 C  t4 Z9 Z; q% ]! h>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值: [$ Y8 Z. D9 C5 s
>> y=eval(y)
, c* s7 M" `8 L3 _5 x0 w& R2 S) b- U1 c4 o) G; ^" a
y =
* D5 d) g6 I! c3 c# y
, h4 c2 R' |  d' L0 ], [   1.1932 + 0.0000i
, \1 a- }& g; b+ \, p4 B% F9 @这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。& j' ]6 v" u- ~3 l
syms a b c d e;
. l# H- |8 _7 n5 i9 U7 R8 h+ s0 Wy=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
% s. C, G$ K. \( R$ o! T& Xx=200:0.1:400;%取值范围,步长
6 w: V; p, {0 h9 m4 ^( Wa=1;b=1;c=1;d=1;e=1;%带入参数
. z6 N, y, y. }# ky=eval(y);
+ B4 X/ Y4 {- tplot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~- O& |8 ^' {2 }" m( g
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-10-24 17:47 , Processed in 0.010223 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

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