转载来自：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('去噪后的图像');