基于CORDIC算法的正余弦信号发生器

（1）CORDIC算法的圆周系统之旋转模式原理简介

CORDIC算法是一种兼容速度快、精度高和方便硬件资源实现的综合方法。它使用多次预定角度的不断旋转来逼近初值角度的近似逼近方法。CORDIC算法有3种旋转系统：圆周系统、线性系统和双曲系统。每种系统下都有两种工作模式：旋转模式和向量模式。其每个系统每种模式下的应用情况如下：
表一：CORDIC算法的3种系统和2种不同模式应用情况

	旋转模式	向量模式
圆周系统	计算sin、cos	计算tan-1
线性系统	计算乘法	计算除法
双曲系统	计算cosh、sinh	计算tanh-1

CORDIC算法圆周系统的旋转模式，可以实现正弦函数和余弦函数的运算。相应的迭代运算如式所示：

其中,

经过n(n趋近于无穷)次旋转，得到的最终结果为

由上式可知，当初始化角度x0=1/An，y0=0时，其中An=1.64676可得出：

按上述步骤迭代，即可在x，y的输出值中得到输入值z0的正弦余弦值。但该cordic算法的局限性在于输入值度数范围为[-99.8829,99.8829]，但在设计的场合中需要目标旋转角度可覆盖整个周期，则需要在本cordic算法的基础上按三角函数的对称性进行预处理。

（2）旋转模式下CORDIC算法的FPGA实现

下图为设计的CORDIC算法的输入输出端口，i_valid信号为输入有效信号，当输入信号有效后，将该信号拉高一个时钟周期后，模块开始对zin中的数据计算正弦余弦值；x_cosz，y_sinz分别表示向量的余弦值和正弦值；o_valid信号则是当输出有效后，拉高一个时钟周期。完整CORDIC算法代码见评论区链接。

（3） CORDIC算法误差分析

在上述的cordic算法的输入中，zi的输入位宽都是16bit，为了分析cordic算法计算误差情况，对cordic计算出的正弦余弦值与实际的值作对比，误差散点图如下图所示：

由上图可知，运用cordic算法计算出的结果中，正弦余弦值误差范围均在10的-4次方以内，该CORDIC算法计算精度达到要求。

（4）基于CORDIC算法的DDS实现

顶层文件的输入输出端口如下图所示：

其中clk为运算的主频时钟，必须采样100Mhz的时钟输入；ad_flag表示adc采集信号有效的标志位，即每在ad_flag的上升沿，锁相放大器就读取一次采集的数据，保证了参考信号的采样频率与ADC信号采样频率一致；16bit的angular_freguency输入信号为要生成的正弦余弦信号的角频率，其小数位位宽为12bit；输出的正弦余弦值均为18bit 小数位位宽16bit。
modelsim仿真分析图如下：