harw.m

来自「以上是我大学毕业设计」· M 代码 · 共 13 行

function y=harw(x,N)
Fs=fft(x,N);
for j=-1:(log(N/2)/log(2))
    if(j==-1)a(1)=ifft(Fs(1:1),1)
    else 
        a((2^j)+1:(2^(j+1)))=ifft(Fs((2^j)+1:(2^(j+1))),(2^(j+1))-(2^j));
    end
end
%a(N/2+1)=ifft(Fs(N/2+1),1);
for j=(N/2+2):N
   % a(j)=conj(a(N+2-j));
end
y=a;