program_18_05_special.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

%采样频率
fs=10000;
nfft=10240;
%轴承滚动体故障信号
fid=fopen('bearingroll2.dat','r');%故障
N=2048;
xdata=fread(fid,N,'int16');
fclose(fid);
xdata=(xdata-mean(xdata))/std(xdata,1);
%功率谱
figure(1);
Y=abs(fft(xdata,nfft));
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');