您好,欢迎来到刀刀网。
搜索
您的当前位置:首页Matlab在DDS系统仿真中的应用

Matlab在DDS系统仿真中的应用

来源:刀刀网
维普资讯 http://www.cqvip.com 计算机与现代化 2006年第6期 JISUANJI YU XIANDAIHUA 总第130期 文章编号:1006-2475(2006)06-0018.04 Matlab在DDS系统仿真中的应用 J能、、、 藏 、、 (广东省教育学院计算机科学系,广东广州510303) 摘要:介绍了直接数字频率合成器(Direct Digital frequency Synthesizer,DDS)的基本结构和工作原理,并在此基础上着重论 述了利用Matlab软件建立DDS的理想和杂散仿真模型的方法。通过建立仿真模型,可以更方便有效地对DDS杂散分布 特点进行分析。 关键词:M'At ̄ah;直接数字频率合成;DDS;系统仿真 中图分类号:TP391,I 文献标识码:A Application of Matlab in Simulation of DDS System XIONG Yah (Department of Computer Science,Guangdong Education Insittute,Guangzhou 510303,China)  ̄:This paper first introduces the structure and fuldanmaal principles of all ideal direct digital frequency synthesizer.Then a method of simulatin8 ideal and spur model of DDS啦Matlab is discussed.Using these simulation models,the efect of the spurs on DDS system can be删yzed lllO/'e egectively. Key words:Matlab;direct dig frequency synthesis;DDS;system simulation 0 引 言 Maflab已经发展成为多学科、多种工作平台的功能强 大的大型软件,同其他的仿真软件相比,Matlab采用 频率控制是现代通信中很重要的一环,而频率控 了数据类型和面向对象编程技术,推出了全新的符号 制的重要研究内容就是要获取宽带、快速、精细、杂散 计算工具包,提供了各种“友善”的交互式操作工具, 小的频率控制信号。直接数字频率合成(DDS)是近 是当今广泛为人们采用的控制系统数字仿真应用软 年来得到迅速发展的一种新的频率合成方法,它从相 件。本文用Maflab软件建立DDS的理想和杂散仿真 位概念出发,以全数字技术直接合成所需要波形,具 模型,并利用DDS仿真模型对其杂散幅度大小和分 有输出频段宽、频率转换速度快、频率分辨率高、输出 布进行了分析,为DDS新结构和新算法的研究搭建 相位连续等一系列优异性能。但是由于DDS数字化 了良好的仿真平台。 的固有特点,决定了其输出频谱杂散较大,这成为 DDS进一步发展的瓶颈,因此杂散分析与降低杂散的 1 DDS的基本工作原理 设计成为DDS研究的重点。通过对DDS的频谱进行 DBS技术是一种把一系列数字量形式的信号通过 计算机仿真,可以找到DDS的杂散分布规律,寻求有 DAC转换成溪拟量形式的信号的合成技术。基本模型如 效的杂散抑制技术,这无论是在理论分析上,还是在 图1所示,主要由时钟频率源 、相位累加器、波形存储器 实际应用中都有极其重要的意义。 Madab是美国MathWolks公司自20世纪80年代 (R0M)、数 豫 器(D/A)、以及低通滤波器(try)组成。 中期推出的数学软件,优秀的数值计算能力和卓越的 在时钟脉冲 的控制下,相位累加器对频率控制字K进 数据可视化能力使其很快在数学软件中脱颖而出。 行累加,使得相位累加器输出的相位序列在每次累加时, 有一相位增量△ = (1 )。相位累加器输出的相位序 收稿日期:2006-02.27 作者简介:熊燕(1970-),女,江西南昌人,广东省教育学院计算机科学系助教,硕士,研究方向:微波移动通信,信号处理。 维普资讯 http://www.cqvip.com 20O6年第6期 熊燕:Matlab在DDS系统仿真中的应用 19 列寻址ROM进行相码.幅码变换,输出离散的幅度编码; 该幅度编码经D/A后得到的阶梯电压。最后经低通滤波 器u 平滑,得到由频率控制字K决定的连续变化的输 出波形。 图l 正弦输出的DDS原理框图 2加S系统仿真模型的建立 2.1理想模型的建立 理想DDS是DDS的理想化数学模型,所谓的理 想参数DDS具备三个条件: (1)没有相位舍位误差,即相位累加器的N位输出都用 于ROM寻址。(2)没有幅度量化误差,即ROM输出值用无限 长的码来表示。(3)DAC的分辨率无穷大、并且DAC具有理 想的数模转换特性。 对理想化的DDS而言,设 为DDS的参考时钟 信号,K为频率控制字,N为相位累加器的位数,fo为 DDS的输出频率,它们之间的关系可以表示为:fo= ×K/2N。 尽管在计算机上数据精度的范围同样受限,但鉴 于DDS系统中数据的精度通常远小于计算机,因此 在计算机上仍可用高精度的浮点数来模拟理想状态。 设在Mathb仿真中,accumulator_y表示相位累加器的 输出信号,mm.y表示ROM查询表的输出信号,则理 想化的DDS中各组成部分可以采用以下方法来用 Matlab进行仿真建模: 1.系统时钟。 在DDS系统中,各部分的时序由系统时钟来驱 动,因此各部分的输出均可看作离散时间序列的函 数,在Matlab中可以用序列n={0,1,2,…}来模拟。 2.相位累加器。 相位累加器在参考时钟的作用下,进行线性相位累 加,当相位累加器累加满量时就会产生一次溢出,完成一 个周期性的动作,因此在相位累加器的输出端可以得到的 波形序列是:n・ 蜘 ,其中n代表的是相位序列中的第 n个时刻,K为频率控制字,N为相位累加器的位数。由上 面分析可知,在 中相位累加器可以用语句 accu ̄,latnr_y=Il (n*K, )来模拟o 3.ROM查询表。 理想状态下,ROM的输出值具有无限长的精度, 所以ROM输出端可以得到理想DDS输出的波形样 本序列为: s(n)=cos( l(II) 在Matlab仿真中,无限精度的正弦值可采用浮点 数来近似,因此ROM表可以采用语句rom_y=COS (accumulator_y (2 pi/2^N))来模拟。 4.DAC转换。 由于在理想状态下,DAC具有理想的数模转换 特性,因此,这部分不产生波形变化。 综上所述,我们可以编写以下Matlab仿真程序, 实现在理想状态下的DDS波形合成过程。 elear K=5;N=10;%频率控制字K为5,相位累加器的位数N为l0 n=1:TN;%时间序列取整数个周期 aeeumuhtor_y=mod(n*K, N);%相位累加器输出序列 rom.y=cos(accumulator_y (2 pi/2^N));%ROM查询表 输出信号序列 [Pxx,W]=periodogram(rom.y,[], onesided ,1024);%求 ROM查询表输出信号的频谱 ifgure; psdplot(Pxx/max(Pxx),W,”,”,”)%画出DDS的频谱图 d on%加坐标线 2.2 杂散模型的建立 实际DDS并不满足理想DDS的3个条件。由于 DDS的全数字化结构,因而实际DDS存在着固有的 误差。其主要误差来源有三个: (1)相位截断误差:在实际DDS电路中。为了取得一定的 频率分辨率,通常相位累加器的位数N取得很大,如N=24, 32,48等,但受体积和成本的,ROM的容量却远小于2 , 因此寻址ROM时,累加器输出相位序列的低B位就被舍去, 而只用其输出的高M(M=N—B)位去寻址,这样就不可避免 地引入误差。(2)幅度量化误差:ROM中存储着正弦波样点的幅 值编码,任何—个呻哥度值要用无限长的比特流才能精确表示,而实 际中ROM的输出位D是个有限值,这样就引入了幅度量化误差。 (3)DAC转换误差:DAC的有限分辨率、非线性特性、瞬间毛 刺及转换速率等非理想转换特性会影响DDS输出频谱的纯 度,产生杂散分量。 1.相位截断情况下仿真模型的建立。 假定DAC具有理想的性能,并且不考虑幅度量 化误差,仅对相位截断给DDS频谱带来的影响进行 分析。设用N位累加器输出高M位对ROM寻址,舍 位B=N—M,频率控制字为K,时钟频率为 。截取高 M位后输出的相位序列 。(n)为: n)=中(n)一( ̄(n)mod2 )=nKmod2N—nK mod2 (1) 根据式(1),相位累加器除了得到理想波形序列 13・K MOD2N外,还包含由相位截断引起的相位误差 序列:£p(n)= (n)一 p(n)=nKmod2B。令 accumulator_y、elTor分别表示理想波形序列和相位误 差序列,则在Matlab中理想波形序列和相位误差序列 维普资讯 http://www.cqvip.com 2O 计算机与现代化 2OO6年第6期 可以表示为:ac.cu ̄nlator_y Il (n*K,2 N),eriDr= Il (n*K,2^B)。令trun表示ROM查询表中查询的 相位序列,则根据式(1)该序列在Matlab中可以表示 为:tirm=c0s(2*pi*(accum-Iotor_y-error)/2 N)o这 样在相位截断情况下DDS的Matlab仿真程序如下: clear K=IO0;N=12;B=5;%频率控制字K为5,相位累加器 的位数N为1O,相位截断的位数为5 n=1:TN%时间序列取整数个周期 accu ̄r.y=mod(n*K,2‘N);%相位累加器输出的理 想序列 哪=mDd(n*K,2旧);%相位漯加器输出的相位 序列 咖=c06(2*pi*(accmnula ̄.y-ertor)/2 N);%相位截 断情况下ROM查询表输出信号序列 [Pxx,W]= 0d ( ,[], one ̄ided ,lou);%求 ROM查询表输出信号的频谱 p6dIplot(P ,W, , , )%画出DDS的频谱图 g甜%加坐标线 2.幅度量化情况下仿真模型的建立。 波形ROM中正弦波幅度码是量化后固化进去 的,因此存在量化误差。量化一般采用舍入方式,量 化阶为q=2-(D-1)(D为量化位数,即ROM字长)。 在Matlda中,量化操作可以通过quaatiz函数来 完成。使用该函数首先要设定量化区间和码本。因 为正弦函数的函数值介于一l到l之间,且量化阶为 q=2一(D一¨,根据舍入量化的基本原理,可以在一l+ l/2D和l—l/2D之间,每隔1/2(D-t)取为一个量化区 间,取每个区间的中值作为对应的量化码本值,在区 间(一ao,一l+1/2D)和[1—1/2D,∞),则分别取一l 和l作为量化码本值。则计算量化区间和量化码本 的Matlab程序为: pat=一1+1/(T(D一1)):1/(T(D一1)):1—1/(T(D一1)); co&hook=一1+l/( D):l/( (D—1)):l一1/(TD) 当量化区间和码本计算好后,作为参数传递给 蛐tiz函数,以仿真幅度量化的过程。这样在幅度量 化情况下DDS的Matlab仿真程序如下: clear K=4;N=8;D=3%频率控制字K为5,相位累加器的位 数N为lO,量化阶数为5 n=1:TN%时间序列取整数个周期 acc ̄r.y mt(n*w ̄1.k,Taccumulmr.n);%相位累 加器输出的信号序列 arm.y=cos(accmnulmor.y (2 p acom,l ̄or.n));% ROM查询表输出信号序列 pat=一1+1/(2^(D一1)):1/(2 (D一1)):1—1/(2 (D— 1))%设定量化区间 codd)00k=一1+1/(TD):1/(T(D一1)):1—1/(2^D)%设 定量化码本值 [index,quanta]=quantiz(mm_y,pat,code.book)%量化函数 进行量化 [ ,W]=periodov..(trunk,[],’onesided ,1024);%求量 化后输出信号的频谱 psdplot(Pxx,W, , , )%画出DDS的频谱图. 3.DAC非线性误差。 DDS杂散信号的另一个主要杂散来源是D/A的 非理想特性产生的,尤其是D/A的动态特性的好坏 对杂散的影响很大。理论上,D/A为一线性器件,它 将所输入代表幅度的数字信号按照一定的比例还原 成为模拟信号,可是由于实际工艺,很难制造出 这样理想的、完全线性的D/A,实际的D/A器件都带 有不同程度的非线性。但由于DAC的非线性误差通 常和具体电路器件有关,所以在本文中我们假设DAC 是理想状态,仿真中没有考虑DAC的非线性误差。 3加S系统仿真结果分析 一善j,Pe奄一h 墨1l鼻l点 3.1 理想状态下加S输出信号的计算机仿真 图2是在Mathb仿真中理想状态下DDS输出波 形的频谱图,图中所有谱线的幅度对主谱幅度进行归 一化处理,横坐标对7c进行了归一化处理。显示的是 (0,fd2)的频谱。从图中可以看出理想状态下DDS 输出波形的频谱纯度很高。杂散分量很小。 { l l .。竹 钠 M ’_’ ̄‘1. 一一 J . m Fn ̄'my(X md/目唧k) 图2理想状态下DDS输出频谱图 3.2相位截断情况下加S输出信号的计算机仿真 图3是用Matlab仿真软件得到的相位截断情况 下DDS的频谱图。这里不考虑幅度量化和D/A转换 器对频谱的影响。 仿真中取频率控制字K=3,累加器的位数N=10,累加器 舍去的位数B=3,输出主谱频率 =等= 。主谱与最强 杂散幅度之比应为: ‘ 5)da=ZOlg[一 1 ≥6‘N—B’曲 ‘ o— U k . 3)d8=42dB 在仿真实验中得出最强杂散幅度为一41.855dB,与理论 计算值非常接近。 维普资讯 http://www.cqvip.com 2OO6年第6期 一崔j重宣兽1暑善-王 熊燕:1VIAtlAh在DDS系统仿真中的应用 一蕾j,P_宴_I E口1量丹-王 一蕾重 亩一督兽1鼻l 2l 相位截断误差分析 。埘埘坷 硼埘图4和图5中N=1 。 0,K=65都相同。埔坷 硼埘 不同的是幅度量化 的比特数D。图4和图5中的最大杂散幅度分别为:一75. l 721dB、一81.643dB。可以看出。当D越大。DDS的杂散功率越 小。而且量化比特数D每增加一位,信噪比将提高约6dB。 h- — ' , m 、 y ,^^ nt_ 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 l 4结束语 l ̄Iom,Jli,.,dl ̄Le'qeney(×_,ld/m,,ple) 图3 N=10,K=3,B=3情况下DDS输出频谱图 本文根据直接数字频率合成器的基本结构和工 3.3幅度量化情况下加S输出信号的计算机仿真 作原理,利用Madab软件建立了DDS的理想和杂散 本文对幅度量化情况下的DIYS频谱进行了 仿真模型,并利用所建立的仿真模型分析了DDS的 Mallab软件仿真实验,在没有考虑相位截断和DAC 杂散幅度大小和分布,找出DDS的杂散分布规律。 非线性的情况下实验结果如图4、图5所示。 本文的工作为DIYS新结构和新算法的研究搭建了良 曩度量化误差分析 好的仿真平台。 参考文献: [1]王育红,张厥盛.直接数字频率合成器的谱质研究[J]. 西安电子科技大学学报。1996。23(1):19—26. 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 l №m -ed Fl。qcI呵(×_HI/Ⅻ叩le) [2]李衍忠。等.DDS谱质分析及其杂散抑制研究综述[J]. 图4 N=10,K=65。B=9t1/况下DDS输出频谱图 I曩度量化误差分析 现代雷达,2O00(8)。33—38. } } [3]汤汉屏.直接数字频率合成杂散抑制方法的研究[J]. 电讯技术,2O00(5):l0一l4. ’ [4]王晓音。聂裕平,庞伟正.DDS输出频谱杂散的抑制[J]. 电子对抗技术。2003,18(6):25—28. I山 IIIil lIIIl 1 [5]黄文梅。熊桂林,杨勇.信号分析与处理——N ̄ltab语言 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 l 及应用[M].长沙:国防科技大学出版社,2000. №Ⅲ・I酬Fleq∞叮(× _I/l_呷Ie) [6]刘敏,魏玲.N-t!-h通信仿真与应用[M].北京:国防工 图5 N=10,K=65,B=10情况下DDS输出频谱图 业出版社,2001. 米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米米 (上接第l7页) 【7] G llle ̄。A K A l’暗La1 An 0p血lal all;0 m for nnn:llal [3]P Q娜 .0I嚏iI all莹or for 0II in蚺l卜 exclusion in c0咖pIller璀蚋 [J].0瑚mm 0l1s 0f the 伽棚 伽衄 吐er联 [J].0硼叩Ilter G嘲mluIlic砒 l坞。 ACM。1981。24(1):9-17. 1991,14(1O):627-632. [8] Mer-db。Ofiveim I ̄uto.Yet AI Di曲 l吐ed [4] K m .A缸ee.b88ed for diE bIlled M1血】Ial E】【du ∞ 血Im[Z].Bn ̄ih l/niver ̄de Fe ̄al 0n[J].AcM Tram.On c0n lter sy蛐e髓,1989,7 de^liIlas ̄1"i8。1996. (1):61-77. [9] 刘丹。刘心松,等.基于读写特征的分布式互斥算法[J]. [5]L 呷ort.11le 0n印 锄[J].JoL ̄ml 0fthe 电子学报。2004。32(2):326—329. ACM,1986,33(2):313-326. [10]王昕哗。叶慧娟,李新.分布式系统设计中的互斥问题 [6]M M 【姗.An al80ritlm for 0II in [J].海军工程大学学报,2003。15(5):1115 112. 灯舶缸ali叠ed 8y鼬咖略[J].ACMTram.饥0m s)耐衄皓, [11] 邢雁.关于分布式互斥系统进程互斥算法的研究[J]. l985,3(2):l4孓159. 自然科学报,21103。19(1):43—46. 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- gamedaodao.com 版权所有 湘ICP备2022005869号-6

违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务