CORDIC算法的matlab实现

大笨象

0 o+ B8 I& f: T/ f1 e
不知道有木有人研究这个，比较偏向于硬件描述。
1 J( S0 v2 Z: S2 e; L今天先发一个sine cose函数的实现。
1 @: R& K6 g1 _9 m有兴趣的可以看一下，m文件见附件。
$ j# u* L+ X1 {. W4 ?$ rfunction v=f_cordic(beta,n)
7 ?6 m9 w# i1 U# ?if beta < -pi/2 || beta > pi/2
    if beta < 0
" o- e2 v$ j% y6 X% m" D- ?        v=f_cordic(beta + pi,n);
    else
4 m: n. |6 I" {* F" D6 g        v=f_cordic(beta - pi,n);
+ w, j- `, d; _& P    end
    v = -v;% flip the sign for second or third quadrant
    return
end
' @) E1 X' }" o: T" T/ ~% Initialization of tables of constants used by CORDIC
3 k2 V) g( R2 |" {7 F2 r. n* L$ O% need a table of arctangents of negative powers of two, in radians:
: |8 j1 e# @/ g# `4 t% angles = atan(2.^-(0:27));
- N) l+ z" b: @angles =  [  ...
    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
1 t6 z, X5 R! I# C  R    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
! s3 j# ^, Z+ ^) `, r' z1 |( q( u' |    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 ...
, Y7 u, n$ O( m% g3 |* g    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
2 J  e9 x- c( H6 ?    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
: v' E6 C" a4 O( [& M    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
0 }# P' @6 ^' F& pKn = Kvalues(min(n, length(Kvalues)));
2 c1 x' i# R* f+ ^- B/ j
% Initialize loop variables:
4 T# w7 p. O" d, rv = [1;0]; % start with 2-vector cosine and sine of zero
( T) z6 L5 D. t  ^- xpoweroftwo = 1;
% Iterations
for j = 1:n;
    if beta < 0
% T- w* l$ V, A# ?7 @        sigma = -1;
; o4 d, o% N' p- _/ u( O- l3 Q    else
0 q* L, Y5 Z  ]* j        sigma = 1;
    end
    % update the angle from table, or eventually by just dividing by two
& K; |; E: E, }    if j <= length(angles)
4 e% z, H) u  D        angle = angles(j);
" l, R$ [7 c  }    else
        angle = angle/2;
% g" e4 Z" O9 x    end
& F0 B" K0 j/ F    factor = sigma * poweroftwo;
) f5 {. B6 e& o0 A    R = [1, -factor; factor, 1];
    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
    poweroftwo = poweroftwo / 2;
( Z- I" v+ k6 K4 F' dend
' B2 u  N! p2 j# q1 {1 }5 _% Adjust length of output vector to be [cos(beta), sin(beta)]:
' ]* n* o7 o, ]0 e) k* \2 Uv = v * Kn;
+ m( e# ?) [! B& ]4 v5 breturn
3 B, }1 T9 W$ M8 J# W
  c* M! G& U8 j9 [$ u使用20位寄存器，仿真波形如下：
8 |+ h2 t; Y0 q% K$ s
. q) J7 ~* Q6 M: e& F! f, `# [

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

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

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010