画相空间轨迹图

clc
clear all
close all
load('NAF200.mat');
[h,l]=size(a)
% namelist = dir('C:\Users\王浩\Desktop\af200\*.mat');
% len = length(namelist);
for i=1:h
%      fileName = ['C:\Users\王浩\Desktop\400张测试集图片的原数据\NAF\' num2str(i) '.mat'];
%     L = load(fileName);
%     filename = namelist(i).name
%     L = load(filename);
%     L.b=L.b';
        X=a(i,:);
%     [p,q] = size(af_test);
%     for j = 1:p
%         X = L.X(j,:)
        
%         X=normalize(X)
        m = 3;
        maxLags =100
        tau = AutoCorrelation(X,maxLags);
        [xn,dn] = PhaSpaRecon(X,tau,m);
        x = xn(1,:);
        y = xn(2,:);
        z = xn(3,:);
        u = (x+y+z)/sqrt(3);
        v = (x+y-2*z)/sqrt(6);
        w = (x-y)/sqrt(2);
        figure('visible','off');
        plot3(u,w,v);
        view(90,0);
        set(gca,'XColor','white') 
        set(gca,'YColor','white')
        set(gca,'ZColor','white')
        f = getframe(gcf)
        imwrite(f.cdata,['C:\Users\王浩\Desktop\NAF\',int2str(i),'in',int2str(i),'.png']);
%         clear r;
end