FFT快速傅里叶变换Matlab自编程序.
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
% fft function;
%
%
function y=fft_new(x);
M=length(x);
n=log2(M);
K=M/2;
even_no=zeros(1,K);
odd_no=zeros(1,K);
f_valve=zeros(1,M);
for i=0:K-1
for x1=0:K-1
even_no(i+1)=even_no(i+1)+( x(2*x1+1)*exp(-j*2*pi*x1*i/K ) )/K;
odd_no(i+1)=odd_no(i+1)+( x(2*x1+1+1)*exp( -j*2*pi*x1*i/K) )/K;
end
end
for i=0:K-1
f_valve(i+1)=(even_no(i+1)+odd_no(i+1)*exp(-j*2*pi*i/M))/2;
end
for i=0:K-1
f_valve(i+K+1)=(even_no(i+1)-odd_no(i+1)*exp(-j*2*pi*i/M))/2;
end
y=f_valve.*M;
程序运行结果