%?裝載原始圖像
I=imread（‘D:\MATLAB\toolbox\images\imdata\cameraman.tif‘）;
nbc=size（I1）;
X?=?im2double（I）;
%?產(chǎn)生噪聲圖像
init=2055415866;
randn（‘seed‘init）;
x=X+randn（size（X））/10;
%?使用?sym4?執(zhí)行圖像的2層小波分解
wname=‘sym4‘;lev=2;
[cl]=wavedec2（xlevwname）;
sigma_s=0.054779;
%?圖像降噪時(shí)，使用wbmpen?函數(shù)選擇閾值主要有ddencmp，thselect，wbmpen，wdcbm四種
alpha=2;
thr_s=wbmpen（clsigma_salpha）;
%?使用軟閾值和保存的低頻信號(hào)，進(jìn)行圖像降噪。閾值去噪函數(shù)主要有wden，wdencmp，wpdencmp，whresh，wpthcoef，wthcoef2六種。
keepapp=1;
xds=wdencmp（‘gbl‘xwnamelevthr_s‘s‘keepapp）;
sigma_h=0.062818;
thr_h=wbmpen（clsigma_halpha）;
%?使用硬閾值和保存的低頻信號(hào)，進(jìn)行圖像降噪
xdh=wdencmp（‘gbl‘xwnamelevthr_h‘h‘keepapp）;
%?畫出原始圖像和降噪后的圖像
figure（1）;
subplot（221）;imshow（I）;title（‘原始圖像‘）;
subplot（222）;imshow（x）;title（‘噪聲圖像‘）;
subplot（223）;imshow（xds）;title（‘軟閾值降噪圖像‘）;
subplot（224）;imshow（xdh）;title（‘硬閾值降噪圖像‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????983??2014-04-06?11:33??jiangzao.m

-----------?---------??----------?-----??----

??????????????????983????????????????????1