小波的应用

clear all;
close all;
I=imread('E:\锦州世园会\psb8.jpg');
X=rgb2gray(I);
[height width]=size(X);
X=im2double(X);
A=imnoise(X,'gaussian',0,0.01); 
wname='sym4';
n=3;
[c,s]=wavedec2(A,n,wname);
for i=1:3
step(i)=s((i+1),1)*s((i+1),2);
end
num(1,1)=s(1,1)*s(1,2)+1; %获取各层各高频分量在c向量中的坐标
num(1,2)=num(1,1)+s(2,1)*s(2,2);
num(1,3)=num(1,2)+s(2,1)*s(2,2);
num(2,1)=num(1,3)+s(2,1)*s(2,2);
num(2,2)=num(2,1)+s(3,1)*s(3,2);
num(2,3)=num(2,2)+s(3,1)*s(3,2);
num(3,1)=num(2,3)+s(3,1)*s(3,2);
num(3,2)=num(3,1)+s(4,1)*s(4,2);
num(3,3)=num(3,2)+s(4,1)*s(4,2);
m=0.02;
C=c;
for i=1:3
[H,V,D]=detcoef2('a',c,s,i);%提取第i层各高频系数
B=[H V D];
[L,T]=size(B);
for k=1:L
for w=1:T