函数 Y=zjfft(x) %????FFT?????????? 抽动; x1=换单（x）； %????????? N=长度（x1）； %？？？？？？ M=log2(N); %？？？？？？ Y=零（1，N）； %???????N??? B=零（1，N）； %------------????????????-------------- Y=x1; %-------------------??----------------------
对于 i=1:M B=Y; 对于 r=1:N %B=Y; zf=((-1).^floor((r-1)/(2.^(i-1)))); %？？？ W=exp(-j*2*pi*mod(r-1,2.^(i-1))/2.^i); %？？？？ Y(r)=zf*B(r)*W.^((zf-1)/(-2))+W.^((zf+1)/2)*B(r+zf*2.^ (i-1)); %
