数学建模社区-数学中国
标题:
急救!matlab编程?
[打印本页]
作者:
lx.lee
时间:
2009-3-31 00:22
标题:
急救!matlab编程?
Lorenz 模型的状态方程表示为:
' F) }0 C2 R, k* `* l
X1(t)=-βX1(t)+X2(t)X3(t)
% B; h/ s A$ p! z" @: y9 a
X2(t)=-σX2(t)+σX3(t)
! \7 j2 z8 @$ Y/ h7 }" {
X3(t)=-X1(t)X2(t)+ρX2(t)-X3(t)
6 _7 K; Q9 u$ o3 O* ]
取σ=10,ρ=28,β=8/3,初值为X1(0)=X2(0)=0,X3=ε,ε为一个小常
7 A, v- z3 J9 k" U, @( z: N
数,记ε=10e-10,
4 G" d& K. d0 [. J7 \/ ^! f
9 l4 s# A: z7 O! F6 S4 M
这个怎么用微分方程解并画出时间和空间曲线
作者:
sxp
时间:
2009-3-31 14:32
方程组两边同时对 t 求导,得到关于x1',x2',x3'的线形方程组(把x1,x2,x3当参量),解得x1',x2',x3'关于x1,x2,x3的表达式(通过符号矩阵运算)。结果形式为[x1',x2',x3']'=A*[x1,x2,x3]';然后用龙格库塔法求出数值解并画出时间和空间曲线
作者:
lx.lee
时间:
2009-3-31 15:49
非常感谢!
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5