X=imread（‘lena.bmp‘）;
subplot（231）;
imshow（X）;
title（‘原始圖像‘）;
Xn=imnoise（X‘gaussian‘00.01）;
subplot（232）;imshow（Xn）;
title（‘含噪圖像‘）;
[thrsorhkeepapp]=ddencmp（‘den‘‘wv‘double（Xn））;
Xd?=?wdencmp（‘gbl‘double（Xn）‘coif5‘2thrsorhkeepapp）;
subplot（234）;imshow（uint8（Xd））;
title（‘默認閾值去噪圖像‘）;

[cl]=wavedec2（double（Xn）2‘coif5‘）;
alpha=1.2;
det1=detcoef2（‘compact‘cl1）;
sigma=median（abs（det1））/0.6745;
thr=wbmpen（clsigmaalpha）;
keepapp=1;
Xb=wdencmp（‘gbl‘cl‘coif5‘2thr‘s‘keepapp）;
subplot（235）;imshow（uint8（Xb））;
title（‘全局閾值去噪圖像‘）;

Xw=wiener2（Xn[5?5]）;
subplot（236）;imshow（