专利名称：脉动信号的功率最大频率点周期的估计方法
技术领域：
本发明涉及一种脉动信号的功率最大频率点周期的估计方法，具体涉及ー种数据处理方法来估计脉动信号的功率最大频率点对应周期的方法。
背景技术：
脉动信号是传感器以及检测仪表中应用较多的信号形式。对于稳定脉动信号的周期測量与估计已有相当成熟的技术与方法。但当传感器输出的脉动信号周期未知且混有各类干扰脉冲时，往往对脉动信号的正真周期值进行測量与估计都比较困难。如涡街流量计中15mnT600_ 口径传感器输出的涡街脉动信号周期就有约0. 27mS 2S范围变化(对应频率
0.2Hf3.7kHz)。同时，流量计在管道振动等干扰又会在信号上混入许多干扰脉沖。这是在干扰脉冲环境下对脉动信号周期进行测量和估计的典型问题。目前，已有的脉动信号周期的估计方法主要有两类一是采用在线性滤波器下直接对周期捕捉值进行数字滤波方法，如一般涡街流量计的涡街脉动信号周期的估计方法， 主要有美国专利5372046、5429001和5576497等；ニ是先对脉动信号进行A/D转换，再采用数字化频谱FFT方法来分析功率最大频率点的周期值，如一些涡街流量计的新的信号处理方法，典型的有中国专利99101209. 7，或采用与第一类方法相结合的方法，如美国专利 US6480812B1、美国专利US6212975B1、欧洲专利645258、中国专利200610029218. 5和中国专利201110069492. 6等，这类方法有需要解决数字化频谱分析的分辨率和实时性难以兼顾的技术问题，或有用了更多的硬件与软件而不利于实现实际的低功耗仪表。
发明内容
本发明的目的是针对已有技术存在的缺点，提供一种脉动信号的功率最大频率点周期的估计方法，其特点是直接根据脉动信号周期的捕捉值数据，采用数据中数值相差最小的数据组数据的平均值A作为脉动信号的功率最大频率点周期的估计值。这样既降低了对脉动信号的功率谱估计的计算量，又提高了对脉动信号周期值估计的精度和正确性。为达到上述目的，本发明的构思是由ー个周期捕捉器对信号Xl的脉动周期进行捕捉得到捕捉值，数据处理器从n个捕捉值数据中取出数值相差最小的m个数据，并把这m 个数据的平均值A作为信号Xl功率最大频率点所对应周期的估计值。根据上述发明的构思，本发明采用以下技术方案
一种脉动信号功率最大频率点周期的估计方法，其特征在于采用的脉动信号信号測量与处理系统脉动信号Xl连接到ー个周期捕捉単元后，连接到一个数据处理单元；周期捕捉单元将信号Xl脉动周期的捕捉值D送给数据处理单元，数据处理単元把连续送入的n个捕捉值D形成一个数据序列C= {c (I)，c (2),…，c (n)}，并在数据序列C中捜索出m个数值相差最小的数据组G，n;再求出数据组G中m个数据的平均值A作为信号Xl功率
最大频率点周期P的估计值；这里，n越大得到的周期估计值实时性越差，m越大时周期估计值的精度越差。
上述脉动信号功率最大频率点周期的估计方法，其特征在于所述的周期捕捉単元可由ー个信号整形単元和一个跳变计时单元组成；信号整形単元将所述的信号Xl整形成脉冲信号X2送到跳变计时単元，跳变计时単元对脉冲信号X2的电平跳变周期进行计时，将毎次的计时值作为信号X2脉冲周期的捕捉值D送入数据处理器単元。上述脉动信号功率最大频率点周期的估计方法，其特征在于所述的数据处理単元把周期捕捉单元每一次新捕捉值D都作为数据队列C中新的c (n)数据，再把原数据队列C 中的c(i)数据变成c(i-l)数据，i=2，3，…n，即随周期捕捉单元每一次送入新的捕捉值D 后使数据处理单元产生新数据队列C= {c (I)，c (2)，…，c (n)}，数据处理器単元根据每一次新数据队列C捜索出ー个新数据组G，再计算出新数据组G中m个数据的平均值A作为信号 Xl功率最大频率点周期P的估计值；同时可用平均值A为周期来产生输出脉冲W。本发明与现有技术相比，具有如下显而易见的突出实质性特点和显著优点
不用一般的A/D转换方法来得到信号的数字化数据，而直接用信号脉动周期的捕捉值
作为数据，采用在捕捉数据中取出数值相差最小的数据组G，并把数据组G的平均值作为信号功率最大频率点对应周期的估计值。本发明不但降低了对脉动信号功率谱估计的计算量，同时提高了对脉动信号周期值估计的实时性和精度。
图I是本发明的脉动信号信号測量与处理系统结构框图。图2是图I中数据处理単元2在得到捕捉值D后对功率最大频率点周期的估计过程不意图。图3是本发明的ー个实施例结构图。
具体实施例方式本发明的优选实施例结合
如下如
实施例一
參见图I和图2，本脉动信号功率最大频率点周期的估计方法，对于ー个具有功率最大频率点对应为周期P的脉动信号XI，所述的信号Xl连接到ー个周期捕捉単元1，周期捕捉単元I将信号Xl脉动周期的捕捉值D送给数据处理单元2，数据处理単元2把连续送入的 n个捕捉值D形成一个数据序列C= {c (I)，c (2),…，c (n)}，并在数据序列C中捜索出m个数值相差最小的数据组G，n > m > 2 ;再求出数据组G中m个数据的平均值A作为信号Xl 功率最大频率点周期P的估计值；这里，n越大得到的估计值实时性越差，m越大时得出的估计值精度越差。通常对于不同脉动信号处理的实时性和精度要求可以选择不同的n和m数值。在一般脉动频率在几十到几千赫兹信号的处理时可以选择64 > n > 8，m的值可以选择为n 值的一半或更小。实施例ニ
本实施例与实施例一基本相同，特别之处如下
上述脉动信号的功率最大频率点周期的估计方法，所述的周期捕捉単元I可由信号整形単元I. I和跳变计时単元I. 2组成；信号整形単元I. I将所述的信号Xl整形成脉冲信号X2送到跳变计时単元I. 2，跳变计时単元I. 2对脉冲信号X2的电平跳变周期进行计时，将毎次的计时值作为信号X2脉冲周期的捕捉值D送入数据处理器単元2。上述脉动信号的功率最大频率点周期的估计方法，所述的数据处理単元2把周期捕捉単元I每一次新捕捉值D都作为数据队列C中新的c (n)数据，再把原数据队列C中的 c(i)数据变成c(i-l)数据，i=2，3，…n，即随周期捕捉単元I每一次送入的新捕捉值D后使数据处理単元2产生新数据队列C= {c (I)，c (2),…，c (n)}，数据处理器単元2根据每ー 次新数据队列C捜索出ー个新数据组G，再计算出新数据组G中m个数据的平均值A作为信号Xl功率最大频率点周期P的估计值；同时可用平均值A为周期来产生输出脉冲W。这里，数据队列C就是随着ー个新捕捉值D进行滚动，同时随着每ー个新捕捉值D 产生ー个新数据组G和对应的平均值A。即随着ー个新捕捉值D滚动得到ー个信号Xl功率最大频率点周期P的估计值。实施例三
參见图3，本实施例与实施例ニ基本相同，特别之处是信号整形単元I. I采用了通常的迟滞比较器电路。其中电阻R2和Rl是通过对X2和Vref的分压形成了对信号Xl的高低迟滞比较电平。这里，计时单元I. 2和数据处理单元2是采用一个单片机MSP430F2350芯片来承担。计时单元I. 2是单片机MSP430F2350芯片内的TBO跳变捕捉与计时功能，信号X2在芯片的26脚TBO输入，信号X2的每个脉冲周期值D通过对TBO脚上的跳变捕捉与计时来得到；数据处理単元2是通过单片机MSP430F2350的程序按图2的信号处理过程来完成，信号 X2的每个脉冲周期值D可直接从芯片内TBO的捕捉数据得到，通过如图2的信号处理程序得出信号Xl功率最大频率点周期P的估计值A，并用A为周期来产生输出脉冲W。
权利要求
1.一种脉动信号功率最大频率点周期的估计方法，其特征在于采用的脉动信号信号测量与处理系统脉动信号Xi连接到一个周期捕捉单元(I)后，连接到一个数据处理单元(2);周期捕捉单元(I)将信号Xl脉动周期的捕捉值D送给数据处理单元(2);所述的数据处理单元(2)把连续送入的η个捕捉值D形成一个数据序列C={c(l)，c(2)，"％c(n)}，并在数据序列C中搜索出m个数值相差最小的数据组G，n ^ m ^ 2 ;再求出数据组G中m个数据的平均值A作为信号Xl功率最大频率点周期P的估计值；这里，η越大得到的估计值实时性越差，m越大时得出的估计值精度越差。
2.根据权利要求I所述的脉动信号功率最大频率点周期的估计方法，其特征在于所述的周期捕捉单元(I)可由一个信号整形单元(I. D和一个跳变计时单元(I. 2)组成；信号整形单元(I. I)将所述的信号Xl整形成脉冲信号X2送到跳变计时单元(I. 2)，跳变计时单元(I. 2)对脉冲信号X2的电平跳变周期进行计时，将每次的计时值作为信号X2脉冲周期的捕捉值D送入数据处理器单元(2)。
3.根据权利要求I所述的脉动信号功率最大频率点周期的估计方法，其特征在于所述的数据处理单元(2)把周期捕捉单元(I)每一次新捕捉值D都作为数据队列C中新的c (η) 数据，再把原数据队列C中的c(i)数据变成c(i-l)数据，i=2，3，…n，即随周期捕捉单元 (I)每一次送入的新捕捉值D后使数据处理单元(2)产生新数据队列C={c(l)，c(2)，… ，c (η)}，数据处理器单元(2)根据每一次新数据队列C搜索出一个新数据组G，再计算出新数据组G中m个数据的平均值A作为信号Xl功率最大频率点周期P的估计值；同时可用平均值A为周期来产生输出脉冲W。
全文摘要
本发明涉及一种脉动信号功率最大频率点周期的估计方法，采用对信号脉动周期的捕捉值数据，直接在捕捉值数据中计算出数值相差最小的数据组的平均值，并把此平均值作为脉动信号功率最大频率点周期的估计值。这样既降低了对脉动信号功率谱估计的计算量，又提高了对脉动信号周期值估计的实时性和精度。
文档编号H03K3/64GK102594306SQ20121006622
公开日2012年7月18日 申请日期2012年3月14日 优先权日2012年3月14日
发明者李斌, 沈天飞, 钟立源, 陈洁, 陈琳 申请人:上海大学