转载

混沌图像---杜芬的凶罩

杜芬振子Duffing Oscillator是一个描写强迫振动的振动子,它是非线性数学中的一个经典方程。至少它有什么用，我也不清楚，感觉它与物理上的弹簧,单摆,双摆有点关系。而我只关心它的图形结构，杜芬振子看上去是双环路。很多混沌吸引子都是双涡旋结构，如洛伦兹的蝴蝶，陈氏吸引子的飞蛾，蔡氏电路的漩涡等，而杜芬振子的形状看着很像凶罩。

这里使用自己定义语法的脚本代码生成混沌图像.相关软件参见:YChaos生成混沌图像.如果你对数学生成图形图像感兴趣,欢迎加入QQ交流群: 367752815

[1]

脚本代码：

[ScriptLines] s=s+t u=y v=a*cos(c*s) - b*y + x - x*x*x x=x+u*t y=y+v*t  [Variables] a=2.078602 b=0.100000 c=0.490400 s=0.000000 t=0.001000 x=1.000000 y=1.000000

图像:

混沌图像---杜芬的凶罩

这图像看上去并不好看,因为随着迭代的增加,杜芬最终到达一个闭合曲线上.

[2]

于是我对其做了些修改,加上随机扰动.

脚本代码：

[ScriptLines]

s=s+t

u=y+rand2(-r,r)

v=a*cos(c*s) - b*y + x - x*x*x

x=x+u*t

y=y+v*t

[Variables]

a=2.090000

b=0.100000

c=0.500000

r=1.000000

s=0.000000

t=0.001000

x=1.000000

y=1.000000

图像:

混沌图像---杜芬的凶罩

修改扰动的范围,可生成不同的图像:

混沌图像---杜芬的凶罩

[3]

前面的图像并不像凶罩,那么下面将其变得扁一些:

[ScriptLines] s=s+t u=y/h+rand2(-r,r) v=a*cos(c*s) - b*y/h + x - x*x*x x=x+u*t y=y+v*t*h  [Variables] a=2.090000 b=0.100000 c=0.500000 h=0.500000 r=1.000000 s=0.000000 t=0.001000 x=1.000000 y=1.000000