洛伦兹混沌模型.pdfVIP

LORENZ MODEL 黎颖 物理学院 121120053 流程图： 0.001 源程序： T=0.001 N=10000 运行结果图： X-Z 图，横坐标X纵坐标Z ，R=0.5，Y0=1 ，Z0=1，X0分别为1和3 （如无特别说明，以下均取Y0=1 ，Z0=1，图示方框中“1,0,5”代表 X与Z随时间变化图，黑线X,红线Z X0=1 ，R=0.5”,“3,0.5”代表X0=3 ，R=0.5，下同） X,Z随时间变化收敛于一点 当X有较大变化时X,Z随时间收敛于同一点 R=2 X,Z随时间变化收敛于一点 当X有较大变化时X,Z随时间收敛于同一点 R=10 X,Z随时间变化收敛于一点 当X有较大变化时X,Z随时间收敛于同一点 R=15 X,Z随时间变化收敛于一点 当X有较大变化时X,Z随时间收敛于同一点 X ，Z随时间变化不收敛 R=25 1,25对应的图形太小 X 的初值有极小变化时随时间变化X ，Z与变化前值将差距放大， 此图中X0=1时值与X0=1.00001时值相差太大以至于显示X0=1 的地方 很小，几乎看不到，但在运行数据中可以清楚看到偏离的放大 以下分别是R=30,50,100,200的X-Z 图，他们都不是收敛的 探讨R与收敛的关系： 通过如下程序来找到让X,Z 随时间收敛的R取值范围 我对收敛的定义是：当X ，Y ，Z有一定偏离时，随着时间的发展，最终X,Y,Z趋于 一个较固定的点（讨论在THEGMA=10,B=8/3条件下进行） 程序编写思路： 1，对同一个R分别取X0=1 ，Y0=1 ，Z0=1和X0=10 ，Y0=10 ，Z0=10， X0=1 ，Y0=1 ，Z0=1经过较长时间后得到X1 ，Z1 X0=10 ，Y0=10 ，Z0=10经过较长时间后得到X2 ，Z2 当ABS(X1-X2)0.0001且ABS(Z1-Z2)0.0001时认为此R处收敛，输出X2 ，Y2 否则认为不收敛，不输出X,Y值 2，对R从0到100，取0.1为补偿循环 3，X0=10 ，Y0=10 ，Z0=10条件太严格，放松条件在X=10 ，Y0=1 ，Z0=1处得到 较宽松条件下的收敛值 源程序如下： 低要求收敛（X0=10 ，Y0=1 ，Z0=1 ）： 此图展示了通过我编写的程序在给定条件下找到的所有收敛点（即X,Z收敛于此点） 此图横坐标为R，纵坐标为X和Z 的值，黑线X ，红线Z， 没有X,Z对应处即不收敛的R范围 图中不便读数，且易把X,Z较小处误认为是不收敛，下面贴出不收敛处及附近数值 从左到右三幅图，每幅图中有三列数据，第一列为R，第二列为X ，第三列为Y ， X ，Y无值处不收敛 高要求收敛（X0=10 ，Y0=10 ，Z0=10 ）： 与低收敛要求相比减少了一部分R区间 图1 图二 图三 图四 虽然这样做是有一定问题的，比如收敛的定义 就不太严格，但是我们至少可以看出，当R值较大时， 尤其是是大于二十几时表现出不收敛，混沌的特性 最后，贴上这里面最漂亮的一张图 Z R值 X R=1，2，3,4,5,6,时，由（1,1,1）随时间演变的轨迹图