数学建模社区-数学中国

标题: 急救！matlab编程？ [打印本页]

---

作者: lx.lee    时间: 2009-3-31 00:22
标题: 急救！matlab编程？
Lorenz 模型的状态方程表示为：
7 u) O4 b8 ^9 T7 p0 |* s6 l: y: EX1(t)=-βX1(t)+X2(t)X3(t)
X2(t)=-σX2(t)+σX3(t)
X3(t)=-X1(t)X2(t)+ρX2(t)-X3(t)
取σ=10,ρ=28,β=8/3,初值为X1(0)=X2(0)=0,X3=ε,ε为一个小常
' Z& }5 }% @- W8 F- Z数，记ε=10e-10，
" P. |0 S. J' p" G7 N6 B9 A
- m  @0 L9 X& S这个怎么用微分方程解并画出时间和空间曲线

---

作者: 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