CORDIC算法的matlab实现

大笨象

不知道有木有人研究这个，比较偏向于硬件描述。
+ w+ L9 W. W2 Y8 Z( j3 u今天先发一个sine cose函数的实现。
# u. C* v2 [8 G$ f3 p有兴趣的可以看一下，m文件见附件。
# S! h( r) v5 Dfunction v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
, C3 D2 g: I$ Q: z5 f0 O9 F# {    if beta < 0
        v=f_cordic(beta + pi,n);
    else
! G  ^% z; m) I" E        v=f_cordic(beta - pi,n);
" a1 N3 C7 C$ u( W: [, C& u' N    end
+ k/ [: q% f" g2 k, ]9 J) @1 T+ V+ k    v = -v;% flip the sign for second or third quadrant
    return
& d1 F3 s8 e2 y) \7 T4 u6 Zend
4 O5 Q4 Q! X9 c5 i0 e. `' K% Initialization of tables of constants used by CORDIC
+ a6 E- D8 L: W% need a table of arctangents of negative powers of two, in radians:
, X0 l4 x! W9 i% angles = atan(2.^-(0:27));
& ^  C& l+ k2 a: u+ g* t5 [angles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
: s6 W, z) H0 E2 z% Y/ `    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
/ p  t- w! ~2 u0 o    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
9 `/ W1 N' T8 Z    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
: @6 C5 p. [9 V( W, ~( ]Kvalues = [ ...
) E( A4 Z1 j# O! d1 I3 H! H    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
3 R9 Q* A3 m: F    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
" ^! Z- n  z- A  L: U9 Q    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
1 v/ O9 Y0 ?, P% D+ f7 v, A( o, C    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));

# u$ L1 X2 S6 S9 ~1 x) r5 d% Initialize loop variables:
$ G: m& D! i% |8 c6 wv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
, \( W  X  A0 C  r% Iterations
for j = 1:n;
    if beta < 0
( _# b+ v/ ^% t        sigma = -1;
; |5 \6 l6 v! A+ }- ~    else
6 ~+ [  ?! i! f2 F) t, w" C4 X        sigma = 1;
* x8 ]. F. Z% \7 ]' j( W# y    end
    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
        angle = angles(j);
    else
        angle = angle/2;
    end
    factor = sigma * poweroftwo;
3 L! n& ~9 D( R4 [8 W" i7 k    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
9 |$ Z' Z8 {+ j/ e8 R" s    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
' \2 j; J/ Y9 N' y8 Vend
% Adjust length of output vector to be [cos(beta), sin(beta)]:
* W, c2 x( |5 }! lv = v * Kn;
return

' L9 T+ A! [! @使用20位寄存器，仿真波形如下：
, b" j" C1 m& ~* `% N1 ~7 M2 c
; x( `$ n  U4 O2 O3 ?

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

6 c/ J3 D0 d! k# m! U) x9 J这个，就是那个。感兴趣的话你还是自己搜搜吧。

大笨象

赛才 发表于 2011-6-9 15:01
$ [! f* o* `5 N) b我表示很给力

) O- |" Z6 a; B0 Y- I; ]. d* O  _给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010