CORDIC算法的matlab实现

大笨象

; |& f1 N# }+ ?, X
不知道有木有人研究这个，比较偏向于硬件描述。
# l5 N( _) f  [$ g今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
6 g( D' z4 @6 d- Z" M* C9 L- N    if beta < 0
        v=f_cordic(beta + pi,n);
" Y7 H) G2 d8 W5 O    else
- l) r/ {+ l% H! f        v=f_cordic(beta - pi,n);
8 {" c, w& s+ M# N. [    end
    v = -v;% flip the sign for second or third quadrant
    return
end
: R3 O  y. O: R- E" T/ k% Initialization of tables of constants used by CORDIC
" N- G- T0 y* {% need a table of arctangents of negative powers of two, in radians:
5 Y& Q, }3 Q2 x% D2 ~% angles = atan(2.^-(0:27));
$ E- y1 M6 Z; O) d6 G. w5 langles =  [  ...
7 w* r, e4 F! g& `    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
4 J. b5 Z. X: }9 \7 Z    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
0 `/ `3 t# y9 e( D) G    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
4 Q! u& V( M1 r+ A6 K" V) I7 k    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
7 E2 y9 v" A; R$ M4 v6 _    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
/ x& I( k5 n& d1 v3 }5 Z% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
2 y$ b4 U; {; y# u4 Y  E6 rKvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
4 o# R* v: Y+ j8 m- i    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
- U0 [7 {/ H6 T7 g' d! K; o    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
$ E. q3 Z' O8 i8 Z4 U$ A    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));
3 }& U2 `2 O) S& M. c/ h5 ~
% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
% Iterations
- X( t8 j& x9 Q7 g- \" ?' @; afor j = 1:n;
    if beta < 0
        sigma = -1;
# w, N* c: D& ]- }) D    else
, Z1 K- D# w( H  p& s8 j. @        sigma = 1;
* a7 a* S" v5 }7 |7 M- q7 {6 L    end
+ q! m* {0 V$ z$ t, x- e    % update the angle from table, or eventually by just dividing by two
4 u' x# J9 [# z+ S0 @& {: W+ z    if j <= length(angles)
$ b& j' j& X" \/ R. D* k" n        angle = angles(j);
0 T6 J+ b% j6 v0 `3 n' k- e0 g6 J    else
        angle = angle/2;
    end
    factor = sigma * poweroftwo;
4 k& S) _/ y! {6 ~8 e7 w4 |& N    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
8 o1 b; i  T2 x' h9 v; r7 C" w    beta = beta - sigma * angle; % update the remaining angle
5 k4 u, e7 ]$ U$ Q% ?6 w) d    poweroftwo = poweroftwo / 2;
/ h% \3 I& E7 c* K* Bend
6 H* y6 Z, i5 z3 k* b% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
return

/ ~3 x4 M' m; O- P/ v+ Y使用20位寄存器，仿真波形如下：
( ?6 N2 f/ G; i7 P8 R

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
6 H& f: s7 `; \% S# B- J这是什么算法？第一次听说

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

大笨象

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

7 j0 A. C# |5 Z. ]6 x' n" m给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010