dyadpartition.sci

小波分解源代码

function dp = dyadpartition(n)

//  Copyright Aldo I Maalouf
J=ceil(log2(n));
m=n;
for j=J-1:-1:0
  if modulo(m,2) == 0 then
     dps(j+1)=m/2;
     m=m/2;
  else
    dps(j+1)=(m-1)/2;
    m=(m+1)/2;
  end
end
dp = cumsum([1;dps]); 
endfunction