cordic算法是一种易于硬件实现的求基本函数的算法,将复杂运算转化为加法、减法和移位操作的组合,避免硬件上不易实现的复杂运算,比如除法,开方等,另外可以避免查找表方法带来的占用较多存储资源的问题。关于cordic的原理很多资料已经介绍的比较详细,这儿不会详细推导,这里只关注如何应用。
一、三角函数运算
1、已知角度值,计算sin、cos
********************代码*****************
import numpy as np
def cordic_for_sincos(alpha,coridc_num):
x0=0.607253
y0=0
z0=alpha
xlast=x0
ylast=y0
zlast=z0
for i in
range(coridc_num):
sign=-1
if zlast>=0:
sign=1
xnow=xlast-ylast*(1/2)**i*sign
ynow=ylast+xlast*(1/2)**i*sign
znow=zlast-np.arctan((1/2)**i)*sign
xlast=xnow
ylast=ynow
zlast=znow
result_cos=xlast
result_sin=ylast
return
result_cos,result_sin
********************代码*****************
算法的大致原理是:在单位圆上,从x轴开始旋转,旋转到角度alpha,此时横坐标对应余弦,纵坐标对应正弦。
仿真代码很简单,其中反正切值根据迭代次数i提前算出写入表格,x0是迭代初值,这个初值是cordic增益的倒数,cordic增益根据迭代次数唯一确定。迭代次数越多,则计算越精确,需要的中间位宽也越多,当然cordic增益值也不同。
2、已知sin、cos值,