zao3

来自「数字图像处理」· 代码 · 共 23 行

load woman
init=2055615866
randn('seed',init);
x=X+18*randn(size(X));
nbc=size(map,1);
wname='sym4';lev=3;
[c,s]=wavedec2(x,lev,wname);
alpha=3;m=3.5*prod(s(1,:));
[thr,nkeep]=wdcbm2(c,s,alpha,m)
sorh='s';
xd=wdencmp('lvd',c,s,wname,lev,thr,sorh);
colormap(pink(nbc));
figure(1);
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('消噪后的信号');