CN109358698B - 一种基于复合频率控制字的直接数字频率合成方法及装置 - Google Patents
一种基于复合频率控制字的直接数字频率合成方法及装置 Download PDFInfo
- Publication number
- CN109358698B CN109358698B CN201811002531.9A CN201811002531A CN109358698B CN 109358698 B CN109358698 B CN 109358698B CN 201811002531 A CN201811002531 A CN 201811002531A CN 109358698 B CN109358698 B CN 109358698B
- Authority
- CN
- China
- Prior art keywords
- frequency
- value
- module
- analog signal
- digital
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 38
- 238000001308 synthesis method Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 50
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 44
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 44
- 238000006243 chemical reaction Methods 0.000 claims abstract description 34
- 238000001914 filtration Methods 0.000 claims abstract description 30
- 230000003321 amplification Effects 0.000 claims description 22
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 22
- 230000007547 defect Effects 0.000 abstract description 3
- 230000002194 synthesizing effect Effects 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 229910052701 rubidium Inorganic materials 0.000 description 3
- IGLNJRXAVVLDKE-UHFFFAOYSA-N rubidium atom Chemical compound [Rb] IGLNJRXAVVLDKE-UHFFFAOYSA-N 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F1/00—Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
- G06F1/02—Digital function generators
- G06F1/022—Waveform generators, i.e. devices for generating periodical functions of time, e.g. direct digital synthesizers
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
本发明公开了一种基于复合频率控制字的直接数字频率合成方法和装置,根据原始频率Fsysclk以及目标合成频率Fout,获得复合频率控制字FTW={X,A,B},根据原始频率Fsysclk对所述数值X进行递归求和,获得求和输出值,每获得B个求和输出值与数值A相加,再进行截取,获得最终相位地址值;将所述的最终相位地址值进行相幅转换,获得幅值数字信号;将所述的幅值数字信号转换为模拟信号,对所述的模拟信号进行滤波后再放大,输出具有目标合成频率Fout的模拟信号;本发明提供的方法和装置克服了传统DDS的相位累加器模数为固定值,频率控制字也必须为正整数,只能合成有限个精确频率信号的缺陷,在任意频点获得了更高的输出准确度。
Description
技术领域
本发明涉及频率合成方法及装置,具体涉及一种基于复合频率控制字的直接数字频率合成方法及装置。
背景技术
频率合成是指由单个或若干个具有较高稳定度和精确度的参考频率信号源通过频率域的线性运算,从而生成具有大致相同稳定度和精确度的其他离散频率的过程。
频率合成器的实现方法大致有三类:直接模拟频率合成方法、锁相频率合成方法和直接数字频率合成方法。直接数字频率合成(DDS)方法从相位的概念出发进行频率合成,采用了数字采样存储技术。由于DDS是一个开环***,无任何反馈环节,因而频率转换时间极短。此外该方法可全数字化实现、便于集成、体积小、重量轻。频率合成器设计中采用DDS方法在当前的电子测量测试领域已经逐渐成为一种主流的做法。
现有技术中的DDS信号产生方法以AT89C51单片机作为微处理器,通过串行端口向DDS数字芯片AD9910发送配置参数。AD9910根据配置参数给定的工作模式和数据,由内部电路计算出其内核工作参数并进行频率合成,然后输出数字DDS信号至其内部集成的14位DAC进行模数转换,最后输出模拟DDS信号。采用专用DDS数字芯片实现频率合成器时,由于专用DDS数字芯片的技术局限性,致使该类频率合成器存在不可避免的性能缺陷:专用DDS数字芯片的相位累加器位数固定,缺乏灵活性;相位累加器位数较小并且合成频率分辨率较低、频点较少,不能保证频率合成器输出的准确度。
现有技术中采用FPGA实现DDS信号发生器主要由FPGA数字电路构造的相位累加器、波形存储器和数字乘法器构成。相位累加器以频率控制字循环累加得到相码,频率控制字由主控单元配置,频率控制字的大小决定输出频率值。波形存储器存储一个周期输出波形的样点。相位累加器的相码作为波形存储器的读地址,从中取出波形样点构成数字波形。波形存储器由FPGA内部的ROM实现。然后经由一个数字乘法器对波形样点乘上一个幅度补偿系数,最后送DAC模块将数字波形转化为模拟波形输出。该方法可以通过主控单元配置浮动补偿系数以及DAC的参考电压值实现幅度灵活精确可控,但其DDS本质与AD9910专用DDS芯片一致,都是采用相位累加器以时钟速率递归求和频率控制字,借助查找表的方式合成频率,使得合成频率分辨率较低、频点较少,不能保证频率合成器输出的准确度。
传统DDS方法的原理如图1所示,相位累加器(位宽为C)输出地址值,该值会线性增大,直到累加器的最大值2C并翻转。具体电路实施时为了节省资源,通常会舍去一定位数,保留到P位(仅使用MSB),以便减小角度振幅转换模块的复杂程度。累加器的输出序列范围0到2P-1映射到单位圆的一圈上,即将0到2P-1范围的二进制值映射到0到2C范围内的弧度角,完成角度至幅度的转换。由转换器输出的Y位数字幅度数字信号再经过一个Y位的DAC芯片转换为模拟信号,最后经过低通滤波器得到所需要的正弦信号。
现有技术中的直接数字频率合成方法,其输出频率公式:
其中,Fout为目标合成频率,Fsysclk为原始频率,一般地,原始频率为时钟模块的采样频率。2C为相位累加器的模数。FTW'为频率控制字,其值为正整数,通常情况下FTW'<2C-1,否则会合成一个奈奎斯特镜像频率。输出频率因此限制在一定的范围之内:
由式II可知,DDS相位累加器模数2C决定DDS的频率分辨率,同时也决定了可能的输出频率个数,即2C-1个。由于相位累加器位宽C为固定整数值,频率控制字又不能为小数,所以输出频率不能是任意值。由式I可见,输出频率和***时钟必须满足下列关系:
针对等于采样速率某一整数约数的输出频率Fout=Fsysclk/Q,Q为整数,代入式III可得到频率控制字FTW':
因为FTW'和Q必须均为整数,所以唯一满足式IV的Q值可以表示为2K,在实际应用中频率控制字可以表示为:
其中C、K皆为整数。
那么,在具体的DDS实现时,存在两种情况:
1、当频率控制字FTW'满足式V时,即FTW'为一个整数,DDS可以输出精确频率。例如,Fsysclk=100MHz,C为32,期望获得Fout=6.25MHz,那么Q=24,由式V计算可得的频率控制字是整数,即FTW'=268435456,此时的相位累加器在累加完成一个循环后地址涵盖了0到2C范围,下一次循环回到初始值,可以输出精确频率值。
2、当频率控制字FTW'不满足式V时,即FTW'为一个非整数,实际累加时又只能取整数,一个周期以后不能回到初始点,所以在一定的范围内每个周期的采样点不同,致使波形幅值不同,波形不稳定。同时,累加一个周期后有相位损失,周期增大,频率减小,出现误差。例如Fsysclk=100MHz,C为32,期望获得Fout=10MHz,那么Q=10,由式V求得FTW'=429496729.6,即表示这种情况下传统DDS方法不能产生准确的10MHz。实际应用中一般将频率控制字四舍五入取整,即FTW'=429496730,此时合成的实际频率为10.00000000931322574615478515625MHz,比期望值大了约9.31mHz。在某些高精度的测试测量领域,如果频率合成器作为标准频率源来使用,这种mHz级别的误差是无法忍受的。
发明内容
本发明的目的在于提供一种基于复合频率控制字的直接数字频率合成方法及装置,用以解决现有技术中直接数字频率合成方法合成频率分辨率较低、频点较少,准确度不高等问题。
为了实现上述任务,本发明采用以下技术方案:
一种基于复合频率控制字的直接数字频率合成方法,用于根据原始频率Fsysclk获得目标合成频率Fout的模拟信号,原始频率Fsysclk的单位为MHz,目标合成频率Fout的单位为MHz,所述的方法包括:
步骤1:获得一组复合频率控制字FTW={X,A,B},其中数值X、数值A、数值B均为整数,且X>0,A≥0,B>0;
当A>0时,A与B为互质数;
当A=0时,B为任意正整数;
步骤2:以原始频率Fsysclk循环执行N'次步骤2.1-2.5,获得一个周期的目标合成频率Fout的模拟信号,N'为将约至M=1时N的值,且N'为正整数,对于第n次执行步骤2.1-2.5,n={1,2,…,N'};
所述的步骤2.1-2.5按照以下步骤进行:
步骤2.1:将初始数值与数值X相加,获得求和中间值,其中当n=1时,初始数值为0;
步骤2.2:判断是否为整数,若为整数,则将求和中间值与数值A相加,获得最终求和值,否则直接将求和中间值作为最终求和值;将所述的最终求和值作为初始相位地址值,再对所述的初始相位地址值进行截取,获得最终相位地址值;
步骤2.3:将所述的最终相位地址值进行相幅转换,获得幅值数字信号;
步骤2.5:若n≤N',则将本次循环时由步骤2.2获得的最终求和值作为初始数值,令n=n+1,执行步骤2.1;否则停止。
本发明具有如下技术特征:
进一步地,在所述的步骤2.2中,对初始相位地址值进行截取时,将初始相位地址值的低k位进行删除,获得最终相位地址值,k为正整数。
进一步地,所述的步骤2.4中,对所述的模拟信号进行滤波时,采用七阶椭圆低通滤波的方法。
本发明还保护一种基于复合频率控制字的直接数字频率合成装置,用于实现以上所述的基于复合频率控制字的直接数字频率合成方法,所述的装置包括控制字获得模块和模拟信号获得模块;
所述的控制字获得模块与模拟信号获得模块连接,用于获得一组复合频率控制字FTW={X,A,B},其中数值X、数值A、数值B均为整数,且X>0,A≥0,B>0;
当A>0时,A与B为互质数;
当A=0时,B为任意正整数;
所述的模拟信号获得模块,用于以原始频率Fsysclk循环执行N'次步骤2.1-2.5,获得一个周期的目标合成频率Fout的模拟信号,N'为将约至M=1时N的值,且N'为正整数,所述的模拟信号获得模块包括相位累加器模块、地址生成器模块,相幅转换器模块以及数模转换、滤波及放大模块;
所述的相位累加器模块与所述的地址生成器模块连接,用于将初始数值与数值X相加,获得求和中间值,其中当n=1时,初始数值为0,所述相位累加器模块的位宽为C;
所述的地址生成器模块与所述的相幅转换器模块连接,用于判断是否为整数,若为整数,则将求和中间值与数值A相加,获得最终求和值,否则将求和中间值作为最终求和值;将所述的最终求和值作为初始相位地址值,再对所述的初始相位地址值进行截取,获得最终相位地址值;
所述的相幅转换器模块与所述的数模转换、滤波及放大模块连接,用于对所述的最终相位地址值进行相幅转换,获得幅值数字信号;
进一步地,所述的装置还包括参考时钟模块和辅助累加器模块;
所述的参考时钟模块与所述的相位累加器模块连接,用于提供原始频率Fsysclk;
所述的辅助累加器模块与所述的地址生成器模块连接,用于存储复合频率控制字中的数值A与数值B。
进一步地,所述的数模转换、滤波及放大模块包括数模转换模块,滤波器模块以及放大驱动模块;
所述的数模转换模块与所述的滤波器模块连接,用于将所述的幅值数字信号转换为模拟信号;
所述的滤波器模块与所述的放大驱动模块连接,用于滤除数模转换模块输出的模拟信号中的杂波和其他干扰信号;
进一步地,所述的模拟信号获得模块由FPGA模块以及数模转换、滤波及放大模块组成。
本发明与现有技术相比具有以下技术特点:
1、本发明提出的基于复合频率控制字的直接数字频率合成方法中复合频率控制字灵活多变,克服了传统DDS的相位累加器模数为固定值,频率控制字也必须为正整数,只能合成有限个精确频率信号的缺陷,在任意频点获得了更高的输出准确度;
2、本发明提出的基于复合频率控制字的直接数字频率合成装置,充分利用了FPGA在数字电路设计中灵活性以及在并行运算处理体系结构设计中的优势,采用FPGA实现采用复合频率控制字的DDS方法,使***装置具备更灵活的体系结构和更快的处理速度。
附图说明
图1为传统DDS方法的硬件实现原理示意图;
图2为本发明提供的DDS方法的流程图;
图3为本发明的一个实施例中提供的用于实现本发明提供的DDS方法的硬件原理示意图;
图4为本发明提供的测试平台结构示意图;
图5为本发明提供的DDS方法与传统DDS方法在输出1.024MHz的模拟信号时的相对频率差示意图;
图6为本发明提供的DDS方法与传统DDS方法在输出2.048MHz的模拟信号时的相对频率差示意图;
图7为本发明提供的DDS方法与传统DDS方法在输出5MHz的模拟信号时的相对频率差示意图;
图8为本发明提供的DDS方法与传统DDS方法在输出10MHz的模拟信号时的相对频率差示意图。
具体实施方式
以下是发明人提供的具体实施例,以对本发明的技术方案作进一步解释说明。
实施例一:
由于传统DDS的相位累加器模数为固定值,频率控制字也必须为正整数,导致传统DDS方法只能合成有限个数的精确频率信号,因此本实施例给出一种基于复合频率控制字的直接数字频率合成方法,根据原始频率Fsysclk合成输出目标合成频率Fout的模拟信号。一般来说,原始频率Fsysclk也就是参考时钟频率,由高稳晶振或原子钟输出的低频标准频率信号倍频获得,用于方法实现时控制处理步骤和节奏,也称为采样频率。
如图2所示,所述的方法包括:
步骤1:获得复合频率控制字FTW={X,A,B},其中数值X、数值A、数值B均为整数,且X>0,A≥0,B>0;
当A>0时,A与B为互质数;
当A=0时,B为任意正整数;
其中,C为相位累加器位宽参数,其取值可以是16、32等。
在本步骤中,采用复合频率控制字的DDS方法输出频率计算公式:
目标合成频率Fout和原始频率Fsysclk确定了M和N的关系,即等于约至最小项的分数原始频率Fsysclk为了符合奈奎斯特采样要求,与传统的DDS方法相比较,N不要求是2的幂的形式,它可以是任意正整数。
式VI可进一步描述为:
其中,M、N和相位累加器的模数(2C)通过长除法确定了与X、Y的关系。X通过取整运算获取,Y为余数,即:
进一步将Y/N化简得到:
式IX中,A和B为整数且互为质数。采用复合频率控制字的DDS方法中,A和B送辅助累加器进行运算处理。
将约至M=1时,得到相应的N,即N'是地址溢出一次经过的实际采样次数或者实际采样次数的附近值。假设N′为一个地址溢出一次实际采样次数,则N′等于N或者N的附近值,即一个周期中N′次采样完成一次从0到地址最大值的累加,当N′次采样后溢出进行下一次循环,如果每个周期的所有的采样点都是重合的点,则A为0,B可以为任意正整数,得到FTW={X,0,B},采用复合频率控制字的精确DDS可以精确输出频率。
如果每个周期的采样点不能重合,即A≠0时,则不能精确输出期望的频率,此时:
如果化简式IX可得到B<N′,且B值是一个小于10的数字,那么复合频率控制字为FTW={X,A,B},一个周期采样后的余数分成均匀的B个等份,均匀的***到整个采样过程当中,每次***A值,使一个周期内的采样覆盖0到2C范围,下一次循环的所有采样点与第一次采样点重合,保证每一个周期的采样点都完全相同,DDS可输出准确的频率。
如果化简式IX不满足B<N′,即Y/N除不尽的情况下,那么将式IX近似化简到B为小于10的一个整数,同时得到一个A值,此时复合频率控制字为FTW′={X,A,B},可保证DDS合成频率更加接近期望值。
在本实施例中,原始频率Fsysclk=100MHz,C=32,采用复合频率控制字的DDS方法合成频率为1.024MHz的信号,复合频率控制字各分量的详细计算过程如下:
Y=M2C-XN=32×232-43980465×3125=347;
可得到复合频率控制字:FTW={43980465,1,9}。
在本实施例中,原始频率Fsysclk=100MHz,C=32,利用本发明的方法合成频率2.048MHz的复合频率控制字为:FTW={87960930,2,9},即X=87960930,A=2,B=9。
在本实施例中,原始频率Fsysclk=100MHz,C=32,利用本发明的方法合成频率5MHz的复合频率控制字为FTW={214748364,4,5},即X=214748364,A=4,B=5。
在本实施例中,原始频率Fsysclk=100MHz,C=32,利用本发明的方法合成频率10MHz的复合频率控制字为:FTW′={429496729,3,5},即X=429496729,A=3,B=5。
在本步骤中,在原始频率Fsysclk下循环执行N'次步骤2.1-2.5,也就是说,在原始频率的作用下每来一个脉冲就执行一次步骤2.1-2.5,下一次脉冲来时,以上次获得的结果作为初始值再循环进行,脉冲是由原始频率进行控制的,当原始频率为1Hz,就是指每秒都会有一个脉冲信号,当有这个脉冲信号的时候,执行一次步骤2.1-2.5。
因此在本发明中原始频率Fsysclk的作用包括:计算数值X、保持整个方法中所有步骤的统一,利用原始频率作为时钟频率,控制循环的间隔时间,同步电路中时钟的基础频率,控制所有步骤的统一性。
在本步骤中,每执行1次步骤2.1-2.5会输出一段模拟信号,由于模拟信号周期性的正弦波的形式,因此每执行1次步骤2.1-2.5输出的是段正弦波的模拟信号,当原始频率为100MHz时,每获得一个周期的10MHz模拟信号,需要循环执行10次步骤2.1-2.5;每获得一个周期的20MHz模拟信号,需要循环执行5次步骤2.1-2.5;每获得一个周期的5MHz模拟信号,需要循环执行20次步骤2.1-2.5。
所述的步骤2.1-2.5按照以下步骤进行:
步骤2.1:将初始数值与数值X相加,获得求和中间值,其中当n=1时,初始数值为0;
在本实施例中,目标合成频率Fout=10MHz,当n=1时,X=429496729,求和中间值为429496729;n=2时,求和中间值为2*X=858993458。
步骤2.2:判断是否为整数,若为整数,则将求和中间值与数值A相加,获得最终求和值,否则将求和中间值作为最终求和值;将所述的最终求和值作为初始相位地址值,再对所述的初始相位地址值进行截取,获得最终相位地址值;
在本步骤中,当初始相位地址值的位数超过地址生成器的最大位宽时,地址生成器会自动溢出清零,即将初始相位地址值清0,也就是说,将最终求和值清0作为初始数值进入下一次循环中。
可选地,在所述的步骤2.2中,对初始相位地址值进行截取时,将初始相位地址值的低k位进行删除,获得最终相位地址值,k为正整数。
在本实施例中,目标合成频率Fout=10MHz,n=2时,初始相位地址值为858993458;n=3时,初始相位地址值为1288490187。n=2和n=3时的初始相位地址进行相同低四位的截断,即k=4,则n=2时,最终相位地址值为85899;n=3时,最终相位地址值为128849。
若初始相位地址值的位数小于k,则将初始相位地址值补0。
步骤2.3:将所述的最终相位地址值进行相幅转换,获得幅值数字信号;
可选地,所述的步骤2.4中,对所述的模拟信号进行滤波时,采用七阶椭圆低通滤波的方法。
步骤2.5:若n≤N',则将本次循环时由步骤2.2获得的最终求和值作为初始数值,令n=n+1,执行步骤2.1,否则停止。
经过N'次循环执行步骤2.1-2.5后,获得了一个周期的目标合成频率Fout的模拟信号,循环的次数根据目标合成频率的模拟信号的长度确定,若需要的模拟信号长度包括多个周期,那么就可以循环多个N'次的步骤2.1-2.5,获得符合长度要求的目标合成频率的模拟信号。
实施例二:
本实施例给出一种基于复合频率控制字的直接数字频率合成装置,用于实现实施例一中所述的基于复合频率控制字的直接数字频率合成方法,如图3所示,所述的装置包括控制字获得模块和模拟信号获得模块。
所述的控制字获得模块与相位累加器模块连接,用于获得一组复合频率控制字FTW={X,A,B},其中数值X、数值A、数值B均为整数,且X>0,A≥0,B>0;
当A>0时,A与B为互质数;
当A=0时,B为任意正整数。
所述的模拟信号获得模块,用于以原始频率Fsysclk循环执行N'次步骤2.1-2.5,获得一个周期的目标合成频率Fout的模拟信号,N'为将约至M=1时N的值,且N'为正整数,所述的模拟信号获得模块包括相位累加器模块、地址生成器模块,相幅转换器模块以及数模转换、滤波及放大模块;
所述的相位累加器模块与所述的地址生成器模块连接,用于将初始数值与数值X相加,获得求和中间值,其中当n=1时,初始数值为0,所述相位累加器模块的位宽为C;
所述的地址生成器模块与所述的相幅转换器模块连接,用于判断是否为整数,若为整数,则将求和中间值与数值A相加,获得最终求和值,否则将求和中间值作为最终求和值;将所述的最终求和值作为初始相位地址值,再对所述的初始相位地址值进行截取,获得最终相位地址值;
所述的相幅转换器模块与所述的数模转换、滤波及放大模块连接,用于对所述的最终相位地址值进行相幅转换,获得幅值数字信号;
可选地,所述的装置还包括参考时钟模块,辅助累加器模块;
所述的参考时钟模块与所述的相位累加器模块连接,用于提供原始频率Fsysclk;
在本实施例中,参考时钟模块采用SRS公司的铷原子钟PRS10。铷原子钟PRS10的输出频率为10MHz正弦信号(VPP=1.41V),经过倍频到100MHz作为原始频率。
所述的辅助累加器模块与所述的地址生成器模块连接,用于存储复合频率控制字中的数值A与数值B;
可选地,所述的数模转换、滤波及放大模块包括数模转换模块,滤波器模块以及放大驱动模块;
所述的数模转换模块与所述的滤波器模块连接,用于将所述的幅值数据转换为模拟信号;
在本实施例中,数模转换模块Maxim公司的16位并行输入DAC芯片MAX5885,将EP1C12Q240I7输出的数字波形转化为模拟波形输出。
所述的滤波器模块与所述的放大驱动模块连接,用于滤除数模转换模块输出的模拟信号中的杂波和其他干扰信号;
在本实施例中,如图3所示,LPF为低通滤波器,为了滤除掉DAC模块输出波形中的杂波和其他干扰信号,采用了七阶椭圆滤波器滤波,提高合成波形的质量。
所述的放大驱动模块用于提高直接数字频率合成装置的驱动能力。
在本实施例中,采用TI公司的超低噪声集成运放芯片OPA847设计实现,使本发明装置具备较大的驱动能力,50Ω负载条件下,输出合成频率信号的幅度不小于3Vpp。
作为一种优选的实施方式,所述的模拟信号获得模块由FPGA模块以及数模转换、滤波及放大模块组成。
在本实施例中,模拟信号获得模块由FPGA模块加上数模转换、滤波及放大模块组成,FPGA模块选用Altera公司的飓风系列的EP1C12Q240I7,在EP1C12Q240I7中设计实现了相位累加器、辅助累加器、地址生成器、相幅转换器等数字逻辑电路,在100MHz原始频率的驱动下并行实现采用复合频率控制字的DDS方法;并且,控制字获得模块的功能也可以由FPGA模块进行实现。
实施例三
在本实施例中,对本发明提供的一种基于复合频率控制字的直接数字频率合成方法及装置的合成准确度进行测试,通过对比本发明的方法和装置与传统DDS方法和装置合成频率的准确度,测试平台如图4所示,图4中标准频率源用于提供原始频率即采样频率,分配放大器将原始频率分别送至DDS的FPGA平台以及通用频率计数器,DDS的PFGA实现平台可以是本发明提供的用于实现本发明提供的DDS方法的PFGA平台,也可以是传统的DDS方法的FPGA实现平台,通用频率计数器的作用是获得相对频率差,TimeLab软件和Stable32软件用于绘制测量结果图,测量结果如图5-8所示,图5-8中横坐标表示测量的平均时间,纵坐标表示相对频率差,在同样采用SRS公司的铷原子钟PRS10作为时钟参考的前提下,本发明合成频率相对频率偏差约为2.00E-12相对于传统DDS方法(相对频率偏差约为2.00E-9)降低了大约3个数量级,可证明本发明具备更高的频率输出准确度。
Claims (7)
1.一种基于复合频率控制字的直接数字频率合成方法,用于根据原始频率Fsysclk获得一个周期的目标合成频率Fout的模拟信号,原始频率Fsysclk的单位为MHz,目标合成频率Fout的单位为MHz,其特征在于,所述的方法包括:
步骤1:获得一组复合频率控制字FTW={X,A,B},其中数值X、数值A、数值B均为整数,且X>0,A≥0,B>0;
当A=0时,B为任意正整数;
步骤2:以原始频率Fsysclk循环执行N'次步骤2.1-2.5,获得一个周期的目标合成频率Fout的模拟信号,N'为将约至M=1时N的值,且N'为正整数,对于第n次执行步骤2.1-2.5,n={1,2,…,N'};
所述的步骤2.1-2.5按照以下步骤进行:
步骤2.1:将初始数值与数值X相加,获得求和中间值,其中当n=1时,初始数值为0;
步骤2.2:判断是否为整数,若为整数,则将求和中间值与数值A相加,获得最终求和值,否则直接将求和中间值作为最终求和值;将所述的最终求和值作为初始相位地址值,再对所述的初始相位地址值进行截取,获得最终相位地址值;
步骤2.3:将所述的最终相位地址值进行相幅转换,获得幅值数字信号;
步骤2.5:若n≤N',则将本次循环时由步骤2.2获得的最终求和值作为初始数值,令n=n+1,执行步骤2.1;否则停止。
2.如权利要求1所述基于复合频率控制字的直接数字频率合成方法,其特征在于,在所述的步骤2.2中,对初始相位地址值进行截取时,将初始相位地址值的低k位进行删除,获得最终相位地址值,k为正整数。
3.如权利要求1所述基于复合频率控制字的直接数字频率合成方法,其特征在于,所述的步骤2.4中,对所述的模拟信号进行滤波时,采用七阶椭圆低通滤波的方法。
4.一种基于复合频率控制字的直接数字频率合成装置,其特征在于,用于实现权利要求1-3任一项权利要求所述的基于复合频率控制字的直接数字频率合成方法,所述的装置包括控制字获得模块和模拟信号获得模块;
所述的控制字获得模块与模拟信号获得模块连接,用于获得一组复合频率控制字FTW={X,A,B},其中数值X、数值A、数值B均为整数,且X>0,A≥0,B>0;
当A=0时,B为任意正整数;
所述的模拟信号获得模块,用于以原始频率Fsysclk循环执行N'次步骤2.1-2.5,获得一个周期的目标合成频率Fout的模拟信号,N'为将约至M=1时N的值,且N'为正整数,所述的模拟信号获得模块包括相位累加器模块、地址生成器模块,相幅转换器模块以及数模转换、滤波及放大模块;
所述的相位累加器模块与所述的地址生成器模块连接,用于将初始数值与数值X相加,获得求和中间值,其中当n=1时,初始数值为0,所述相位累加器模块的位宽为C;
所述的地址生成器模块与所述的相幅转换器模块连接,用于判断是否为整数,若为整数,则将求和中间值与数值A相加,获得最终求和值,否则将求和中间值作为最终求和值;将所述的最终求和值作为初始相位地址值,再对所述的初始相位地址值进行截取,获得最终相位地址值;
所述的相幅转换器模块与所述的数模转换、滤波及放大模块连接,用于对所述的最终相位地址值进行相幅转换,获得幅值数字信号;
5.如权利要求4所述的基于复合频率控制字的直接数字频率合成装置,其特征在于,所述的装置还包括参考时钟模块和辅助累加器模块;
所述的参考时钟模块与所述的相位累加器模块连接,用于提供原始频率Fsysclk;
所述的辅助累加器模块与所述的地址生成器模块连接,用于存储复合频率控制字中的数值A与数值B。
7.如权利要求6所述的基于复合频率控制字的直接数字频率合成装置,其特征在于,所述的模拟信号获得模块由FPGA模块以及数模转换、滤波及放大模块组成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811002531.9A CN109358698B (zh) | 2018-08-30 | 2018-08-30 | 一种基于复合频率控制字的直接数字频率合成方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811002531.9A CN109358698B (zh) | 2018-08-30 | 2018-08-30 | 一种基于复合频率控制字的直接数字频率合成方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109358698A CN109358698A (zh) | 2019-02-19 |
CN109358698B true CN109358698B (zh) | 2020-03-20 |
Family
ID=65350162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811002531.9A Expired - Fee Related CN109358698B (zh) | 2018-08-30 | 2018-08-30 | 一种基于复合频率控制字的直接数字频率合成方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109358698B (zh) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7606849B2 (en) * | 2004-05-10 | 2009-10-20 | Advantest Corporation | Method and apparatus for improving the frequency resolution of a direct digital synthesizer |
US9547327B2 (en) * | 2006-11-09 | 2017-01-17 | Alan N. Willson, Jr. | Excess-fours processing in direct digital synthesizer implementations |
JP4369497B2 (ja) * | 2007-05-17 | 2009-11-18 | 株式会社日立製作所 | ダイレクト・デジタル・シンセサイザおよびそれを用いた核磁気共鳴装置 |
CN102571036B (zh) * | 2011-12-07 | 2013-02-27 | 中国电子科技集团公司第十研究所 | 用dds产生任意精准频率的方法 |
CN104300978A (zh) * | 2014-10-13 | 2015-01-21 | 西安电子工程研究所 | 一种线性度可控的线性调频信号产生方法 |
CN104753530B (zh) * | 2015-04-16 | 2017-07-04 | 电子科技大学 | Dds中相位修正及非均匀相幅转换方法及装置 |
CN104753502B (zh) * | 2015-04-23 | 2017-08-29 | 成都理工大学 | 一种基于fpga的dds信号发生器及其实现方法 |
CN105846819A (zh) * | 2016-03-23 | 2016-08-10 | 上海航天测控通信研究所 | 一种基于fpga的直接数字频率合成方法及合成器 |
-
2018
- 2018-08-30 CN CN201811002531.9A patent/CN109358698B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN109358698A (zh) | 2019-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1873612A1 (en) | Phase-coherent signal generator | |
CN1677870B (zh) | 利用谐波抵消的线性补偿电路 | |
US3882403A (en) | Digital frequency synthesizer | |
US5774084A (en) | Method and apparatus for translating digital data into an analog signal | |
JPH0548003B2 (zh) | ||
JPS5931897B2 (ja) | 周波数合成装置 | |
US5673212A (en) | Method and apparatus for numerically controlled oscillator with partitioned phase accumulator | |
CN103117731A (zh) | 一种基于fpga的正弦信号发生器实现方法及实现装置 | |
JPH03253108A (ja) | ダイレクト・デジタル・シンセサイザー及び信号発生方法 | |
US20050146360A1 (en) | Multi-stage numeric counter oscillator | |
US4185247A (en) | Means for reducing spurious frequencies in a direct frequency synthesizer | |
EP0712211B1 (en) | Phase modulation having individual placed edges | |
KR101341138B1 (ko) | 에일리어싱된 주파수에 대한 위상 고정 | |
CN109358698B (zh) | 一种基于复合频率控制字的直接数字频率合成方法及装置 | |
CN208752198U (zh) | 一种带反馈校正的多通道超声波任意波形信号发生器 | |
CN113985361A (zh) | 一种改善dds输出信号频谱质量的***及方法 | |
CN110768665B (zh) | 一种二倍时钟采样速率的dds信号扫频源*** | |
CN112073064A (zh) | 基于自适应校正的宽带信号生成装置及方法 | |
US7209937B2 (en) | Method and apparatus for generation of arbitrary mono-cycle waveforms | |
Tang et al. | A suppressing method for spur caused by amplitude quantization in DDS | |
CN107943204B (zh) | 数字频率合成方法以及装置 | |
KR20000035114A (ko) | 직접 디지털 신시사이저 | |
Shan et al. | Design and implementation of a FPGA-based direct digital synthesizer | |
CN108923770A (zh) | 一种三角波信号产生方法 | |
CN111399587B (zh) | 一种基于梯形波的超声脉冲激励信号合成方法及其*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200320 |
|
CF01 | Termination of patent right due to non-payment of annual fee |