CORDIC算法的matlab实现

大笨象

% N1 L6 e9 f# ?& [
0 Z8 z0 `* [4 w" r+ A不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
7 S4 P$ v- ~; i7 E% ^: t$ W7 j    if beta < 0
) @$ T; h" \, y  r        v=f_cordic(beta + pi,n);
- e9 A; P3 O2 x' l7 N9 V1 N    else
        v=f_cordic(beta - pi,n);
    end
6 G2 `: I" t# \: P    v = -v;% flip the sign for second or third quadrant
" |- p3 k& `: j1 B5 I3 w( L    return
5 \( {# b' P- {9 u" Y. I4 n0 G, rend
% Initialization of tables of constants used by CORDIC
) e+ B) G9 U. E% l5 }% need a table of arctangents of negative powers of two, in radians:
' S4 u8 u5 B8 Z9 I' ^; o$ ?; R% angles = atan(2.^-(0:27));
angles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
/ }9 X. N/ p* p& u/ S2 T+ ~    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
( A( R5 W0 M0 P  y    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
& d: ]5 d, Z2 y    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
6 T7 t, q! Z9 V9 I; X    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
2 }1 h  z& ?5 @/ J+ V3 P! \  l( T2 M% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
0 l4 `: V6 b- b* l    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
1 g' m" X# O' D    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
3 [' m) s" L0 u) ~3 D7 ^9 n3 Y    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
/ Q8 S( T, r6 e; Z# t* G* t0 `/ O    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
: w& {1 R9 X7 L' wKn = Kvalues(min(n, length(Kvalues)));
. l* B- i- V0 h4 |: i6 s6 Y! M
6 l9 _4 M1 J% I) v% Initialize loop variables:
4 P: v/ \( X2 p! |v = [1;0]; % start with 2-vector cosine and sine of zero
% K1 U, O& O3 i# |) ^% g6 }poweroftwo = 1;
6 a8 t1 J3 V# X- i% Iterations
2 S3 h7 u& k& c4 zfor j = 1:n;
0 n* s9 S* q: f) ~! F1 T$ u! G    if beta < 0
/ w8 C  M) q, c+ X& U5 l        sigma = -1;
    else
        sigma = 1;
    end
  [% E( v, Z# {7 f! [    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
        angle = angles(j);
& n7 [3 S. W6 e: R6 m3 H- L6 k    else
. ^9 H0 Y5 W+ \5 d$ Y! G5 W3 `        angle = angle/2;
/ Y% W. I& H- t# N+ H    end
9 `7 ?$ Q- _6 ?0 M! }    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;
end
% Adjust length of output vector to be [cos(beta), sin(beta)]:
7 f0 W  a, m# F$ U& Z& f, t  nv = v * Kn;
return
) |- p/ m, \, M+ Y) a6 N
8 x3 ~' [+ f/ D; @% V使用20位寄存器，仿真波形如下：

1 K- ?# l& M3 _2 E

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
5 Y: O1 j) ~* L8 Q5 [) ]  J+ k这是什么算法？第一次听说

, J$ O+ ~  e5 b7 U( h2 A. O" u这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

赛才 发表于 2011-6-9 15:01
0 J* b- O, @% N9 ^5 H1 Y3 S我表示很给力

2 e5 p- V, C. w: S8 P) G给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010