现代频率合成技术课程设计

1、基于M A T L A B的D D S设计与仿真摘 要：利用matlab仿真工具建立数字频率合成器DDS勺仿真模型，便于我们直截了当地 了解DDS勺工作原理和各部分模块的功能，而且便于我们分析DDS勺工作性能和各种参数指标。实验背景随着技术和器件水平的提高，称之为直接数字式频率合成器（DDS新的频率合成技 术得到飞速的发展。DDS在相对带宽、频率转换时间、相位连续性、正交输出、高分辨力以及集 成化等一系列指标方面，已远远超过了传统勺频率合成器所达到勺水平，完成了频率合成技术 的又一次飞跃。DDS与传统的DS和IS 一起构成了现代频率合成技术体系，将频率合成技术推 向了一个新的阶段。DDS勺原理

2、：数字频率合成是从相位概念出发直接合成所需波形的一种新技术，它采用一个恒定的输入 参数时钟，通过数据处理的方式产生频率、相位可调的输出信号。DDS系统由相位累加器、波形ROM、D/A转换器和低通滤波器构成。它具有频率分辨率高、频率切换时相位连 续等优点。DDS是继直接合成技术和锁相环式频率合成技术之后的第三代频率合成技术。他的工作原 理是基于相位与幅度的对应关系，通过改变频率控制字（K）来改变相位累加器（位数为N的相位累加速度，然后在固定时钟的控制下取样，取样得到的相位值（去取相位累加器 的高M位）通过相位幅度转换得到与相位值对应的幅度序列，幅度序列通过数模 转换及低通滤 波得到正弦波输出。下

3、图为DDS勺原理图。图1 DDS原理框图其中，K为频率控制字，为基准时钟频率，N为相位累加器的字长，D为ROh数据位及D/A转换器的字长。相位累加器在基准时钟的控制字下以步长K做累加，把相加后的结果送至相位累加器的输入，相位累加器一方面在上一时钟周期作用后产生的新的想位数据反 馈到自身的输入端，在下一个时钟的作用下继续与频率控制数据K相加，另一方面将这个值作为取样地址输出，送入正弦查找表ROM作为波形ROM勺地址，对波形ROM进行寻址。波形ROM俞出D位的幅度码S(n)经D/A转换器变成阶梯波S(t)，再经过低通滤波器平滑后就可以得到合成的信号波形。合成的信号波形取决于波形ROM中存放的幅度码

4、，因此用DD列以产生任意波形。(1)频率预置与调节K被称之为频率控制字，也叫相位增量DDS方程为，为输出频率，为时钟频率，当K= 1时，DDS输出最低频率(也即频率分辨率)为/，而DDS的最大输出频率由 Nyquist采样定理决定，即/2，也就是说K的最大值为-1。因此只要N足够 大，DDS可以得到很 细的频率间隔。要改变DDS的输出频率，只要改变频率控制字K即可(2)累加器图2累加器框图相位累加器由N位加法器和N位寄存器级联构成，来一个时钟脉冲，加法器将控制 字K与 寄存器输出的累加相位相加，再将相加后的结果送入寄存器的数据输入端。寄存 器将加法器在上 一个时钟作用后产生的相位数据反馈到加法

5、器的输入端，使加法器在下一个时钟作用下继续与频 率控制字进行相加。这样，相位累加器在时钟的作用下进行相 位累加。当相位累加器加满时产生 一次溢出，完成一个周期性的动作。(3)波形存储器用累加器输出的数据作为波形存储器的取样地址，进行波形的相位一幅值转化，即可在给定的时间上确定波形的抽样幅值。N位的寻址ROM相当于把的正弦信号离 散成具有 个样本值的序列，若波形ROMt D位数据，则2N个样值的幅值以D为二进制 数值固化 在ROM中，按照地址的不同可以输出相应正弦信号的幅值。相位一幅度变化原理图如下所示：图3相位一幅度变化原理图(4) D/A转换器D/A转换器的作用是把合成的正弦波数字量转换为模

6、拟量。正弦幅度量化序列S(n)经过D/A转换后变成了包络为正弦的阶梯波S (t )。需要注意的是频率合器对D/A转换器的分辨率有一定要求，D/A转换器的分辨率越高，合成的正弦波S(t)台阶数就越多，输出波形的精度也就越高。(5)低通滤波器，必须在两边土处的非谐波分量，幅值包络为辛格函数。因此为了去除主频D/A转化器的输出端接入截止频率为/2的低通滤波器。一、DDS性能DDS的频率合成原理及实现技术与传统的直接合成DS的锁相合成完全不同，在性能上也很独特。相对宽度当频率控制字k=1时，最低输出频率为=/M,式中M=2N当累加器字长N很大时，最低输出频率达Hz，mHz量级都是不困难的，可认为DDS

7、勺最低合成频率接近于零频。DDS的最高输出频率受限于时钟频率 和采样定理，=/2。在实际应用中，考虑到输出滤波器的非理想特性，一般采用=X 40%这样的DDS的相对带宽为=MX40%= X 40%频率分辨率DDS的最小频率步进量就是它的最低输出频率，即二二-二一。也可以采用十进制的相位累加器，那么M=10N可见只要累加器有足够的字长，实现非常精密的分辨率也没有多 大的苦难，正像前面介绍的一样，可达Hz、mHz、甚至PHz频率步进量。是传统频率合成技术所无可达到的。频率转换时间DDS的频率转换时间近似认为是即时的，这是因为它的相位序列在时间上是离散的。在频率 控制字K改变以后，要经一个时钟周期之

8、后才能按新的相位增量累加，所以可以说它的频率转换 时间就是频率控制字的传输时间，即一个时钟周期。目前，集成DDS产品的频率转换时间可达10ns的量级。这是常用锁相频率合成所无法达到的。频率转换时的相位连续性当频率控制字从K1变为K2之后，它是在已有的累加相位nK1S之上，再每次累加K2S,相位 函数的曲线是连续的，只是在改变频率瞬间其斜率发生了突变，因而保持了输出信号相位的连续 性。这一点对利用相位信息的那些系统很重要。相位连续可避免信息的丢失，相位不连续会导致频谱的扩散，不利于频谱资源的有效利用可输出正交信号有些应用场合要用到正交信号输出，即同时输出S 1(t)=sin(2 n t)和S2(

9、t)=cos(2 n t) 在DDS中，只要分别在两个ROM中存储 和两个函数表，即可同时输出正交信号，实现框图如下图所示。图4可输出正交信号的DDS框图可输出任意波形若在ROM中存储其他所需的波形函数表，DDS即可输出相应的周期性的波形，因此，更新 ROM中的数据，使DDS俞出方波、三角波、锯齿波等等。调制性能由于DDS是全数字的，用频率控制字K可直接调整输出信号的频率与相位，所以很易于在 DDS1实现数字调制和调相，很多DDS”品都具有数字调制功能。噪声和杂散因为DDS是数字技术，先构成离散信号，再变换成模拟信号输出，因而噪声与杂散的存在 是必然的。这是我们要特别关注的，以下对影响DDS俞

10、出的杂散来源进行分析。杂散分类及其影响相位截断对输出信号频谱的影响在实际应用中，为了提高频率分辨率，相位累加器的位数N尽量做得大，这就要求R0的容量很大。当N=32寸，就需要RO的容量达4GB这在实现上是很困难的。一般都是用相位 序列的高A位寻址ROM舍去相位序列的低B二N-A位。这就引入了相位截断误差。以下的分析设K 和 互质。当它们有公约数时，可以先化简，然后可归于以下模型或理想信 号模型：(n)是周期为入的阶梯波，其中，k=K mod(n)可以看作是对周期为入的锯齿波e(t)的采样，即表示相位截断误差信号(t)锯齿波e(t)的频谱为：，如图5所示。虚线表示锯齿波e(t)，实 线图5(n)

11、与e(t)关系其中，R=K mod信号波形序列为：nn 因为(n )，所以有:为阶梯波的梯度。经过ROM相位一幅度转换后输出的n7t-7t -由相位舍位引起的输出误差信号为：n把 S( n)看作对连续信号的采样：nns nt)的频谱为-o)-o),e(t)的频谱E如上述，所以e(t)nt)的频谱为二者的频域卷积:误差信号经过D/A后的模拟信号表达式为隔(t)= AS( t)h(t),所以最终输出的误差信 号频谱为鸟(=ASH ), H对信号频谱的影响是一个SINC函数的幅度调制，对谱线的位置分布没有影响，从信号的频谱结构角度可以把它带来的影响忽略。由上式可以得 出，相位截断效应带来的频谱杂散位

12、于：,n可见，有相位截断的DDS输出信号频谱杂散分量十分丰富，并且有大量的频率值落在O,fo之间，从理论上也无法将其完全滤掉。由的表达式还可以发现，它和频率控制字K有关，不同的输出频率其频谱结构也会不同，这就为提高信号质量增加了困难。幅度量化误差对输出信号频谱的影响ROM数据位宽的有限使得DDS输出的信号为阶梯波。这将对输出信号引入幅度量化误差(n)。设ROM数据位宽为D,不考虑相位截断效应，幅度量化误差表示为：由取整函数in t( )的性质可知，误差函数(n)是一个位于区间,上的随机序列，概率密度服从均匀分布，幅度量化对信号频谱的影响可以归于白噪声。但是当 K 和 成整数比例关系时，由噪声信

13、号表达式可知，多个周期的噪声具有重复性，因此噪 声分布也具有一定的周期性。D/A非理想特性及参考时钟相位噪声对输出信号频谱的影响D/A的非理想特性主要包括动态非线性、静态非线性、有限分辨率及内部闪烁噪声等。D/A 的非理想特性难于建模，不同的器件性能各异，只能根据具体的器件参数分别考虑。而且D/A 对输出信号频谱的影响跟相位截断效应和幅度量化误差比起来是很小的。只有在对信号频谱质 量要求非常高的应用中才根据具体情况对这个问题进行深入的研究。这也是为什么在DDS技术发展的二十多年间，对这个问题的研究不是很多。随着电子制作工艺的迅猛发展，D/A期间的非理想特性逐渐得到改善，一般可以通过选用高性能的

14、D/A器件来满足应用的要求。参考时钟的相位噪声会传递到输出信号中去。DDS实际上可以看作一个分频器。从理论上来说，输出信号的相位噪声会对参考时钟的相位噪声有dB的改善。在对DDS信号频谱作理论分析时，这些都不是主要的研究对象，但是在实际应用中，尤其是在对信号质量有苛刻要求的场合，这些因素的影响必须引起重视。MATLA仿真和结果分析4.1仿真参数的设置本次设计的要求为输出频率为1MHz相位累加器的位数N=10，并调N和频率控制字K的值分析相位截断误差等各类误差的变化。根据理想DDS勺频谱图(如图6所示)可知杂散分量对输出信号质量影响最为严重。当 趋近于 时，也趋向于，两者很难区分，信号质量无法保

15、证。为了使低通滤波器有效地滤除杂散，一般小于。因此在此我们分别取、12MHz、16MHz对结果作比较。由公式算得取整后的K分别为128、85、64。图6理想DD输出的频结构4.2 仿真结果4.2.1 K参数的影响研究为了直观的比较相位截断前后的波形时域和频域上的差别，以下用 A=4来模拟R0寻址的相位序列的高4位。在不同时钟频率和频率控制字K下阶段前后的时域和频域图如下图所示。(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)图7仿真结果从图7中可以看到，当输出频率和累加器位数N一定时，频率控制字越小波形更加平滑，这是因为当N位一定时，ROM中的存储的幅度值位个，累加器每次

16、累加步长为K,因此一周期内R0 口的采样数据数为floor()个(floor为Matlab求整函数，求向零靠近的整数)，由此可见K越小，一周期门内采样值越多，波形完整性越强，这结论跟实验结果完全符合。再次我们可以看到当累加器低N-M位被截取后，波形的完整性进一步恶化，这是因为当K 一定时，ROM实际被采样到数据数为floor()个，显然比没被截断前少，导致波形平滑度降低。以上实验得出的时域波形中，更加直观的看出相位阶段 带来的额外 量化误差。4.2.2 N参数的影响研究为了分析相位累加器位数N寸输出波形的影响，我们分别取N=8、10、12作为一组比较系。以下为本次实验结果。(M=4 ,K=63

17、,0.80.60.40.2-0.2-0.4-0.6-0.8x 10(1)图8不同N下的DDS俞出由图8（ 1）,（ 5）可以看出，当N比较小时，输出波形有由明显的低频分量，而且，因此引起明显的波形更接近于三角波，这是因为N=8时，一个周期内采样数据数为量化误差。以上信号是未通过低通滤波其处理的，因此由（6）更直观地看到从D/A出来的模拟信号是在DD原理中介绍的阶梯波S（t）。从频谱的分析来看，相位阶段前后的频谱由明显的杂散分量的加入，而且，随N的减小而增多。综上，在实际DD设计中，相位累加器的位数N要尽可能的取大，以此降低幅度量 化带来的误差。但N越大需要的R0存储空间越大，实际中很难实现。因

18、此R0压缩技术和 截断误差的降低技术是DD设计中的关键问题所在。423功率谱分析以下图中（1）、（ 3）分别为单1如N=12寸的有量化误差的信号功率谱图。当N从10 变到12时，则信号功率与量化噪声总功率之比下降约12dB。这是因为当ROM采用。位二进 制数保存正弦函数值时，量化误差为：其中,R :x表示对x做最靠近x的取整运算。显然，与S（n）有相同的序列周,因此幅度量化误差在频谱中没有引入新的杂散成分,而是表现为均匀的噪声基底。通常在一个周期内被认为是在间均匀分布的噪声，则由量化引起的信噪比为：由上式可见，量化位数D每增加一位，则SNR将提高6dB(1)(2)(3)(4)图9 DDS杂散功

19、率谱从图9的(2),( 4)中的频域图比较中可以看到，处理截断后的误差信号也成周期出现，对此在截断误差的表达式n中，已知为周期入的阶梯波，其中为K和的最大公约数。由数字信号处理理论易知的谱线以为周期，在区间(0,)内的谱线由根谱线所组成，这个结论跟实验结果吻合。附录1)DD实现及分析相位截断前后的波形及频谱的代码：clear all;N=12;%累加器的位数；K=63; %fix(2m)*0.222);%频率控制字，即累加的步长；M=6;%截取累加器的高10位；D=8;%8bit DACg=gcd(KQN);pe=2AN/(g);n=1:pe;pp=pe+1; %频谱分析点数add_y=mod

20、(n*K,2AN); %累加器的输出表达式；error=mod( n*K,2A(N-M);rom_x=si n(2*pi*add_y/(2AN);rom_y=si n(2*pi*(add_y-error)/(2AN);Fs=16*10A6;%米样频率%此时的相位增量为2*pi*K/MTs=1/Fs;%米样频率t=(0:pe-根据需要输出的信号点数，取得时间1)*Ts;%t=(0:Nd-1)*Ts; % plot(t,rom_x);xlabel(t/s);ylabel(A);title(N=，num2str(N),截断前的时域图);grid on;Fy=abs(fft(rom_x,4*pe);对

21、信号进行傅里叶变换len_Fy=length(Fy); %此处的数字角频率f=(0:len_Fy)/len_Fy*Fs;(W=2*pi*f/fs )%输出信号的频率，将数字频率换成模拟频率 转*次”、十count二floor(len_Fy/2); %频谱显示的点数figure %plot(f(1:count),Fy(1:count) );显示多幅图像grid on;xlabel(f/Hz);ylabel(A);title(N=，num2str(N),截断前的频域图);figureplot(t,rom_y);xlabel(t/s);ylabel(A);title(N=，num2str(N),截断

22、后的时域图);grid on;Fy二abs(fft(rom_y,4*pe); %对信号进行傅里叶变换len_Fy=length(Fy); %此处的数字角频率f=(0:len_Fy)/len_Fy*Fs%输 出信号的 频率，将 数字频率转 换成模拟 频率(W=2*pi*f/fs ) count=floor(len_Fy/2); % figure % plot(f(1:count),Fy(1:count); grid on;xlabel(f/Hz); ylabel(A);频谱显示的点数title(N=,num2str(N),显示多幅图像2)分析相位截断误差和幅度量化误差的功率谱的代码clearall;N=10; %累加器的位数；K=63; %(20)*0.222);%频率截断后的频域图);控制字，即累加的步长；M=4;%截取累加器的高10位；D=8;%8bit DACg=gcd(KQN);pe=2AN/(g);n=1:pe;fc=16*10A6;pp=pe+1; %频谱分析点数add_y=mod(n*K,2AN); %累加器的