CORDIC算法的matlab实现

大笨象

4 E3 a; i, n! j6 q不知道有木有人研究这个，比较偏向于硬件描述。
! W6 A" ]7 ~/ m  t& J$ J/ D今天先发一个sine cose函数的实现。
' b5 G0 T  A1 J0 ]! J" i有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
! n- F  c$ u8 Y7 Jif beta < -pi/2 || beta > pi/2
2 t$ z1 ]4 R. E( \: l; g% q$ w& A: G    if beta < 0
- D4 X& h' `" v& a6 V: T        v=f_cordic(beta + pi,n);
# R4 o' V! j9 w% X: y5 X4 }; N/ v( j9 u    else
# x* k3 a9 W* f0 J, e, v        v=f_cordic(beta - pi,n);
" d1 }& a( u0 k5 K( f    end
    v = -v;% flip the sign for second or third quadrant
9 `# }5 C* A3 J, r  [  o  ]7 @    return
end
% Initialization of tables of constants used by CORDIC
0 s; Y9 T" l) J: `0 D! L% need a table of arctangents of negative powers of two, in radians:
0 R" ]. e0 C4 i. ~' i+ J1 x% angles = atan(2.^-(0:27));
+ b# n' F9 |. z9 bangles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
+ u& I! L+ w( g5 y+ `; q    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
" L. O+ Q0 S; F+ c# j    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    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 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
7 q) N& \9 y7 ]    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
! j  v' N- O$ x9 D* p: G/ }    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
' x) Z& d$ h$ L9 k1 S0 ?    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
/ D& m& V1 N4 wKn = Kvalues(min(n, length(Kvalues)));

2 E' y$ X0 R9 l3 x% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
( K5 ^1 B" ^7 y) P- W' I2 J0 s7 U% Iterations
+ o5 T- K! e; u% O# d5 \5 z. @" N$ Nfor j = 1:n;
; n' n  b! j0 g* T8 S    if beta < 0
        sigma = -1;
    else
. `9 U3 p, _  j0 \( Q2 S3 m        sigma = 1;
    end
6 I% E% W: B- x' g1 x1 [    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
% ?8 d7 C& h7 p8 F, U        angle = angles(j);
0 V5 X/ w6 T& z; n" h    else
8 k% g# A9 [# B) q* D* C        angle = angle/2;
1 l& g8 q1 z. m/ z1 C+ m  D    end
" s% N3 k, f1 z& O/ Z( o( d    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
6 I6 K) v/ |  }6 O    beta = beta - sigma * angle; % update the remaining angle
" w0 A2 E  M* H    poweroftwo = poweroftwo / 2;
end
2 W$ D' c: z! p3 d, Q% G! M3 h1 R% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
% U) d1 N' o' h4 @2 Y- Preturn

* p4 |* U; m9 [/ o) _- P* Q( g  ?使用20位寄存器，仿真波形如下：

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

赛才 发表于 2011-6-9 15:01
7 ~* I% i) F9 F4 b9 l我表示很给力

* `, _, F& s! M给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010