CORDIC算法的matlab实现

大笨象

. h# h9 h7 g* \/ N7 i不知道有木有人研究这个，比较偏向于硬件描述。
* S' e( w0 C  f3 w& @& c9 @  a$ A今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
+ r0 P/ Q$ a* |- M6 y  P    if beta < 0
# U+ \4 V% L6 D! V# Z        v=f_cordic(beta + pi,n);
. j) c6 w' ]5 \( u$ n    else
4 {( o) U1 d! j9 {. ^4 i7 e        v=f_cordic(beta - pi,n);
+ K/ k: j5 q& n7 N% H    end
    v = -v;% flip the sign for second or third quadrant
    return
end
% Initialization of tables of constants used by CORDIC
- n  \) [2 ~  {! }$ [% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
$ W8 e( w# [8 N( C0 Xangles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
! ^7 q4 [7 B: r: E& M2 e9 y    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
4 y' [0 ~, c* V# I, R9 \4 I    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
9 `8 H$ R" z' F' y$ f) Z# T    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
6 K; L, n" @4 t( t4 \    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
6 c" m) c8 d3 v& E0 {- o    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
* l" h% D9 h3 h3 W$ S" `3 K) nKvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
, O, i. h( J  j8 v5 e    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
$ ~: g# C& B/ I+ E    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
5 q* @  A: R" l3 F2 A. w% R9 P    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
% [( d0 z: n# J. I' j1 I7 w( Z    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

. K+ h' _$ O) L, @0 z  U* \# Y2 q% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
. ^0 t/ u: }4 Hpoweroftwo = 1;
% Iterations
3 p+ t! f3 s  F: |# e1 N- bfor j = 1:n;
    if beta < 0
( n. ~/ x% l" Q; _        sigma = -1;
    else
        sigma = 1;
    end
    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
$ O2 @* C( e$ F" B: _  M9 \( p* \3 N        angle = angles(j);
4 D5 A( n7 c; x6 o  W1 x9 j    else
        angle = angle/2;
% H/ s, r; {; i4 I    end
    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
! C4 t% K' s9 h8 F- a    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
9 O3 c5 z, c7 T* c6 j9 V$ c+ i4 W- Rend
% Adjust length of output vector to be [cos(beta), sin(beta)]:
1 g' W% W: o" d. e' U1 u( ]6 Rv = v * Kn;
return
4 k* i, W- j" ^# w
使用20位寄存器，仿真波形如下：
* z  g5 `* Q- g$ p* D- b
2 `/ d9 W# w/ X9 d

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

赛才 发表于 2011-6-9 15:01
5 f# e6 L. L6 H* L" W) Q( O我表示很给力

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010