matlab小波去噪

matlab的小波去噪，里面附有多种小波去噪方法，有注释，程序

function thr = Birge_Massart（c s）
%用Birge-Massart策略估计阈值
% 给定分解层数j，对j+1及更高层，所有系数保留
% 对第i层（1<=i<=j），保留绝对值最大的n_i个系数，n_i由下式确定
%   n_i = M/（j+2-i）^alpha
% 其中M和alpha为经验系数，缺省情况下M取第一层分解后的系数长度L（1）
% 一般情况下，M满足 L（1）<=M<=2L（1）
% alpha的取值因用途不同，压缩情况下一般取alpha=1.5，降噪一般取alpha=3

%设定alpha、M和j的值
alpha = 3;
M = prod（s（1 :））;
j = size（s 1） - 2;

%计算每层的阈值
for i = 1:j
    n = round（ M/（j+2-i）^alpha ）;                   %每层保留的系数个数
    D = detcoef2（‘compact‘ c s i）;               %提取第i层的所有细节系数 
    if n >= size（D）
        thr（: i） = zeros（3 1）;
    else
        D = sort（abs（D））;                           %排序
        thr（i） = D（end-n）;                          %设置第i层的阈值
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-06-11 12:23  小波去噪\
     文件         844  2007-07-09 15:05  小波去噪\Birge_Massart.m
     文件         767  2007-07-09 15:55  小波去噪\den1.m
     文件        1113  2007-07-09 15:55  小波去噪\den1_10.m
     文件        1600  2007-07-09 15:54  小波去噪\den1_5_1.m
     文件         796  2007-07-09 15:55  小波去噪\den1_9.m
     文件         756  2007-07-09 15:55  小波去噪\den2.m
     文件         856  2007-06-07 08:28  小波去噪\den3.m
     文件        2808  2007-07-09 15:57  小波去噪\den4.m
     文件        4380  2018-06-11 12:31  小波去噪\denoise0615.m
     文件        5379  2007-07-09 15:22  小波去噪\denoise0701.m
     文件         492  2007-06-28 19:33  小波去噪\Donoho.m
     文件      150923  2000-10-07 15:18  小波去噪\lena.png
     文件      476195  2000-10-07 15:05  小波去噪\lena_color.png
     文件         492  2007-07-09 15:18  小波去噪\meanfilter.m
     文件         581  2007-07-09 15:18  小波去噪\medfilter.m
     文件         189  2007-07-02 09:22  小波去噪\MSE.m
     文件         281  2007-07-09 15:52  小波去噪\MSE_color.m
     文件        1051  2007-07-01 09:57  小波去噪\noise.m
     文件          97  2018-06-11 12:30  小波去噪\PSNR.m
     文件         113  2007-07-02 09:29  小波去噪\PSNR_color.m
     文件        1754  2007-06-28 21:27  小波去噪\wdenoise.m
     文件        1778  2007-07-09 16:03  小波去噪\文件清单及说明.txt