mysvd00.asv

matlab 时频分析上面程序第八章的故障诊断的小程序

%%%%%此文件用于进行相空间重构，并绘制相空间特征值及特征值斜率的变化趋势

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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;