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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程
9 Z  d6 O* j1 F* j/ t
* E' E3 G/ J4 Q5 m. C+ i>> clear
; ~5 Y3 }8 p: b+ h>> syms a b c d e;
! d  L- q: o! k- r4 w>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')) F! T) \4 e3 P3 Q/ O

  Y+ l( e/ X% m9 g8 W' K! N+ Ry =
4 A2 J' ^+ h2 h) T # n( n, E/ I& @8 U1 D9 Z
(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))
, i+ g  `3 w. M8 S  q( z
  X. ?$ _+ L4 X>> a=1;b=1;c=1;d=1;e=1;%若常数已知4 j  i) {! N+ d% Q! [; S

) R( d+ e0 k5 i, T( Q>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值$ {  X1 _% M/ H7 h4 r( m
>> y=eval(y)
: f  u' U- k9 E) s3 G  N# B0 [/ c
% g8 |/ L3 l! m5 qy =0 p3 e6 d) n2 L4 R  y8 _

4 a( L; J, [/ W1 G, A8 K/ R% Q   1.1932 + 0.0000i
1 I! U. o' t  \/ b# Q) U* e这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
# p. x- j& L' |4 E6 Qsyms a b c d e;
9 t8 C/ g, N( n  d. P; oy=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
5 J5 `7 _# N( F3 Q: W2 h3 m# K7 Ux=200:0.1:400;%取值范围,步长
' v' s; v! y" qa=1;b=1;c=1;d=1;e=1;%带入参数
! @& a# |7 W3 W5 _2 T" `$ Y! ky=eval(y);
6 h$ L# G  I* ]" r- G& w) D6 Yplot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~4 T, D. `* ~6 g9 q# }- h
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-11-1 07:29 , Processed in 0.009860 second(s), 6 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

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