%%%%%此文件用于進行相空間重構，并繪制相空間特征值及特征值斜率的變化趨勢

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all

%load mydata4xx1;
%B=xx1(3,:);
%B=0.2*xx1(1,:)+0.1*xx1(2,:)+xx1(3,:);
%load mydatayy;
%B=yy(3,:);

%load zh600_3_2_Waveform;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% load xie720-3Mpa-5_Waveform;
% x=a8(5001:11000);
% t=a0(5001:11000);

load bbb.m;

x=a8(5001:11000);
t=a0(5001:11000);
B=x-mean(x);
%load mydata2yy;
%B=YY;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%load zh480_2_4_Waveform;
%x=a6(1:5000);
%B=x-mean(x);


%$$load mydataxx1;
%$$B=0.5*xx1(1,:)+0.3*xx1(2,:)+xx1(3,:);
%$$t=0:0.002:9.9998;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%@@load mytime;
%@@load myxiedata41;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%@@B=xiedata41(4,:);
T=length(B);
tau=19;     %%% 設置時間延遲
n=60;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Y=reconstitution(B,T,n,tau); %%%% 重構吸引子矩陣，將樣本長度為T的信號B重構為n維，時間延遲為tau的吸引子矩陣
[U E V]=svd(Y);   %%%奇異值分解
%break;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:n
   e(i)=E(i,i);   %%%保存特征值
end
  
for i=1:(n-1)
     xielv(i)=e(i)-e(i+1);    %%%計算特征值比值
end

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   
  figure (2);
  subplot(212);
  plot(e(1:n));    %%%特征值變化趨勢
  
  subplot(211);
  plot(t,B);%%%特征值斜率變化趨勢
  
  
  break;
  hunhe=e;
  save mydatahunhesvd hunhe;