CORDIC算法的matlab实现

大笨象

不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
/ C3 L7 G! d; x4 j7 \3 {; u! c    if beta < 0
        v=f_cordic(beta + pi,n);
% Z$ j4 @8 l3 o0 ^9 h# I    else
; ?' U. p3 D6 L" z) q/ N8 i        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
) ?; K" {. v: |2 e, W1 r    return
end
7 ?1 g2 w7 h4 q4 n6 S' |% Initialization of tables of constants used by CORDIC
% need a table of arctangents of negative powers of two, in radians:
# H  u1 a7 z* P% E6 A/ R$ \; ~% angles = atan(2.^-(0:27));
& J' I2 _1 K& E/ L% Q1 Zangles =  [  ...
. b" G+ T2 f- E    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
" X, O( f3 I9 ~7 p- ?- ~& N6 _    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
; T' {8 K' [$ Y& c* X% f( p+ [    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
0 n. l, v0 ^1 a6 `* Y1 b9 {9 N    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
6 d# F6 t+ v0 k! S    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
  v( i& q9 v3 [2 v- K' V" W5 l# ^% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
1 W7 `: Z$ `, rKvalues = [ ...
) R7 J5 T1 R) ?1 x- t5 o% L    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
, h  k4 S6 D9 U) M4 N9 i    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

% Initialize loop variables:
4 p2 t! ^% V; r7 _5 v( d: E1 Iv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
9 O9 u3 `1 Y/ @+ Z7 h! z% Iterations
for j = 1:n;
! X  C! b2 B- P( I1 J    if beta < 0
        sigma = -1;
/ j" l/ I; r' E    else
        sigma = 1;
    end
: t: L$ p: b) S1 {$ Y    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
        angle = angles(j);
    else
+ V. H. Q5 l' e& A4 v2 L0 w        angle = angle/2;
    end
    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
) y& j! B# h1 v    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
end
# Z3 k! w% v, _% Adjust length of output vector to be [cos(beta), sin(beta)]:
) Z" v2 t' y" r8 Hv = v * Kn;
! i% H1 [2 }$ j+ ]/ b' W  W. z$ Zreturn

3 H! N; B; _7 B# ]7 G0 ?使用20位寄存器，仿真波形如下：

. m9 {: s/ z9 I  K% }

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

  V( X, W! M! ?6 A3 q% k* T' N这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

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

3 i# q' Y3 X, y0 I! \给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010