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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all;
clc;
%load myshuju02
%B=yy(2,:);
%load data1;
%B=a1(1:4238);

%load mydata4xx1;
%B=0.2*xx1(1,:)+0.1*xx1(2,:)+xx1(3,:);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%load zh600_3_2_Waveform;
load xie720-3Mpa-5_Waveform;


B=a8(5001:11000);
%load mydata2yy;
%B=YY;

t=a0(5001:11000);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%load mydataxx1;
%B=0.5*xx1(1,:)+0.3*xx1(2,:)+xx1(3,:);
%t=0:0.002:9.9998;
%T=length(B);

%break
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%$$load mytime;
%$$load myxiedata41;
%$$B=xiedata41(4,:);
T=length(B);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% 設置基本參數
T0=1;    %%%%設定重構相空間矩陣保留特診征值維數
T1=30;
n=60;       %%% 重構相空間的的行數
tau=19;     %%% 設置時間延遲

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%重構相空間
ee=e;
for i=1:n
    if i>=T0 & i<=T1
       ee(i)=e(i);           
    else
       ee(i)=0;
    end
end
%%%%%%%%%%%%%%%%%%%%%
E=diag(ee(1:n),0);
m=length(Y(1,:));%%%% 重構時間序列的長度
EE=zeros(n,m);
for i=1:m   %%%恢復特征值矩陣   
    if i<=n
    EE(i,i)=E(i,i);
    end
end
    
 AA=U*EE*V';   %%% 利用SVD反變換進行相空間重構
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 YY=myresignal(AA,tau);  %%%對重構的矩陣通過平均對應元素恢復提取的特征信號
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %fs=5681.32;   %%xie chi
 %N=5000;          %采樣點數
 %x=YY(1:5000);
 %x=x-mean(x);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%xf=fft(x,N)/N;   % 進行傅立葉變換
%xf=fftshift(xf);   %雙邊復數譜
%df=fs/N;          %頻率分辨率Hz
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%繪制雙邊幅值譜
%f=(-N/2:N/2-1)*df;    %頻域序列

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   
  figure (7);
  %subplot(221);
  %plot(e(1:n));    %%%特征值變化趨勢
  %plot(H(1:n));    %%%特征譜變化趨勢
  
  %subplot(222);
  %plot(xielv(1:n-1));%%%特征值斜率變化趨勢
  
  subplot(211);
  plot(t,B(1:T));%%%繪制特征值分解前的原始信號
  
  subplot(212);
  plot(t,YY(1:T));%%%繪制特征值分解后的重構信號
  
  save mydatashice-a8 YY;

    break;
  subplot(313);
  plot(f,abs(xf));
  xlabel('f/Hz');
  
  save mydatayy10-10 YY;