数学建模社区-数学中国

标题: 急救!matlab编程? [打印本页]

作者: lx.lee    时间: 2009-3-31 00:22
标题: 急救!matlab编程?
Lorenz 模型的状态方程表示为: 6 O4 i8 C. [6 B" v
X1(t)=-βX1(t)+X2(t)X3(t)
, Q2 d7 n# U! q! K% w' nX2(t)=-σX2(t)+σX3(t)
. B" O/ b; k7 ~X3(t)=-X1(t)X2(t)+ρX2(t)-X3(t)
: z7 x( `0 r3 Z! k取σ=10,ρ=28,β=8/3,初值为X1(0)=X2(0)=0,X3=ε,ε为一个小常 ) y& K2 C3 s/ D( @3 A$ }! W7 ]6 `
数,记ε=10e-10, % j# X0 h; h: S: I5 r6 t2 z' x8 ?
. k6 f( Q. \4 |8 Q( z3 W, H3 n
这个怎么用微分方程解并画出时间和空间曲线
作者: 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