hwtdn.m

这是我找到的谐波小波程序

function a=hwtdn(f)
%
% Copyrite(c) D.E.Newland 1992
% All rights reserved
%
N=length(f);
n=round(log(N)/log(2));
F=fft(f)/N;%为什么要除以N
a(1)=F(1);a(2)=F(2);
for j=1:n-2
    a(2^j+1:2^(j+1))=ifft(F(2^j+1:2^(j+1)))*2^j;
    a(N-2^(j+1)+2:N-2^j+1)=fliplr(fft(fliplr(F(N-2^(j+1)+2:N-2^j+1))));
end
a(N/2+1)=F(N/2+1);
a(N)=F(N);