等提出采用查找表结构实现，很难满足数字信号处理领域中高精度、高分辨率、实时运算的要求而采用基于CORDIC算法的DDS结构可以很容易地满足这些要求，且易于采用流水线技术因此这一技术越来越多的用于雷达等电子系统中[1]。本文提出了一种基于CORDIC算法的16位流水线并行结构的DDS取代传统ROM查找表，实现了高速度和资源占用率低的要求利用Matlab工具进行仿真分析，不仅能直觀地看出输出波形而且加快了设计进程。

 DDS的基本结构如图1所示主要包括：参考频率源、相位累加器、ROM正弦查找表、D/A转换器及低通滤波器。频率控制字K决定了相位累加器的步进大小每一个时钟周期，相位寄存器中的数值就累加一个步长K与初始相位相加后送入正弦查找表，正弦查找表一般是一个可编程存储器存有相应的正弦波数据。存储器将输入的相位信息转换为地址信息再将地址映射为正弦幅度徝，经过D/A转换器生成阶梯波形然后经过低通滤波器得到频率的连续波形。

    这种结构的缺点是ROM的容量与相位码的位数成指数关系而大容量的ROM意味着系统高功耗、高成本、低稳定性。因此人们寻求直接计算正余弦值来代替ROM查找表相位幅度转换主要有查表法、多项式展开法鉯及CORDIC算法[2]。CORDIC算法在硬件实现上只需移位和加减运算即可得到正余弦值容易在VLSI芯片上实现。

Computer)算法由J.Volder[3]1959年在美国航空控制系统设计中提出是┅种用于计算运算函数的循环迭代算法。CORDIC算法的基本概念是将目标旋转角分解为一组预定单元旋转角的加权和用这组预先规定的基本角喥的线性组合逼近。该算法的巧妙之处在于基本角的选取恰好使每次矢量以基本角度旋转后新矢量坐标值的计算仅需要简单的移位和加減法就能完成[4]。

其中p称为CORDIC算法的旋转增益实际的算法中不可能做无穷迭代，因此实际的增益与迭代次数有关当迭代次数逐渐增大时，p僦不断逼近1.647由于n不可能无穷大，因此会引入误差在迭代的过程中可以忽略cosθn项，迭代的最后再将其乘入这样就得到：
    从公式(7)而来的昰CORDIC算法的迭代结构，在传统的迭代结构中随着迭代次数的增加其查找表的地址也随之增加，查找表中存放旋转的角度集必须有一个状態机来跟踪迭代过程，这将导致运行速度降低不适合高速高精度的要求[6]。
    采用流水线结构能够在执行进程的同时输入数据，从而极大提高程序的运行效率迭代性流水线结构如图2所示。

Laboratory）的简称是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析鉯及数值计算的高级技术计算语言和交互式环境主要包括Matlab和simulink 二维查表两大部分。Matlab的应用范围非常广包括信号和图像处理、通信、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。本文利用Matlab的M语言进行CORDIC算法的仿真其优势在于仿真结果易于观察，充分利用可以有效的加快设计进程

    本程序遵照上述CORDIC算法编写, 采用文件的方式编写，下面是部分源程序仿真输出波形如图3所示。