Cordic 算法的原理介绍

转载于http://www.cnblogs.com/touchblue/p/3535968.html

cordic 算法知道正弦和余弦值，求反正切，即角度。

采用用不断的旋转求出对应的正弦余弦值，是一种近似求解发。

旋转的角度很讲求，每次旋转的角度必须使得 正切值近似等于 1/(2^N)。旋转的目的是让Y轴趋近与0。把每次旋转的角度累加，即得到旋转的角度和即为正切值。

比如Y轴旋转45度，则值减小1/2;

   再旋转26.56505°，再减少1/4;

   再旋转角度14.03624º，再减少1/8; 依次减少1/16， 1/32......，最后Y轴的值无限小，趋近于0 。

比如X=1， Y=1，的角度，角度是45°。经过一次旋转，要使得Y=0，这个角度必须是45°。

 

 

 

 

 

 

如上图

如图中，直角坐标系中点（X0，Y0）逆时钟旋转角度θ，变换成坐标（X1，Y1）,那么用X0，Y0,以及θ的三角函数，如果表示X1，Y1呢？

请想象，如果坐标也旋转角度θ，那么X1，Y1的坐标依然是(X0，Y0)。接着往下看：

看完以上这副图，就该明白这个等式了：

x1= x0cos θ − y0sin θ
y1= x0sin θ + y0cos θ

再把这个式子化成正切函数。

 

 

Cordic 算法的思想是通过迭代的方法，不断的旋转特定的角度(这个特定的角度就是使得 Y为上次的1/2)，使得累计旋转的角度的和无限接近某一设定的角度，

每次旋转的角度的θ = arctan( 1/(2^n) )；

具体迭代如下表：Z0 =30°，Y0=0，X0 = 0.6073

输入30°，经过9次迭代后， Z0 = 0，Y0=0.5006， X0 = 0.8657

x '(i +1)= ( x 'i− y 'i（σi）2−i)
y '(i +1)= ( x 'i（σi）2−i+ y 'i)

(当 i = 0)
x '1= 0.607 − 0 ⋅ (+1 ) ⋅ 1 = 0.607

y '1= 0.607 ⋅ (+1 ) ⋅ 1 + 0 = 0.607

 通过Cordic算法后，得到y9=0.5006 (=sin(30°))

                                 x9=0.8657 (=cos(30°))

所以也可以用cordic算法求出正切值的。

 

或者求反正切值：

计算公式：