iharw.m

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

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