Normalshrink阈值去噪matlab代码

%im 输入滤波对象
%levels为NSCT最大尺度

% function result=Normalshrink（imlevels）
% [mn]=size（im）;
% L_im=length（im）;
% belt=（log（L_im/levels））^0.5;
% delt=median（abs（im（:）））/0.6745;
% delt_temp=std（im（:））;
% t_temp=belt*delt/delt_temp;
% result=zeros（mn‘double‘）;
% for i=1:m
%     for j=1:n
%         ww=t_temp/（abs（im（ij））*exp（abs（im（ij））-t_temp））;
%         if abs（im（ij））%             result（ij）=0;
%         elseif im（ij）>t_temp
%             result（ij）=im（ij）-ww*t_temp;
%         else
%             result（ij）=im（ij）+ww*t_temp;
%         end
%     end
% end
function result=Normalshrink（imlevels）
[mn]=size（im）;
L_im=length（im（:））;
belt=（log（L