f14_3.m

《MATLAB时频分析技术及其应用》源程序代码

%采样频率
fs=10000;
%齿轮点蚀声音信号
fid=fopen('GEAR2.dat','r');%故障
N=1024;
xdata=fread(fid,N,'int16');
fclose(fid);
xdata=(xdata-mean(xdata))/std(xdata,1);
%时域波形
figure(1);
plot(1:N,xdata);
xlabel('时间 t/n');
ylabel('电压 V/v');

%db5小波进行4层分解
%一维小波分解
[c,l] = wavedec(xdata,4,'db5');
a4 = wrcoef('a',c,l,'db5',4);
%重构第1-4层细节系数
d4 = wrcoef('d',c,l,'db5',4);
d3 = wrcoef('d',c,l,'db5',3);
d2 = wrcoef('d',c,l,'db5',2);
d1 = wrcoef('d',c,l,'db5',1);

%显示细节号
figure(2)
subplot(5,1,1);
plot(d1,'LineWidth',2);
ylabel('d1');
subplot(5,1,2);
plot(d2,'LineWidth',2);
ylabel('d2');
subplot(5,1,3);
plot(d3,'LineWidth',2);
ylabel('d3');
subplot(5,1,4);
plot(d4,'LineWidth',2);
ylabel('d4');
xlabel('时间 t/s');
subplot(5,1,5);
plot(a4,'LineWidth',2);
ylabel('a4');
xlabel('时间 t/s');