一种超低频振动参数量值的快速测量方法及装置.pdf

上传人:li****8 文档编号:4289962 上传时间:2018-09-13 格式:PDF 页数:9 大小:423.96KB
返回 下载 相关 举报
摘要
申请专利号:

CN201210149961.X

申请日:

2012.05.14

公开号:

CN102661782A

公开日:

2012.09.12

当前法律状态:

授权

有效性:

有权

法律详情:

授权|||实质审查的生效IPC(主分类):G01H 11/06申请日:20120514|||公开

IPC分类号:

G01H11/06

主分类号:

G01H11/06

申请人:

中国航空工业集团公司北京长城计量测试技术研究所

发明人:

梁志国; 武腾飞; 张大鹏; 李新良; 张大治

地址:

100095 北京市海淀区温泉镇环山村

优先权:

专利代理机构:

代理人:

PDF下载: PDF下载
内容摘要

本发明涉及一种超低频振动参数量值的快速测量方法及装置,属于力学测量技术领域。测量方法是利用小于一个波形周期的局部波形获得振动信号的幅度值和频率,从而实现超低频振动的快速测量;测量装置包括振动传感设备、数据采集系统和计算机数据处理设备。本发明的方法及装置,以虚拟仪器方式实现超低频正弦振动信号幅度值和频率值的快速测量,可以仅仅使用不足一个周期的部分波形获得其振动幅度、频率等量值,达到获得低频振动快速测量与实时刷新数据的效果,避免了传统振动测量方法的必需采样多个波形周期后才能进行幅度和频率的测量,本发明可用于制作超低频振动测量系统、测量仪器设备等,并可以用于超低频振动参数的计量校准。

权利要求书

1.一种超低频振动参数量值的快速测量方法,其特征在于,包括以下步骤:1)利用振动传感设备将表征超低频振动的位移、速度或者加速度物理信号转变成电信号波形;2)利用数据采集系统,将步骤1)所述电信号波形转换成表征振动信号波形的波形数据序列yi,其中i表示数据序列中采样点序号且i=1,2,…,n,n为数据序列的长度,再将数据序列yi送给计算机数据处理设备;3)计算机数据处理设备对输入的数字序列进行正弦波拟合,获得拟合正弦振动波形的幅度角频率和初始相位具体为:3.1对正弦振动波形信号进行数学描述,为其中,A为正弦振动波形的幅值,f为正弦振动波形的频率,为正弦振动波形的初始相位;对步骤2)中输入到计算机数据处理设备的波形数据序列yi进行数学描述,为yi=y(ti)=y((i-1)·Δτ),i=1,2,…,n其中Δτ为采样时间间隔且Δτ=1/v;3.2计算机数据处理设备对波形数据序列yi(i=1,2,…,n)进行正弦波形四参数拟合,获得拟合信号:其中,为正弦振动波形的幅度拟合值,为正弦振动波形的角频率拟合值,为正弦振动波形的初始相位拟合值,为正弦振动波形的直流分量值拟合值;正弦振动波形频率拟合值为: f ^ = ω ^ · v 2 π ]]>4)根据步骤3)的拟合处理结果,以虚拟仪器方式显示超低频振动参数。2.根据权利要求1所述的一种超低频振动参数量值的快速测量方法,其特征在于,所述步骤3.2中进行正弦波形四参数拟合的过程为:(1)设定初始参数,具体包括:待测的正弦波频率预期值为f0,角频率预期值ω0=2πf0/v,波形采集序列所含信号不足一个周期且个数为p,波形采集的数据序列占用时间长度为τ,则f0≤1/τ,选取另一个足够小的正数因子q使得被估计的正弦频率f0>q/τ,因而f0∈[q/τ,2/τ];(2)设定拟合迭代停止条件中收敛参数为he,he是一个接近于0的足够小的正数;(3)从已知采样时刻t1,t2,...,tn得到的正弦波采集样本为y1,y2,...,yn,用计点法获得信号波形占用时间长度为τ=(n-1)/v,确定目标频率f0的存在区间[q/τ,2/τ];(4)确定迭代左边界频率fL==q/τ和迭代左边界角频率ωL=2πfL/v以及迭代右边界频率fR=2/τ和迭代右边界角频率ωR=2πfR/v,令中值角频率ωM=(ωR+ωL)/2;(5)在迭代左边界角频率ωL、迭代右边界角频率ωR和中值角频率ωM上分别利用频率已知的三参数正弦波形拟合公式计算各自的拟合残差ρ(ωL)、ρ(ωR)和ρ(ωM);(6)若ρ(ωL)<η·ρ(ωM),则令迭代右边界角频率ωR=ωM,迭代左边界角频率ωL不变,重复执行步骤(5)~步骤(6),其中η为判据因子;(7)若ρ(ωL)≥η·ρ(ωM),则必有ωR<2ω0,确定左边界角频率为ωL,右边界角频率ωR,按照优选法原则,选取两个中值角频率分别为ωM=ωL+0.618×(ωR-ωL)和ωT=ωR-0.618×(ωR-ωL);(8)在左边界角频率ωL上执行频率已知的三参数正弦曲线拟合获得拟合幅度AL、拟合相位拟合直流分量DL、拟合残差ρL,在右边界角频率ωR上执行频率已知的三参数正弦曲线拟合获得拟合幅度AR、拟合相位拟合直流分量DR、拟合残差ρR,在中值角频率ωM上执行频率已知的三参数正弦曲线拟合获得拟合幅度AM、拟合相位拟合直流分量DM、拟合残差ρM,在中值角频率ωT上执行频率已知的三参数正弦曲线拟合获得拟合幅度AT、拟合相位拟合直流分量DT、拟合残差ρT;(9)若ρM<ρT,则最小拟合残差ρ=ρM,有ω0∈[ωT,ωR],参量更新为ωL=ωT,ωT=ωM,ωM=ωL+0.618×(ωR-ωL);若ρM≥ρT,则最小拟合残差ρ=ρT,有ω0∈[ωL,ωM],参量更新为ωR=ωM,ωM=ωT,ωT=ωR-0.618×(ωR-ωL);(10)判定若|(ρM-ρT)/ρT|≥he,则重复步骤(5)~步骤(10);若|(ρM-ρT)/ρT|<he,则停止迭代,此时又分为两种情况:若最小拟合残差ρ=ρT,获得四参数拟合正弦曲线参数为ρ,拟合过程结束;若最小拟合残差ρ=ρM,获得四参数拟合正弦曲线参数为ρ,拟合过程结束。3.根据权利要求1所述的一种超低频振动参数量值的快速测量方法,其特征在于,所述步骤2)中数据序列的长度n以及数据采样速率v根据被测振动信号频率预估值f设定,设定原则为保证每个波形周期内要多于500个采样点并且所采样的数据少于振动信号的2个波形周期,采样的数据的取值范围取为0.1~2个波形周期数。4.根据权利要求2所述的一种超低频振动参数量值的快速测量方法,其特征在于,所述步骤(2)中he的优选取值范围为1×10-40~1×10-20。5.根据权利要求2所述的一种超低频振动参数量值的快速测量方法,其特征在于,所述步骤(6)中η的优选取值范围为:1≤η≤1.5。6.一种超低频振动参数量值的快速测量装置,其特征在于,包括振动传感设备、数据采集系统和计算机数据处理设备,其中待测的表征超低频振动的位移、速度或者加速度的物理信号进入振动传感设备后转换成电信号波形输出给数据采集系统,数据采集系统再将电信号波形转换为波形数字序列后送入计算机数据处理设备,计算机数据处理设备再将输入的数字序列进行处理,得到超低频振动的幅值和频率测量结果并对结果进行输出。7.根据权利要求6所述的一种超低频振动参数量值的快速测量装置,其特征在于,所述振动传感设备采用光电测量方案,即振动传感设备由激光干涉仪、信号采集处理电路和波形输出电路组成。8.根据权利要求6所述的一种超低频振动参数量值的快速测量装置,其特征在于,所述振动传感设备采用机电测量方案,即振动传感设备由振动传感器和波形输出电路组成。

说明书

一种超低频振动参数量值的快速测量方法及装置

技术领域

本发明涉及一种超低频振动参数量值的快速测量方法及装置,属于力学测
量技术领域。

背景技术

振动是自然界里一个极为基本的物理现象,振动测量技术是力学计量测试
领域的基本技术之一,意义重大且影响深远。很多应用场合的测量需要获取速
度、位移、加速度量值,这类传感器大多是通过振动测量的方式进行校准和标
定的,力学领域里影响深远的材料疲劳特性试验、飞机等大型装备的定寿命试
验等均是通过振动及其测量方式实现的。其它如地震监测、石油勘探、核试验
监测等很多场合需要使用振动测量手段。

描述振动的物理量通常使用位移、速度、加速度三种之一进行,三者既有
区别又互相联系,并有确定函数关系,速度是位移的微分变量,而加速度则是
速度的微分变量。振动测量所关心的最终结果主要是振动幅度、振动频率以及
振动相位延迟,其它参数如振动波形失真、振动漂移等则属于辅助参数。

振动测量中,数学表述均使用正弦波表示,其幅度、频率等参数量值的确
定通常需要多个振动波形周期,在绝大多数情况下,都能获得很好的效果,但
是在面对超低频振动测量时,表现出测量响应时间过长的难题。例如,用于地
震监测技术中的超低频振动台,其振动周期为1000s,一个周期波形产生需要相
当长的时间,按照传统的振动测量方法,在这样长的时间里,测量系统无法获
取有关振动的参数信息,仅仅可测量获得一段振动波形而已。这也导致很难快
速获得超低频振动的幅度和频率等信息,振动幅度和频率等参数量值的测量值
的实时刷新很难实现。

发明内容

本发明的目的是为了克服现有技术的局陷,提出一种超低频振动参数量值
的快速测量方法及装置,以虚拟仪器的方式在超低频正弦振动测量过程中实现
振动幅度和频率等参数量值的快速测量,其核心思路是使用少于一个波形周期
的正弦振动波形(例如十分之一个波形周期)获取振动信号的幅度和频率等参
数量值。

本发明是通过以下技术方案实现的。

本发明的一种超低频振动参数量值的快速测量方法,具体包括以下步骤:

1)利用振动传感设备将表征超低频振动的位移、速度或者加速度物理信号
转变成电信号波形;

2)利用数据采集系统,将步骤1)所述电信号波形转换成表征振动信号波
形的波形数据序列yi,其中i表示数据序列中采样点序号且i=1,2,…,n,n为数据
序列的长度,再将数据序列yi送给计算机数据处理设备;

3)计算机数据处理设备对输入的数字序列进行正弦波拟合,获得拟合正弦
振动波形的幅度角频率和初始相位具体为:

3.1对正弦振动波形信号进行数学描述,为


其中,A为正弦振动波形的幅值,f为正弦振动波形的频率,为正弦振动波
形的初始相位;

对步骤2)中输入到计算机数据处理设备的波形数据序列yi进行数学描述,

yi=y(ti)=y((i-1)·Δτ),i=1,2,…,n

其中Δτ为采样时间间隔且Δτ=1/v;

3.2计算机数据处理设备对波形数据序列yi(i=1,2,…,n)进行正弦波形四
参数拟合,获得拟合信号:


其中,为正弦振动波形的幅度拟合值,为正弦振动波形的角频率拟合值,
为正弦振动波形的初始相位拟合值,为正弦振动波形的直流分量值拟合值;

正弦振动波形频率拟合值为:

f ^ = ω ^ · v 2 π ]]>

4)根据步骤3)的拟合处理结果,以虚拟仪器方式显示超低频振动参数。

上述步骤3.2中进行正弦波形四参数拟合的过程为:

(1)设定初始参数,具体包括:待测的正弦波频率预期值为f0,角频率预期
值ω0=2πf0/v,波形采集序列所含信号不足一个周期且个数为p,波形采集的数据
序列占用时间长度为τ,则f0≤1/τ,选取另一个足够小的正数因子q,使得被估计
的正弦频率f0>q/τ,因而f0∈[q/τ,2/τ];

(2)设定拟合迭代停止条件中收敛参数为he,he是一个接近于0的足够小的
正数;

(3)从已知采样时刻t1,t2,...,tn得到的正弦波采集样本为y1,y2,...,yn,用计点
法获得信号波形占用时间长度为τ=(n-1)/v,确定目标频率f0的存在区间[q/τ,2/τ];

(4)确定迭代左边界频率fL=q/τ和迭代左边界角频率ωL=2πfL/v以及迭代右边
界频率fR=2/τ和迭代右边界角频率ωR=2πfR/v,令中值角频率ωM=(ωR+ωL)/2;

(5)在迭代左边界角频率ωL、迭代右边界角频率ωR和中值角频率ωM上分别
利用频率已知的三参数正弦波形拟合公式计算各自的拟合残差ρ(ωL)、ρ(ωR)和
ρ(ωM);

(6)若ρ(ωL)<η·ρ(ωM),则令迭代右边界角频率ωR=ωM,迭代左边界角频率ωL
不变,重复执行步骤(5)~步骤(6),其中η为判据因子;

(7)若ρ(ωL)≥η·ρ(ωM),则必有ωR<2ω0,确定左边界角频率为ωL,右边界角频
率ωR,按照优选法原则,选取两个中值角频率分别为

ωM=ωL+0.618×(ωR-ωL)


ωT=ωR-0.618×(ωR-ωL);

(8)在左边界角频率ωL上执行频率已知的三参数正弦曲线拟合获得拟合幅
度AL、拟合相位拟合直流分量DL、拟合残差ρL,在右边界角频率ωR上执行
频率已知的三参数正弦曲线拟合获得拟合幅度AR、拟合相位拟合直流分量
DR、拟合残差ρR,在中值角频率ωM上执行频率已知的三参数正弦曲线拟合获得
拟合幅度AM、拟合相位拟合直流分量DM、拟合残差ρM,在中值角频率ωT
上执行频率已知的三参数正弦曲线拟合获得拟合幅度AT、拟合相位拟合直
流分量DT、拟合残差ρT;

(9)若ρM<ρT,则最小拟合残差ρ=ρM,有ω0∈[ωT,ωR],参量更新为ωL=ωT,
ωT=ωM,ωM=ωL+0.618×(ωR-ωL);若ρM≥ρT,则最小拟合残差ρ=ρT,有
ω0∈[ωL,ωM],参量更新为ωR=ωM,ωM=ωT,ωT=ωR-0.618×(ωR-ωL);

(10)判定若|(ρM-ρT)/ρT|≥he,则重复步骤(5)~步骤(10);若
|(ρM-ρT)/ρT|<he,则停止迭代,此时又分为两种情况:若最小拟合残差ρ=ρT,
获得四参数拟合正弦曲线参数为ρ,拟合过程结
束;若最小拟合残差ρ=ρM,获得四参数拟合正弦曲线参数为
ρ,拟合过程结束。

一种超低频振动参数量值的快速测量装置,用于实现上述超低频振动参数
量值的快速测量方法,包括振动传感设备、数据采集系统和计算机数据处理设
备,其中待测的表征超低频振动的位移、速度或者加速度物理信号进入振动传
感设备后转换成电信号波形输出给数据采集系统,数据采集系统再将电信号波
形转换为波形数字序列后送入计算机数据处理设备,计算机数据处理设备再将
输入的数字序列进行处理,得到超低频振动的参数量值也即幅值和频率的测量
结果并对结果进行输出。

有益效果

本发明的方法及装置,以虚拟仪器方式实现超低频正弦振动信号幅度值和
频率值的快速测量,可以仅仅使用不足一个周期的部分波形获得其振动幅度、
频率等量值,例如在仅有十分之一个波形周期下即可以实现振动参量的测量,
达到获得低频振动快速测量与实时刷新数据的效果。本发明具有快速测量超低
频振动量值的特点,可用于制作超低频振动测量系统、测量仪器设备等,并可
以用于超低频振动参数的计量校准。

附图说明

图1为本发明采用光电测量方案时快速测量装置的结构示意图;

图2为本发明采用机电测量方案时快速测量装置的结构示意图。

具体实施方式

下面结合附图和实施例对本发明做进一步说明。

实施例

一种超低频振动参数量值的快速测量方法,具体包括以下步骤:

1)利用振动传感设备将表征超低频振动的位移、速度或者加速度物理信号
转变成电信号波形;

2)利用数据采集系统,将步骤1)所述电信号波形转换成表征振动信号波
形的波形数据序列yi,其中i表示数据序列中采样点序号且i=1,2,…,n,n为数据
序列的长度,再将数据序列yi送给计算机数据处理设备;

数据序列的长度n以及数据采样速率v根据被测振动信号频率预估值f设定,
设定原则为保证每个波形周期内要多于500个采样点并且所采样的数据少于振
动信号的2个波形周期,采样的数据的取值范围取为0.1~2个波形周期数;

3)计算机数据处理设备对输入的数字序列进行正弦波拟合,获得拟合正弦
振动波形的幅度角频率和初始相位具体为:

3.1对正弦振动波形信号进行数学描述,为


其中,A为正弦振动波形的幅值,f为正弦振动波形的频率,为正弦振动波
形的初始相位;

对步骤2)中输入到计算机数据处理设备的波形数据序列yi进行数学描述,

yi=y(ti)=y((i-1)·Δτ),i=1,2,…,n

其中Δτ为采样时间间隔且Δτ=1/v;

3.2计算机数据处理设备对波形数据序列yi(i=1,2,…,n)进行正弦波形四
参数拟合,获得拟合信号:


其中,为正弦振动波形的幅度拟合值,为正弦振动波形的角频率拟合值,
为正弦振动波形的初始相位拟合值,为正弦振动波形的直流分量值拟合值;

正弦振动波形频率拟合值为:

f ^ = ω ^ · v 2 π ]]>

4)根据步骤3)的拟合处理结果,以虚拟仪器方式显示超低频振动参数。

上述步骤3.2中进行正弦波形四参数拟合的过程为:

(1)设定初始参数,具体包括:待测的正弦波频率预期值为f0,则ω0=2πf0/v,
波形采集序列所含信号不足一个周期且个数为p,波形采集序列占用时间长度为
τ,则f0≤1/τ,选取另一个足够小的正数因子q(此处选取q=1×10-5),使得被估
计的正弦频率f0>q/τ,因而f0∈[q/τ,2/τ];

(2)设定拟合迭代停止条件中收敛参数为一个接近于0的足够小的正数he,
此处选取he=1×10-20;

(3)从已知采样时刻t1,t2,...,tn得到的正弦波采集样本为y1,y2,...,yn,用计点
法获得信号波形占用时间长度为τ=(n-1)/v,确定目标频率f0的存在区间[q/τ,2/τ];

(4)确定迭代左边界频率fL=q/τ和迭代左边界角频率ωL=2πfL/v以及迭代右边
界频率fR=2/τ和迭代右边界角频率ωR=2πfR/v,令中值角频率ωM=(ωR+ωL)/2;

(5)在迭代左边界角频率ωL、迭代右边界角频率ωR和中值角频率ωM上分别
利用频率已知的三参数正弦波形拟合公式计算各自的拟合残差ρ(ωL)、ρ(ωR)和
ρ(ωM);

(6)若ρ(ωL)<η·ρ(ωM),则令迭代右边界角频率ωR=ωM,迭代左边界角频率ωL
不变,重复执行步骤(5)~步骤(6),其中η为判据因子,η的取值范围为1~1.5;

(7)若ρ(ωL)≥η·ρ(ωM),则必有ωR<2ω0,确定左边界角频率为ωL,右边界角频
率ωR,按照优选法原则,选取两个中值角频率分别为

ωM=ωL+0.618×(ωR-ωL)


ωT=ωR-0.618×(ωR-ωL);

(8)在左边界角频率ωL上执行频率已知的三参数正弦曲线拟合获得拟合幅
度AL、拟合相位拟合直流分量DL、拟合残差ρL,在右边界角频率ωR上执行
频率已知的三参数正弦曲线拟合获得拟合幅度AR、拟合相位拟合直流分量
DR、拟合残差ρR,在中值角频率ωM上执行频率已知的三参数正弦曲线拟合获得
拟合幅度AM、拟合相位拟合直流分量DM、拟合残差ρM,在中值角频率ωT
上执行频率已知的三参数正弦曲线拟合获得拟合幅度AT、拟合相位拟合直
流分量DT、拟合残差ρT;

(9)若ρM<ρT,则最小拟合残差ρ=ρM,有ω0∈[ωT,ωR],参量更新为ωL=ωT,
ωT=ωM,ωM=ωL+0.618×(ωR-ωL);若ρM≥ρT,则最小拟合残差ρ=ρT,有
ω0∈[ωL,ωM],参量更新为ωR=ωM,ωM=ωT,ωT=ωR-0.618×(ωR-ωL);

(10)判定若|(ρM-ρT)/ρT|≥he,则重复步骤(5)~步骤(10);若
|(ρM-ρT)/ρT|<he,则停止迭代,此时又分为两种情况:若最小拟合残差ρ=ρT,
获得四参数拟合正弦曲线参数为ρ,拟合过程结
束;若最小拟合残差ρ=ρM,获得四参数拟合正弦曲线参数为
ρ,拟合过程结束。

一种超低频振动参数量值的快速测量装置,用于实现上述超低频振动参数
量值的快速测量方法,包括振动传感设备、数据采集系统和计算机数据处理设
备,其中待测的表征超低频振动的位移、速度或者加速度物理信号进入振动传
感设备后转换成电信号波形输出给数据采集系统,数据采集系统再将电信号波
形转换为波形数字序列后送入计算机数据处理设备,计算机数据处理设备再将
输入的数字序列进行处理,得到超低频振动的幅值和频率测量结果并对结果进
行输出。

所述振动传感设备可采用如下两种技术方案:一种是光电测量方案,即振
动传感设备由激光干涉仪、信号采集处理电路和波形输出电路组成,此时超低
频振动参数量值的快速测量装置如图1所示;第二种是机电测量方案,即振动传
感设备由振动传感器和波形输出电路组成,此时超低频振动参数量值的快速测
量装置如图2所示。

目前用于地震监测所用的超低频振动台,其输出频率下限值可以达到
1mHz,一个波形周期达到1000s,约17分钟,使用通常的振动测量方法进行正常
测量,17分钟以内不会有任何响应,更不能进行快速测量其相位差值。将本发
明的超低频振动的快速测量方法及装置用于解决该问题,通过波形测量手段,
在仅仅采集十分之一个周期波形条件下,在2分钟之内即可以快速获得测量结
果,在噪声信号幅度比为2%的情况下,可以获得幅度估计误差优于±3%、频率
估计误差优于±2%、相位估计误差优于±0.6°的效果。

以上所述为本发明的较佳实施例而已,本发明不应该局限于该实施例和附
图所公开的内容。凡是不脱离本发明所公开的精神下完成的等效或修改,都落
入本发明保护的范围。

一种超低频振动参数量值的快速测量方法及装置.pdf_第1页
第1页 / 共9页
一种超低频振动参数量值的快速测量方法及装置.pdf_第2页
第2页 / 共9页
一种超低频振动参数量值的快速测量方法及装置.pdf_第3页
第3页 / 共9页
点击查看更多>>
资源描述

《一种超低频振动参数量值的快速测量方法及装置.pdf》由会员分享,可在线阅读,更多相关《一种超低频振动参数量值的快速测量方法及装置.pdf(9页珍藏版)》请在专利查询网上搜索。

1、(10)申请公布号 CN 102661782 A (43)申请公布日 2012.09.12 C N 1 0 2 6 6 1 7 8 2 A *CN102661782A* (21)申请号 201210149961.X (22)申请日 2012.05.14 G01H 11/06(2006.01) (71)申请人中国航空工业集团公司北京长城计 量测试技术研究所 地址 100095 北京市海淀区温泉镇环山村 (72)发明人梁志国 武腾飞 张大鹏 李新良 张大治 (54) 发明名称 一种超低频振动参数量值的快速测量方法及 装置 (57) 摘要 本发明涉及一种超低频振动参数量值的快速 测量方法及装置,属于。

2、力学测量技术领域。测量方 法是利用小于一个波形周期的局部波形获得振动 信号的幅度值和频率,从而实现超低频振动的快 速测量;测量装置包括振动传感设备、数据采集 系统和计算机数据处理设备。本发明的方法及装 置,以虚拟仪器方式实现超低频正弦振动信号幅 度值和频率值的快速测量,可以仅仅使用不足一 个周期的部分波形获得其振动幅度、频率等量值, 达到获得低频振动快速测量与实时刷新数据的效 果,避免了传统振动测量方法的必需采样多个波 形周期后才能进行幅度和频率的测量,本发明可 用于制作超低频振动测量系统、测量仪器设备等, 并可以用于超低频振动参数的计量校准。 (51)Int.Cl. 权利要求书2页 说明书5。

3、页 附图1页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 5 页 附图 1 页 1/2页 2 1.一种超低频振动参数量值的快速测量方法,其特征在于,包括以下步骤: 1)利用振动传感设备将表征超低频振动的位移、速度或者加速度物理信号转变成电信 号波形; 2)利用数据采集系统,将步骤1)所述电信号波形转换成表征振动信号波形的波形数 据序列y i ,其中i表示数据序列中采样点序号且i1,2,n,n为数据序列的长度,再 将数据序列y i 送给计算机数据处理设备; 3)计算机数据处理设备对输入的数字序列进行正弦波拟合,获得拟合正弦振动波形的 幅度角频率和初始相。

4、位具体为: 3.1对正弦振动波形信号进行数学描述,为 其中,A为正弦振动波形的幅值,f为正弦振动波形的频率,为正弦振动波形的初始相 位; 对步骤2)中输入到计算机数据处理设备的波形数据序列yi进行数学描述,为 y i y(t i )y(i-1),i1,2,n 其中为采样时间间隔且1/v; 3.2计算机数据处理设备对波形数据序列y i (i1,2,n)进行正弦波形四参数拟 合,获得拟合信号: 其中,为正弦振动波形的幅度拟合值,为正弦振动波形的角频率拟合值,为正弦振 动波形的初始相位拟合值,为正弦振动波形的直流分量值拟合值; 正弦振动波形频率拟合值为: 4)根据步骤3)的拟合处理结果,以虚拟仪器方。

5、式显示超低频振动参数。 2.根据权利要求1所述的一种超低频振动参数量值的快速测量方法,其特征在于,所 述步骤3.2中进行正弦波形四参数拟合的过程为: (1)设定初始参数,具体包括:待测的正弦波频率预期值为f 0 ,角频率预期值 0 2f 0 /v,波形采集序列所含信号不足一个周期且个数为p,波形采集的数据序列占用时间 长度为,则f 0 1/,选取另一个足够小的正数因子q使得被估计的正弦频率f 0 q/ ,因而f 0 q/,2/; (2)设定拟合迭代停止条件中收敛参数为h e ,h e 是一个接近于0的足够小的正数; (3)从已知采样时刻t 1 ,t 2 ,.,t n 得到的正弦波采集样本为y 。

6、1 ,y 2 ,.,y n ,用计点法 获得信号波形占用时间长度为(n-1)/v,确定目标频率f 0 的存在区间q/,2/; (4)确定迭代左边界频率f L q/和迭代左边界角频率 L 2f L /v以及迭代右 边界频率f R 2/和迭代右边界角频率 R 2f R /v,令中值角频率 M ( R + L )/2; (5)在迭代左边界角频率 L 、迭代右边界角频率 R 和中值角频率 M 上分别利用频 率已知的三参数正弦波形拟合公式计算各自的拟合残差( L )、( R )和( M ); (6)若( L )( M ),则令迭代右边界角频率 R M ,迭代左边界角频率 权 利 要 求 书CN 1026。

7、61782 A 2/2页 3 L 不变,重复执行步骤(5)步骤(6),其中为判据因子; (7)若( L )( M ),则必有 R 2 0 ,确定左边界角频率为 L ,右边界角 频率 R ,按照优选法原则,选取两个中值角频率分别为 M L +0.618( R - L ) 和 T R -0.618( R - L ); (8)在左边界角频率 L 上执行频率已知的三参数正弦曲线拟合获得拟合幅度A L 、拟合 相位拟合直流分量D L 、拟合残差 L ,在右边界角频率 R 上执行频率已知的三参数正弦 曲线拟合获得拟合幅度A R 、拟合相位拟合直流分量D R 、拟合残差 R ,在中值角频率 M 上执行频率已。

8、知的三参数正弦曲线拟合获得拟合幅度A M 、拟合相位拟合直流分量D M 、拟 合残差 M ,在中值角频率 T 上执行频率已知的三参数正弦曲线拟合获得拟合幅度A T 、拟 合相位拟合直流分量D T 、拟合残差 T ; (9)若 M T ,则最小拟合残差 M ,有 0 T , R ,参量更新为 L T , T M , M L +0.618( R - L );若 M T ,则最小拟合残差 T ,有 0 L , M ,参量更新为 R M , M T , T R -0.618( R - L ); (10)判定若|( M - T )/ T |h e ,则重复步骤(5)步骤(10);若|( M - T )/。

9、 T | h e ,则停止迭代,此时又分为两种情况:若最小拟合残差 T ,获得四参数拟合正弦曲 线参数为,拟合过程结束;若最小拟合残差 M ,获得 四参数拟合正弦曲线参数为,拟合过程结束。 3.根据权利要求1所述的一种超低频振动参数量值的快速测量方法,其特征在于,所 述步骤2)中数据序列的长度n以及数据采样速率v根据被测振动信号频率预估值f设定, 设定原则为保证每个波形周期内要多于500个采样点并且所采样的数据少于振动信号的2 个波形周期,采样的数据的取值范围取为0.12个波形周期数。 4.根据权利要求2所述的一种超低频振动参数量值的快速测量方法,其特征在于,所 述步骤(2)中h e 的优选取。

10、值范围为110 -40 110 -20 。 5.根据权利要求2所述的一种超低频振动参数量值的快速测量方法,其特征在于,所 述步骤(6)中的优选取值范围为:11.5。 6.一种超低频振动参数量值的快速测量装置,其特征在于,包括振动传感设备、数据采 集系统和计算机数据处理设备,其中待测的表征超低频振动的位移、速度或者加速度的物 理信号进入振动传感设备后转换成电信号波形输出给数据采集系统,数据采集系统再将电 信号波形转换为波形数字序列后送入计算机数据处理设备,计算机数据处理设备再将输入 的数字序列进行处理,得到超低频振动的幅值和频率测量结果并对结果进行输出。 7.根据权利要求6所述的一种超低频振动参。

11、数量值的快速测量装置,其特征在于,所 述振动传感设备采用光电测量方案,即振动传感设备由激光干涉仪、信号采集处理电路和 波形输出电路组成。 8.根据权利要求6所述的一种超低频振动参数量值的快速测量装置,其特征在于,所 述振动传感设备采用机电测量方案,即振动传感设备由振动传感器和波形输出电路组成。 权 利 要 求 书CN 102661782 A 1/5页 4 一种超低频振动参数量值的快速测量方法及装置 技术领域 0001 本发明涉及一种超低频振动参数量值的快速测量方法及装置,属于力学测量技术 领域。 背景技术 0002 振动是自然界里一个极为基本的物理现象,振动测量技术是力学计量测试领域的 基本技。

12、术之一,意义重大且影响深远。很多应用场合的测量需要获取速度、位移、加速度量 值,这类传感器大多是通过振动测量的方式进行校准和标定的,力学领域里影响深远的材 料疲劳特性试验、飞机等大型装备的定寿命试验等均是通过振动及其测量方式实现的。其 它如地震监测、石油勘探、核试验监测等很多场合需要使用振动测量手段。 0003 描述振动的物理量通常使用位移、速度、加速度三种之一进行,三者既有区别又互 相联系,并有确定函数关系,速度是位移的微分变量,而加速度则是速度的微分变量。振动 测量所关心的最终结果主要是振动幅度、振动频率以及振动相位延迟,其它参数如振动波 形失真、振动漂移等则属于辅助参数。 0004 振动。

13、测量中,数学表述均使用正弦波表示,其幅度、频率等参数量值的确定通常需 要多个振动波形周期,在绝大多数情况下,都能获得很好的效果,但是在面对超低频振动测 量时,表现出测量响应时间过长的难题。例如,用于地震监测技术中的超低频振动台,其振 动周期为1000s,一个周期波形产生需要相当长的时间,按照传统的振动测量方法,在这样 长的时间里,测量系统无法获取有关振动的参数信息,仅仅可测量获得一段振动波形而已。 这也导致很难快速获得超低频振动的幅度和频率等信息,振动幅度和频率等参数量值的测 量值的实时刷新很难实现。 发明内容 0005 本发明的目的是为了克服现有技术的局陷,提出一种超低频振动参数量值的快速 。

14、测量方法及装置,以虚拟仪器的方式在超低频正弦振动测量过程中实现振动幅度和频率等 参数量值的快速测量,其核心思路是使用少于一个波形周期的正弦振动波形(例如十分之 一个波形周期)获取振动信号的幅度和频率等参数量值。 0006 本发明是通过以下技术方案实现的。 0007 本发明的一种超低频振动参数量值的快速测量方法,具体包括以下步骤: 0008 1)利用振动传感设备将表征超低频振动的位移、速度或者加速度物理信号转变成 电信号波形; 0009 2)利用数据采集系统,将步骤1)所述电信号波形转换成表征振动信号波形的波 形数据序列y i ,其中i表示数据序列中采样点序号且i1,2,n,n为数据序列的长度,。

15、 再将数据序列y i 送给计算机数据处理设备; 0010 3)计算机数据处理设备对输入的数字序列进行正弦波拟合,获得拟合正弦振动波 形的幅度角频率和初始相位具体为: 说 明 书CN 102661782 A 2/5页 5 0011 3.1对正弦振动波形信号进行数学描述,为 0012 0013 其中,A为正弦振动波形的幅值,f为正弦振动波形的频率,为正弦振动波形的初 始相位; 0014 对步骤2)中输入到计算机数据处理设备的波形数据序列y i 进行数学描述,为 0015 y i y(t i )y(i-1),i1,2,n 0016 其中为采样时间间隔且1/v; 0017 3.2计算机数据处理设备对波。

16、形数据序列y i (i1,2,n)进行正弦波形四参 数拟合,获得拟合信号: 0018 0019 其中,为正弦振动波形的幅度拟合值,为正弦振动波形的角频率拟合值,为正 弦振动波形的初始相位拟合值,为正弦振动波形的直流分量值拟合值; 0020 正弦振动波形频率拟合值为: 0021 0022 4)根据步骤3)的拟合处理结果,以虚拟仪器方式显示超低频振动参数。 0023 上述步骤3.2中进行正弦波形四参数拟合的过程为: 0024 (1)设定初始参数,具体包括:待测的正弦波频率预期值为f 0 ,角频率预期值 0 2f 0 /v,波形采集序列所含信号不足一个周期且个数为p,波形采集的数据序列占用时间 长度。

17、为,则f 0 1/,选取另一个足够小的正数因子q,使得被估计的正弦频率f 0 q/ ,因而f 0 q/,2/; 0025 (2)设定拟合迭代停止条件中收敛参数为h e ,h e 是一个接近于0的足够小的正数; 0026 (3)从已知采样时刻t 1 ,t 2 ,.,t n 得到的正弦波采集样本为y 1 ,y 2 ,.,y n ,用计 点法获得信号波形占用时间长度为(n-1)/v,确定目标频率f 0 的存在区间q/,2/ ; 0027 (4)确定迭代左边界频率f L q/和迭代左边界角频率 L 2f L /v以及 迭代右边界频率f R 2/和迭代右边界角频率 R 2f R /v,令中值角频率 M 。

18、( R + L )/2; 0028 (5)在迭代左边界角频率 L 、迭代右边界角频率 R 和中值角频率 M 上分别利 用频率已知的三参数正弦波形拟合公式计算各自的拟合残差( L )、( R )和( M ); 0029 (6)若( L )( M ),则令迭代右边界角频率 R M ,迭代左边界角 频率 L 不变,重复执行步骤(5)步骤(6),其中为判据因子; 0030 (7)若( L )( M ),则必有 R 2 0 ,确定左边界角频率为 L ,右边 界角频率 R ,按照优选法原则,选取两个中值角频率分别为 0031 M L +0.618( R - L ) 0032 和 0033 T R -0.6。

19、18( R - L ); 说 明 书CN 102661782 A 3/5页 6 0034 (8)在左边界角频率 L 上执行频率已知的三参数正弦曲线拟合获得拟合幅度A L 、 拟合相位拟合直流分量D L 、拟合残差 L ,在右边界角频率 R 上执行频率已知的三参 数正弦曲线拟合获得拟合幅度A R 、拟合相位拟合直流分量D R 、拟合残差 R ,在中值角频 率 M 上执行频率已知的三参数正弦曲线拟合获得拟合幅度A M 、拟合相位拟合直流分量 D M 、拟合残差 M ,在中值角频率 T 上执行频率已知的三参数正弦曲线拟合获得拟合幅度 A T 、拟合相位拟合直流分量D T 、拟合残差 T ; 0035。

20、 (9)若 M T ,则最小拟合残差 M ,有 0 T , R ,参量更新为 L T , T M , M L +0.618( R - L );若 M T ,则最小拟合残差 T ,有 0 L , M ,参量更新为 R M , M T , T R -0.618( R - L ); 0036 (10)判定若|( M - T )/ T |h e ,则重复步骤(5)步骤(10);若|( M - T )/ T |h e ,则停止迭代,此时又分为两种情况:若最小拟合残差 T ,获得四参数拟合正 弦曲线参数为,拟合过程结束;若最小拟合残差 M , 获得四参数拟合正弦曲线参数为,拟合过程结束。 0037 一种超。

21、低频振动参数量值的快速测量装置,用于实现上述超低频振动参数量值的 快速测量方法,包括振动传感设备、数据采集系统和计算机数据处理设备,其中待测的表征 超低频振动的位移、速度或者加速度物理信号进入振动传感设备后转换成电信号波形输出 给数据采集系统,数据采集系统再将电信号波形转换为波形数字序列后送入计算机数据处 理设备,计算机数据处理设备再将输入的数字序列进行处理,得到超低频振动的参数量值 也即幅值和频率的测量结果并对结果进行输出。 0038 有益效果 0039 本发明的方法及装置,以虚拟仪器方式实现超低频正弦振动信号幅度值和频率值 的快速测量,可以仅仅使用不足一个周期的部分波形获得其振动幅度、频率。

22、等量值,例如在 仅有十分之一个波形周期下即可以实现振动参量的测量,达到获得低频振动快速测量与实 时刷新数据的效果。本发明具有快速测量超低频振动量值的特点,可用于制作超低频振动 测量系统、测量仪器设备等,并可以用于超低频振动参数的计量校准。 附图说明 0040 图1为本发明采用光电测量方案时快速测量装置的结构示意图; 0041 图2为本发明采用机电测量方案时快速测量装置的结构示意图。 具体实施方式 0042 下面结合附图和实施例对本发明做进一步说明。 0043 实施例 0044 一种超低频振动参数量值的快速测量方法,具体包括以下步骤: 0045 1)利用振动传感设备将表征超低频振动的位移、速度或。

23、者加速度物理信号转变成 电信号波形; 0046 2)利用数据采集系统,将步骤1)所述电信号波形转换成表征振动信号波形的波 形数据序列y i ,其中i表示数据序列中采样点序号且i1,2,n,n为数据序列的长度, 再将数据序列y i 送给计算机数据处理设备; 说 明 书CN 102661782 A 4/5页 7 0047 数据序列的长度n以及数据采样速率v根据被测振动信号频率预估值f设定,设 定原则为保证每个波形周期内要多于500个采样点并且所采样的数据少于振动信号的2个 波形周期,采样的数据的取值范围取为0.12个波形周期数; 0048 3)计算机数据处理设备对输入的数字序列进行正弦波拟合,获得。

24、拟合正弦振动波 形的幅度角频率和初始相位具体为: 0049 3.1对正弦振动波形信号进行数学描述,为 0050 0051 其中,A为正弦振动波形的幅值,f为正弦振动波形的频率,为正弦振动波形的初 始相位; 0052 对步骤2)中输入到计算机数据处理设备的波形数据序列y i 进行数学描述,为 0053 y i y(t i )y(i-1),i1,2,n 0054 其中为采样时间间隔且1/v; 0055 3.2计算机数据处理设备对波形数据序列y i (i1,2,n)进行正弦波形四参 数拟合,获得拟合信号: 0056 0057 其中,为正弦振动波形的幅度拟合值,为正弦振动波形的角频率拟合值,为正 弦振。

25、动波形的初始相位拟合值,为正弦振动波形的直流分量值拟合值; 0058 正弦振动波形频率拟合值为: 0059 0060 4)根据步骤3)的拟合处理结果,以虚拟仪器方式显示超低频振动参数。 0061 上述步骤3.2中进行正弦波形四参数拟合的过程为: 0062 (1)设定初始参数,具体包括:待测的正弦波频率预期值为f 0 ,则 0 2f 0 / v,波形采集序列所含信号不足一个周期且个数为p,波形采集序列占用时间长度为,则 f 0 1/,选取另一个足够小的正数因子q(此处选取q110 -5 ),使得被估计的正弦频 率f 0 q/,因而f 0 q/,2/; 0063 (2)设定拟合迭代停止条件中收敛参。

26、数为一个接近于0的足够小的正数h e ,此处 选取h e 110 -20 ; 0064 (3)从已知采样时刻t 1 ,t 2 ,.,t n 得到的正弦波采集样本为y 1 ,y 2 ,.,y n ,用计 点法获得信号波形占用时间长度为(n-1)/v,确定目标频率f 0 的存在区间q/,2/ ; 0065 (4)确定迭代左边界频率f L q/和迭代左边界角频率 L 2f L /v以及 迭代右边界频率f R 2/和迭代右边界角频率 R 2f R /v,令中值角频率 M ( R + L )/2; 0066 (5)在迭代左边界角频率 L 、迭代右边界角频率 R 和中值角频率 M 上分别利 用频率已知的三。

27、参数正弦波形拟合公式计算各自的拟合残差( L )、( R )和( M ); 0067 (6)若( L )( M ),则令迭代右边界角频率 R M ,迭代左边界角频 说 明 书CN 102661782 A 5/5页 8 率 L 不变,重复执行步骤(5)步骤(6),其中为判据因子,的取值范围为11.5; 0068 (7)若( L )( M ),则必有 R 2 0 ,确定左边界角频率为 L ,右边 界角频率 R ,按照优选法原则,选取两个中值角频率分别为 0069 M L +0.618( R - L ) 0070 和 0071 T R -0.618( R - L ); 0072 (8)在左边界角频率。

28、 L 上执行频率已知的三参数正弦曲线拟合获得拟合幅度A L 、 拟合相位拟合直流分量D L 、拟合残差 L ,在右边界角频率 R 上执行频率已知的三参 数正弦曲线拟合获得拟合幅度A R 、拟合相位拟合直流分量D R 、拟合残差 R ,在中值角频 率 M 上执行频率已知的三参数正弦曲线拟合获得拟合幅度A M 、拟合相位拟合直流分量 D M 、拟合残差 M ,在中值角频率 T 上执行频率已知的三参数正弦曲线拟合获得拟合幅度 A T 、拟合相位拟合直流分量D T 、拟合残差 T ; 0073 (9)若 M T ,则最小拟合残差 M ,有 0 T , R ,参量更新为 L T , T M , M L 。

29、+0.618( R - L );若 M T ,则最小拟合残差 T ,有 0 L , M ,参量更新为 R M , M T , T R -0.618( R - L ); 0074 (10)判定若|( M - T )/ T |h e ,则重复步骤(5)步骤(10);若|( M - T )/ T |h e ,则停止迭代,此时又分为两种情况:若最小拟合残差 T ,获得四参数拟合正 弦曲线参数为,拟合过程结束;若最小拟合残差 M , 获得四参数拟合正弦曲线参数为,拟合过程结束。 0075 一种超低频振动参数量值的快速测量装置,用于实现上述超低频振动参数量值的 快速测量方法,包括振动传感设备、数据采集系统。

30、和计算机数据处理设备,其中待测的表征 超低频振动的位移、速度或者加速度物理信号进入振动传感设备后转换成电信号波形输出 给数据采集系统,数据采集系统再将电信号波形转换为波形数字序列后送入计算机数据处 理设备,计算机数据处理设备再将输入的数字序列进行处理,得到超低频振动的幅值和频 率测量结果并对结果进行输出。 0076 所述振动传感设备可采用如下两种技术方案:一种是光电测量方案,即振动传感 设备由激光干涉仪、信号采集处理电路和波形输出电路组成,此时超低频振动参数量值的 快速测量装置如图1所示;第二种是机电测量方案,即振动传感设备由振动传感器和波形 输出电路组成,此时超低频振动参数量值的快速测量装置。

31、如图2所示。 0077 目前用于地震监测所用的超低频振动台,其输出频率下限值可以达到1mHz,一个 波形周期达到1000s,约17分钟,使用通常的振动测量方法进行正常测量,17分钟以内不会 有任何响应,更不能进行快速测量其相位差值。将本发明的超低频振动的快速测量方法及 装置用于解决该问题,通过波形测量手段,在仅仅采集十分之一个周期波形条件下,在2分 钟之内即可以快速获得测量结果,在噪声信号幅度比为2的情况下,可以获得幅度估计误 差优于3、频率估计误差优于2、相位估计误差优于0.6的效果。 0078 以上所述为本发明的较佳实施例而已,本发明不应该局限于该实施例和附图所公 开的内容。凡是不脱离本发明所公开的精神下完成的等效或修改,都落入本发明保护的范 围。 说 明 书CN 102661782 A 1/1页 9 图1 图2 说 明 书 附 图CN 102661782 A 。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 物理 > 测量;测试


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1