马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
时间序列相空间重构程序
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
上面的程序运行重构的只是2维图,如何重构三维图形,请高手指教!!谢谢!! |