CORDIC算法的matlab实现

大笨象

* A2 w/ B! R6 b& n4 L; E( p
不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
, S4 `3 `: b6 t9 x: \" M有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
5 A9 m& [" i( i, m( b: v- W    if beta < 0
        v=f_cordic(beta + pi,n);
, e9 y) q( s) p8 H- O0 i1 Y    else
        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
    return
  x! z( i; s8 `9 Y3 oend
% Initialization of tables of constants used by CORDIC
: _/ H' y  R- l0 ^6 x% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
angles =  [  ...
4 _: e- u$ H; t9 ~0 {    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
5 B8 {6 e! l5 z    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
6 {  l: @* y& v5 r    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
, ?3 A: }$ D1 u. j! Z- B! \: h* I    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
  j& b' X) S9 h9 NKvalues = [ ...
( J- \" B* j6 K; m. d$ j    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
6 z6 \6 U/ j3 d& U    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
2 v  U+ z7 s5 C8 K6 w4 f- A& }- e    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
1 T+ }/ ?; e9 K) z0 }( Y" Z& j    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));
# O# m# g* r" U2 ]) S8 N9 J
1 L% \$ e1 J0 e. `% @7 N6 r1 l! s% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
% k; B2 q% U8 @% Iterations
) A/ k( I) \$ M( T( Qfor j = 1:n;
0 g$ {! y+ B( y# A6 l$ ^, Y    if beta < 0
        sigma = -1;
, P* U) ?& ?& r, C2 D    else
        sigma = 1;
    end
) h' `' w. Y0 P# k+ v! Q8 e    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
        angle = angles(j);
    else
1 u" L4 g2 B/ l& u8 m        angle = angle/2;
! l: t- j/ p/ A& [    end
    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
' I# D7 [- N9 P0 X. E1 \: \& E" vend
, Z0 m  o- L) a4 z7 G6 s1 z, W# L% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
8 `# M: z( N( T2 X- jreturn
3 h. E: K* }, w7 D; n  [/ S2 u
% |! |2 i  N1 h使用20位寄存器，仿真波形如下：
9 D8 M6 f6 f% G8 v& X' Y2 p

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

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

) y) B5 g% I4 m2 {给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010