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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程
% E% c" G; }1 O- }
  e( S" h8 m0 @5 ?# A>> clear
% c9 }$ g+ m# N3 `>> syms a b c d e;
, ^/ B! v; G0 P3 K' e" O3 X- _% E1 l  T>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
# E1 P: S) S+ R4 }3 f# R ; m1 F2 a- f( w
y =; U8 y: W" m" r

  y$ t: @+ R8 z  J(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))
- y0 h% W4 g& e! ]# N* o2 F' R' A
, p; N; A; Q: S+ {>> a=1;b=1;c=1;d=1;e=1;%若常数已知& N1 m& U; W- Q, @! C8 |" F
8 I* P7 V3 m1 i/ c3 H
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值+ \; {) l$ x" t* w9 B/ _) K
>> y=eval(y) , F7 l& A6 ~4 A
9 ^; H5 k  b  t- C, h
y =
4 Z- M% F2 _$ ?8 d8 T
) ?' o. i/ N# H7 q6 o# [1 Z" n. @1 o   1.1932 + 0.0000i
( ^8 z  p4 ]0 g# n4 b. e这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
! ^4 @# _5 G4 b6 X5 E+ P0 L' |: a; Nsyms a b c d e;
4 D9 u; w- h) ry=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
) V+ Q/ f% d7 @" C( dx=200:0.1:400;%取值范围,步长  ?9 ~; w: [5 y+ p9 O5 j- J
a=1;b=1;c=1;d=1;e=1;%带入参数
$ W' z; l/ W* O- s; G8 b& Ly=eval(y);
& N6 c7 c, u+ Q5 i$ eplot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~& G! [. Z' H1 L( d2 e* @( K
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-11-1 07:31 , Processed in 0.028218 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

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