数字信号处理自己编写的FFT程序
function[FX]=myfft(x)
N0=length(x);
ifN0==1
FX=x;
elseN0>1
N=2^ceil((log(N0))/(log(2)));
fori=0:(N/2)-1
x1(i+1)=x(2*(i+1));
x2(i+1)=x(2*(i+1)-1);
end
FX1=myfft(x1);
FX2=myfft(x2);
WN=exp(-j*2*pi/N);
fori=0:(N/2)-1
FX(i+1)=FX1(
下载地址
用户评论
