waveletdenoising 小波变换应用

小波变换应用，小波软阈值及硬阈值去噪算法及对此算法的改进算法。

%硬阈值去噪
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