?? main_largestlyapunov_rosenstein2.m
字號(hào):
% 計(jì)算混沌時(shí)間序列 Lyapunov 指數(shù) - Logistic 序列% 計(jì)算 Henon 映射 x(n+1) = 1 - a * x(n)^2 + y(n); y(n+1) = b * x(n)clcclear allclose alldisp('********** Logistic ***********')%-----------------------------------------------------------------% 產(chǎn)生 Logistic 混沌時(shí)間序列lambda = 4;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; % 時(shí)延d = 2; % 嵌入維tmax = 30; % 最大離散步進(jìn)時(shí)間p = 10; % 序列平均周期%-----------------------------------------------------------------% 調(diào)用計(jì)算函數(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)
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -