CORDIC算法的matlab实现

大笨象

( Y: E) G+ B& N+ Q/ W不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
. w" {4 Y$ U* p$ |& m1 nfunction v=f_cordic(beta,n)
+ s0 C+ ~1 N* ~: m2 F' d" Yif beta < -pi/2 || beta > pi/2
" ~& U0 `" r( b* t( a    if beta < 0
' M  O' {9 G4 b) v% ~+ \        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
$ x1 s# E, W# c! a' H. R  n    return
end
% Initialization of tables of constants used by CORDIC
% need a table of arctangents of negative powers of two, in radians:
' \- O: w7 z% t$ b) |5 E: f% angles = atan(2.^-(0:27));
angles =  [  ...
/ w" I% B  Q7 M/ |- s4 Z    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
) U, A8 F4 Y9 K0 t1 d    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
3 B* l5 q3 i, |3 \/ K    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
/ [$ w: f# e5 p5 U    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
* u& f) `5 @; L9 ?0 P    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
0 ^# M( O* U- X. t- m: `6 S: C% q% 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 ...
. N$ [3 ?5 f2 a( g  w    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
: i5 C% P8 }" T    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
" V0 U7 z2 `" {% ]' }9 m    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
0 X1 _! E/ w) f8 lKn = Kvalues(min(n, length(Kvalues)));
' M. {' [: |1 l4 u# {* j6 P* r
8 G) [3 g# j" x% W9 V% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
  O# C$ a4 ]  _poweroftwo = 1;
% Iterations
& B1 U8 j% \6 F' ~! `for j = 1:n;
, u$ F. A& ]& c% H6 f4 u    if beta < 0
        sigma = -1;
    else
        sigma = 1;
    end
* q4 U. G9 {8 M" s# W- z    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
# T+ g8 i: E0 m        angle = angles(j);
    else
: d! c* S4 G7 |9 T6 d) S9 K5 q" K        angle = angle/2;
    end
1 j) T6 p" K3 e( K8 j$ F    factor = sigma * poweroftwo;
5 C/ h4 d/ t. p1 t  g) N, `2 o    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
& z' j1 l, M; _" a    beta = beta - sigma * angle; % update the remaining angle
" l: s: X8 n. A    poweroftwo = poweroftwo / 2;
end
% Adjust length of output vector to be [cos(beta), sin(beta)]:
: K2 H/ ^( Y0 W. ~( Wv = v * Kn;
1 q4 J0 i5 ^# a, Creturn
0 s$ f) V& u* k, Z' m0 ^
! _  R; d* N( v使用20位寄存器，仿真波形如下：
2 [) o: N' O9 W0 h5 o) {

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

赛才 发表于 2011-6-9 15:01
# U* f: E$ u# A  X. F我表示很给力

: `2 ]& N! n7 r给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010