研究了基于CORDIC算法的反正弦函数运算器的FPGA实现.DOC

快速精确的反正弦函数运算在现代工程中应用广泛。为了提高反正弦函数的精度和计算能力， 研究了基于CORD IC算法的反正弦函数运算器的FPGA 实现， 并通过改进算法减小了误差， 使误差精度达到10--4 数量级。并在X ilinx 的XC5VSX50T芯片上验证， 结果表明该运算器的数据处理速度快， 精度较高， 适用于高速大数据量的数据信号处理领域。 在数字信号处理领域， 高速高精度的反正弦函数发生器有着广泛的应用。目前在FPGA 上实现反正弦函数运算器采用的主要方法是查找表法和泰勒公式展开法。查找表法所需要的存储单元随着精度的增加或输入值范围的增大而成指数增加； 泰勒公式展开法将函数简化成一系列的乘法和加法运算，但是在FPGA 上实现乘法运算既耗时又占用大量资源。本文设计了基于CORDIC算法的反正切函数计算模块，并且根据IEEE-754单精度浮点数据格式对输入输出数据进行处理，实现了高精度的浮点反正切函数的计算。 1 反正切函数实现原理 CORDIC（Coordinate Rotation Digital Computer）算法即坐标旋转数字计算方法，是J.D.Volder1于1959年首次提出，主要用于三角函数、双曲线、指数、对数的计算。该算法通过基本的加和移位运算代替乘法运算，使得矢量的旋转和定向的计算不再需要三角函数、乘法、开方、反三角、指数等函数。 CORDIC算法有旋转模式和向量模式两种计算模式。旋转模式可以用来计算一个输入角的正弦、余弦，向量模式可以计算给定向量的角度和长度。 CORDIC算法的基本迭代公式为：   从上式可以看出，CORDIC算法在向量模式可以计算出给定向量（X,Y）的长度和角度，即从平面坐标到极坐标的变换。 2 数据格式转换接口模块 本文设计的反正切函数硬件模块输入为IEEE-754单精度浮点数据，而模块内部迭代使用的是定点整型数据，因此需要进行转换。 在图1的输入数据转换接口示意图中，X、Y为输入的IEEE-754浮点数据格式，输入范围是（-∞，+∞），经过接口模块转换为整型定点数据Xn、Yn,其表示范围是[-1 +1]. 矢量（X,Y）在平面坐标系中的角度为arctan（Y/X），它只与Y和X的比值有关，与Y和X的实际长度无关。用这个性质可以以X和Y中绝对值最大的值作为归一化数值，将X和Y的范围重新映射在[-1 +1]之间，实现（X,Y）到（Xn,Yn）的转换。 计算结束后输出结果Z.Z是32位定点整型数据，且232被定义为2π，将其规格化为IEEE-754格式的过程如图2所示。在对Z进行规格化之前，需要进行前导零检测，以确定规格化时尾数左移的位数和指数位的大小，前导零的检测硬件上可以用casex语句实现。 3 整体设计以及仿真综合 浮点反正切函数的硬件结构包括了三个主要部分，即浮点数据格式转换接口模块、CORDIC内核计算模块和浮点输出数据转换接口模块，如图3所示。 采用Quartus II对设计进行FPGA综合，FPGA芯片选择EP2C70F896C6,硬件环境为Altera公司的DE2-70平台，总共需要1 522个逻辑单元，占用芯片资源的2%,最高工作频率为100 MHz. 4 Nios II中反正切函数的自定义指令实现 反正切函数与Nios II CPU的接口采用multi-cycle custom instruction,dataa和datab为输入数据，result为结果输出，要从C语言中直接调用自定义指令，需要一个宏定义接口。可以在system.h文件中找到自定义指令的宏定义，如： #define ALT_CI_CORDIC_ATAN2_N 0 #define ALT_CI_CORDIC_ATAN2（A,B） __builtin_custom_inii （ALT_CI_CORDIC_ATAN2_N,（A），（B）） 为了正确调用自定义指令，在主程序中重新做以下宏定义： #define ATAN2（A,B） __builtin_custom_fnff（ALT_CI_CORDIC_ATAN2_N,（A），（B）） 与system.h文件中系统自动生成的宏定义不同之处在于将宏定义的接口说明由__builtin_custom_inii改为了__builtin_custom_fnff.系统自动生成的宏定义默认输入输出皆为整型数据，改为__builtin_custom_fnff就是通知系统这是一条输入输出都是单精度浮点数据类型的用户自定义指令。这样由CPU调用时就不会出现数据类型