第4章cordic算法原理及实现.pptVIP

第4章 CORDIC算法原理及实现 何宾 2009.09 本章概述 本章介绍了CORDIC算法的原理，着重介绍了的三个 坐标系及其两种模式下的CORDIC的实现原理，以及迭 代的算法的实现方法。 在此基础上，详细介绍了CORDIC算法在FPGA上的 实现方法及实现过程，并对其性能进行了详细的讨论。 CORDIC简介 坐标旋转数字计算机(Coordinate Rotation Digital Computer，CORDIC)算法可以追溯到1957年由J.Volder发 表的一篇文章。 在上个世纪五十年代，在大型实际的计算机中的实 行移位相加受到了当时技术上的限制，所以使用CORDIC 变得非常必要。到了七十年代，惠普公司和其他公司生产 了手持计算器，许多计算器使用一个内部CORDIC单元来 计算所有的三角函数(那时求一个角度的正切值需要延迟 大约1秒中)。 CORDIC简介 二十世纪八十年代，随着高速度乘法器与带有大存 储量的通用处理器的出现，CORDIC算法变得无关紧要了。 然而在二十一世纪的今天，对于FPGA来说， CORDIC一定是在数字信号处理应用中(比如：多输入多输 出（MIMO），波束形成以及其他自适应系统)计算三角函 数的备选技术。 CORDIC算法原理--圆坐标系旋转原理 如图4.1，在xy坐标平面内将点(x1,y1)旋转θ角度到点 (x2,y2)。其关系用下式表示： （4.1） 这被称为是平面旋转、向量旋转或者线性(矩阵)代数 中的Givens旋转。上面的方程组同样可写成矩阵向量形式： CORDIC算法原理--圆坐标系旋转原理 上面的方程组同样可写成矩阵向量形式： CORDIC算法原理--圆坐标系旋转原理 如果去除项cosθ，得到伪旋转方程式: （4.5） 如图4.2所示，即旋转的角度是正确的，但是x 与y 的 值增加倍cos`θ(由于cos`θ1，所以模值变大)。 CORDIC算法原理--圆坐标系旋转原理 CORDIC算法原理--圆坐标系旋转原理 并不能通过适当的数学方法去除cosθ项 ,然而随后 发现去除cosθ项可以简化坐标平面旋转的计算操作。 CORDIC方法的核心是伪旋转角度，其中tanθ=2-i。 故方程可表示为： （4.6） 表4.1给出用CORDIC算法中每个迭代(i)的旋转角度 (精确到9位小数) CORDIC算法原理--圆坐标系旋转原理 表4.1用CORDIC算法中每个迭代(i)的旋转角度 CORDIC算法原理--圆坐标系旋转原理 在这里，把变换改成了迭代算法。将各种可能的旋转 角度加以限制，使得对任意角度的旋转能够通过一系列连 续小角度的旋转迭代i来完成。旋转角度遵循法 则： ，遵循这样的法则，乘以正切项变成了移位 操作。 前几次迭代的形式为：第1次迭代旋转45o，第2次迭代 旋转26.6o，第3次迭代旋转14o等。 CORDIC算法原理--圆坐标系旋转原理 很显然，每次旋转的方向都影响到最终要旋转的累积 角度。在–99.7≤θ≤99.7的范围内的任意角度都可以旋 转。满足法则的所有角度的总和为99.7。对于该范围之外 的角度，可使用三角恒等式转化成该范围内的角度。当 然，角分辨率的数据位数与最终的精度有关。 cos45xcos26.5xcos14.03xcos7.125…xcos0.0139 =0 1?0.607252941=1.6467602。因此，在13次旋转以 后，为