%采樣頻率
fs=5000;
nfft=10240;
%尖銳碰摩故障信號(hào)
%fid=fopen('rotor0.dat','r');%正常狀態(tài)
%fid=fopen('rotor1.dat','r');%早期尖銳觸碰
fid=fopen('d329_6_1.dat','r');%晚期平鈍型觸碰
N=2048;
x=fscanf(fid,'%f',N);
fclose(fid);
xdata=(x-mean(x))/std(x,1);

figure(1);
%時(shí)域圖
subplot(2,2,1);
plot(1:N,x);
xlabel('時(shí)間 t/s');
ylabel('幅值 A/W');
%功率譜
Y=abs(fft(xdata,nfft));
subplot(2,2,2);
plot((0:nfft/2-1)/nfft*fs,Y(1:nfft/2));
xlabel('頻率 f/Hz');
ylabel('功率譜 P/W');

%5層小波包分解
T=wpdec(xdata,5,'db2');
%重構(gòu)低頻信號(hào)
y1=wprcoef(T,[5,1]);
%y1的波形
figure(2);
subplot(2,2,1);
plot(1:N,y1);
xlabel('時(shí)間 t/n');
ylabel('電壓 V/v');
%y1的功率譜
Y1=abs(fft(y1,nfft));
subplot(2,2,2);
plot((0:nfft/2-1)/nfft*fs,Y1(1:nfft/2));
xlabel('頻率 f/Hz');
ylabel('功率譜 P/W');