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

谁会编程解二阶微分方程

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

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程4 @( R2 L2 e2 x: N4 ^8 q

9 \3 _/ f: u. b>> clear
( `' U6 O: {- O- [4 `' [>> syms a b c d e;
$ |2 R7 O2 H0 K>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
1 a* @7 [- ]# X$ ~  G # F4 H/ ?* ]# F( K3 U. }$ w
y =
# Y& _9 |, G+ I9 }0 y
: m, e0 O4 x  N(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))
9 r: X7 G6 T9 `4 E
& {  s! q2 _) Z6 C. K& U>> a=1;b=1;c=1;d=1;e=1;%若常数已知& b! C" \7 R7 q7 }# O2 n2 b# R" A
& [- m  {: ]' G9 N9 z# G+ d4 x
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值7 M( W6 B- Q! F  [
>> y=eval(y) ! W* K8 s3 G  P8 a+ n
  I, f/ }/ |  G4 F) l; b7 R9 Y( g  F. U
y =6 E$ a' r! t7 ^7 A

& t# N" s" l3 V2 Q  l, z6 B  k, L5 G   1.1932 + 0.0000i
& V5 V0 H3 d* ^6 ?1 D: F这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
# i7 T9 J1 [3 `4 l3 Tsyms a b c d e;
) L" W8 H( `+ q- l8 `) P5 }( ]y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');, k0 {3 b' ]$ F4 b7 ?' T  ]: w0 [
x=200:0.1:400;%取值范围,步长6 Y+ f5 q# ^. K: H
a=1;b=1;c=1;d=1;e=1;%带入参数4 y) o1 a# J4 `% U+ ^* H- M) a
y=eval(y);/ l2 C' x3 E5 Q3 s$ O+ I
plot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~" h$ y% Q4 [' d  |
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-8-1 18:08 , Processed in 0.058774 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

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