现代频率合成技术课程设计-基于MATLAB的DDS设计与仿真

基于MATLAB的DDS设计与仿真摘要：利用matlab仿真工具建立数字频率合成器DDS的仿真模型，便于我们直截了当地了解DDS的工作原理和各部分模块的功能，而且便于我们分析DDS的工作性能和各种参数指标。实验背景随着技术和器件水平的提高，称之为直接数字式频率合成器（DDS）新的频率合成技术得到飞速的发展。DDS在相对带宽﹑频率转换时间﹑相位连续性﹑正交输出﹑高分辨力以及集成化等一系列指标方面，已远远超过了传统的频率合成器所达到的水平，完成了频率合成技术的又一次飞跃。DDS与传统的DS和IS一起构成了现代频率合成技术体系，将频率合成技术推向了一个新的阶段。DDS的原理：数字频率合成是从相位概念出发直接合成所需波形的一种新技术，它采用一个恒定的输入参数时钟，通过数据处理的方式产生频率﹑相位可调的输出信号。DDS系统由相位累加器﹑波形ROM﹑D/A转换器和低通滤波器构成。它具有频率分辨率高﹑频率切换时相位连续等优点。DDS是继直接合成技术和锁相环式频率合成技术之后的第三代频率合成技术。他的工作原理是基于相位与幅度的对应关系，通过改变频率控制字（K）来改变相位累加器（位数为N）的相位累加速度，然后在固定时钟的控制下取样，取样得到的相位值（去取相位累加器的高M位）通过相位幅度转换得到与相位值对应的幅度序列，幅度序列通过数模转换及低通滤波得到正弦波输出。下图为DDS的原理图。图1DDS原理框图其中，K为频率控制字，fc为基准时钟频率，N为相位累加器的字长，D为ROM数据位及D/A转换器的字长。相位累加器在基准时钟fc的控制字下以步长K做累加，把相加后的结果送至相位累加器的输入，相位累加器一方面在上一时钟周期作用后产生的新的想位数据反馈到自身的输入端，在下一个时钟的作用下继续与频率控制数据K相加，另一方面将这个值作为取样地址输出，送入正弦查找表ROM，作为波形ROM的地址，对波形ROM进行寻址。波形ROM输出D位的幅度码S(n)经D/A转换器变成阶梯波S(t频率预置与调节K被称之为频率控制字，也叫相位增量。DDS方程为f0=fcK/2N,f0为输出频率，f累加器图2累加器框图相位累加器由N位加法器和N位寄存器级联构成，来一个时钟脉冲fc波形存储器用累加器输出的数据作为波形存储器的取样地址，进行波形的相位—幅值转化，即可在给定的时间上确定波形的抽样幅值。N位的寻址ROM相当于把0o~360o的正弦信号离散成具有2N图3相位—幅度变化原理图D/A转换器D/A转换器的作用是把合成的正弦波数字量转换为模拟量。正弦幅度量化序列S

（n）经过D/A转换后变成了包络为正弦的阶梯波S（t）。需要注意的是频率合器对D/A转换器的分辨率有一定要求，D/A转换器的分辨率越高，合成的正弦波S(t)台阶数就越多，输出波形的精度也就越高。低通滤波器对D/A输出阶梯波S（t）进行频谱分析，可知S（t）中除主频fo外，还存在分布在fc，2fc……两边±fo处的非谐波分量，幅值包络为辛格函数。因此为了去除主频DDS性能DDS的频率合成原理及实现技术与传统的直接合成DS的锁相合成完全不同，在性能上也很独特。相对宽度当频率控制字K=1时，最低输出频率为fomax=fc/M，式中M=2N，当累加器字长N很大时，最低输出频率达Hz，DDS的最高输出频率受限于时钟频率fc和采样定理，fomax=fc/2。在实际应用中，考虑到输出滤波器的非理想特性，一般采用fomax=fc×40%。这样的DDS的相对带宽为fomax频率分辨率DDS的最小频率步进量就是它的最低输出频率，即Δfo=fomin=fcM=fc2N频率转换时间DDS的频率转换时间近似认为是即时的，这是因为它的相位序列在时间上是离散的。在频率控制字K改变以后，要经一个时钟周期之后才能按新的相位增量累加，所以可以说它的频率转换时间就是频率控制字的传输时间，即一个时钟周期Tc=1频率转换时的相位连续性当频率控制字从K1变为K2之后，它是在已有的累加相位nK1δ之上，再每次累加K2δ，相位函数的曲线是连续的，只是在改变频率瞬间其斜率发生了突变，因而保持了输出信号相位的连续性。这一点对利用相位信息的那些系统很重要。相位连续可避免信息的丢失，相位不连续会导致频谱的扩散，不利于频谱资源的有效利用。可输出正交信号有些应用场合要用到正交信号输出，即同时输出S1(t)=sin(2πfot)和S2(t)=cos(2πfo在DDS中，只要分别在两个ROM中存储sinθ和cos图4可输出正交信号的DDS框图可输出任意波形若在ROM中存储其他所需的波形函数表，DDS即可输出相应的周期性的波形，因此，更新ROM中的数据，使DDS输出方波﹑三角波﹑锯齿波等等。调制性能由于DDS是全数字的，用频率控制字K可直接调整输出信号的频率与相位，所以很易于在DDS上实现数字调制和调相，很多DDS产品都具有数字调制功能。噪声和杂散因为DDS是数字技术，先构成离散信号，再变换成模拟信号输出，因而噪声与杂散的存在是必然的。这是我们要特别关注的，以下对影响DDS输出的杂散来源进行分析。杂散分类及其影响相位截断对输出信号频谱的影响在实际应用中，为了提高频率分辨率，相位累加器的位数N尽量做得大，这就要求ROM的容量很大。当N=32时，就需要ROM的容量达4GB，这在实现上是很困难的。一般都是用相位序列的高A位寻址ROM，舍去相位序列的低B=N-A位。这就引入了相位截断误差。以下的分析设K和2BεεP(n)是周期为λ=2B(k,2B)的阶梯波，其中，k=Kmod2B。εP(n)可以看作是对周期为λTc的锯齿波e(t)的采样，即εPt图5εP(n)锯齿波e(t)的频谱为：E(ω)其中，ωx=R2Bωc,R=Kmod2BS因为εP(n)<<2S(n)≈由相位舍位引起的输出误差信号为：?S把ΔS(n)看作对连续信号的采样：?Ssin(2π?0t)的频谱为δ(ω-ω0)-δ(ω+ω0),e(t)的频谱E(ω)如上述，所以e(t)Sin(2π?0t)的频谱为二者的频域卷积：φ误差信号ΔS（t）的频谱为?Sω?S误差信号经过D/A后的模拟信号表达式为ΔS0(t)=ΔS(t)h(t)，所以最终输出的误差信号频谱为ΔS0(ω)=ΔS(ω)H(ω)，H(ω)对信号频谱的影响是一个SINC函数的幅度调制，对谱线的位置分布没有影响，从信号的频谱结构角度可以把它带来的影响忽略。由上式可以得出，相位截断效应带来的频谱杂散位于：?=mFc可见，有相位截断的DDS输出信号频谱杂散分量十分丰富，并且有大量的频率值落在[0,f0]之间，从理论上也无法将其完全滤掉。由ωx的表达式还可以发现，它和频率控制字K幅度量化误差对输出信号频谱的影响ROM数据位宽的有限使得DDS输出的信号为阶梯波。这将对输出信号引入幅度量化误差εP(n)。设ROM数据位宽为Dε由取整函数int(·)的性质可知，误差函数εP(n)是一个位于区间[-12·2B,12·2B]D/A非理想特性及参考时钟相位噪声对输出信号频谱的影响D/A的非理想特性主要包括动态非线性、静态非线性、有限分辨率及内部闪烁噪声等。D/A的非理想特性难于建模，不同的器件性能各异，只能根据具体的器件参数分别考虑。而且D/A对输出信号频谱的影响跟相位截断效应和幅度量化误差比起来是很小的。只有在对信号频谱质量要求非常高的应用中才根据具体情况对这个问题进行深入的研究。这也是为什么在DDS技术发展的二十多年间，对这个问题的研究不是很多。随着电子制作工艺的迅猛发展，D/A期间的非理想特性逐渐得到改善，一般可以通过选用高性能的D/A器件来满足应用的要求。参考时钟的相位噪声会传递到输出信号中去。DDS实际上可以看作一个分频器。从理论上来说，输出信号的相位噪声会对参考时钟的相位噪声有20logFc/f0MATLAB仿真和结果分析仿真参数的设置本次设计的要求为输出频率为1MHz，相位累加器的位数N=10，并调N和频率控制字K的值分析相位截断误差等各类误差的变化。根据理想DDS的频谱图（如图6所示）可知杂散分量Fc-f0对输出信号质量影响最为严重。当f0趋近于Fc/2时，Fc-f0也趋向于Fc/2，两者很难区分，信号质量无法保证。为了使低通滤波器有效地滤除杂散，f0一般小于图6理想DDS输出的频结构仿真结果K参数的影响研究为了直观的比较相位截断前后的波形时域和频域上的差别，以下用A=4来模拟ROM寻址的相位序列的高4位。在不同时钟频率Fc（1）（2）（3）（4）（5）（6）（7）（8）（9）（10）（11）图7仿真结果从图7中可以看到，当输出频率和累加器位数N一定时，频率控制字越小波形更加平滑，这是因为当N位一定时，ROM中的存储的幅度值位2N个，累加器每次累加步长为K，因此一周期内ROM的采样数据数为floor(2N/K)个（floor为Matlab求整函数，求向零靠近的整数），由此可见K越小，一周期门内采样值越多，波形完整性越强，这结论跟实验结果完全符合。再次我们可以看到当累加器低N-M位被截取后，波形的完整性进一步恶化，这是因为当K一定时，ROM实际被采样到数据数为floor(N参数的影响研究为了分析相位累加器位数N对输出波形的影响，我们分别取N=8、10、12作为一组比较系。以下为本次实验结果。（M=4，K=63,Fc=16MHz(1)(2)(3)(4)(5)(6)(7)(8)图8不同N下的DDS输出由图8（1），（5）可以看出，当N比较小时，输出波形有由明显的低频分量，而且波形更接近于三角波，这是因为N=8时，一个周期内采样数据数为2NK从频谱的分析来看，相位阶段前后的频谱由明显的杂散分量的加入，而且，随N的减小而增多。综上，在实际DDS设计中，相位累加器的位数N要尽可能的取大，以此降低幅度量化带来的误差。但N越大需要的ROM存储空间越大，实际中很难实现。因此ROM压缩技术和截断误差的降低技术是DDS设计中的关键问题所在。功率谱分析以下图中（1）、（3）分别为N=10和N=12时的有量化误差的信号功率谱图。当N从10变到12时，则信号功率与量化噪声总功率之比下降约12dB。这是因为当ROM采用D位二进制数保存正弦函数值时,量化误差为:?其中,R〔x〕表示对x做最靠近x的取整运算。显然,?Mn与S(n)有相同的序列周期2N/2N,K,因此幅度量化误差在频谱中没有引入新的杂散成分,而是表现为均匀的噪声基底。通常在一个周期内SNR由上式可见,量化位数D每增加一位,则SNR将提高6dB。（1）（2）（3）（4）图9DDS杂散功率谱从图9的（2），（4）中的频域图比较中可以看到，处理截断后的误差信号也成周期出现，对此在截断误差的表达式

?Sn=2π2NεPnsin2πnK2N中，已知εPn为周期λ=2附录DDS实现及分析相位截断前后的波形及频谱的代码：clearall;N=12;%累加器的位数；K=63;%fix((2^N)*;%频率控制字，即累加的步长；M=6;%截取累加器的高10位；D=8;%8bitDACg=gcd(K,2^N);pe=2^N/(g);n=1:pe;pp=pe+1;%频谱分析点数add_y=mod(n*K,2^N);%累加器的输出表达式；error=mod(n*K,2^(N-M));rom_x=sin(2*pi*add_y/(2^N));rom_y=sin(2*pi*(add_y-error)/(2^N));Fs=16*10^6;%采样频率%此时的相位增量为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),',','截断前的时域图']);gridon;Fy=abs(fft(rom_x,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));gridon;xlabel('f/Hz');ylabel('A');title(['N=',num2str(N),',','截断前的频域图']);figureplot(t,rom_y);xlabel('t/s');ylabel('A');title(['N=',num2str(N),',','截断后的时域图']);gridon;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));gridon;xlabel('f/Hz');ylabel('A');title(['N=',num2str(N),',','截断后的频域图']);2）分析相位截断误差和幅度量化误差的功率谱的代码clearall;N=10;%累加器的位数；K=63;%((2^N)*;%频率控制字，即累加的步长；M=4;%截取累加器的高10位；D=8;%8bitDACg=gcd(K,2^N);pe=2^N/(g);n=1:pe;fc=16*10^6;pp=pe+1;%频谱分析点数add_y=mod(n*K,2^N);%累加器的输出表达式；error=mod(n*K,2^(N-M));rom_y=cos(2*pi*(add_y-error)/(2^N));%相位到幅度映射的输出；pat=-1+1/(2^(D-1)):1/(2^(D-1)):1-1/(2^(D-1));%设定量化区间codebook=-1+1/(2^D):1/(2^(D-1