load noisbloc
x = noisbloc;
t = wpdec(x,3,'sym4');
plot(t);
%在運行上面命令出現(xiàn)的圖形窗口中將Node Label從Depth-postion切換到
%Index并單擊索引值為0的結(jié)點，即根結(jié)點，如圖4-9所示
%全局閾值處理
t1 = t;
sorh = 'h'; 
%硬閾值
thr = wthrmngr('wp1ddenoGBL','penalhi',t); 
%設(shè)置消噪閾值
cfs = read(t,'data');
cfs = wthresh(cfs,sorh,thr);
t1 = write(t1,'data',cfs);
figure(2);
plot(t1)
%在運行上面命令后出現(xiàn)的圖形窗口中再將Node Label從Depth-postion切換到
%Index并單擊索引值為0的結(jié)點，即根結(jié)點，如圖4-10所示
%各個結(jié)點分別進行閾值處理
t2 = t;
sorh = 's';
thr(1) = wthrmngr('wp1ddenoGBL','penalhi',t);
thr(2) = wthrmngr('wp1ddenoGBL','sqtwologswn',t);
tn = leaves(t);
for k=1:length(tn)
node = tn(k);
cfs = read(t,'data',node);
numthr = rem(node,2)+1;
cfs = wthresh(cfs,sorh,thr(numthr));
t2 = write(t2,'data',node,cfs);
end
figure(3);
plot(t2)
%在運行上面命令后出現(xiàn)的圖形窗口中再將Node Label從Depth-postion切換到Index并單擊
%索引值為0的結(jié)點，即根結(jié)點，如圖4-11所示