Cordic 算法的原理介绍

admin 2017-6-24 3960


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算法求出正切值的。

 

或者求反正切值:

计算公式:

 


最新回复 [0]
返回