CN104391174B - 一种基于数字信号处理器平台的微分测频方法及系统 - Google Patents

一种基于数字信号处理器平台的微分测频方法及系统 Download PDF

Info

Publication number
CN104391174B
CN104391174B CN201410676742.6A CN201410676742A CN104391174B CN 104391174 B CN104391174 B CN 104391174B CN 201410676742 A CN201410676742 A CN 201410676742A CN 104391174 B CN104391174 B CN 104391174B
Authority
CN
China
Prior art keywords
differentiator
module
input
output end
signal
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.)
Active
Application number
CN201410676742.6A
Other languages
English (en)
Other versions
CN104391174A (zh
Inventor
夏国明
林晨
施芹
裘安萍
苏岩
丁衡高
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201410676742.6A priority Critical patent/CN104391174B/zh
Publication of CN104391174A publication Critical patent/CN104391174A/zh
Application granted granted Critical
Publication of CN104391174B publication Critical patent/CN104391174B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

本发明提供一种基于数字信号处理器平台的微分测频方法,包括以下步骤:步骤1,提供微分测频系统;步骤2,带通滤波器对输入的谐振信号进行滤波;步骤3,裂相器将滤波后的信号分成两路信号,一路信号进行希尔伯特变换,一路信号进行低通滤波;步骤4,第一微分器模块与第二微分器模块将步骤2中的两路信号分别进行两次微分;步骤5,延迟步骤3中的两路信号和步骤4中第一次微分后的两路信号;步骤6,利用交叉相乘法或绝对值法得到信号的频率值的平方;步骤7,低通滤波器对频率值的平方进行滤波后输出,并将输出信号反馈于带通滤波器。本发明提供的微分测频方法稳定性好、抗噪能力强、精度高、灵敏度高、易于实现。

Description

一种基于数字信号处理器平台的微分测频方法及系统
技术领域
本发明涉及数字信号处理与硬件编程领域,特别涉及一种对信号微分测频的算法在数字信号处理器平台上的实现。
背景技术
随着我国生产和科学技术事业的发展,许多电参量的测量方案、测量结果都与频率有着非常密切的关系,因此频率的测量也尤为重要。近年来,以频率作为输出信号的MEMS谐振式惯性器件的精度都在不断提高,因此相匹配的高精度低噪声的频率测量技术也成为了各国专家研究对象。
目前国内外主要的频率测量电路的基本原理有三种:(1)通过对被测信号进行傅里叶变换测量周期信号的频率,这种方法具有良好的抗噪性,但在傅里叶变换过程中要进行时域的截断,这会造成频谱的泄露,导致测量结果的失真,精准度较差。(2)在标准的闸门时间内对被测频率信号的周期个数进行计数而得出被测频率值,这种方法测量方便、读数直接,但是计数所导致的量化误差,限制测频精度的提高;尽管有改进的多周期同步法和延迟链法对量化误差进行限制,但是多周期同步法是以牺牲系统带宽为代价,而延迟链法对延迟单元的精度要求极高,实现难度较大,而且存在延迟链长度分布不均以及延迟抖动等问题,使其实际精度远低于理论值。(3)基于锁相环的测频电路,利用压控振荡器(VCO)控制信号与频率成正比的特性实现频率识别;锁相环测频法优点在于灵敏度与精度高,但是锁相环是一个闭环系统,如果频率信号变化过大过快,会出现失锁现象,锁相环测频不具有良好的稳定性。2013年,南京理工大学的夏国明提出了一种sigma-delta PLL频率测量电路及方法(夏国明,裘安萍,施芹,石然,苏岩,丁衡高。sigma-delta PLL频率测量电路及方法。公开号:CN103487648A),该方法分辨率高、易于实现,在数字部分抗干扰性好,但是,模拟部分信号会受到外界的干扰,降低精度,且sigma-delta不具有较高的稳定性。
综上所述,目前普遍使用的频率测量方法无法满足精度与稳定性的双方面要求,且无法消除信号中包含的倍频项,难以适应新MEMS谐振式惯性器件的高精度频率读取的要求。
发明内容
为了克服现有技术的不足,本发明提供一种稳定性好、抗噪能力强、精度高、灵敏度高、易于实现的一种基于数字信号处理器平台利用微分算法的可消除倍频项影响的高精度MEMS谐振式惯性器件频率测量方法及系统。
一种基于数字信号处理器平台的微分测频方法,包括以下步骤:
步骤1,提供微分测频系统,包括带通滤波器、裂相器、第一微分器模块、第二微分器模块、第一延迟模块、第二延迟模块、第三延迟模块、第四延迟模块、交叉相乘法或绝对值法推导频率算法模块、低通滤波器;
步骤2,带通滤波器对MEMS惯性器件的谐振信号进行滤波,得到信号Asin(ωt);
步骤3,裂相器将滤波后的信号分成两路信号,一路信号进行希尔伯特变换,一路信号进行低通滤波;
步骤4,第一微分器模块与第二微分器模块将步骤2中的两路信号分别进行一次微分和两次微分;
步骤5,第一延迟模块、第二延迟模块、第三延迟模块、第四延迟模块分别延迟步骤3中的两路信号和步骤4中第一次微分后的两路信号;
步骤6,交叉相乘法或绝对法值推导频率算法模块利用交叉相乘法或绝对值法得到信号的频率值的平方;
步骤7,低通滤波器对频率值的平方进行滤波后输出,并将输出信号反馈于带通滤波器。
一种基于数字信号处理器平台的微分测频系统,包括带通滤波器、裂相器、第一低通滤波器、希尔伯特变换器、第一微分器模块、第二微分器模块、第一延迟模块、第二延迟模块、第三延迟模块、第四延迟模块、交叉相乘法或绝对值法推导频率算法模块、第二低通滤波器;所述带通滤波器以反馈回来的输出信号为控制量,选择不同的参数,调节带通滤波器的带宽,滤除输入信号中的高倍频项;裂相器包括第一低通滤波器和希尔伯特变换器,带通滤波器输出的信号分成两路信号分别给第一低通滤波器和希尔伯特变换器进行低通滤波和希尔伯特变换,输出两路正交的正弦信号;第一微分器模块与第二微分器模块将低通滤波器、希尔伯特变换输出的两路信号分别进行一次微分和二次微分;第一延迟模块、第二延迟模块、第三延迟模块、第四延迟模块将低通滤波后的信号、希尔伯特变换后的信号、两个正交的正弦信号、两路一次微分后信号延迟,消除两个正交的正弦信号、两路一次微分后信号和两路二次微分后信号之间的延迟和超前的时序关系;交叉相乘法或绝对值法推导频率算法模块,将消除延迟和超前的时序关系后的信号通过平方或绝对值推导算法,计算得出频率的平方后输出;第二低通滤波器,对交叉相乘法或绝对值法推导频率算法模块输出的信号进行低通滤波后输出和反馈给带通滤波器。
与现有技术相比,本发明的显著优点为:(1)系统为开环,因此稳定性很好;(2)全系统都在数字信号处理器内实现,没有模拟电路,不易受干扰,抗噪能力强;(3)设置了可调参数的带通滤波器,带宽能随着测量频率的变化而变化,解决了输入信号包含高倍频项使频率测量精度降低的问题,同时也降低到了噪声干扰,提高了测频的精度以及灵敏度;(4)系统协调应用了裂相器,微分器以及交叉相乘法或绝对法值推导频率算法,消除由于希尔伯特变换器与低通滤波器造成的信号随频率变化而产的生幅值波动,同时也避免了极小数除极小数造成的大误差。在很大程度上提高了系统的精度与灵敏度。(5)系统结构简单易于实现。
下面结合附图及具体实施例对本发明作进一步详细说明。
附图说明
图1为本发明的系统框图;
图2为本发明的方法流程图;
图3是本发明的交叉相乘法或绝对值法推导频率算法模块的交叉相乘法结构示意图;
图4是本发明的交叉相乘法或绝对值法推导频率算法模块绝对值相加法结构示意图;
图5是本发明的交叉相乘法或绝对值法推导频率算法模块的绝对值比较法结构示意图;
图6是本发明的FIR滤波器的波纹现象。
具体实施方式
结合图2,一种基于数字信号处理器平台的微分测频方法,包括以下步骤:
步骤1,提供上述微分测频系统,包括带通滤波器100、裂相器200、第一微分器模块300A、第二微分器模块300B、第一延迟模块410A、第二延迟模块410B、第三延迟模块420A、第四延迟模块420B、交叉相乘法或绝对值法推导频率算法模块400、低通滤波器500;
步骤2,MEMS惯性器件的谐振信号输入到参数可调带通滤波器100的输入端口,参数可调带通滤波器100根据控制端接收到的来自输出端的反馈,首先利用算法对反馈回来的值进行判断,看系统是否处于以下三种情况:(1)系统处于刚开机的初始状态,输出端没有数据输出;(2)惯性器件损坏输出异常;(3)惯性器件的工作环境超出设计指标。若为以上任意一种情况,对信号进行在设计的测量范围以内的全通滤波。否则根据控制端接收到的信号调整FIR滤波器参数,降低带通滤波器的通带宽度,使通带宽度上限小于反馈回来的频率的1.5倍且大于反馈回来的频率的1.3倍,下限小于反馈回来的信号的0.8倍且大于反馈回来信号的0.6倍,若超过默认的全通滤波器的范围,上限或下限以全通滤波器的范围为标准,从而抑制频带外的噪声对振荡信号相位噪声的影响。用设计好的滤波器对信号进行滤波后,将信号Asin(ωt)输出。
步骤3,信号Asin(ωt)输入进裂相器200,裂相器200由希尔伯特变换器210与低通滤波器220构成,信号分成两路:一路进行希尔伯特变换,输出C=-b×Acos(ωt);一路进行低通滤波,输出S=a×Asin(ωt)。S与C信号中的参数a与b是关于ω的函数,其值会随着ω的变化而在数值1周围产生小幅的波动,最大幅度小于10%,这是由于FIR滤波器的波纹现象造成的,如图6所示,如果不加处理它会影响测频精度,后面的步骤中会通过算法将其影响消除。
步骤4,将信号S与C分别输入到第一微分器模块300A与第二微分器模块300B,经第一次微分后输出信号经两次微分后输出信号微分器是由根据前后数据用数学算法推导得出当前的数据点的导数值,输出数据会存在滞后现象。
步骤5,第一延迟模块410A、第二延迟模块410B、第三延迟模块420A、第四延迟模块420B分别延迟步骤3中的两路信号和步骤4中第一次微分后的两路信号。第一延迟模块(410A)延迟信号S,第二延迟模块(410B)延迟信号C,第三延迟模块(420A)延迟信号第四延迟模块(420B)延迟信号
步骤6,在步骤5中,利用第一延迟模块410A、第二延迟模块410B将信号S、C、之间的延迟和超前的时序关系消除后,保证信号的同步关系,利用绝对值比较法公式计算推导得到信号的频率值的平方ω2,其公式为(其中使用了第一全波整流模块420A、第二全波整流模块420B、第三全波整流模块420C、第四全波整流模块420D,数值比较模块430,除法器440):
步骤7,将除法器的输出信号ω2输入到低通滤波器500。低通滤波器500是根据加速度计的输入环境与要求而设计,用于降低数字信号处理器与外部干扰而造成的高频噪声,由于输出信号ω2是与惯性器件的敏感量成正比,低通截止频率可以根据惯性器件的带宽来设置,抑制高频信号的干扰。
作为本发明的一种改进,在步骤5中,利用第一延迟模块410A、第二延迟模块410B将信号S、C、之间的延迟和超前的时序关系消除后,保证信号的同步关系,利用绝对值相加法公式计算推导得到信号的频率值的平方ω2,其公式为(其中使用了第一全波整流模块420A、第二全波整流模块420B、第三全波整流模块420C、第四全波整流模块420D,第一加法器430A、第二加法器430B,除法器440):
如果直接使用公式
就会出现零除以零的情况,造成除法器的出错。同时当分子与分母的数字接近于零时,这时噪声对于输出结果的影响就会很大,导致输出数据的大幅波动,而绝对值算法利用裂相器避免了这些情况的发生,同时消除由于希尔伯特变换器210与低通滤波器220造成的信号随频率变化而产的生幅值波动,提高了测频的精度与稳定性。
公式③和公式②的方法比较,公式③判断模块代替了公式②加法模块,使结构更加简单,功耗更低,但是精度也稍有下降。
作为本发明的另一种改进,在步骤5中,利用第一延迟模块410A、第二延迟模块410B、第三延迟模块420A、第四延迟模块420B将信号S、C、 之间的延迟和超前的时序关系消除后,保证信号的同步关系;在步骤6中利用交叉相乘公式计算推导得到信号的频率值的平方ω2,其公式为((其中使用了第一乘法器430A、第二乘法器430B、第三乘法器430C、第四乘法器430D,第一减法器440A、第二减法器440B,除法器450):
算法消除了由于希尔伯特变换器210与低通滤波器220造成的信号随频率变化而产的生幅值波动,提高了测频的精度与稳定性。
经过测试,利用公式①进行测频的精度与稳定性最高,利用公式②次之,利用公式③相比精度与稳定性较低。
结合图1,一种基于数字信号处理器平台的微分测频系统,包括带通滤波器100、裂相器200、第一低通滤波器210、希尔伯特变换器220、第一微分器模块300A、第二微分器模块300B、第一延迟模块410A、第二延迟模块410B、第三延迟模块420A、第四延迟模块420B、交叉相乘法或绝对值法推导频率算法模块400、第二低通滤波器500;
所述带通滤波器100以反馈回来的输出信号为控制量,选择不同的参数,调节带通滤波器的带宽,滤除输入信号中的高倍频项;
裂相器200将带通滤波器100输出的信号分成两个正交的正弦信号后;
第一低通滤波器210将裂相器200分出的一路信号进行低通滤波;
希尔伯特变换器220将裂相器200分出的另一路信号进行希尔伯特变换;
第一微分器模块300A与第二微分器模块300B将低通滤波器210、希尔伯特变换220输出的两路信号分别进行一次微分和二次微分;
第一延迟模块410A、第二延迟模块410B、第三延迟模块420A、第四延迟模块420B将低通滤波后的信号、希尔伯特变换后的信号、两个正交的正弦信号、两路一次微分后信号延迟,消除两个正交的正弦信号、两路一次微分后信号和两路二次微分后信号之间的延迟和超前的时序关系;
交叉相乘法或绝对值法推导频率算法模块400,将消除延迟和超前的时序关系后的信号通过平方或绝对值推导算法,计算得出频率的平方后输出;第二低通滤波器500,对交叉相乘法或绝对值法推导频率算法模块400输出的信号进行低通滤波后输出和反馈给带通滤波器100。
实施例一
MEMS惯性器件与带通滤波器100输入端连接;
带通滤波器100输出端与裂相器200输入端连接;
裂相器200的两个输出端分别连接低通滤波器210和希尔伯特变换器220的输入端;
低通滤波器210输出端分别与第一微分器模块300A输入端和第一延迟模块410A连接;
第一微分器模块300A包括第一微分器310A和第二微分器320A,第一微分器310A输入端与低通滤波器210输出端连接,第一微分器310A输出端分别连接第二微分器320A和第三延迟模块420A的输入端,第二微分器320A输出端与交叉相乘法或绝对值法推导频率算法模块400连接;
希尔伯特变换器220输出端分别与第二微分器模块300B输入端和第二延迟模块410B输入端连接;
第二微分器模块300B包括第三微分器310B和第四微分器320B,第三微分器310B输出端分别与第四微分器320B和第四延迟模块420B输入端连接,第四微分器320B输出端与交叉相乘法或绝对值法推导频率算法模块400连接;
交叉相乘法或绝对值法推导频率算法模块400包括第一乘法器430A、第二乘法器430B、第三乘法器430C、第四乘法器430D,第一减法器440A、第二减法器440B,除法器450;所述第一乘法器430A输入端分别与第一延迟模块410A和第四延迟模块420B连接,输出端与第一减法器440A输入端连接;第二乘法器430B输入端分别与第二微分器320A和第四延迟模块420B连接,输出端与第二减法器440B输入端连接;第三乘法器430C输入端分别与第二延迟模块410B、第三延迟模块420A连接,输出端与与第一减法器440A输入端连接;第四乘法器430D输入端分别与第四微分器320B、第三延迟模块420A连接,输出端与第二减法器440B输入端连接;第一减法器440A、第二减法器440B的输出端分别与除法器450输入端连接;除法器450输出端与第二低通滤波器500输入端连接。
第二低通滤波器500输出端与带通滤波器100连接。
实施例二
MEMS惯性器件与带通滤波器100输入端连接;
带通滤波器100输出端与裂相器200输入端连接;
裂相器200的两个输出端分别连接低通滤波器210和希尔伯特变换器220的输入端;
低通滤波器210输出端分别与第一微分器模块300A输入端和第一延迟模块410A连接;
第一微分器模块300A包括第一微分器310A和第二微分器320A,第一微分器310A输入端与低通滤波器210输出端连接,第一微分器310A输出端连接第二微分器320A输入端,第二微分器320A输出端与交叉相乘法或绝对值法推导频率算法模块400连接;
希尔伯特变换器220输出端分别与第二微分器模块300B输入端和第二延迟模块410B输入端连接;
第二微分器模块300B包括第三微分器310B和第四微分器320B,第三微分器310B输出端与第四微分器320B输入端连接,第四微分器320B输出端与交叉相乘法或绝对值法推导频率算法模块400连接;
交叉相乘法或绝对值法推导频率算法模块400包括第一全波整流模块420A、第二全波整流模块420B、第三全波整流模块420C、第四全波整流模块420D、第一加法器430A、第二加法器430B、除法器440;第一全波整流模块420A输入端与第一延迟模块410A输出端连接,输出端与第二加法器430B输入端连接;第二全波整流模块420B输入端与第二延迟模块410B输出端连接,输出端与第二加法器430B输入端连接;第三全波整流模块420C输入端与第二微分器320A输出端连接,输出端与第一加法器430A输入端连接;第四全波整流模块420D输入端与第四微分器320B输出端连接,输出端与第一加法器430A输入端连接;第一加法器430A、第二加法器430B输出端分别与除法器440输入端连接;除法器450输出端与第二低通滤波器500输入端连接;
第二低通滤波器500输出端与带通滤波器100连接。
实施例三
MEMS惯性器件与带通滤波器100输入端连接;
带通滤波器100输出端与裂相器200输入端连接;
裂相器200的两个输出端分别连接低通滤波器210和希尔伯特变换器220的输入端;
低通滤波器210输出端分别与第一微分器模块300A输入端和第一延迟模块410A连接;
第一微分器模块300A包括第一微分器310A和第二微分器320A,第一微分器310A输入端与低通滤波器210输出端连接,第一微分器310A输出端连接第二微分器320A输入端,第二微分器320A输出端与交叉相乘法或绝对值法推导频率算法模块400连接;
希尔伯特变换器220输出端分别与第二微分器模块300B输入端和第二延迟模块410B输入端连接;
第二微分器模块300B包括第三微分器310B和第四微分器320B,第三微分器310B输出端与第四微分器320B输入端连接,第四微分器320B输出端与交叉相乘法或绝对值法推导频率算法模块400连接;
交叉相乘法或绝对值法推导频率算法模块400包括第一全波整流模块420A、第二全波整流模块420B、第三全波整流模块420C、第四全波整流模块420D、数值比较模块430、除法器440;第一全波整流模块420A输入端与第一延迟模块410A输出端连接,输出端分别与除法器440、数值比较模块430输入端连接;第二全波整流模块420B输入端与第二延迟模块410B输出端连接,输出端分别与除法器440、数值比较模块430输入端连接;第三全波整流模块420C输入端与第四微分器320B输出端连接,输出端与除法器440输入端连接;第四全波整流模块420D输入端与第二微分器320A输出端连接,输出端与除法器440输入端连接;数值比较模块430输出端与除法器440连接,除法器440输出端与第二低通滤波器500输入端连接;
第二低通滤波器500输出端与带通滤波器100连接。

Claims (10)

1.一种基于数字信号处理器平台的微分测频方法,其特征在于,包括以下步骤:
步骤1,提供微分测频系统,包括带通滤波器(100)、裂相器(200)、第一微分器模块(300A)、第二微分器模块(300B)、第一延迟模块(410A)、第二延迟模块(410B)、第三延迟模块(420A)、第四延迟模块(420B)、交叉相乘法或绝对值法推导频率算法模块(400)、低通滤波器(500);
步骤2,带通滤波器(100)对MEMS惯性器件的谐振信号进行滤波,得到信号Asin(ωt);
步骤3,裂相器(200)将滤波后的信号分成两路信号,一路信号进行希尔伯特变换,一路信号进行低通滤波;
步骤4,第一微分器模块(300A)将步骤2中的一路信号分别进行一次微分和两次微分;第二微分器模块(300B)将步骤2中的另一路信号分别进行一次微分和两次微分;
步骤5,第一延迟模块(410A)、第二延迟模块(410B)、第三延迟模块(420A)、第四延迟模块(420B)分别延迟步骤3中的两路信号和步骤4中第一次微分后的两路信号;
步骤6,交叉相乘法或绝对值法推导频率算法模块(400)利用交叉相乘法或绝对值法得到信号的频率值的平方;
步骤7,低通滤波器(500)对频率值的平方进行滤波后输出,并将输出信号反馈于带通滤波器(100)。
2.根据权利要求1所述的基于数字信号处理器平台的微分测频方法,其特征在于,步骤2具体过程为:
带通滤波器(100)根据控制端接收到的来自输出端的反馈,判断系统是否处于以下三种情况:系统处于刚开机的初始状态,输出端没有数据输出;惯性器件损坏输出异常;惯性器件的工作环境超出设计指标;
若满足以上任意一种情况,对谐振信号进行在测量范围以内的全通滤波;否则根据控制端接收到的信号调整滤波器参数,降低带通滤波器的通带宽度,抑制频带外的噪声对振荡信号相位噪声的影响。
3.根据权利要求1所述的基于数字信号处理器平台的微分测频方法,其特征在于,步骤3中,
一路Asin(ωt)信号进行希尔伯特变换,输出C=-b×Acos(ωt);
一路Asin(ωt)信号进行低通滤波,输出S=a×Asin(ωt);
S与C信号中的参数a与b是关于ω的函数,其值会随着ω的变化而在数值1周围产生波动。
4.根据权利要求1所述的基于数字信号处理器平台的微分测频方法,其特征在于,步骤4中,
将信号S与C分别输入到第一微分器模块(300A)与第二微分器模块(300B),经第一次微分后输出信号经两次微分后输出信号
5.根据权利要求1所述的基于数字信号处理器平台的微分测频方法,其特征在于,步骤5的具体过程为:
第一延迟模块(410A)延迟信号S;
第二延迟模块(410B)延迟信号C;
第三延迟模块(420A)延迟信号
第四延迟模块(420B)延迟信号
用第一延迟模块(410A)、第二延迟模块(410B)、第三延迟模块(420A)、第四延迟模块(420B)将信号S、C、之间的延迟和超前的时序关系消除。
6.根据权利要求5所述的基于数字信号处理器平台的微分测频方法,其特征在于,步骤6中利用公式①得到信号Asin(ωt)的频率值的平方ω2
或步骤6中利用绝对值相加法得到信号Asin(ωt)的频率值的平方ω2,见公式②
或步骤6中利用绝对值比较法得到信号Asin(ωt)的频率值的平方ω2,见公式③
7.一种基于数字信号处理器平台的微分测频系统,其特征在于,包括带通滤波器(100)、裂相器(200)、第一低通滤波器(210)、希尔伯特变换器(220)、第一微分器模块(300A)、第二微分器模块(300B)、第一延迟模块(410A)、第二延迟模块(410B)、第三延迟模块(420A)、第四延迟模块(420B)、交叉相乘法或绝对值法推导频率算法模块(400)、第二低通滤波器(500);
所述带通滤波器(100)以反馈回来的输出信号为控制量,选择不同的参数,调节带通滤波器的带宽,滤除输入信号中的高倍频项;
裂相器(200)包括第一低通滤波器(210)和希尔伯特变换器(220),带通滤波器(100)输出的信号分成两路信号分别给第一低通滤波器(210)和希尔伯特变换器(220)进行低通滤波和希尔伯特变换,输出两路正交的正弦信号;
第一微分器模块(300A)与第二微分器模块(300B)将低通滤波器(210)、希尔伯特变换(220)输出的两路信号分别进行一次微分和二次微分;
第一延迟模块(410A)、第二延迟模块(410B)、第三延迟模块(420A)、第四延迟模块(420B)将低通滤波后的信号、希尔伯特变换后的信号、两个正交的正弦信号、两路一次微分后信号延迟,消除两个正交的正弦信号、两路一次微分后信号和两路二次微分后信号之间的延迟和超前的时序关系;
交叉相乘法或绝对值法推导频率算法模块(400),将消除延迟和超前的时序关系后的信号通过平方或绝对值推导算法,计算得出频率的平方后输出;
第二低通滤波器(500),对交叉相乘法或绝对值法推导频率算法模块(400)输出的信号进行低通滤波后输出和反馈给带通滤波器(100)。
8.根据权利要求7所述的基于数字信号处理器平台的微分测频系统,其特征在于,
MEMS惯性器件与带通滤波器(100)输入端连接;
带通滤波器(100)输出端与裂相器(200)输入端连接;
裂相器(200)的两个输出端分别连接低通滤波器(210)和希尔伯特变换器(220)的输入端;
低通滤波器(210)输出端分别与第一微分器模块(300A)输入端和第一延迟模块(410A)连接;
第一微分器模块(300A)包括第一微分器(310A)和第二微分器(320A),第一微分器(310A)输入端与低通滤波器(210)输出端连接,第一微分器(310A)输出端分别连接第二微分器(320A)和第三延迟模块(420A)的输入端,第二微分器(320A)输出端与交叉相乘法或绝对值法推导频率算法模块(400)连接;
希尔伯特变换器(220)输出端分别与第二微分器模块(300B)输入端和第二延迟模块(410B)输入端连接;
第二微分器模块(300B)包括第三微分器(310B)和第四微分器(320B),第三微分器(310B)输出端分别与第四微分器(320B)和第四延迟模块(420B)输入端连接,第四微分器(320B)输出端与交叉相乘法或绝对值法推导频率算法模块(400)连接;
交叉相乘法或绝对值法推导频率算法模块(400)包括第一乘法器(430A)、第二乘法器(430B)、第三乘法器(430C)、第四乘法器(430D),第一减法器(440A)、第二减法器(440B),除法器(450);所述第一乘法器(430A)输入端分别与第一延迟模块(410A)和第四延迟模块(420B)连接,输出端与第一减法器(440A)输入端连接;第二乘法器(430B)输入端分别与第二微分器(320A)和第四延迟模块(420B)连接,输出端与第二减法器(440B)输入端连接;第三乘法器(430C)输入端分别与第二延迟模块(410B)、第三延迟模块(420A)连接,输出端与与第一减法器(440A)输入端连接;第四乘法器(430D)输入端分别与第四微分器(320B)、第三延迟模块(420A)连接,输出端与第二减法器(440B)输入端连接;第一减法器(440A)、第二减法器(440B)的输出端分别与除法器(450)输入端连接;除法器(450)输出端与第二低通滤波器(500)输入端连接;
第二低通滤波器(500)输出端与带通滤波器(100)连接。
9.根据权利要求7所述的基于数字信号处理器平台的微分测频系统,其特征在于,
MEMS惯性器件与带通滤波器(100)输入端连接;
带通滤波器(100)输出端与裂相器(200)输入端连接;
裂相器(200)的两个输出端分别连接低通滤波器(210)和希尔伯特变换器(220)的输入端;
低通滤波器(210)输出端分别与第一微分器模块(300A)输入端和第一延迟模块(410A)连接;
第一微分器模块(300A)包括第一微分器(310A)和第二微分器(320A),第一微分器(310A)输入端与低通滤波器(210)输出端连接,第一微分器(310A)输出端连接第二微分器(320A)输入端,第二微分器(320A)输出端与交叉相乘法或绝对值法推导频率算法模块(400)连接;
希尔伯特变换器(220)输出端分别与第二微分器模块(300B)输入端和第二延迟模块(410B)输入端连接;
第二微分器模块(300B)包括第三微分器(310B)和第四微分器(320B),第三微分器(310B)输出端与第四微分器(320B)输入端连接,第四微分器(320B)输出端与交叉相乘法或绝对值法推导频率算法模块(400)连接;
交叉相乘法或绝对值法推导频率算法模块(400)包括第一全波整流模块(420A)、第二全波整流模块(420B)、第三全波整流模块(420C)、第四全波整流模块(420D)、第一加法器(430A)、第二加法器(430B)、除法器(440);第一全波整流模块(420A)输入端与第一延迟模块(410A)输出端连接,输出端与第二加法器(430B)输入端连接;第二全波整流模块(420B)输入端与第二延迟模块(410B)输出端连接,输出端与第二加法器(430B)输入端连接;第三全波整流模块(420C)输入端与第二微分器(320A)输出端连接,输出端与第一加法器(430A)输入端连接;第四全波整流模块(420D)输入端与第四微分器(320B)输出端连接,输出端与第一加法器(430A)输入端连接;第一加法器(430A)、第二加法器(430B)输出端分别与除法器(440)输入端连接;除法器(450)输出端与第二低通滤波器(500)输入端连接;
第二低通滤波器(500)输出端与带通滤波器(100)连接。
10.根据权利要求7所述的基于数字信号处理器平台的微分测频系统,其特征在于,
MEMS惯性器件与带通滤波器(100)输入端连接;
带通滤波器(100)输出端与裂相器(200)输入端连接;
裂相器(200)的两个输出端分别连接低通滤波器(210)和希尔伯特变换器(220)的输入端;
低通滤波器(210)输出端分别与第一微分器模块(300A)输入端和第一延迟模块(410A)连接;
第一微分器模块(300A)包括第一微分器(310A)和第二微分器(320A),第一微分器(310A)输入端与低通滤波器(210)输出端连接,第一微分器(310A)输出端连接第二微分器(320A)输入端,第二微分器(320A)输出端与交叉相乘法或绝对值法推导频率算法模块(400)连接;
希尔伯特变换器(220)输出端分别与第二微分器模块(300B)输入端和第二延迟模块(410B)输入端连接;
第二微分器模块(300B)包括第三微分器(310B)和第四微分器(320B),第三微分器(310B)输出端与第四微分器(320B)输入端连接,第四微分器(320B)输出端与交叉相乘法或绝对值法推导频率算法模块(400)连接;
交叉相乘法或绝对值法推导频率算法模块(400)包括第一全波整流模块(420A)、第二全波整流模块(420B)、第三全波整流模块(420C)、第四全波整流模块(420D)、数值比较模块(430)、除法器(440);第一全波整流模块(420A)输入端与第一延迟模块(410A)输出端连接,输出端分别与除法器(440)、数值比较模块(430)输入端连接;第二全波整流模块(420B)输入端与第二延迟模块(410B)输出端连接,输出端分别与除法器(440)、数值比较模块(430)输入端连接;第三全波整流模块(420C)输入端与第四微分器(320B)输出端连接,输出端与除法器(440)输入端连接;第四全波整流模块(420D)输入端与第二微分器(320A)输出端连接,输出端与除法器(440)输入端连接;数值比较模块(430)输出端与除法器(440)连接,除法器(440)输出端与第二低通滤波器(500)输入端连接;
第二低通滤波器(500)输出端与带通滤波器(100)连接。
CN201410676742.6A 2014-11-21 2014-11-21 一种基于数字信号处理器平台的微分测频方法及系统 Active CN104391174B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410676742.6A CN104391174B (zh) 2014-11-21 2014-11-21 一种基于数字信号处理器平台的微分测频方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410676742.6A CN104391174B (zh) 2014-11-21 2014-11-21 一种基于数字信号处理器平台的微分测频方法及系统

Publications (2)

Publication Number Publication Date
CN104391174A CN104391174A (zh) 2015-03-04
CN104391174B true CN104391174B (zh) 2017-09-01

Family

ID=52609099

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410676742.6A Active CN104391174B (zh) 2014-11-21 2014-11-21 一种基于数字信号处理器平台的微分测频方法及系统

Country Status (1)

Country Link
CN (1) CN104391174B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104777355B (zh) * 2015-04-21 2018-11-13 上海昕石电子科技有限公司 一种快速检测频率的方法和装置
CN107748308B (zh) * 2017-11-16 2019-09-24 哈尔滨理工大学 一种基于相位判断的电网输电质量评价系统
CN107846017B (zh) * 2017-11-16 2019-08-27 哈尔滨理工大学 基于相位补偿幅值相除的电网输电质量评价方法
CN109917154B (zh) * 2019-03-13 2021-04-06 东南大学 一种硅微谐振式加速度计带宽测试系统及方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204330882U (zh) * 2014-11-21 2015-05-13 南京理工大学 一种基于数字信号处理器平台的微分测频系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001111386A (ja) * 1999-10-04 2001-04-20 Nippon Columbia Co Ltd デジタル信号処理装置
CN1821802A (zh) * 2006-03-16 2006-08-23 江苏金智科技股份有限公司 电力电源的高精度测频测幅模块及方法
CN201540331U (zh) * 2009-10-20 2010-08-04 西安瑞日电子发展有限公司 多通道高精度同步测频装置
CN103487648B (zh) * 2013-09-25 2015-12-02 南京理工大学 sigma-delta PLL频率测量电路及方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204330882U (zh) * 2014-11-21 2015-05-13 南京理工大学 一种基于数字信号处理器平台的微分测频系统

Also Published As

Publication number Publication date
CN104391174A (zh) 2015-03-04

Similar Documents

Publication Publication Date Title
CN104391174B (zh) 一种基于数字信号处理器平台的微分测频方法及系统
CN107786201B (zh) 一种基于锁频环的二阶广义积分器结构及锁相环同步方法
CN206259920U (zh) 一种基于双二阶广义积分器锁相环
CN203014748U (zh) 微机械陀螺仪闭环驱动自动增益控制电路
CN106230446B (zh) 一种抗干扰电路、移动终端及抑制谐波干扰的方法
CN101187682B (zh) 一种在线自适应谐波电流检测系统
Karimi-Ghartemani et al. Comparison of two methods for addressing DC component in phase-locked loop (PLL) systems
CN103684351B (zh) 具有pid指数因子的自适应惯性滤波方法
CN105048995B (zh) 巴特沃斯数字滤波器及利用其实现频率自适应的滤波方法
CN103472302A (zh) 用单相光伏并网逆变器检测电网电压相位的方法
CN204330882U (zh) 一种基于数字信号处理器平台的微分测频系统
CN109358228A (zh) 基于双增强型锁相环的电网电压正负序分量实时估计方法
CN110601678A (zh) 一种iir滤波器实现零相位的方法及其装置
CN106303879A (zh) 一种基于时域分析的检测装置及检测方法
CN112217509B (zh) 基于简化梳状滤波器的主动降噪锁频环装置及方法
CN111624408B (zh) 基于共生多元泛函计算的电力线路对地电容实时测量方法
CN111308892B (zh) 一种基于模型观测的对传感器噪声抑制比可调的闭环控制方法
CN103312324B (zh) 短波段信号的生成方法及系统
Tan et al. Model linearization analysis for three-phase unbalanced phase-locked loop techniques
CN109951186A (zh) 一种基于数字锁相环的扫频信号源
CN108933978A (zh) 一种降低噪声的方法、装置、终端和计算机存储介质
CN107831661A (zh) 一种高精度工业机器人机械振动信号跟踪方法
CN203368443U (zh) 短波段信号的生成系统
CN103346874B (zh) 一种基于dds的数字通信时钟同步系统
CN106841868B (zh) 一种微波光纤延迟线电特性的测量装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant