hareresholddenoising.m

对信号进行小波分解

load S20F20.DAT
x=S20F20;
y=x(:,2);

x=[];
k=[];
N=length(y);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  4.硬消噪采用cycle_spinning技术
%  累加量
z5=zeros(N,1);
%  平移变换频移法
for i=1:N;
    z=circshift(y.',i-1).';  %  源信号右平移
    [z1,z2]=lwt(z,'db4');    %  小波正变换
    z2=zeros(1,N/2);         %  高频分量全部为零（主要噪声，硬消噪）
    z3=ilwt(z1,z2,'db4');    %  小波反变换  
    z4=circshift(z3.',-(i-1)).';  %  变换后信号左平移  
    z5=z5+z4/N;
    %z5=z4;  %  平均  ?????????
end;
    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  5.显示
error=norm(y-z5)/norm(y);  % 相对误差
figure(1);
subplot(3,1,1)
plot(y);
legend('original signal');
subplot(3,1,2);
plot(z5);
legend('denoised signal');
subplot(3,1,3);
%plot(1:6001,y,1:6001,z5,'r:');
compare=y-z5;
plot(compare);
legend('compare');
%plot(compare);
%legend('copmpare');