- 1、本文档共13页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
混沌与非线性科学作业5670718120980刘艳作业5
《混沌与非线性科学》作业5、6、7
0718120980 刘艳
作业5
1964 年,法国天文学家伊侬(Henon)从研究球状星团中提取出Henon 映射:
⎧ 2
x 1 by ax
+ −
⎨ i+1 i i
⎩ y i+1 xi
( )
a,b 为参数,现在取b 0.3 ,a 0.5,1.0和1.4 ,三种情况讨论 x, y 的运动轨
迹。
解:采用Matlab 程序实现:
function f=h(x1,y1,a,b)
x(1)=x1;
y(1)=y1;
for i=1:5000
x(i+1)=1+b*y(i)-a*(x(i))^2;
y(i+1)=x(i);
end
plot(x,y,b.);
disp(y)
( ) ( ) a 0.5,b 0.3
Case1. 给定初始值为 x 0 ,y 0 0,0 , 时,所得Henon 映射的运
动轨迹为:
图1a Henon 映射运动轨迹
稳定值
(x4 ,y 4 )
(x 2 , y 2 )
(x 3 , y 3 )
(x 5 , y 5 )
(x0 ,y 0 )
(x 1 , y 1 ) 稳定值
图1b 迭代过程细节图
( ) ( )
由图1b 可以看出,从初始值开始迭代分为两个岔道,由x , y 、x ,y 、
1 1 3 3
( ) ( ) ( )
x5 , y 5 等下标为奇数的点构成一条分支,而下标为偶数的点 x 2 ,y 2 、x 4 ,y 4 、
(x 6 ,y 6 )等构成另外一条分支。且最终都收敛于各自的稳定值。具体的数值有:
x 0 x 1 x 0.5 x 1.175 x 0.4597
⎧ 0 ⎧ 1 ⎧ 2 ⎧ 3 ⎧ 4
⎨ ⎨ ⎨ ⎨
文档评论(0)