%硬閾值去噪
clear?all;
close?all;
X=imread（‘woman.tif‘）;
[height?width]=size（X）;
X=im2double（X）;
A=imnoise（X‘gaussian‘00.01）;?
wname=‘sym4‘;
n=3;
[cs]=wavedec2（Anwname）;

for?i=1:3
step（i）=s（（i+1）1）*s（（i+1）2）;
end
num（11）=s（11）*s（12）+1;?%獲取各層各高頻分量在c向量中的坐標
num（12）=num（11）+s（21）*s（22）;
num（13）=num（12）+s（21）*s（22）;

num（21）=num（13）+s（21）*s（22）;
num（22）=num（21）+s（31）*s（32）;
num（23）=num（22）+s（31）*s（32）;

num（31）=num（23）+s（31）*s（32）;
num（32）=num（31）+s（41）*s（42）;
num（33）=num（32）+s（41）*s（42）;

m=0.02;
C=c;
for?i=1:3
[HVD]=detcoef2（‘a‘csi）;%提取第i層各高頻系數
B=[H?V?D];
[LT]=size（B）;
for?k=1:L
for?w=1:T
sigma=median（abs（B（kw）））/0.6745;
end
end
th=sigma*sqrt（2*（log（L*T）/log（2.71828）））;%計算閾值
ch=c（1num（4-i1）:num（4-i3）+step（4

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1401??2009-11-23?19:10??小波閾值去噪\ha.m

?????文件???????1425??2009-11-23?19:09??小波閾值去噪\so.m

?????文件???????1546??2009-11-23?19:09??小波閾值去噪\soft_ha.m

?????目錄??????????0??2010-11-25?20:29??小波閾值去噪

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

?????????????????4372????????????????????4