CORDIC算法的matlab实现

大笨象

6 Y4 {  [, x& i- A0 j
不知道有木有人研究这个，比较偏向于硬件描述。
; I. `, A* e# a今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
( G" N& s2 K5 O4 }function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
    if beta < 0
        v=f_cordic(beta + pi,n);
) G* B8 S( F; r    else
        v=f_cordic(beta - pi,n);
, F( N& R  P: K0 \8 g* `+ ?9 }    end
    v = -v;% flip the sign for second or third quadrant
, X% p  I2 O, i+ r% x& q" r    return
end
7 @4 ^+ X( i5 ^+ U% Initialization of tables of constants used by CORDIC
9 S  U& c# a- a  k% need a table of arctangents of negative powers of two, in radians:
5 ]" h% w: L' N4 e3 h2 l% angles = atan(2.^-(0:27));
angles =  [  ...
) H: [: r$ w* @& V- x. R0 j    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
1 l" l  [. K3 Y4 _    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
$ d5 m% u# k( Z+ ~    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
* q* j  V$ H: |5 j- }+ _) Q: _6 a% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
& f3 e* M( u5 `0 m7 GKvalues = [ ...
    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 ...
- c6 V/ V' Y, H8 N1 [  ]: {    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
% a5 J7 g" _/ r; v/ Z4 z+ d( o! l    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
  d% t3 w' U# k% z: }Kn = Kvalues(min(n, length(Kvalues)));
" q3 }2 x! u. L2 ]3 b* L& K
% Initialize loop variables:
5 s5 I  n7 n0 X- u" yv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
% Iterations
1 b3 S0 P7 E' Y# z6 Q6 [. h! y7 ^for j = 1:n;
    if beta < 0
        sigma = -1;
7 T7 v. s3 {! n  _: R' ]/ o    else
4 v! a( ?* o8 Q) ]0 D6 [# c& {+ d5 s        sigma = 1;
    end
  ]6 c' e: J& p3 ?! Y    % update the angle from table, or eventually by just dividing by two
  ]& K) {8 B! L- j' i8 }! B" _( E    if j <= length(angles)
2 _% d, `: O$ c0 Z: S0 o4 a        angle = angles(j);
: ?& e) i+ s4 \) d! `, u. l    else
        angle = angle/2;
    end
    factor = sigma * poweroftwo;
! C, ^% \; j+ j; b6 `    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
3 ~9 G% ^9 I! d- M: @9 x3 S5 f+ oend
% Adjust length of output vector to be [cos(beta), sin(beta)]:
; i7 H' b* C* X' Q+ @. ~v = v * Kn;
- ]% w& s8 |* L8 [return

& _% i2 P% }1 A& d, j" E* Q& |4 g! ~1 }使用20位寄存器，仿真波形如下：
/ U' E' D9 E$ ]" n( y

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

& a' I$ Y/ O5 i这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

赛才 发表于 2011-6-9 15:01
1 D# ?" I. F+ y我表示很给力

+ U/ a( |3 \" f, c6 {8 W( n给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010