《一种超低频振动参数量值的快速测量方法及装置.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 。