twiddles.m

来自「ADSP-TS101S and ADSP-TS201S Real and Com」· M 代码 · 共 13 行

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);