基于小波变换的图像去噪算法及代码

毕业设计的题目是小波去噪，参考了一些网上的资料，自己加以补充和修改，希望有用。

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

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

     文件        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

     文件       4356  2007-07-09 15:58  毕设\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

     文件        101  2007-05-23 15:27  毕设\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

     目录          0  2009-09-15 16:45  毕设

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

               653225                    23