第十次作业---洛伦兹模型

---

2013301020024 谢章权

---

• 摘要

  　　研究流体力学的洛伦兹模型问题。

• 背景内容

  　　理想流体力学方程为
  

  $$\frac{dx}{dt}=\delta(y-z)$$
  $$\frac{dy}{dt}=-xz+rx-y$$
  $$\frac{dz}{dt}=xy-bz$$

• 正文

　　　１．模拟洛伦兹模型

　取值为　

$$\delta=10$$ 　　 $$b=\frac{8}{3}$$
当取不同的r值时，模拟结果如下

为进一步观察，我们取r=25
则此时x-z平面模拟图如下
！x-z
关于vpython动态模拟代码如下
[vpython](computationalphysics-2013301020024/10th/Lorenz_model_x_z.py)
- 总结
使用vpython动态模拟洛伦兹模型时发现，随着运算次数的增加，所描绘的曲线越来越 不平滑，这应是vpython自带运算处理算法所致。