CORDIC算法的matlab实现

大笨象

% C8 _% a, _" Y3 x3 S8 h2 v: H
不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
  X* V4 K% E  _7 E8 efunction v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
    if beta < 0
        v=f_cordic(beta + pi,n);
  h; n# G$ V4 R4 s7 L4 b0 {! q) L    else
        v=f_cordic(beta - pi,n);
    end
' J- J* {- e+ M% z    v = -v;% flip the sign for second or third quadrant
    return
: e' N2 G0 W) j3 g2 Gend
7 a; S- [4 `$ A) W% Initialization of tables of constants used by CORDIC
1 y3 `+ A, M' t; Z: ^, Z% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
& B! J: q) r0 E: uangles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
( q  ~& S7 |9 l( Q) W    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
: p" K: ~; t- [# n    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
% k! R% _6 W2 ~! ~8 l    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
2 _8 m- E# d1 ~    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
5 l2 N3 \! O! \! \9 e3 z% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
/ p, v: r  J) f" d: `    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
9 G  r0 I/ d% r    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

* l' u0 Y9 u" d% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
$ c, B" {  F. l% l% Iterations
: c3 k2 _: V  x7 w* Z+ ?+ pfor j = 1:n;
    if beta < 0
& c0 F" s& A. U, T% e  z6 B3 r* z        sigma = -1;
    else
# S& t  F0 g4 s" \        sigma = 1;
    end
5 W; c8 X. \, ~; M( a6 T4 g3 Q    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
* c2 p, t1 o2 L8 q4 g$ Z! ^        angle = angles(j);
    else
5 {1 z9 U: z7 R9 ]+ _2 ]0 V        angle = angle/2;
    end
    factor = sigma * poweroftwo;
* C% g: e: }+ O) k' d" B9 q1 _$ E    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
  Q; G! i3 [% ^/ \% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
" d" s3 I6 E+ P) P' o2 O: X2 hreturn

+ l' N5 y& S* R. W7 K6 W% ]" y) ]使用20位寄存器，仿真波形如下：

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

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

; y5 d8 G0 t, A% `5 E  s8 ^给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010