function [xn,dn] = PhaSpaRecon(s,tau,m)% 混沌序列的相空間重構 (phase space reconstruction)% [xn,dn] = PhaSpaRecon(s,tau,m)% 輸入參數：    s          混沌序列%               tau        重構時延%               m          重構維數% 輸出參數：    xn         相空間中的點序列(每一列為相空間中一個點)%               dn         一步預測的目標len = length(s);if (len-1-(m-1)*tau < 1)    disp('err: delay time or the embedding dimension is too large!')    xn = [];    dn = [];else    xn = zeros(m,len-1-(m-1)*tau);    for i = 1:m        xn(i,:) = s(1+(i-1)*tau : len-1-(m-i)*tau);   % 相空間重構，每一列為一個點     end    dn = s(2+(m-1)*tau : end);    % 預測的目標end