CORDIC算法的matlab实现

大笨象

. _) ~7 K3 o5 r7 i! u: O) _, g不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
" ?9 O: e" a- a8 o+ _if beta < -pi/2 || beta > pi/2
/ c. ~* {+ m' ^6 w' u& q  L4 G( D    if beta < 0
2 y' P" f0 e- Q( B* M" E2 l4 @        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
0 ]* }$ r1 f/ h, P; o% W& ~, e    end
; \; R  Z7 \  c( X& m0 \5 n    v = -v;% flip the sign for second or third quadrant
    return
5 D0 k0 [: h6 P& v2 f" d" m8 Dend
% Initialization of tables of constants used by CORDIC
% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
+ j1 C* u$ x5 j& O9 d0 fangles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
0 F" _+ E. [, G. N% {  s6 f+ _    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
$ A8 t; |. v6 J' K* Z    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
. Z3 E: u" Y  Q7 w' }) X    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
9 n' f6 A7 {" a- w% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
2 q3 ^& Z/ L$ m* yKvalues = [ ...
2 Y& F& z$ ?2 [    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
8 V$ d5 D4 T% f$ r  G3 U$ R    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
7 e0 f6 O8 Q) K: W" kKn = Kvalues(min(n, length(Kvalues)));
% L/ U, o& Q. I5 y5 l# P" k( l
* n  p0 z: ^/ p# m- b% Initialize loop variables:
9 G) e# T  \. ?& P  ~8 _( W( lv = [1;0]; % start with 2-vector cosine and sine of zero
' U% i. I3 L& r4 O$ C' npoweroftwo = 1;
% Iterations
for j = 1:n;
2 ]  g: B+ `  F/ ~% b% W    if beta < 0
        sigma = -1;
, e6 L: I* A4 _3 P/ U  U    else
7 h  O  o. T9 c+ A( E        sigma = 1;
, P- K' d; l% Z6 d" K: Y    end
    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
        angle = angles(j);
    else
        angle = angle/2;
    end
    factor = sigma * poweroftwo;
1 d! t8 Y! [4 H) y/ `$ w; U4 x3 ^# C- ?    R = [1, -factor; factor, 1];
* m" ?9 |  p9 N# E0 i  m9 [! i    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
' S* _  _4 Z) h2 I    poweroftwo = poweroftwo / 2;
7 q% `  V( ]3 w: ~3 cend
0 v- Z9 @' }8 g% Adjust length of output vector to be [cos(beta), sin(beta)]:
& U; v2 D: `, ?v = v * Kn;
return
  F7 s! R- s2 k" ?2 g) i
2 T  r: B+ F& O# R' r- X+ C8 V6 T) w使用20位寄存器，仿真波形如下：

' R" Z1 ]3 Z3 P+ g

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

赛才 发表于 2011-6-9 15:01
7 G) f& l$ Y+ [: X/ a% z0 D# _我表示很给力

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010