例程18-7.m

《小波分析理论与Matlab7实现例程》一书里面各章节的源码

%实验数据
%fid=fopen('rotor00.dat','r');%正常状态
fid=fopen('rotor01.dat','r');%0.35mm
%fid=fopen('rotor02.dat','r');%0.77mm
%fid=fopen('rotor03.dat','r');%1.05mm
%fid=fopen('rotor04.dat','r');%1.35mm
N=2048;
fs=5000;
nfft=10240;
x=fscanf(fid,'%f',N);
fclose(fid);
xdata=(x-mean(x))/std(x,1);

figure(1);
%时域图
subplot(2,2,1);
plot((1:N)/fs,x);
xlabel('时间 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');
%重构低频信号
y1=wprcoef(T,[5,1]);
%y1的波形
figure(2);
subplot(2,2,1);
plot(1:N,y1);
xlabel('时间 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');