CORDIC算法的matlab实现

大笨象

7 ~2 L3 H, {  i: _1 M! [
不知道有木有人研究这个，比较偏向于硬件描述。
! X% s( Z. o: y今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
$ H1 M) `* I8 n1 S3 p" h3 Q1 x    if beta < 0
+ b- t+ U$ Y! b0 }8 ^  r        v=f_cordic(beta + pi,n);
8 B* X4 j  G, y' j8 U! q    else
& b8 z' i8 P. M4 \3 c* D% h        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
    return
end
3 a+ r/ H" H& v& m. y% 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));
angles =  [  ...
) Q! C/ c/ }& P8 M( |* D    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
' }. }% v& h; @    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
7 t! k3 E; v2 ~4 P6 r3 N& E    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
- h3 i+ Q# G1 v& i; Q    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
! J! _0 e, a& c& Y    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
$ k( s/ h4 i5 s5 @9 }( q    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
- ?# E2 @1 M0 s& X% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
5 ~8 I6 U. {2 M( @    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 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
* \4 c' U: r6 l! s+ i' QKn = Kvalues(min(n, length(Kvalues)));
3 H0 a- ?+ ^' @
% Initialize loop variables:
2 h3 J/ u4 a1 J' K% Bv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
% Iterations
4 u1 N- e. \/ P% b% cfor j = 1:n;
    if beta < 0
" X! z. ]0 ~, x$ Y: z9 B        sigma = -1;
    else
1 x$ w% ]' n) G        sigma = 1;
# O0 {) S6 \. p9 L* q    end
, s0 ^% a9 U% k, |* h    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
        angle = angles(j);
. e9 ^1 E, O' P( A4 K9 J+ l    else
        angle = angle/2;
4 \/ q3 q' k' j! L/ P! F9 r" n8 g    end
: P& Y' w% c1 G7 C+ w6 V    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
" p8 Y' h  j8 W- e1 o4 t  r7 n    v = R * v; % 2-by-2 mtrix multiply
, W; N/ S, A; T    beta = beta - sigma * angle; % update the remaining angle
  U& d5 |- Z* ^3 |$ e/ o+ E    poweroftwo = poweroftwo / 2;
3 o' ^& A+ F' O1 Zend
2 A6 f0 m6 D- _7 j1 C# c: ]7 \  p% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
) z% G5 @. @* J  o/ Zreturn

使用20位寄存器，仿真波形如下：

7 i  H" ~, C) c5 z

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

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

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010