dfft.m

基于MATLAB的fft电力谐波分析 matlab代码

function [y,N]=dfft(x)

for m=1:inf  %找x的长度对应的2的最低幂次
    if 2^m>=length(x);break;end
end 

N=2^m;

if length(x<N)
  x1=[x,zeros(1,N-length(x))];%补零
end

y=daoxu(x1,N); %调用倒序排列函数

for mm=1:m % 做m级蝶形运算，每一级都有2^(mm-1)个蝶形结 
Nz=2^mm;u=1;  
WN=exp(-i*2*pi/Nz); 
for j=1:Nz/2 % 蝶形运算过程,在进行第mm级运算时需要2^(mm-1)个 蝶形 
for k=j:Nz:N % 蝶形运算的跨越间隔为Nz=2^mm 
kp=k+Nz/2;  %k,kp为碟形运算的对应单元下标
t=y(kp)*u; % 蝶形运算的乘积项 
y(kp)=y(k)-t;
y(k)=y(k)+t; 
end 
u=u*WN; 
end 
end