function [xn,dn] = PhaSpaRecon(s,tau,m,T)% 混沌序列的相空間重構(gòu) (phase space reconstruction)% [xn, dn, xn_cols] = PhaSpaRecon(s, tau, m)% 輸入?yún)?shù)：    s          混沌序列(列向量)%               tau        重構(gòu)時(shí)延%               m          重構(gòu)維數(shù)%               T          直接預(yù)測(cè)步數(shù)% 輸出參數(shù)：    xn         相空間中的點(diǎn)序列(每一列為一個(gè)點(diǎn))%               dn         一步預(yù)測(cè)的目標(biāo)(行向量)[rows,cols] = size(s);if (rows>cols)    len = rows;    s = s';else    len = cols;endif (nargin < 4)    T = 1;endif (nargout==1)    if (len-(m-1)*tau < 1)        disp('err: delay time or the embedding dimension is too large!')        xn = [];    else        xn = zeros(m,len-(m-1)*tau);        for i = 1:m            xn(i,:) = s(1+(i-1)*tau : len-(m-i)*tau);   % 相空間重構(gòu)，每一行為一個(gè)點(diǎn)         end    end    elseif (nargout==2)        if (len-T-(m-1)*tau < 1)        disp('err: delay time or the embedding dimension is too large!')        xn = [];        dn = [];    else        xn = zeros(m,len-T-(m-1)*tau);        for i = 1:m            xn(i,:) = s(1+(i-1)*tau : len-T-(m-i)*tau);   % 相空間重構(gòu)，每一行為一個(gè)點(diǎn)         end        dn = s(1+T+(m-1)*tau : end);    % 預(yù)測(cè)的目標(biāo)    end    end