转载来自:http://blog.csdn.net/bluecol/article/details/46390105
%装载原始图像 img=imread('ddd.png'); x = img ; %x包含原始图像 %产生噪声图像 init=2055615866; rand('seed',init); noise = 12*randn(size(x)); x=uint8(double(x)+noise); subplot(1,3,1),image(mat2gray(x)); title('含噪图像'); %使用wdencmp进行图像去噪 %得到x的小波分解结构 n=2; w='db5'; [c,l]=wavedec2(x,n,w); %小波系数的阈值处理 [thr,sorh,keepapp]=ddencmp('den','wv',x); [xd,cxd,lxd,perf0,perfl2]=wdencmp('gbl',c,l,w,thr,'s',l); %画出去噪后的图象 subplot(1,2),image(mat2gray(xd)); title('去噪后的图像');