CORDIC算法的matlab实现

大笨象

不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
2 G+ {5 d  f! i: P/ n  y    if beta < 0
3 ~5 v# g* t8 D1 ]& Q        v=f_cordic(beta + pi,n);
5 J7 T. H- v, A5 \3 B( t3 o    else
6 |# e$ c$ E& J        v=f_cordic(beta - pi,n);
1 M' e% A1 X, X- L' j! j6 d    end
    v = -v;% flip the sign for second or third quadrant
    return
end
% Initialization of tables of constants used by CORDIC
5 _8 }6 v$ p. b0 q5 a% s6 x- s% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
angles =  [  ...
* o% b3 |5 l1 R3 n6 p    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
7 k* U' u' j% o: w    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
6 [9 [1 ~* M  ]$ b! @0 V    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
3 ]3 t  r* l, Q& P- F+ G  q* r    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
7 |  c0 f8 l  J7 }% z    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
+ o/ J2 ^! s& x6 J, T/ h9 z% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
/ z) B( a9 A# R/ Z7 }2 F8 m    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
( Q5 l* P- x+ e: K    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
/ X* q+ l! Q' Y% H8 C    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
  D( q/ Q- c: M0 J4 l: L7 B$ ]    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

% Initialize loop variables:
0 T  i! ~$ Y  I- R# z- K" hv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
+ _" J" H9 }+ _0 t9 I; _% Iterations
+ ]2 d* z0 L# J" ufor j = 1:n;
5 h9 ?' B' J7 t    if beta < 0
        sigma = -1;
    else
$ ~4 S" I* s! ~& `        sigma = 1;
    end
    % update the angle from table, or eventually by just dividing by two
: M3 s# k+ R3 `0 l( _  G7 s    if j <= length(angles)
! y- A, ]5 t6 t! w& C        angle = angles(j);
' M: A+ ^* _7 m' ~    else
        angle = angle/2;
    end
7 u4 i4 q4 X% H2 W3 {) K    factor = sigma * poweroftwo;
' Y' @# K) A+ h    R = [1, -factor; factor, 1];
* E5 b! a4 F2 n1 A/ V    v = R * v; % 2-by-2 mtrix multiply
* T2 x+ d* U( Q    beta = beta - sigma * angle; % update the remaining angle
" E. p: Q' t, m    poweroftwo = poweroftwo / 2;
$ {" q  T) _$ s. Zend
% Adjust length of output vector to be [cos(beta), sin(beta)]:
. i) g% e: h" N% Q  sv = v * Kn;
return

使用20位寄存器，仿真波形如下：
- Q- j. ~! J1 u8 z$ j0 U4 J% ^8 A
/ y: a9 r3 x" o$ f0 o! v( T$ ?

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
这是什么算法？第一次听说

+ i9 E2 P# A$ w这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

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

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010