load woman;
init=2055615866;
randn('seed',init);
x=X+20*randn(size(X));
[thr,sorh,keepapp]=ddencmp('den','wv',x);
xd=wdencmp('gbl',x,'sym4',2,thr,sorh,keepapp);
colormap(pink(255));
sm=size(map,1);
subplot(2,2,1),
image(wcodemat(X,sm));
subplot(2,2,2),
image(wcodemat(x,sm));
subplot(2,2,3),
image(wcodemat(xd,sm));
%%%=================================================
[c,s]=wavedec2(x,3,'sym4');
[thr3,nkeep3]=wdcbm2(c,s,1.5,2.7*prod(s(1,:)));
[xd3,cxd,sxd,perf0,perf12]=wdencmp('lvd',c,s,'sym4',3,thr3,'s');
subplot(2,2,4),image(wcodemat(xd3,sm));