twiddles.m

ADI TS系列的DSP FFT源代码

function [id] = twiddles(name_file, N)

   for k=1:N
       twid(k)=complex(cos(2*pi*(k-1)/N), -sin(2*pi*(k-1)/N));
   end;

       fin = fopen(name_file, 'wt');
           for k=1:N/2
               fprintf(fin, '%e\n%e\n', real(twid(k)), imag(twid(k)));
           end;
       fclose(fin);