【资料图】

1、function xn=myfft(x)N=length(x);M=log2(N);xtmp=zeros(1,N);value=zeros(1,M);for i=0:N-1 repr=i; for t=1:1:M repr=bitshift(i,1-t); value(t)=bitand(repr,1); end pos=0; for k=1:1:M pos=pos+value(k)*2^(M-k); end xtmp(pos+1)=x(i+1);endfor i=1:M deepth=2^(i-1); width=2^(M-i); for t=1:2^i:N for k=1:deepth tmp=xtmp(t+k-1); wn=width*(k-1); xtmp(t+k-1)=tmp+exp(-j*2*pi*wn/N)*xtmp(t+k+deepth-1); xtmp(t+k+deepth-1)=tmp-exp(-j*2*pi*wn/N)*xtmp(t+k+deepth-1); end endendxn=xtmp;。

本文就为大家分享到这里，希望小伙伴们会喜欢。