基于流水线结构的cordic算法的实现汇编.docxVIP

基于流水线结构的cordic算法的实现摘要系统在处理数据的时候，一个指令周期含有多个时钟脉冲，每个脉冲周期由不同的部件完成不同的操作。非流水线结构是指一个指令周期完成以后再接受下一条处理数据的指令；而流水线结构，每个时钟脉冲都接受下一条处理数据的指令，只是不同的部件做不同的事情，就象生产线流水操作一样，并不是等一个或一批产品做完，再接受下一批生产命令，而是每个工序完成以后，立即接受下一批生产任务。这样提高了系统处理数据的速度。随着超大规模集成电路(Very Large Scale Integrated circuits , VLSI)技术的飞速发展,经常需要用硬件快速和精确地进行三角函数值的计算,而坐标旋转算法(Coordinate Rotational Digital Computer, CORDIC)能够将多种难以用硬件电路直接实现的复杂的三角函数运算分解为统一的加减、移位操作,极大地降低了硬件设计的复杂性。在现代信号处理中,经常会遇到三角函数、超越函数和坐标转化等问题。传统的实现方法有查找表多项式展开等方法，但是这些方法在精度、速度、简单性和效率方面往往不能兼顾。CORDIC 算法则可以很好地兼顾这几方面的要求。 CORDIC算法只使用移位和加减运算，硬件实现简单。使用流水线结构，每级CORDIC使用独立的单元，这样使运算速度非常快。当流水线填满之后，每个时钟周期就会得到一个结果。从CORDIC算法的基本原理出发，讨论其工作过程以及旋转角的覆盖范围，在此基础上，给出具有流水线结构的FPGA实现方法以及增益因子的大小与流水线级数的确定关系，给出了verilog实现算法，在Quartus6.0调试与仿真，验证采用FPGA实现的CORDIC算法的有效性。关键词： CORDIC算法 旋转 迭代 verilog语言 流水线CORDIC算法原理CORDIC，即坐标旋转算法（Coordinate Rotations Digital Computr），是由Voider J 等人于1959 年在设计美国航空导航控制系统的过程中提出来的算法。其基本思想是通过一系列固定的、与运算基数相关的角度不断偏摆以逼近所需的旋转角度。由于这些固定的角度与计算基数有关，运算只有移位和加减。这在FPGA实现上就有了很大的优势，可以大大的节约硬件资源。CORDIC算法包括圆周系统、线性系统、双曲系统，利用三种系统可以实现多种超越函数，再利用函数之间的关系就可以实现更广泛的函数运算。该算法可以把某个坐标在直角坐标系(x, y)和极坐标系之间进行坐标变换。1971 年Walter J提出了统一CORDIC 算法，引入参数m 将CORDIC 算法实现的三种迭代模式圆周、线性和双曲线变换都包括了进来。统一于一个表达式下，形成目前所用到的CORDIC 算法最基本的数学基础。如图1-1，初始角为向量旋转角度后得到新的向量，根据简单的几何关系有：经过整理式（1）和式（2），可以得到：从上述分析可知，假设向量经过i 次旋转，每次旋转的角度记作,每一次旋转的角度方向记作(?1或+1) ,这样经过n次旋转后，总共转过的角度为我们可将去除，就得到了伪旋转公式如下：即旋转所得到的（y,x)较真实的值变大倍，但是这样可以简化cordic过程。正是通过每次旋转的角度的累积来接近所需要旋转的角度，这样就需要一个因子来判断下一步的旋转向什么方向，即应该加上还是应该减去新旋转的角度。所以，上述旋转方程应该添加一个旋转因子，即：在这里，引入一个角度累加器，来判别旋转因子下面写出经过推广的CORDIC 算法：其中m =1对应图（1）的圆周旋转运算，而m = ?1或0分别对应双曲旋转运算和线性旋转运算，表示旋转的方向，=-1表示顺时针旋转，=1表示逆时针旋转。为了表示方便，对式（5）中的做强制约束，使其等于，这样带来的好处是就转换成了硬件中的移位运算。式（6）和式（5）等价于以下的迭代公式：使结果的旋转称为旋转模式，使的旋转称为向量模式，旋转过程会带来系数的变换，为了校正，引入校正系数K。经过计算，当n趋近于无穷大时，，，最终可以通过乘以1/K来矫正旋转得到的数据。根据推导，式（9）的n次迭代可以得到如下结果（m=1，旋转模式）给定初始条件，迭代结果为所以，将所需产生的角度值作为输入，迭代结果输出和就是需要的三角函数值。采用的迭代方程组为式（9）。设计分析与实验过程因为角度测定范围要覆盖，确定数据格式为 ***.****,****，****，*共计16位，最高位是符号位，而后是两个整数位，低13位为小数位。按照上面提及的迭代序列：0,1,2，…n-1，所覆盖的角度只有，而一般都要求覆盖，所以我们采用增加迭代的方法来扩大角度覆盖范围，即增加两个i = 0的迭代，则移位序列为0