CORDIC算法的matlab实现

大笨象

不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
) U5 ?7 K) f3 t, Tfunction v=f_cordic(beta,n)
8 r1 b- H/ g1 }( I0 B! Xif beta < -pi/2 || beta > pi/2
" V; t$ J- H9 S2 e3 K4 k  o; B    if beta < 0
: S$ u; a0 h  B8 w9 o' j2 u, }! N* w        v=f_cordic(beta + pi,n);
    else
4 L% T3 J9 d$ y1 P2 x2 y        v=f_cordic(beta - pi,n);
$ f0 x+ s( {5 N' v' M5 [0 h    end
    v = -v;% flip the sign for second or third quadrant
    return
end
% Initialization of tables of constants used by CORDIC
! D3 R1 }# Z; Z+ b  I! L% need a table of arctangents of negative powers of two, in radians:
1 {8 C, ^  n& ^; l. [% angles = atan(2.^-(0:27));
# i& ^! p( U5 b7 P6 {0 q1 k* ?angles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
" R" F$ E+ ]/ I" K6 X9 w8 ]    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
5 S$ a9 H6 Y) @% l) S% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
( F1 X' p* ^. U) b4 m; a  P  x9 qKvalues = [ ...
. D2 h( h. _* [& f: c1 [; P  V    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
9 |9 h6 S: ^1 [    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
# j. j; Y' O4 ]1 d$ O5 E    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
) x& I0 P3 E) t! B  w' D/ r) s    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
$ H( i% G; \. S/ K    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
, \/ V) x* A5 Q9 \: c3 z4 O    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
' L* a# f0 {$ z3 \8 P1 |- n' QKn = Kvalues(min(n, length(Kvalues)));

* H: T! f  X9 a& J% Initialize loop variables:
$ B# l/ ]9 Z8 Xv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
% Iterations
4 q* @  _; K2 M! Lfor j = 1:n;
# J8 O, ?- b% w5 U7 ], @  k; g    if beta < 0
        sigma = -1;
    else
8 M" \, l7 o% s% Y( n6 t        sigma = 1;
    end
6 s7 t  Y) w. t, q    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
        angle = angles(j);
    else
        angle = angle/2;
- b" b6 @8 n$ Z- M0 o4 _6 j5 O    end
    factor = sigma * poweroftwo;
6 c& Z# C* Y7 E0 V    R = [1, -factor; factor, 1];
1 I( X* t5 b* V& j) q" U4 |# S5 [8 v! W    v = R * v; % 2-by-2 mtrix multiply
- g) B. \4 e$ d" k    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
4 |. q; z! n' ]4 x  B( N, `end
% Adjust length of output vector to be [cos(beta), sin(beta)]:
6 X: h3 X% |( c% R& G/ W9 jv = v * Kn;
return
& I( {/ e3 q; y2 q; D2 _
  {; f$ `' r. g9 D9 |- B5 N  X使用20位寄存器，仿真波形如下：
9 E) R4 b9 I) y8 h

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
$ _$ I3 F! l, V7 c: p这是什么算法？第一次听说

- V* x: n2 i, H这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

赛才 发表于 2011-6-9 15:01
# n1 {8 z& e2 M/ b0 U1 ~  x9 I1 a我表示很给力

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010