CORDIC算法的matlab实现

大笨象

6 J# \4 {  V: ?1 ^( h
不知道有木有人研究这个，比较偏向于硬件描述。
+ ~$ |+ w5 N# l0 w今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
+ g8 U1 T; V% H' `# c5 {1 ^function v=f_cordic(beta,n)
9 G# N+ X9 `* e9 Mif beta < -pi/2 || beta > pi/2
    if beta < 0
        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
    end
    v = -v;% flip the sign for second or third quadrant
' x" B# f  R& r4 c1 ?4 |    return
end
; o0 q0 m+ @% U5 m' A" t" H% Initialization of tables of constants used by CORDIC
% N; e( m; {' T1 ^% ?! m3 k  M/ y/ z% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
angles =  [  ...
, Z# Q* j) F4 Q# W8 u    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
$ X8 [/ [: s3 r% L    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
! k0 @. }9 Q! ^8 ~  D    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
$ j- A  u" P7 |4 A: c; P    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
9 {$ X' e+ {" Z$ h. t* L8 R    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 ...
    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
  w5 l* `! L- [; A8 s$ P    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));
, X9 M4 I7 f/ ~- I+ a! y% c
% Initialize loop variables:
7 o/ E$ Y& q7 p  g/ w/ Zv = [1;0]; % start with 2-vector cosine and sine of zero
poweroftwo = 1;
% Iterations
0 O2 v9 ?- K# L9 gfor j = 1:n;
+ b; `# W$ [) P0 H7 M, D% z, X    if beta < 0
8 N! Y% X, b$ x  \: j" U" z$ Q- q. p        sigma = -1;
( a$ ?$ J$ j: l+ j2 i    else
8 M0 p  y* H& m" y        sigma = 1;
* X, m* }0 d% p* {: S4 ?7 }    end
    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
4 b3 I  c5 I! l/ n: g: N9 P; T! g        angle = angles(j);
    else
! t) C; L& j- Z; F) s1 @        angle = angle/2;
" L, u3 c/ D/ o- c* d6 U, ]( B2 C    end
    factor = sigma * poweroftwo;
    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
" J# M; e# y6 V8 c, d& x    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
9 ^  E4 k, E- i6 E1 L( s6 nend
( B% y6 ^( c8 @3 f% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
4 n/ h- R  O, G/ mreturn

使用20位寄存器，仿真波形如下：
* ]; _3 D" }5 D

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
% }) k) b4 Y; c8 `' Q/ x# R这是什么算法？第一次听说

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

大笨象

赛才 发表于 2011-6-9 15:01
3 {. ?7 c* j+ C$ L我表示很给力

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010