% 計算混沌時間序列 Lyapunov 指數(shù) - Logistic 序列% 計算 Henon 映射 x(n+1) = 1 - a * x(n)^2 + y(n); y(n+1) = b * x(n)clcclear allclose alldisp('********** Logistic ***********')%-----------------------------------------------------------------% 產(chǎn)生 Logistic 混沌時間序列lambda = 3.5;x = 0.1;k1 = 400;                   % 前面的迭代點(diǎn)數(shù)k2 = 10000;                 % 后面的迭代點(diǎn)數(shù)f = zeros(k1+k2,length(lambda));for i = 1:k1+k2    x = lambda .* x .* (1 - x);    f(i,:) = x;endX = f(k1+1:end,:);%-----------------------------------------------------------------% Henon 入口參數(shù)fs = 1;                 % 采樣頻率t = 1;                  % 時延d = 2;                  % 嵌入維tmax = 30;              % 最大離散步進(jìn)時間p = 10;                % 序列平均周期%-----------------------------------------------------------------% 調(diào)用計算函數(shù)Y = Lyapunov_rosenstein_2(X,fs,t,d,tmax,p);I = [0:length(Y)-1]';linear = [2:9]';  % 線性區(qū)域F1 = polyfit(I(linear),Y(linear),1);Y1 = polyval(F1,I(linear),1);%-----------------------------------------------------------------% 結(jié)果顯示figuresubplot(211); plot(I,Y,'k.-'); grid; xlabel('i'); ylabel('y(i)'); hold on;plot(I(linear),Y1,'r-'); hold off;subplot(212); plot(I(1:end-1),diff(Y),'k.-'); grid; xlabel('n'); ylabel('slope');Lyapunov1_2 = F1(1)Lyapunov1_e = F1(1)*log(2)