xiaobolw.m

本源代码主要实现谐波小波的分解及重构的matlab程序。

load('E:/高珍/Test401-200.mat');
t1=double(Data1_Time);
t=(0:length(t1)-1)./5000;
y1=double(Data1_Channel_0____);
plot(t,y1);
Xlabel('时间, t/s');
Ylabel('电压, V');
[c,l]=wavedec(y1,5,'db5');
ca5=appcoef(c,l,'db5',5);
cd5=detcoef(c,l,5);
cd4=detcoef(c,l,4);
cd3=detcoef(c,l,3);
cd2=detcoef(c,l,2);
cd1=detcoef(c,l,1);
figure(2);
subplot(321);plot(ca5);
Ylabel('ca5');
subplot(322);plot(cd5);
Ylabel('cd5');
subplot(323);plot(cd4);
Ylabel('cd4');
subplot(324);plot(cd3);
Ylabel('cd3');
subplot(325);plot(cd2);
Ylabel('cd2');
subplot(326);plot(cd1);
Ylabel('cd1');
%对低频小波系数进行重构
cd5l=length(cd5);
d5=zeros(1,cd5l);
cd4l=length(cd4);
d4=zeros(1,cd4l);
cd3l=length(cd3);
d3=zeros(1,cd3l);
cd2l=length(cd2);
d2=zeros(1,cd2l);
cd1l=length(cd1);
d1=zeros(1,cd1l);
cL=[ca5;d5';d4';d3';d2';d1'];
y1dbL=waverec(cL,l,'db5');
%对高频小波系数进行重构
ca5l=length(ca5);
a5=zeros(1,ca5l);
cd1l=length(cd1);
d1=zeros(1,cd1l);
cd2l=length(cd2);
d2=zeros(1,cd2l);
cd3l=length(cd3);
d3=zeros(1,cd3l);
cd4l=length(cd4);
d4=zeros(1,cd4l);
cH=[a5';cd5;d4';d3';d2';d1'];
y1dbH=waverec(cH,l,'db5');
%对重构系数做频谱分析
Fs=5000;
w1=hanning(length(y1dbL));
[Pxx1,f1]=periodogram(y1dbL,w1,120000,Fs);
figure(3);
plot(f1,Pxx1);
Xlabel('f/Hz');
Ylabel('幅值量N');
w2=hanning(length(y1dbH));
[Pxx2,f2]=periodogram(y1dbH,w2,120000,Fs);
figure(4);
plot(f2,10*log10(Pxx2));
Xlabel('f/Hz');
Ylabel('Power Spectral Density(dB/Hz)');