這個問題的背景是壹個關於混沌的叫做Hénon映射的離散時間動態系統模型。
其中,經典的Hénon圖取a?=?1.4,b?=?0.3,這是題目中給出的數據。
(1)根據需要編寫函數:
功能?【x,?y]?=?henon(x0,?妳呢?n)
x?=?[0;?zeros(N,1)* NaN];
y?=?x;
為了什麽?我?=?1?:?普通
x(i+1)?=?1?+?y(i)?-?1.4*x(i)^2;
y(i+1)?=?0.3?*?x(壹);
end的輸入參數分別為初始值x0、y0和叠代次數,輸出為叠代結果。
(2)調用函數
x0?=?0.1;?y0?=?0.2;
【x,?y]?=?henon(x0,?妳呢?100);
disp([x(end)?Y(end)])結果如何?xk = -1.2493,yk = 0.3805。
(3)調用函數
x0?=?0;?y0?=?0;
【x,?y]?=?henon(x0,?妳呢?30000);
plot(x,y,?'.')
xlabel?x,?伊拉貝爾?y