CORDIC算法的matlab实现

大笨象

( {* q2 h+ Q- w% @* W& w7 V0 P
不知道有木有人研究这个，比较偏向于硬件描述。
2 @6 B4 O; m0 U8 T5 `, f8 J  b+ @% r今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
4 w+ n( d6 i& ?+ ~" u    if beta < 0
        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
6 s: B8 ]+ S* L; b) u3 M6 o0 i    end
3 E6 B0 u: @) G. ?    v = -v;% flip the sign for second or third quadrant
    return
end
8 G! k  a  r+ G/ R% Initialization of tables of constants used by CORDIC
2 Q* }# U  `  e/ j: }% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
% a7 D* q8 r" c) x# @angles =  [  ...
. ~. s2 h8 Z! k! k- v4 D    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
) W" B  e3 K( m% X1 U- K& O    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
4 {3 j! D, a; j# L$ y    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
1 P/ L) e2 @$ V4 m# ?/ z8 y4 _( x% }/ Y    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
4 f) k5 g; O$ x5 b* i% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
9 {, X4 R! m# h; S- w& jKvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
3 ^9 n' X$ e% c& J    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
8 b0 F& A: {  s  k/ _/ \  k: r* u    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
7 p; u# l; `8 x0 y    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
* X6 |+ w) D, n2 I    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
$ u5 |9 Z" B. ~: ^# R0 IKn = Kvalues(min(n, length(Kvalues)));
) p7 `+ K% \5 h8 u* g0 h, R1 F
. |) Y8 c8 H6 q/ [' e  J% Initialize loop variables:
9 Y' d& ~+ o1 v' o( c% ?6 E, |v = [1;0]; % start with 2-vector cosine and sine of zero
4 |1 H$ r: \2 ~4 O  @/ b! ppoweroftwo = 1;
' k6 a& p6 R; f: u. ]* L% Iterations
* y9 r$ v2 J. H( G' Zfor j = 1:n;
3 j! w% N5 {7 M( B    if beta < 0
        sigma = -1;
    else
2 ?: h# O8 K) v- t* X9 n: N9 ]        sigma = 1;
- G+ `- ?9 R" `; D# h1 i    end
# `3 J* K! b5 G" D3 U8 V    % update the angle from table, or eventually by just dividing by two
5 q2 I9 P; |; V* n( H8 ]    if j <= length(angles)
        angle = angles(j);
    else
        angle = angle/2;
    end
    factor = sigma * poweroftwo;
7 b  d6 `% W/ F+ `0 e3 [2 Y    R = [1, -factor; factor, 1];
/ X/ h0 M; h1 C1 Q4 {& F    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
end
5 H# f- V# U7 q' n% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
return

使用20位寄存器，仿真波形如下：
" ]6 O! ]5 _: L6 G
  T6 Y' f5 L+ ~, o: q; i

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

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

% P8 T2 l  p+ ~+ g, Q9 H给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010