相空间重构 matlab 代码

function [xndn] = PhaSpaRecon（staum）
% 混沌序列的相空间重构 （phase space reconstruction）
% [xn dn xn_cols] = PhaSpaRecon（s tau m）
% 输入参数：    s          混沌序列（列向量）
%               tau        重构时延
%               m          重构维数
% 输出参数：    xn         相空间中的点序列（每一列为一个点）
%               dn         一步预测的目标（行向量）

[rowscols] = size（s）;
if （rows>cols）
    len = rows;
    s = s‘;
else
    len = cols;
end

if （nargout==2）
    
    if （len-1-（m-1）*tau < 1）
        disp（‘err: delay time or the embedding dimension is too large!‘）
        xn = [];