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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程6 t) e. }# O5 Z8 z& h/ x" I6 V

. O. \$ ]" R3 T( F>> clear8 H1 |+ P; `9 A
>> syms a b c d e;
/ ?# `" F. W5 R- L: H>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')( o( g+ L: O5 k4 J1 s# Q) V

: c6 x( U- f& w2 h+ ]9 v- h/ }y =
/ |, d4 i! r! {1 R
2 [- l% `7 o# S* v% c0 {, B1 ]6 @: 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))0 w& b" k7 Q) \0 J

. G: G6 n  d8 q  {>> a=1;b=1;c=1;d=1;e=1;%若常数已知, X2 w2 @) a0 F
. E: z: ~4 J# N' [* y7 M0 d
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值; d2 a- d  b" C- C" n
>> y=eval(y) 7 G. Q5 Q$ O/ m! O

1 @* V1 d4 t1 Ry =
: n7 D7 i7 q9 J5 e' Y$ {  g
5 s1 p3 Q) U2 }: ?1 A1 R9 ]$ M   1.1932 + 0.0000i
6 a6 J4 X+ E( ^. S0 T1 x这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
4 a: h  K0 V! }! gsyms a b c d e;
! ?; j5 ^7 K1 ?, }y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
; p! c- Q# |  qx=200:0.1:400;%取值范围,步长' y' u; x. `- h2 q' `
a=1;b=1;c=1;d=1;e=1;%带入参数
: K' M4 @9 V5 c5 L' uy=eval(y);
3 n" Z9 f6 K# c8 {plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~+ A8 S! R9 S9 O. r# J; w9 L
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2025-7-18 10:17 , Processed in 0.133395 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

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