CORDIC算法的matlab实现

大笨象

' a, c) h9 k, X' X
不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
" A  \% P3 A. d0 p0 `3 G+ A/ Ffunction v=f_cordic(beta,n)
$ I4 m" B' W9 K& \6 [( G4 Wif beta < -pi/2 || beta > pi/2
) m7 Z$ W7 v) Q4 z& P    if beta < 0
0 }1 [5 J7 x) |7 S        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
    end
1 Y: }" l6 v$ w1 p' a# |    v = -v;% flip the sign for second or third quadrant
    return
end
% Initialization of tables of constants used by CORDIC
# b# y* q' B, C/ E+ y0 q' o% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
angles =  [  ...
' D* V. e7 O' g5 e4 A/ S  U    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
- b: r: @& r# u$ q9 t3 t    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
1 Y9 N0 G* n4 g  h) w2 \    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
: l/ b" p0 I/ X% m5 I; I' ~2 `    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
. H. G. p9 z  x1 O    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
4 x; G4 O' l8 W    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
! F: i' X+ [. p8 e: L    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
/ B- V' l9 K) l6 rKn = Kvalues(min(n, length(Kvalues)));

% Initialize loop variables:
7 e% k$ I( P; \' K" lv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
  o2 G' F; N& ]; y) d" h' P( x% Iterations
- |! \% G7 U: |( \! g5 Rfor j = 1:n;
    if beta < 0
        sigma = -1;
    else
        sigma = 1;
    end
    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
  b6 K- S; z9 j9 t1 r        angle = angles(j);
/ T/ E" q( L" |( _: I2 l    else
        angle = angle/2;
1 O0 P2 f: K; [# h) s) H    end
    factor = sigma * poweroftwo;
5 h( N) i, U. E# C* l: A) U    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
, M  d. Y7 B: E4 j& e1 p    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
end
% Adjust length of output vector to be [cos(beta), sin(beta)]:
* L& P. [- Q$ r5 m4 Z5 E( P" @+ s) Dv = v * Kn;
5 C1 w5 W0 f( ^5 areturn

+ U- z; N. }( G9 h9 Q- U使用20位寄存器，仿真波形如下：
  z4 W9 }2 \1 C: Z

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

赛才 发表于 2011-6-9 15:01
$ ?8 z# O5 V" r' c& e3 `- u, o我表示很给力

5 O: M; ?' Q* L/ H) V; T9 M给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010