FFT快速傅里叶变换Matlab自编程序

来源：飒榕旅游知识分享网

% 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

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;

程序运行结果

-52.52x 10 Matlab内置函数自编程序Amplitude(V/Hz)1.510.50 61010Frequency(Hz)7108 50Matlab内置函数自编程序 Phase Angle(degree)0-50-100-150 61010Frequency(Hz)7108

因篇幅问题不能全部显示，请点此查看更多更全内容

查看全文