CORDIC算法的matlab实现

大笨象

0 M& C! w7 L8 y% u& L
不知道有木有人研究这个，比较偏向于硬件描述。
% Y/ N; }6 A' V" S$ F2 E今天先发一个sine cose函数的实现。
% V9 r. X0 B, ?" c9 q4 U有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
    if beta < 0
2 u8 [4 O9 l; J        v=f_cordic(beta + pi,n);
    else
) O  x/ C3 W" z        v=f_cordic(beta - pi,n);
$ c- @) o, v! u* Y    end
' s3 K/ n2 @; A9 @    v = -v;% flip the sign for second or third quadrant
    return
end
* Z, f2 \" |+ i2 i1 i6 _) T% Initialization of tables of constants used by CORDIC
9 ]  U9 E; L# y; [5 T- w9 V$ W: S0 S% need a table of arctangents of negative powers of two, in radians:
, m2 G* l1 P- J( u) u4 r* i! t% angles = atan(2.^-(0:27));
' }; {# ~5 s2 C/ A- n( zangles =  [  ...
7 M0 g  J. v  N8 X; x    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
* H9 s5 m1 Y+ ]    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
: @7 ^4 G8 z0 }& z& e$ ~7 f1 y" r    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
% E, w8 t8 E, M: D    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
  t+ N# s6 P1 D% ]  w4 H- G. M) ]; e% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
* R3 P- X& ^3 \7 f1 ?    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 ...
4 O( J9 i( k3 g: l/ @    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

% Initialize loop variables:
7 o, M* u) P$ |# a4 p! Fv = [1;0]; % start with 2-vector cosine and sine of zero
7 y4 V1 c- m, ~* O/ hpoweroftwo = 1;
% X# H. z$ i, F2 V' [' l1 l2 u% Iterations
9 S# F0 j' Z) M% pfor j = 1:n;
+ q4 ^% p' m2 D) b5 l* @: \    if beta < 0
        sigma = -1;
( r. v) A# G4 s$ m* h    else
        sigma = 1;
    end
9 L2 l. z2 N& j5 _: f( e    % update the angle from table, or eventually by just dividing by two
% e# Z2 r+ H/ @9 A    if j <= length(angles)
        angle = angles(j);
    else
        angle = angle/2;
* N& M+ l# x& t    end
# X: H% F5 C% t5 ]3 \) p) I( H) F- e    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
# E. o* k5 M  m; j& ]" K    beta = beta - sigma * angle; % update the remaining angle
) y: {7 I. H3 j$ Z' |* {# A    poweroftwo = poweroftwo / 2;
6 s, U% E+ U& m: nend
! h; e( ?4 a3 z1 S7 ]% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
return

* ^) @4 [( W% p9 j2 m& @  R8 ~使用20位寄存器，仿真波形如下：

& z2 Y$ f( a0 j: z

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
6 S3 g# g+ T2 T这是什么算法？第一次听说

. F5 A. x% @) U) x这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

赛才 发表于 2011-6-9 15:01
) C5 v5 y7 D# p7 Y+ f8 G( C2 r我表示很给力

7 `) w$ p) k9 _4 p$ s  @给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010