% 計(jì)算混沌時(shí)間序列 Lyapunov 指數(shù) - Logistic 序列% 使用平臺(tái) - Matlab6.5 / Matlab7.0% 作者：陸振波，海軍工程大學(xué)% 歡迎同行來(lái)信交流與合作，更多文章與程序下載請(qǐng)?jiān)L問(wèn)我的個(gè)人主頁(yè)% 電子郵件：luzhenbo@yahoo.com.cn% 個(gè)人主頁(yè)：http://luzhenbo.88uu.com.cnclcclearclose all%-----------------------------------------------------------------% 產(chǎn)生 Logistic 混沌時(shí)間序列lambda = 4;k1 = 10000;             % 前面的迭代點(diǎn)數(shù)k2 = 10000;             % 后面的迭代點(diǎn)數(shù)x0 = rand(1);x = zeros(k1+k2,1);for i = 1:k1+k2    x0 = lambda*x0*(1-x0);    x(i) = x0;endx = x(k1+1:end);%-----------------------------------------------------------------% Logistic 入口參數(shù)tau = 1;                % 時(shí)延m = 2;                  % 嵌入維taumax = 30;            % 最大離散步進(jìn)時(shí)間P = 100;                % 序列平均周期fs = 1;                 % 采樣頻率%-----------------------------------------------------------------Y1 = Lyapunov_rosenstein_2(x,tau,m,taumax,P,fs);figuresubplot(211); plot(Y1); grid; xlabel('i'); ylabel('y(i)'); title('Logistic')subplot(212); plot(diff(Y1)); grid; xlabel('n'); ylabel('slope');linear_zone = [1:7]';  % 線性區(qū)域F1 = polyfit(linear_zone,Y1(linear_zone),1);Lyapunov_2 = F1(1)                  % 以 2 為底的指數(shù)Lyapunov_e = Lyapunov_2/1.4427      % 以 e 為底的指數(shù)