m=
tau=
P=N-(m-1)*tau
for ii=1:m
for jj=1:P
data(jj,ii)=x(jj+(ii-1)*tau)
end
end
y1=data(:,1)
y2=data(:,2)
y3=data(:,3)
figure(1)
plot3(y1,y2,y3)
view(2)
xlabel('x'),ylabel('粗梁y'蚂罩)
title('Lorenz重构闷凳闹吸引子图')
function X=reconstitution(data,N,m,tau)%该册扰函数用来重构相空间
% m为岩樱嵌入空间维数
% tau为时间延迟
% data为输入时间序列
% N为时间序州枣旦列长度
% X为输出,是m*n维矩阵
M=N-(m-1)*tau%相空间中点的个数
for j=1:M %相空间重构
for i=1:m
X(i,j)=data((i-1)*tau+j)
end
end
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)