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(