CORDIC算法的matlab实现

大笨象

, M4 ?( n7 J  f$ S* Y: \不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
6 w) J" T. \! X  T9 u5 D有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
    if beta < 0
        v=f_cordic(beta + pi,n);
- x, c  F  ?, u  |. ?$ M: ]    else
        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
& D" X$ a+ u' X& ~& e3 |+ |    return
* \% A: S1 w3 A3 K# Z+ v5 J0 `end
% j. l2 t& H' I% 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));
% U) @% C1 X. ^$ xangles =  [  ...
$ V+ i; r& u/ W# F. b  o    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
( t) S/ R( j: r% X3 R5 i7 o; q    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
' g1 o; L4 m* m5 ?) V2 g& _    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
0 {" `: H; v1 ~! [9 H% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
' e3 Y% S! o& Z7 P0 D  IKvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
) v  V1 ?& A  A1 M6 l+ J    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
5 I) J4 y3 M  W    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
. p4 {1 Q4 x1 d7 Q8 w" |    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

% Initialize loop variables:
1 }# O+ ]; Z. ~9 ^1 wv = [1;0]; % start with 2-vector cosine and sine of zero
  {- L* Y* U/ R% R' Mpoweroftwo = 1;
  s' _" G- ^* _  m4 p& K% Iterations
2 ~! m8 U; T# j+ }( B7 @for j = 1:n;
    if beta < 0
  m3 c6 k' S. y+ u2 }        sigma = -1;
    else
        sigma = 1;
+ ^# s# A4 e( S0 r# R9 G9 K' t    end
    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
( G$ B* G- q- `( |4 D; l        angle = angles(j);
3 X8 Q9 Z3 t' I. W6 c0 u# T  m    else
0 |# z1 b( O. l: C6 o, V        angle = angle/2;
    end
6 m  u4 ]$ [; b+ B, b6 l" K6 c    factor = sigma * poweroftwo;
7 y. x4 A" K, u# q9 z$ P3 n4 e! p    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
; K) K4 L" L0 E    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
end
3 @  u: \4 y) h* |) \% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
return
- T3 W+ T4 ?9 X, c. L2 k( U( ?
使用20位寄存器，仿真波形如下：

. A" T/ a' ?6 W( B8 R+ ?

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
4 I3 k- r& H6 j' `9 ?: ]) M6 U, d8 x这是什么算法？第一次听说

1 Q% `; g1 G* ^% w这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

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

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010