yunr
2009-01-27, 23:07
有一个matlab 画图题困扰me一整天了
问题描述，用以原点为圆心的单位圆为初始条件，进行henon映射，循环10000次。
我自己写的程序如下：
单位圆
alpha=0:pi/50:2*pi;
r=1;
x=r*cos(alpha);
y=r*sin(alpha);
plot(x,y,'-');
axis equal
--------------------------
henon映射
a=1.4;
b=0.3;
x(1)=0.4;
y(1)=0.4;
for i =1:10000
x(i+1)=1-a*x(i)^2+y(i);
y(i+1)=b*x(i);
end
plot (x,y,'.')
问题所在： 这两个程序单独运行时都正确，但是，第二个程序的初始条件是(x1,y1)--(0.4,0.4). 而并非想要的初始条件单位圆。该怎么样把单位圆作为初始条件放在henon映射程序中呢？也就是把0.4换成单位圆。。。 请
筒子们帮个忙拉，多谢奥 :heart:
问题描述，用以原点为圆心的单位圆为初始条件，进行henon映射，循环10000次。
我自己写的程序如下：
单位圆
alpha=0:pi/50:2*pi;
r=1;
x=r*cos(alpha);
y=r*sin(alpha);
plot(x,y,'-');
axis equal
--------------------------
henon映射
a=1.4;
b=0.3;
x(1)=0.4;
y(1)=0.4;
for i =1:10000
x(i+1)=1-a*x(i)^2+y(i);
y(i+1)=b*x(i);
end
plot (x,y,'.')
问题所在： 这两个程序单独运行时都正确，但是，第二个程序的初始条件是(x1,y1)--(0.4,0.4). 而并非想要的初始条件单位圆。该怎么样把单位圆作为初始条件放在henon映射程序中呢？也就是把0.4换成单位圆。。。 请
筒子们帮个忙拉，多谢奥 :heart: