波变换

开发平台：

Matlab

den2.m：源码内容

							function X = den2(x, wname, n, thr, alpha)
% 对软硬阈值函数的改进方法
% w_jk < λ时 直接置零
% w_jk >= λ时 w_jk' = sign(w_jk) * sqrt(|w_jk|^2 - αλ^2)
[C, S] = wavedec2(x, n, wname);                     %对图像进行小波分解
dcoef = C( prod(S(1, :)) + 1 : end);                %高频部分系数
ind = find( abs(dcoef) < thr) + prod(S(1, :));      %小于阈值thr的系数
C(ind) = 0;                                         %   直接置零
ind = find( abs(dcoef) >= thr) + prod(S(1, :));     %大于阈值thr的系数
C(ind) = sign(C(ind)) .* sqrt( C(ind).^2 - alpha * thr^2 );
                                                    %按照公式处理
                                                    
X = waverec2(C, S, wname);                          %重构图像

						