CORDIC算法的matlab实现

大笨象

不知道有木有人研究这个，比较偏向于硬件描述。
今天先发一个sine cose函数的实现。
2 @* k. U' i8 Q# }有兴趣的可以看一下，m文件见附件。
function v=f_cordic(beta,n)
# V& t0 ?/ ]" y& g* [if beta < -pi/2 || beta > pi/2
5 Y3 j) K: K. @3 x" X2 N    if beta < 0
        v=f_cordic(beta + pi,n);
    else
        v=f_cordic(beta - pi,n);
* u2 Q  n7 U# j* u4 }/ [; X! A4 U. a    end
/ ^0 i$ ]1 h5 n0 \    v = -v;% flip the sign for second or third quadrant
3 ?" R* `5 L7 x) T5 ]    return
+ {* l# F! j/ F* F( B. jend
% Initialization of tables of constants used by CORDIC
% need a table of arctangents of negative powers of two, in radians:
( h( {( X8 a) J, ]% angles = atan(2.^-(0:27));
: y2 ^9 }4 K' a) ?9 |  \angles =  [  ...
3 b. v3 |' O2 e5 t5 w, ]    0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
+ ]1 K) f  ]' J7 ]    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
# h! ]! H! x9 n) N; [- e    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
: L8 K2 i8 r% N5 b5 G- ]    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
( b/ _* H7 j7 d, l. b    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    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 ...
9 j, n" o; ~: R) p: M) T    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
! R3 r3 k( ^9 ~0 m* _    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
Kn = Kvalues(min(n, length(Kvalues)));
& f% {  n0 M5 d$ V; I5 m  E
' d; {1 P/ l% p4 E# t% Initialize loop variables:
' F* Q$ y) ~1 V2 s: ]v = [1;0]; % start with 2-vector cosine and sine of zero
" C4 z5 Y) h& A) Y- l8 W9 U- xpoweroftwo = 1;
% Iterations
for j = 1:n;
6 C  w3 G* l' f! D; b3 o4 ?- Z    if beta < 0
        sigma = -1;
    else
        sigma = 1;
    end
    % update the angle from table, or eventually by just dividing by two
; r( u" W' X8 ~" x$ k2 M7 o( R    if j <= length(angles)
        angle = angles(j);
    else
3 i) K9 m+ m3 C  |) t" X        angle = angle/2;
    end
    factor = sigma * poweroftwo;
" P$ F! A$ {6 z    R = [1, -factor; factor, 1];
( Q, F' i. P3 h    v = R * v; % 2-by-2 mtrix multiply
( h% a% g1 f, r/ z& U! c/ I    beta = beta - sigma * angle; % update the remaining angle
& d1 `6 ~, L7 _& g    poweroftwo = poweroftwo / 2;
0 ?: E9 T( M2 `end
5 p% M! v0 E' T! h% Adjust length of output vector to be [cos(beta), sin(beta)]:
v = v * Kn;
return
9 `/ D! p/ K6 o' L1 I3 Q& _
使用20位寄存器，仿真波形如下：
# L+ k( J9 P7 y, y5 n$ h

杨帆

我表示学习了

zhangjichang

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

赛才

我表示很给力

大笨象

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

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

大笨象

赛才 发表于 2011-6-9 15:01
  T# I: a" |0 L& e, `. N我表示很给力

给力就好。呵呵

魅影骑士

顶一下！

yingzhen

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

羅雲琦

有什麽用？求解

jt202010