CN109959391B - 一种用于光纤陀螺动态性能评测的信号同步方法 - Google Patents

一种用于光纤陀螺动态性能评测的信号同步方法 Download PDF

Info

Publication number
CN109959391B
CN109959391B CN201910257472.8A CN201910257472A CN109959391B CN 109959391 B CN109959391 B CN 109959391B CN 201910257472 A CN201910257472 A CN 201910257472A CN 109959391 B CN109959391 B CN 109959391B
Authority
CN
China
Prior art keywords
fiber
signal
optic gyroscope
signals
gyroscope
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
CN201910257472.8A
Other languages
English (en)
Other versions
CN109959391A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201910257472.8A priority Critical patent/CN109959391B/zh
Publication of CN109959391A publication Critical patent/CN109959391A/zh
Application granted granted Critical
Publication of CN109959391B publication Critical patent/CN109959391B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)

Abstract

本发明公开了一种用于光纤陀螺动态性能评测的同步信号采集方法。包括如下步骤:(1)将光纤陀螺固定在角振动台上,利用校准过的高精度编码器输出转速信息;(2)同步信号采集模块实现对光纤陀螺和编码器数据的采集;(3)利用本发明提出的方法计算出延迟,并对光纤陀螺采集程序进行反馈控制;(4)同步信号采集模块将同步的两路信号打包发送至上位机。本发明可以实现对两路数据的同步接收,适用于任意频率的光纤陀螺和角振动台编码器,具有普适性;实验结果显示本方法可以大幅降低两信号间同步误差。本发明可以很好地满足计量领域对多路信号同步性的要求,有广泛的应用前景。

Description

一种用于光纤陀螺动态性能评测的信号同步方法
技术领域
本发明涉及信号同步方法,具体涉及到一种用于光纤陀螺动态性能评测的信号同步方法。
技术背景
光纤陀螺是基于Sagnac效应的角速率传感器,广泛应用于军事、姿态控制、航空航天等领域。随着现代飞航武器对陀螺工作环境的要求越来越高,经常处于振动、颠簸、急转弯等情景,在这种高动态环境下,光纤陀螺的动态性能成为影响精度的重要原因。因此评定并改善动态特性已经成为光纤陀螺研究的一项重要内容。
目前常用的测试方法为利用角振动台提供摇摆,采集角振动台的编码器输出值作为参考信号,与陀螺输出值进行比较,实现对光纤陀螺动态性能的评定。但这两信号并非严格同步,由于光纤陀螺输出信号与编码器信号间存在延迟,这会对光纤陀螺动态性能评测结果造成较大误差。因此,必须使用信号同步方法,来实现陀螺与编码器信号的同步采集。
现有的信号同步方法有互相关函数法和插值法,但互相关函数法在低采样频率下所能达到的精度较低。在1000 Hz采样频率下,互相关函数法无法分辨1ms内的延时,而陀螺和编码器的采样频率通常在1k以下,所以这一方法并不适用。而插值法通过函数拟合近似求得同步时刻的采样值,精度的提高以运算量增加为代价,不适合工业使用,且不同算法的复杂程度和适用范围也不相同,不具有普适性。
基于以上原因,提出了一种用于光纤陀螺动态性能评测的同步信号采集方法,可以实现对两路数据的同步接收,适用于目前工业应用的陀螺和编码器的信号频率,且不失普适性;实验结果显示本方法可以大幅降低两信号间同步误差。本方法使用方便,具有普适性,可以很好地满足计量领域对多路信号同步性的要求,有广泛的应用前景。
发明内容
本发明的目的是针对现有信号同步采集方法的不足,提出一种精确控制光纤陀螺与角振动台编码器输出信号同步的方法,本发明适用于任意采样频率的光纤陀螺及角振动台,可以充分满足目前对多路信号间同步性的要求。
为了实现上述目的,本发明提供如下技术方案:
一种用于光纤陀螺动态性能评测的同步信号采集方法:包括以下步骤:
(1)角振动台以预设频率、振幅做简谐运动;
(2)由于角振动台做简谐运动,因此编码器输出的角度值可以表示为:
s1=A*sin(2π*f*t)
由于两路信号存在延时Δt,光纤陀螺输出的角度值为:
s2=A*sin[2π*f*(t+Δt)]
可以得到两信号同步误差为:
Figure BDA0002014196930000021
其中:s1为编码器输出的角度值,s2为光纤陀螺输出的角度值,A为角振动台预设的振幅,f为角振动台预设振动频率,Δt为两信号间的时间延时量;
可以看出,信号同步误差与振幅、振动频率、信号间延时相关,其中
Figure BDA0002014196930000022
分量中Δt只对该分量造成时序上的移动,对同步误差Δs幅值造成的影响可以忽略,因此可以将
Figure BDA0002014196930000023
视为Δs的幅值量;设同步误差Δs的幅值为A’,可以得到
Figure BDA0002014196930000024
(3)信号同步采集模块以相同的采样频率对编码器及光纤陀螺进行采样,其中编码器输出信号为角度,光纤陀螺输出信号为角速度,由此获得光纤陀螺输出信号系列Dc{m}和编码器输出信号系列Ds{m};
(4)编码器输出信号系列Ds{m}对应的角度值为s1,通过积分计算出光纤陀螺输出信号系列Dc{m}对应的角度值为s2,通过快速傅里叶变换出同步误差Δs的幅值A’,其过程为:
Δs信号可以表示为:Δs(t)=Δs(t+kT)
其中T为函数周期,k=0,1,2,…;
将其用如下的傅立叶级数表征:
Figure BDA0002014196930000031
其中a0、an、bn和cn分别为:
Figure BDA0002014196930000032
其快速傅里叶变换的表达式为
Figure BDA0002014196930000033
其中N为采样点数,i表示虚数,Δs(n)为Δs的离散信号。
快速傅里叶计算结果乘以2/N再求模即可得到基于连续信号傅立叶级数等效的各次谐波幅值,
Figure BDA0002014196930000034
为求振动频率f下的幅值,只需取
Figure BDA0002014196930000035
其中Fs为采样频率;
此时计算出的cn即为A’,最后通过公式:
Figure BDA0002014196930000036
计算出两信号间的时间延时量;
(5)信号同步采集模块根据计算出的延时结果,向光纤陀螺数据采集模块传输相应延时时间量。光纤陀螺采集模块提前或延迟相应时间向光纤陀螺下达发送数据的指令,并接受陀螺返回的相应数据,即可实现角编码器信号与光纤陀螺信号在时域上的精确对准,达到信号同步的目的;
(6)信号同步采集模块将采集到的两路同步信号打包发送至上位机。
与现有技术相比,本发明的优点在于:由于本发明采用了新的信号同步方法,相比于传统信号同步方法具有对采样频率要求低、具有普适性的特点,利用角编码器与光纤陀螺输出的信号进行对比分析,实现对光纤陀螺动态性能进行评测。解决了光纤陀螺动态性能评测技术中对于两路信号高同步性要求的问题,可以大幅提高两信号的同步性,减小同步误差对后续计量分析的影响。本方法使用方便,无需对原有的测试采集电路进行改动;可以有效地解决多路信号的不同步问题,对计量领域有很重大的意义。
附图说明
图1是测试系统图。
图2是本方法流程图。
图3是延时补偿前两路信号差值图。
图4是延时补偿后两路信号差值图。
图中:角振动台1,编码器2,光纤陀螺3,信号同步采集模块4,上位机5。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对发明内容做进一步说明。
测试系统如图1所示,包括角振动台1、光纤陀螺2、编码器3、信号同步采集模块4、上位机5;其中,光纤陀螺2固定在角振动台1表面,其敏感轴与角振动台1转轴平行;信号同步采集模块4同步采集光纤陀螺3与编码器2的输出数据,并根据本发明提出的算法计算出两信号的延时;信号同步采集模块4可以对角光纤陀螺采集模块2进行反馈,从而实现延时补偿;信号同步采集模块4将采集到的数据打包发送至上位机5。
一种用于光纤陀螺动态性能评测的信号同步方法,流程如图2所示,具体包括如下步骤:
(1)角振动台以预设频率、振幅做简谐运动;
(2)由于角振动台做简谐运动,因此编码器输出的角度值可以表示为:
s1=A*sin(2π*f*t)
由于两路信号存在延时Δt,光纤陀螺输出的角度值为:
s2=A*sin[2π*f*(t+Δt)]
可以得到两信号同步误差为:
Figure BDA0002014196930000051
其中:s1为编码器输出的角度值,s2为光纤陀螺输出的角度值,A为角振动台预设的振幅,f为角振动台预设振动频率,Δt为两信号间的时间延时量;
可以看出,信号同步误差与振幅、振动频率、信号间延时相关,其中
Figure BDA0002014196930000052
分量中Δt只对该分量造成时序上的移动,对同步误差Δs幅值造成的影响可以忽略,因此可以将
Figure BDA0002014196930000053
视为Δs的幅值量;设同步误差Δs的幅值为A’,可以得到
Figure BDA0002014196930000054
(3)信号同步采集模块以相同的采样频率对编码器及光纤陀螺进行采样,其中编码器输出信号为角度,光纤陀螺输出信号为角速度,由此获得光纤陀螺输出信号系列Dc{m}和编码器输出信号系列Ds{m};
(4)编码器输出信号系列Ds{m}对应的角度值为s1,通过积分计算出光纤陀螺输出信号系列Dc{m}对应的角度值为s2,通过快速傅里叶变换出同步误差Δs的幅值A’,其过程为:
Δs信号可以表示为:Δs(t)=Δs(t+kT)
其中T为函数周期,k=0,1,2,…;
将其用如下的傅立叶级数表征:
Figure BDA0002014196930000055
其中a0、an、bn和cn分别为:
Figure BDA0002014196930000061
其快速傅里叶变换的表达式为
Figure BDA0002014196930000062
其中N为采样点数,i表示虚数,Δs(n)为Δs的离散信号。
快速傅里叶计算结果乘以2/N再求模即可得到基于连续信号傅立叶级数等效的各次谐波幅值,
Figure BDA0002014196930000063
为求振动频率f下的幅值,只需取
Figure BDA0002014196930000064
其中Fs为采样频率;
此时计算出的cn即为A’,最后通过公式:
Figure BDA0002014196930000065
计算出两信号间的时间延时量;
(5)信号同步采集模块根据计算出的延时结果,向光纤陀螺数据采集模块传输相应延时时间量。光纤陀螺采集模块提前或延迟相应时间向光纤陀螺下达发送数据的指令,并接受陀螺返回的相应数据,即可实现角编码器信号与光纤陀螺信号在时域上的精确对准,达到信号同步的目的;
(6)信号同步采集模块将采集到的两路同步信号打包发送至上位机。
本发明所述的用于光纤陀螺动态性能评测的信号同步方法,可以大幅降低两信号间同步误差;且适用于任意采集频率下的光纤陀螺及编码器。本方法可以实现两路信号同步采集上传,满足光纤陀螺动态性能评定对信号同步性的需求;同样适用于其他应用场景的多路信号同步。
实施例:
下面结合附图对本发明一种实施例做进一步说明。
(1)将光纤陀螺3通过螺钉固定在角振动台1上,光纤陀螺3的敏感轴与角振动台1的转轴平行,接通电源预热5分钟。
(2)在上位机3上输入测试振幅10°,振动频率1Hz,信号采集频率2kHz。信号同步采集模块4启动并向陀螺发送接收指令;角振动台开始运动。
(3)信号同步采集模块4接收光纤陀螺3数据以及编码器2数据,同步误差如图3所示,采集2s数据后,利用最小二乘法拟合出同步误差的幅值A’,通过公式
Figure BDA0002014196930000071
计算出此时的延时时间量Δt=600ms。若直接使用此时收集到的两路信号进行陀螺动态性能评测,高达0.05°的同步误差会对结论造成不可忽视的影响。
(4)根据计算出的结果Δt=600ms,对陀螺进行相应的延时反馈控制;并将收到的两路同步信号打包发送至上位机5。
(5)如图4所示,此时两信号最大误差量为0.002°,若此误差全部由同步误差造成,由上述公式可以求得,此时对应的两信号间延时为25微秒。因此可以认为光纤陀螺与编码器这两路信号间延迟小于25微秒。
(6)结束测试,关闭系统。

Claims (1)

1.一种用于光纤陀螺动态性能评测的同步信号采集方法,其特征在于包括以下步骤:
(1)角振动台以预设频率、振幅做简谐运动;
(2)由于角振动台做简谐运动,因此编码器输出的角度值表示为:
s1=A*sin(2π*f*t)
由于两路信号存在延时,光纤陀螺输出的角度值为:
s2=A*sin[2π*f*(t+Δt)]
得到两信号的同步误差为:
Figure FDA0002599773740000011
其中:s1为编码器输出的角度值,s2为光纤陀螺输出的角度值,A为角振动台预设的振幅,f为角振动台预设振动频率,Δt为两信号间的时间延时量;
信号同步误差与振幅、振动频率、信号间的时间延时量相关,其中
Figure FDA0002599773740000012
分量中Δt只对该分量造成时序上的移动,对同步误差Δs幅值造成的影响忽略不计,因此将
Figure FDA0002599773740000013
视为Δs的幅值量;设同步误差Δs的幅值为A’,得到:
Figure FDA0002599773740000014
(3)信号同步采集模块以相同的采样频率对编码器及光纤陀螺进行采样,其中编码器输出信号为角度,光纤陀螺输出信号为角速度,由此获得光纤陀螺输出信号系列Dc{m}和编码器输出信号系列Ds{m};
(4)编码器输出信号系列Ds{m}对应的角度值为s1,通过积分计算出光纤陀螺输出信号系列Dc{m}对应的角度值为s2,通过快速傅里叶变换计算得出同步误差Δs的幅值A’,其过程为:
Δs信号表示为:Δs(t)=Δs(t+kT)
其中T为函数周期,k=0,1,2,…;
将其用如下的傅里叶级数表征:
Figure FDA0002599773740000015
Figure FDA0002599773740000021
其中a0、an、bn和cn分别为:
Figure FDA0002599773740000022
其快速傅里叶变换的表达式为:
Figure FDA0002599773740000023
其中N为采样点数,i表示虚数,Δs(n)为Δs的离散信号,ΔS(K)为Δs的傅里叶级数的系数;
快速傅里叶计算结果乘以2/N再求模,得到基于连续信号傅里叶级数等效的各次谐波幅值:
Figure FDA0002599773740000024
为求振动频率f下的幅值,只需取
Figure FDA0002599773740000025
其中Fs为采样频率;
此时计算出的cn即为A’,最后通过公式:
Figure FDA0002599773740000026
计算出两信号间的时间延时量;
(5)信号同步采集模块根据计算结果,向光纤陀螺数据采集模块传输相应时间延时量;光纤陀螺数据采集模块提前或延迟相应时间向光纤陀螺下达发送数据的指令,并接受陀螺返回的相应数据,即可实现角编码器信号与光纤陀螺信号在时域上的精确对准,达到信号同步的目的;
(6)信号同步采集模块将采集到的两路同步信号打包发送至上位机。
CN201910257472.8A 2019-04-01 2019-04-01 一种用于光纤陀螺动态性能评测的信号同步方法 Active CN109959391B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910257472.8A CN109959391B (zh) 2019-04-01 2019-04-01 一种用于光纤陀螺动态性能评测的信号同步方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910257472.8A CN109959391B (zh) 2019-04-01 2019-04-01 一种用于光纤陀螺动态性能评测的信号同步方法

Publications (2)

Publication Number Publication Date
CN109959391A CN109959391A (zh) 2019-07-02
CN109959391B true CN109959391B (zh) 2020-11-24

Family

ID=67025523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910257472.8A Active CN109959391B (zh) 2019-04-01 2019-04-01 一种用于光纤陀螺动态性能评测的信号同步方法

Country Status (1)

Country Link
CN (1) CN109959391B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112067892B (zh) * 2020-08-06 2023-10-13 北京自动化控制设备研究所 一种信号转换同步装置及使用其的信号转换同步方法
CN113237491A (zh) * 2021-04-22 2021-08-10 北京航天计量测试技术研究所 一种数字陀螺仪的频率特性测试装置及其测试方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101464150A (zh) * 2009-01-09 2009-06-24 哈尔滨工程大学 一种光纤陀螺动态性能的测试方法
CN105157694A (zh) * 2015-09-23 2015-12-16 浙江大学 一种准确评测光纤陀螺输出延迟时间的装置和方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8149417B2 (en) * 2010-01-27 2012-04-03 Honeywell International Inc. Synchronous radiation hardened fiber optic gyroscope
CN105091835B (zh) * 2014-05-07 2018-07-24 中国航空工业第六一八研究所 一种转台动态角位置误差测量精度提高装置及校正方法
CN104950169B (zh) * 2015-06-19 2017-08-15 浙江大学 一种高速光纤陀螺频率特性的测试方法与系统
CN107064559A (zh) * 2017-02-23 2017-08-18 西安航天精密机电研究所 一种基于角摇摆运动的sins加速度计频率特性测试方法
CN109443390B (zh) * 2018-12-05 2021-01-12 中国电子科技集团公司第十三研究所 一种数字陀螺仪相位延时测试系统
CN109405851B (zh) * 2018-12-14 2020-12-18 浙江大学 基于自准直指零仪和光纤陀螺的测试装置及测试方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101464150A (zh) * 2009-01-09 2009-06-24 哈尔滨工程大学 一种光纤陀螺动态性能的测试方法
CN105157694A (zh) * 2015-09-23 2015-12-16 浙江大学 一种准确评测光纤陀螺输出延迟时间的装置和方法

Also Published As

Publication number Publication date
CN109959391A (zh) 2019-07-02

Similar Documents

Publication Publication Date Title
US5060175A (en) Measurement and control system for scanning sensors
CN109959391B (zh) 一种用于光纤陀螺动态性能评测的信号同步方法
CN111879339A (zh) 一种mems陀螺仪温度误差补偿方法
CN107525524B (zh) 一种基于三轴同步转台的惯性导航系统时延确定方法
CN109323711A (zh) 一种陀螺仪模态反转零位自校正方法及系统
CN108681239B (zh) 一种两轴一体陀螺加速度计解耦伺服控制回路系统及方法
CN111024119B (zh) 一种三轴mems陀螺仪快速标定方法
CN110108299A (zh) 一种硅微机械陀螺仪标度因数在线自校准系统
CN110426033A (zh) 基于松耦合imu阵列导航系统的时间同步算法
CN111624671B (zh) 旋转加速度计重力梯度仪重力梯度解调相位角确定方法及装置
CN100533063C (zh) 光纤陀螺热噪声、散粒噪声及强度噪声分离的方法
CN108469849A (zh) 一种随机角振动控制方法
CN109931959B (zh) 硅微陀螺仪正交误差校正方法
CN112577514A (zh) 一种mems惯性器件的标定方法
Xing et al. Offline calibration for MEMS gyroscope g-sensitivity error coefficients based on the Newton iteration and least square methods
CN105509867A (zh) 一种利用光纤陀螺测量卫星结构微角振动方法
CN107843389B (zh) 用于扭摆系统的冲量测量误差评定方法
CN109579878B (zh) 基于频扫激励信号的惯性元件误差模型快速辨识方法
CN110440779B (zh) 一种高动态mems谐振式陀螺仪力平衡控制方法
CN112304336A (zh) 一种高频角振动转台控制方法
CN109212629A (zh) 旋转加速度计重力梯度仪角运动误差补偿装置及方法
CN110672121B (zh) 一种控制力矩陀螺框架动态响应测试方法及系统
CN109212619A (zh) 一种旋转加速度计重力梯度仪线运动误差补偿装置及方法
Jatiningrum et al. Modelling an angular accelerometer using frequency-response measurements
RU2783617C1 (ru) Способ компенсации движения виброподставки лазерного гироскопа и устройство для его реализации

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