signal_denoising3.m

MATLAB 小波分析高级技术——书中的一些应用程序

load leleccum;
indx = 1:1024; 
x = leleccum(indx);
% 产生含噪信号
init=2055615866; 
randn('seed',init); 
nx = x + 18*randn(size(x)); 
[c,l]=wavedec(nx,3,'db5');
%设置尺度向量
n=[1,2,3];
%设置阀值向量
p=[120,110,100];
%对高频系数进行阀值处理
nc=wthcoef('d',c,l,n,p);
%对修正后的小波分解结构进行重构
rx = waverec(nc,l,'db5');
subplot(2,2,1);
plot(x);
title('原始信号');
subplot(2,2,2);
plot(nx);
title('含噪信号');
subplot(2,2,3)
plot(rx);
title('消噪后的信号');