CORDIC算法的matlab实现

大笨象

: P, X1 r- h8 b/ {1 g% s
不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
9 S, ?. M" `$ B' W9 n* W1 u有兴趣的可以看一下，m文件见附件。
* J, R2 P& Q$ b' Z* ^. A( K5 Yfunction v=f_cordic(beta,n)
  O9 |" O) u4 H& T4 J- Uif beta < -pi/2 || beta > pi/2
" g* ^- N. c0 J$ ~    if beta < 0
        v=f_cordic(beta + pi,n);
. O  r! c9 }* J2 a7 N' W* X/ ^    else
  F' U& `* Y4 |  |        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
9 }) l; H& a. y3 T9 d6 V" [- Y. c) `    return
end
" d0 [% p" f( E6 p% Initialization of tables of constants used by CORDIC
8 n$ d) M* U: Y% need a table of arctangents of negative powers of two, in radians:
7 C& x, Q+ v7 `8 Z: l% angles = atan(2.^-(0:27));
8 c  u0 M) a- l+ ]; N7 Oangles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
& N# `7 L" c6 E: @" J8 B    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
( i+ a8 B/ I2 \3 F: ?# R6 H. b) W    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
6 C5 Y1 W$ R# \" J    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
2 y0 H$ r) O% H) S' f7 n4 R% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
! ?4 I* p# R1 p# d- q    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
+ _9 x7 q' ?" e    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
* e* M' }: o9 f- Q+ L    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
+ ?- o  h, ^0 }' V& e    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

% Initialize loop variables:
+ l& |& f# }4 O( tv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
' w  k" [4 {8 k, A4 A) j" X4 @% Iterations
for j = 1:n;
    if beta < 0
$ j4 T9 O( Z5 Z        sigma = -1;
' S0 v# |/ L# k$ }% o' X" h* Y    else
  n( Z* ^2 J" t$ W7 e$ N        sigma = 1;
8 h& g2 j: }1 ~( H    end
    % update the angle from table, or eventually by just dividing by two
; F. C  p" M' M" D0 K5 }; O    if j <= length(angles)
        angle = angles(j);
/ {4 k- x4 {+ Z) o! l    else
7 W0 U  C; W+ M+ m/ D, a7 r        angle = angle/2;
    end
    factor = sigma * poweroftwo;
' l& c% b2 K2 Y* |8 b* m4 d# p    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
/ Y9 U1 u% L+ X" B3 }0 M    beta = beta - sigma * angle; % update the remaining angle
2 R4 N% e& C/ }2 y+ l* ~, }9 A; s$ y    poweroftwo = poweroftwo / 2;
end
  S" C# p7 o* [4 @  Z% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
return
: R" d; S2 m7 ~5 A% \
使用20位寄存器，仿真波形如下：
  X  f0 i" M1 e* G
- ~5 v. Q0 ?5 O0 a% z  ?: ?

杨帆

我表示学习了

zhangjichang

这是什么算法？第一次听说

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
  l. c0 _, Q) I7 f# `) O$ g这是什么算法？第一次听说

" ~% y5 J' @9 C这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

赛才 发表于 2011-6-9 15:01
我表示很给力

; H* x7 |  {9 `3 \+ y给力就好。呵呵

魅影骑士

顶一下！

yingzhen

费解。。。。。。。。。。。。。。。。

羅雲琦

有什麽用？求解

jt202010