function y=FFT1(x,q)
n=0;
t=q;
%解决补零问题
if (length(x)<q)
q=length(x);
end
while t>1 %计算几点FFT？即2的多少次方
t=t/2;
n=n+1;
end
if(length(x)<2^n)
bak=zeros(1,2^n);
for p=1:q
bak(p)=x(p);
end
x=bak;
end
q=2^n;
%计算FFT
for m=1:n %n是2的n次方等于m m级 最外层循环
for h=1:2^(m-1) %每一级由几个部分组成？
for k1=((h-1)*q/(2^(m-1))+1):((h-1)*q/(2^(m-1))+q/(2^m)) %每个部分的上半部不含虚数的
y(k1)=x(k1)+x(2^(n-m)+k1); %计算式
end
mark=1;
for k2=((h-1)*(q/(2^(m-1)))+q/(2^m)+1):((h-1)*(q/(2^(m-1)))+q/(2^(m-1)))%下半部分含虚数的那部分
w=0:2^(m-1):(q/2-2^(m-1)); %每一级的w参数 其分布规律
y(k2)=(x(k2-2^(n-m))-x(k2))*exp(-j*(2*pi/q)*w(mark)); %计算式
mark=mark+1;
end
x=y;
end
end
%到位序列运算
c=zeros(q,n);
bak2=zeros(1,2^n);
for a=0:(q-1)
d=a;
for b=1:n %二进制的转换
c((a+1),b)=rem(d,2);
d=floor(d/2);
end
g=0;
for e=1:n %倒叙——矩阵角标
g=c((a+1),e)*2^(n-e)+g;
end
bak2(g+1)=y(a+1); %为防止数据交换覆盖，备份
end
y=bak2; %备份后的还原
FFT.rar （1个子文件） 
FFT1.m 2KB
林当时
- 粉丝: 113
- 资源: 1万+
评论0