- load flujet;
- subplot(2,2,1);
- image(X);
- colormap(map);
- axis square;
- init=2055615899;
- randn('seed',init);
- X1=X+20*randn(size(X));
- subplot(2,2,2);
- image(X1);
- colormap(map);
- T=wpdec2(X1,1,'sym2');
- thr=8.342;
- NT=wpthcoef(T,0,'s',thr);
- X2=wprcoef(NT,1);
- subplot(2,2,3);
- image(X2);
- colormap(map);