function [Lyapunov1]=LargestLyapunov(data,m,tau,P)% 用最小數(shù)據(jù)量法計算混沌時間序列 Lyapunov 指數(shù)% tau = 1;                        % 時延% m = 14;                          % 嵌入維% data ;    % 列向量，即 n*1 taumax = 50;                    % 最大離散步進(jìn)時間% P = 2;                        % 序列平均周期% 作者:Adu,武漢大學(xué),adupopo@163.comY = lyapunov_small(data,tau,m,P);figure(1)plot(Y(1:taumax),'-b'); grid; xlabel('i'); ylabel('y(i)');n=input('請輸入要擬合長度n(默認(rèn)為8)=');if isempty(n)    n=8else   nendlinear_zone = [1:n]';          % 線性區(qū)域F = polyfit(linear_zone,Y(linear_zone),1);Lyapunov1 = F(1);yp=polyval(F,1:n+20);hold onplot(1:n+20,yp,'-r')hold off