CORDIC算法的matlab实现

大笨象

1 q. a4 q9 g2 M2 f! K; @
7 y  C& Q" ~1 ?; H$ ~6 f# Z不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
$ u; f) Y- |+ e& p有兴趣的可以看一下，m文件见附件。
) g  [+ O# b9 n+ F) N+ o" ^" X  qfunction v=f_cordic(beta,n)
8 d& t7 e, ~7 V# _if beta < -pi/2 || beta > pi/2
  V- A: Z2 B( r0 @& _) i+ s* Y: a    if beta < 0
; j& @* c/ w; F' n& P+ \& k        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
5 k/ \# b, Y3 Q0 ^0 g5 o    end
    v = -v;% flip the sign for second or third quadrant
    return
7 w: u9 p0 M. `end
: K" _4 V5 U  }) @# H* R% 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));
. s  m# V& C5 x; t' x* Z7 j( n- Q; kangles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
7 P: \7 D! [9 a6 L/ P    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
; W: M" A+ p( u+ P( \! w6 r    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
% D# j+ r! J9 B% C4 o2 `0 S    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
" P! T) h! v  C" y: Q6 A    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
* c7 e$ Y9 {0 a1 z  n9 }/ yKvalues = [ ...
- a3 F9 |; R% p3 l    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
# \9 ]) S- i3 ^. J, f8 g    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
" x: B6 T: x7 x  x0 a. b0 b    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
' j; t# }& D" p0 I& m: N    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
. j/ X+ h; E7 e& `3 TKn = Kvalues(min(n, length(Kvalues)));
7 q) _+ c& J. l' l
% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
" }' j( r, }0 {5 F; J2 l' Y% Iterations
1 |5 \+ W  D/ x2 I" o! e1 D' E" Pfor j = 1:n;
    if beta < 0
; c; z; B- o5 T, o! ^) u3 g        sigma = -1;
: a  @. W2 b3 G9 \2 X' L    else
        sigma = 1;
    end
    % update the angle from table, or eventually by just dividing by two
8 O. r! G" b* Q9 j( d0 @  T    if j <= length(angles)
! q8 }  t6 j4 N        angle = angles(j);
' z, a% a' P4 [- I: Z. V1 N    else
; b- W: b: G7 P1 g2 F& d+ y        angle = angle/2;
    end
$ T% w. _  N) A; l    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
end
# d/ O/ M0 w1 F& T% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
return

% C4 Z: e) r3 ^# [* X- Y使用20位寄存器，仿真波形如下：
4 e/ @0 O0 G( R6 v

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

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

4 S6 H4 @1 Y3 O. M给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010