基于小波阈值去噪

基于小波去噪，里面含有硬阈值去噪，软阈值去噪，半软阈值去噪等，已编译成功 MATLAB程序

close all;
clc;clear all;
% 产生中频信号
t=0:0.1:5*pi;
xn=5*cos（t）+0.95*randn（1length（t））; 
N=length（xn）;
figure
plot（xn）
xlabel（‘采样点数‘）;ylabel（‘幅值/m‘）;
% title（‘含噪正弦信号‘）;
axis（[0N-1010]）;
lev=5;
yn=wden（xn‘rigrsure‘‘h‘‘one‘lev‘db5‘）; % 硬阈值 [XDCXDLXD]=wden（XTPTRSORHSCALN‘wname‘）N为小波分解的层数
%函数wden用于一维信号的自动消噪。X为原始信号，[CL]为信号的小波分解，N为小波分解的层数。*TPTR=‘rigrsure‘，自适应阈值选择使用Stein的无偏风险估计原理。
%SORH是软阈值或硬阈值的选择（分别对应‘s‘和‘h‘）。XD为消噪后的信号，[CXDLXD]为消噪后信号的小波分解结构。SCAL指所使用的阈值是否需要重新调整，*SCAL=‘one‘??不调整；
figure
plot（1:Nyn）
xlabel（‘采样点数‘）;ylabel（‘幅值/m‘）;
% title（‘a.硬阈值处理‘）;
axis（[0N-1010]）;

yn2=wden（xn‘rigrsure‘‘s‘‘one‘lev‘db5‘）; % 软阈值
figure
plot（1:Nyn2）
xlabel（‘采样点数‘）;ylabel（‘幅值/m‘）;
% title（‘b.软阈值处理‘）;
axis（[0N-1010]）;

wname=‘sym5‘;
[cl]=wavedec（xn5wname）;%使用小波基函数 ‘wname‘ 对一维信号 X 进行 N 层分解。N必须是正整数
sigma=wnoisest（cl1）;%估计一维小波的系数的标准偏差 %
Thr=sigma*sqrt（2*log（N））;

car5=wrcoef（‘a‘clwname5）;%对一维信号的分解结构[cl]用指定的小波函数进行重构，当‘type=a‘ 时对信号的低频部分进行重构，此时N可以为0.
yn=car5;
yn2=car5;
yn3=car5;
yn_thr=car5;

for j=1:5
    cdr=wrcoef（‘d‘clwnamej）;   %当‘type=d‘时，对信号的高频部分进行重构，此时N为正整数。   
    cdr3=yasser（cdrThr）;  % Yasser降噪
    yn3=yn3+cdr3;
end
figure
plot（1:Nyn3）
xlabel（‘采样点数‘）;ylabel（‘幅值/m‘）;
% title（‘c.Yasser阈值处理‘）;
axis（[0N-1010]）;

for j=1:5
    cdr=wrcoef（‘d‘clwnamej）;      
    cdr_thr=wthr（cdrThr）;  % 改进阈值降噪
    yn_thr=yn_thr+cdr_thr;
end
figure
plot（1:Nyn_thr）
% title（‘d.改进阈值处理‘）;
xlabel（‘采样点数‘）;ylabel（‘幅值/m‘）;
axis（[0N-1010]）;

% 参数计算
yy = [yn;yn2;yn3;yn_thr];
for j=1:4
    sum=0;sum1=0;sum2=0;
    for k=1:N
        sum=sum+xn（k）^2;
        sum1=sum1+（xn（k）-yy（jk））^2;
        sum2=sum2+yy（jk）^2;
    end
    rm=sqrt（sum1/N）;
    snr=10*log10（sum/sum1）;
    per=sqrt（sum2/N）/sqrt（sum/N）;
    fprintf（‘程序结束:\n RMSE is %fSNR is %fPER is %f.\n‘rmsnrper）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2195  2016-03-07 16:30  小波去噪\denoise.asv

     文件       2195  2016-03-07 16:30  小波去噪\denoise.m

     文件        262  2015-12-22 00:48  小波去噪\hardthr.m

     文件        283  2015-12-22 00:48  小波去噪\softhr.m

     文件        477  2015-12-22 00:48  小波去噪\wthr.m

     文件        305  2015-12-22 00:48  小波去噪\yasser.m

     目录          0  2016-03-14 15:23  小波去噪

----------- ---------  ---------- -----  ----

                 5717                    7