juzhiquzao.m

数字图像处理

 load woman
init=2055615866
randn('seed',init);
x=X+18*randn(size(X));
[c,s]=wavedec2(x,3,'sym4');
n=[1,2,3];
p=[150,120,60];
nc=wthcoef2('d',c,s,n,p,'s');
rx=waverec2(nc,s,'sym4');
nbc=size(map,1);
colormap(pink(nbc));
subplot(2,2,1);
 image(wcodemat(X,nbc));
title('原始信号');
subplot(2,2,2);
image(wcodemat(x,nbc));
title('含噪信号');
subplot(2,2,3);
 image(wcodemat(xd,nbc));
 title('消噪后的信号');