mm_dwt.sci

小波分解源代码

function dwtm = MM_DWT(dwt,threshold);
//  Copyright Aldo I Maalouf
 [n,nscales] = size(dwt);
 D = nscales-1;
 J0 = D-1;
 m = dwt(:,2:nscales); 
 tmp = ModulusMaxima(m,1,1,threshold);
 dwtm(1)=0;
 dwtm(2:nscales)=tmp;
 dwtm(:,1) 	   = ones(n,1);
 endfunction