1、调用wden函数进行一维小波去噪，其调用格式如下： [XD,CXD,雉搽妤粲LXD]=wden(x,tptr,sorh,scal,N,'wname&#泌驾台佐39;) 其中x表示原始信号，tptr表示阀值类型，sorh表示软阀值(s表示)或硬阀值(h表示)，N表示分解层数，wname小波函数类型，XD去除噪声后信号，CXD表示各层分量，LXD表示各层分量对应的长度。示例1对应调用wden函数对应的m文件如下： t=0:0.001:1;%%时间序列%% x=sin(5*pi*t);%%原始信号%% n=randn(size(t));%%表示产生一个噪声信号%% y=x+n;%%含有噪声的信号%% yd=wden(y,'heursure','s','one',3,'sym5'); subplot(3,1,1) %%表示显示三行一列第一个图%%%%设置第一幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(x,'r'); xlabel('n'); ylabel('幅值')； title('原始信号')； axis([0 1000 -5 5]) subplot(3,1,2) %%表示显示三行一列第二个图%%%%设置第二幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(y,'k'); xlabel('n'); ylabel('幅值')； title('含有噪声信号')； axis([0 1000 -5 5]); subplot(3,1,3) %%表示显示三行一列第三个图%%%%设置第三幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(yd,'c'); xlabel('n'); ylabel('幅值')； title('消除噪声后信号')； axis([0 1000 -5 5]);
2、运行示例1对应的m文件程序，结果如图：
3、调用ddencmp函数，wdcbm函数分别得到对应的全局阈值，分层阈值，调用格式分别如下： [thr,sorh,keepapp]=ddencmp('den’,’wv’,x) [thr1,nkeep]=wdcbm(c,l,N);其中den表示去噪，wv小波，x原始输入信号，keepapp=1，表示保留低频系数不变，thr表示全局阈值，sorh表示阈值类型，thr1表示分层阈值
4、调用wdencmp函数对一维或二维信号进行全局阀值消噪，调用格式如下： [XD,CXD,LXD]=wdencmp('gbl’,x,'wname',N,thtr,sorh,keepapp)示例2对应调用wdencmp函数进行全局阈值去噪对应的m文件如下： t=0:0.001:1; x=sin(5*pi*t);%%原始信号%% n=randn(size(t));%%表示产生一个噪声信号%% y=x+n;%%获得全局阈值程序%% [thr,sorh,keepapp]=ddencmp('den','wp',y);%%利用全局阈值进行去噪%% yd=wdencmp(‘gbl',y,'sym4',4,thr,sorh,keepapp); subplot(3,1,1) %%表示显示三行一列第一个图%%%%设置第一幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(x,'r'); xlabel('n'); ylabel('幅值')； title('原始信号')； axis([0 1000 -5 5]) subplot(3,1,2) %%表示显示三行一列第二个图%%%%设置第二幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(y,'k'); xlabel('n'); ylabel('幅值')； title('含有噪声信号')； axis([0 1000 -5 5]); subplot(3,1,3) %%表示显示三行一列第三个图%%%%设置第三幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(yd,'c'); xlabel('n'); ylabel('幅值')； title('利用全局阈值消除噪声后信号')； axis([0 1000 -5 5]);
5、运行示例2对应的m文件，结果如图：
6、调用wdencmp函数对一维或二维信号进行分层阀值消噪，调用格式如下 [XD,CXD,LXD]=wdencmp('lvd',c,l,'wname',N,thtr,sorh,keepapp)示例3对应调用wdencmp函数进行分层阈值去噪对应的m文件如下：t=0:0.001:1; x=sin(5*pi*t);%%原始信号%% n=randn(size(t));%%表示产生一个噪声信号%% y=x+n;%%获得分层阈值%% [c,l]=wavedec(y,4,'sym4'); [thr1,nkeep]=wdcbm(c,l,4)%%利用分层阈值进行去噪%% yd=wdencmp('lvd',c,l,'sym4',4,thr1,'s'); subplot(3,1,1) %%表示显示三行一列第一个图%%%%设置第一幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(x,'r'); xlabel('n'); ylabel('幅值')； title('原始信号')； axis([0 1000 -5 5]) subplot(3,1,2) %%表示显示三行一列第二个图%%%%设置第二幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(y,'k'); xlabel('n'); ylabel('幅值')； title('含有噪声信号')； axis([0 1000 -5 5]); subplot(3,1,3) %%表示显示三行一列第三个图%%%%设置第三幅图像的x轴、y轴、标题信息以及坐标轴范围信息%% plot(yd,'c'); xlabel('n'); ylabel('幅值')； title('利用分层阈值消除噪声后信号')； axis([0 1000 -5 5]);
7、运行示例3对应m文件程度，结果如图：
