CORDIC-based DDFS Architectures - Digital Design of Signal Processing Systems: A Practical Approach - page 509

809859
720633
620162
510013
391906
267683
139283
8710
-122000
-250805
-375697];
tableY = [
65435
195315
322148
443953
558831
664988
760769
844678
915406
971849
1013127
1038596
1047857
1040767
1017437
978229];
cosine = [];
sine = [];
for theta = 0:.01:pi/2
angle =[0:.01:pi/2];
theta = round(theta * 2^(N-1)); % convert theta in Q1.N 1 format
% Seperating bits of theta as b=b[1],b[2],....,b[N]
for k=1:N
b(N+1-k) = rem(theta, 2);
theta = fix(theta/2);
end
% Compute index for M = 4;
index = b(1)*2^3 + b(2)*2^2 + b(3)*2^1 + b(4)*2^0+1;
x(4)=tableX(index);
y(4)=tableY(index);
% Recoding b[M+1], b[M+2], ..., b[n] as r with +1,-1
for k=M+1:N
r(k) = 2*b(k) - 1;
end
% Simplified iterations for modified CORDIC rotations
for k=M+1:N,