samtest286.m

汽轮机组故障诊断中横向裂纹振动的数据

clear
fid=fopen('E:\汽轮机振动论文\date\date\故障数据\286-动静碰磨\rub2-dh.txt','r');
A=fscanf(fid,'%f ');
%for i=1024:8192
 %   A(i)=0;
    
 %end
 fclose(fid);
 fs=200;
 N=length(A);
 nf=0:N/2-1;
xf=nf*fs/N;
subplot(2,1,1) 
plot(A);
y=fft(A);
ylyp=abs(y)*2/1024 ;
subplot(2,1,2) 
plot(xf,ylyp(1:512));