xiaobods.m

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

load('E:/高珍/Test410-200.mat');
t1=double(Data1_Time);
t=(0:length(t1)-1)./1000;
y1=double(Data1_Channel_0____);
plot(t,y1);
Xlabel('时间, t/s');
Ylabel('电压, V');
[c,l]=wavedec(y1,3,'db5');
ca3=appcoef(c,l,'db5',3);
cd3=detcoef(c,l,3);
cd2=detcoef(c,l,2);
cd1=detcoef(c,l,1);
figure(2);
subplot(221);plot(ca3);
Ylabel('ca3');
subplot(222);plot(cd3);
Ylabel('cd3');
subplot(223);plot(cd2);
Ylabel('cd2');
subplot(224);plot(cd1);
Ylabel('cd1');
%对低频小波系数进行重构
cd3l=length(cd3);
d3=zeros(1,cd3l);
cd1l=length(cd1);
d1=zeros(1,cd1l);
cd2l=length(cd2);
d2=zeros(1,cd2l);
cL=[ca3;d3';d2';d1'];
y1dbL=waverec(cL,l,'db5');
%对高频小波系数进行重构
ca3l=length(ca3);
a3=zeros(1,ca3l);
cd1l=length(cd1);
d1=zeros(1,cd1l);
cd2l=length(cd2);
d2=zeros(1,cd2l);
cH=[a3';cd3;d2';d1'];
y1dbH=waverec(cH,l,'db5');
%对重构系数做频谱分析
Fs=1000;
w1=hanning(length(y1dbL));
[Pxx1,f1]=periodogram(y1dbL,w1,51200,Fs);
figure(3);
plot(f1,Pxx1);
Xlabel('f/Hz');
Ylabel('幅值量N');
w2=hanning(length(y1dbH));
[Pxx2,f2]=periodogram(y1dbH,w2,51200,Fs);
figure(4);
plot(f2,10*log10(Pxx2));
Xlabel('f/Hz');
Ylabel('Power Spectral Density(dB/Hz)');