CORDIC算法的matlab实现

大笨象

4 X, M# [' Y5 K) F  o. C( h
; V9 h% s8 p- |, @. W  U不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
if beta < -pi/2 || beta > pi/2
, d8 ]' D' [$ C' m3 z1 B- {0 t    if beta < 0
1 B2 A5 |& \) {2 u: Y2 i2 `# i        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
    end
! N* r. a" r% g  J8 `9 ^9 _; I$ I; [    v = -v;% flip the sign for second or third quadrant
* w# C7 X& q6 L1 Z/ [2 N$ I    return
end
8 h$ w6 Z* {: Y0 }% Initialization of tables of constants used by CORDIC
# N# S: g2 d1 O% need a table of arctangents of negative powers of two, in radians:
% angles = atan(2.^-(0:27));
angles =  [  ...
. x/ G8 m2 Q6 N  J. S3 o% p    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
* G) |3 B! K: D/ p) ^4 ], T+ X    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
) w$ \9 I* X  r5 ]8 L6 {1 v: T' T: Y! ]1 u    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
: l9 s: W5 z: w+ ~* G    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
3 r7 I" _4 z0 I( z7 x, ?+ s& L( \    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
/ R$ C+ U" [* n& P    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
) P% U+ E( ^; k+ Q6 E    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
  @9 r) r7 e; Y3 D+ H8 S& p% and a table of products of reciprocal lengths of vectors [1, 2^-j]:
Kvalues = [ ...
    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
8 ], h! X# C/ v& X7 x: {    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 ...
    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
* E$ E3 u0 T+ U6 VKn = Kvalues(min(n, length(Kvalues)));

0 Z. k& \7 _3 Z$ y6 A% Initialize loop variables:
v = [1;0]; % start with 2-vector cosine and sine of zero
  s' r$ S. {5 j6 Npoweroftwo = 1;
7 x$ l5 ~" m7 ?- K' I, ~4 F% Iterations
for j = 1:n;
# p1 P: R: I5 g5 ^" H, G    if beta < 0
; H' h) n4 D2 M        sigma = -1;
    else
        sigma = 1;
4 d( J1 C$ u4 K    end
0 p* V4 J3 |4 S9 w  A    % update the angle from table, or eventually by just dividing by two
    if j <= length(angles)
' U; k5 L3 P- c" M        angle = angles(j);
    else
        angle = angle/2;
    end
: g/ R! M: B% {6 p    factor = sigma * poweroftwo;
) ^& [, M1 ~0 x    R = [1, -factor; factor, 1];
" n9 v5 ^2 p3 z    v = R * v; % 2-by-2 mtrix multiply
    beta = beta - sigma * angle; % update the remaining angle
8 ?% w' Z5 g$ m. N* A" ~    poweroftwo = poweroftwo / 2;
% r( s+ p# _" X% z# b" ]end
% Adjust length of output vector to be [cos(beta), sin(beta)]:
/ x. ~( S$ m7 S% Dv = v * Kn;
return
+ \; D4 X$ S" l2 k& }
3 u$ @+ ?7 j0 \! v& G- r" c& t- q使用20位寄存器，仿真波形如下：
8 A; F) e* u( E, K

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

zhangjichang 发表于 2011-6-8 22:31
: N% O- [" g+ ^7 F8 s+ p这是什么算法？第一次听说

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

大笨象

赛才 发表于 2011-6-9 15:01
4 x& s8 N- c0 `+ s! j: d( k8 Y我表示很给力

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010