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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程( l  @* h- o- w9 M
( X& W8 f7 E+ J3 n+ T0 n* r
>> clear' e# f( }9 P6 i  N! F
>> syms a b c d e;* @9 t: l- O" E* Y
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')7 C; ?0 M, h- U# ~+ ]: J
% a% w" {5 w) e  e, P* y. R$ M
y =/ f. |: N8 j" x' p4 ]4 |+ K

. L0 v+ j/ q" R9 X' V+ x(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))
  v# X# F: k& w! S# D' [  f% E8 Z
" f; Q; Z0 K5 Q0 b- ?1 A>> a=1;b=1;c=1;d=1;e=1;%若常数已知
. i9 x" W0 C1 V( C8 C! _6 P
: S/ Q9 G4 n$ G1 i1 O* z  {, l2 f>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
0 o3 A. r7 I' i% G>> y=eval(y) 6 J- B; N9 l; f6 t) h- f* N

4 K: v# s9 o* k! ]- Ey =0 b* a2 K! Q. y7 D' x( F
  x0 ~  x, o" R* y
   1.1932 + 0.0000i
( z+ b% e& ]+ o5 D% J  r' M这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
$ B# T+ c' @. @3 ~syms a b c d e;
  h! d! h/ X$ S8 oy=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
7 v. I, H( q4 M7 C& Y; Jx=200:0.1:400;%取值范围,步长1 b" d2 I5 f6 Z# Q% @
a=1;b=1;c=1;d=1;e=1;%带入参数- {9 Q) _3 T* D- v, d' E0 h: u& G
y=eval(y);
9 @! e( N% O- I; [8 Splot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~
) }# ^5 G/ u5 i  Q# Avb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-5-6 05:00 , Processed in 0.052599 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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