CORDIC算法的matlab实现

大笨象

8 a6 \' `3 y% ~% B& O7 R6 `& |
不知道有木有人研究这个，比较偏向于硬件描述。
! r& Z% c8 q4 C8 I0 C: a% E今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
3 ]# L6 [5 E) \7 lif beta < -pi/2 || beta > pi/2
    if beta < 0
* }2 w0 x  p$ H0 n2 o        v=f_cordic(beta + pi,n);
3 N6 Q2 g+ _" p' d7 `# r% U  p    else
+ y, l* ~! s4 f$ x        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
    return
end
$ |; C+ Q: H8 j: w% Initialization of tables of constants used by CORDIC
9 H) x: A: {: n9 w% need a table of arctangents of negative powers of two, in radians:
; p1 N5 c& {) W, Q# \% angles = atan(2.^-(0:27));
* y* M; k$ u% Y1 z* k: Q: `1 Qangles =  [  ...
* W6 @! l# {' {8 v$ J& w    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
0 x1 |" \* Y( D! |, F* i" ]    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
$ D0 }7 g4 e! ^8 @  E  U8 K    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
" p* v& {' x6 Y% H( @* h; t    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
2 G+ U& A9 X8 F0 d9 g    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
8 x5 @6 h  X/ }( D    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
7 Q) }; l, f6 Y9 R6 v* X1 _' _# W    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
& K. u; L- k/ ]1 {- C5 Z    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));
2 Z  B; W3 F/ b' B
% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
% Iterations
for j = 1:n;
    if beta < 0
        sigma = -1;
: ], `+ o4 p, z! \& i    else
  p4 \3 ^& s, F, e0 C: R3 o7 N        sigma = 1;
9 [! p! Y# w% C+ M/ G1 R    end
9 ?( L, w) m/ p$ _7 z    % update the angle from table, or eventually by just dividing by two
' p/ H! j3 C/ d( L. l% k    if j <= length(angles)
        angle = angles(j);
    else
/ W$ n  I% l& L2 A( I2 V9 o        angle = angle/2;
    end
( \" M: _/ U, C: h( Y7 i    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
7 L# T' q- S# s( ]6 T    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
" q2 X. r; J5 X) t& B# c    poweroftwo = poweroftwo / 2;
4 h- t& a" T& t/ wend
1 M) D" _+ S+ M$ j4 w* m  T% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
- J8 r, m4 v9 I, D+ @2 @: _! O7 \return
& w) z* n$ a1 ~- h5 X# n  U( M3 {
使用20位寄存器，仿真波形如下：

8 A2 r0 B2 Q( w7 v0 e7 I

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

1 L* m1 m$ Y' r# j+ \1 w: K$ J3 @6 x2 ]5 @6 X这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

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

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010