test_r2_dit.m

来自「dft algorithm(no fft) (matlab)」· M 代码 · 共 50 行

%------------------------------------------------------------
%   Complete this code that implements
%   Radix-2 DIT Algorithm for 8-point Input Sequence.
%------------------------------------------------------------
xn = [ 0 1 2 3 4 5 6 7 ]
N=length(xn);

%------------------------------------------------------------
%   Bit Reversal for Input Sequence
%------------------------------------------------------------
stage=0
xn = bit_reverse(xn)

%------------------------------------------------------------
%   Stage 1
%------------------------------------------------------------
stage=1
for ii=0:(N/2-1)
   [r1, r0] = dec_to_bin(ii);
   i0 = (N/2)*r1 + (N/4)*r0 + (N/8)*0;
   i1 = (N/2)*r1 + (N/4)*r0 + (N/8)*1;
   NN = N/4;
   WNkn = exp(-j*2*pi/NN) ^ 0;
   [xn(i0+1), xn(i1+1)] = r2_dit_ce(xn(i0+1), xn(i1+1), WNkn);
end
xn

%------------------------------------------------------------
%   Stage 2
%------------------------------------------------------------
stage=2
for ii=0:(N/2-1)
   [r1, r0] = dec_to_bin(ii);
   i0 = ?
   i1 = ?
   NN = ?
   WNkn = ?
   [xn(i0+1), xn(i1+1)] = r2_dit_ce(xn(i0+1), xn(i1+1), WNkn);
end
xn

%------------------------------------------------------------
%   Stage 3
%------------------------------------------------------------
stage=3
for ii=0:(N/2-1)
?
end
xn