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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程$ ?9 P# l: p  {( `9 A

; h9 L( W- A+ }5 i7 ?; {>> clear
% r( U7 E6 Q' Z. x>> syms a b c d e;: G& i& y$ L" D$ A2 J4 p+ v
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')$ W) Y" F9 |. K6 C5 P9 j3 _

. L6 s0 z) ]# H3 Q( p' V, Fy =% g* ~2 }1 m: }2 s4 J4 F

$ L, U- d2 K% Q8 b(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)): w. U- F4 h5 N; D
6 @6 Y- L: x! I5 m7 J  Z' v4 ]7 L
>> a=1;b=1;c=1;d=1;e=1;%若常数已知/ `- [! a( G8 F0 s

$ h$ W% H$ q8 C- m>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
+ x5 G5 {. {1 ?" a# h" a& a* R>> y=eval(y) $ ^8 P- j2 q6 D( x, ^" {

& F; W0 R9 V* @9 \! q# Py =2 A# m6 b% B+ ?

( q* L) H6 \4 v: U   1.1932 + 0.0000i; b5 _  ~& ^6 p1 X$ X
这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
0 e0 m; p; G  A" G- I0 ksyms a b c d e;2 e. D# Z0 c5 y8 W. x
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');  f) E9 U2 C8 t/ a* @5 h% Y$ q
x=200:0.1:400;%取值范围,步长
& [( f; \, D4 ~3 Fa=1;b=1;c=1;d=1;e=1;%带入参数( |! X: N$ `1 a( @, P
y=eval(y);
$ L+ \, }; ~8 t! x7 P+ `plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~
1 R/ A- G4 F8 M$ ~( @* P8 lvb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-10-1 15:31 , Processed in 0.013526 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

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