wfft.m

给出了数字信号处理课程设计的基本程序

function y=WFFT(x,N)
M=log2(N);
long=length(x);
if(long<N)
     x=[x,zeros(1,N-long)];
end
NV2=N/2
NM1=N-1
J=1
for I=1:NM1  %100
  if(I<J) 
    T=x(J)
   x(J)=x(I)
    x(I)=T
  end   
    K=NV2    %5
    
    while(K<J)  %6
       J=J-K
       K=K/2
   end
       J=J+K  %7
end
for m=1:M
    a=2^m;b=(N/a)-1;
    u=exp(-j*2*pi/a);
       for lg=0:b
           c=lg*a;d=c+a/2;e=a/2;
           for k=1:e
               w=u^(k-1);
               T=w*x(c+k+a/2);
               D=x(c+k);
               x(c+k)=D+T;
               x(c+k+a/2)=D-T;
           end
       end
end
y=x