前几天一个同事看到我写的双摆程序后问我:"这吊儿郎当的是什么玩意?",我当时给了他一个白眼说:"没文化,No education.这么高逼格的东西,你竟然说吊儿浪荡!"
上一篇文章是关于“三体”的混沌图像,这一篇就讲与之有些相似的“双摆”。双摆是物理学中的一个概念,依稀记得中学物理讲过单摆。先给下单摆与双摆的定义:
单摆:由一根不可伸长、质量不计的绳子,上端固定,下端系一个质点的装置。
双摆:是一个摆的支点装在另一摆的下部所形成的组合物体。双摆有两个摆角,所以有两个自由度。双摆是多自由度振动系统的最简单的力学模型之一,它也是一种混沌实例。
那么问题来了:双摆这种混沌系统,如果给定初始条件,以后任意时刻的运动状态都可以算出来吗?
答案是:理论上可以,但实际计算中,由于误差,会使运算结果与实际相差越来越大。这一点与三体很像,非线性数学就是这么任性。
下面四个GIF动画图像为双摆的录屏:
这里使用我定义语法的脚本代码生成混沌图像。相关软件参见:YChaos生成混沌图像。如果你对数学生成图形图像感兴趣,欢迎加入QQ交流群: 367752815.图像的生成方法是:使用第二个摆点的运动轨迹生成图像.
先上脚本代码:
[ScriptLines]
s=[static]p
w=[static]q
a=l*l*(m+n)
b=l*n*h*cos(s-w)
c=l*(n*h*sin(s-w)*v*v + (m+n)*g*sin(s))
d=n*h*l*cos(s-w)
e=n*h*h
f=n*h*(-l*sin(s-w)*u*u + g*sin(w))
i=(c*e - b*f)/(b*d - a*e)
j=(a*f - c*d)/(b*d - a*e)
u=u+i*t
v=v+j*t
s=s+u*t
w=w+v*t
x=l*sin(s) + h*sin(w)
y=-l*cos(s) - h*cos(w)
[Variables]
g=9.800000
h=2.000000
l=1.000000
m=1.000000
n=2.000000
p=1.000000
q=2.000000
t=0.000500
u=0.000000
v=0.000000
这里使用的是角度变化实现双摆,参考的资料是:
http://sebug.net/paper/books/scipydoc/double_pendulum.html
看图像:
相关软件:
相关文章: