CN107014367B - 一种优化光纤陀螺动态性能的方法 - Google Patents

一种优化光纤陀螺动态性能的方法 Download PDF

Info

Publication number
CN107014367B
CN107014367B CN201710228060.2A CN201710228060A CN107014367B CN 107014367 B CN107014367 B CN 107014367B CN 201710228060 A CN201710228060 A CN 201710228060A CN 107014367 B CN107014367 B CN 107014367B
Authority
CN
China
Prior art keywords
loop
closed
gyro
gyroscope
fiber
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
CN201710228060.2A
Other languages
English (en)
Other versions
CN107014367A (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.)
707th Research Institute of CSIC
Original Assignee
707th Research Institute of CSIC
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 707th Research Institute of CSIC filed Critical 707th Research Institute of CSIC
Priority to CN201710228060.2A priority Critical patent/CN107014367B/zh
Publication of CN107014367A publication Critical patent/CN107014367A/zh
Application granted granted Critical
Publication of CN107014367B publication Critical patent/CN107014367B/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
    • G01C19/00Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/58Turn-sensitive devices without moving masses
    • G01C19/64Gyrometers using the Sagnac effect, i.e. rotation-induced shifts between counter-rotating electromagnetic beams
    • G01C19/72Gyrometers using the Sagnac effect, i.e. rotation-induced shifts between counter-rotating electromagnetic beams with counter-rotating light beams in a passive ring, e.g. fibre laser gyrometers

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Electromagnetism (AREA)
  • Power Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)

Abstract

本发明涉及一种优化光纤陀螺动态性能的方法,包括依次连接sagnac效应、光电转换、滤波放大、AD转换、余弦响应偏置调制与解调、一次积分数字存储、数字截取以及依次连接二次积分、D/A转换、功率放大、功率放大、相位调制和线圈延迟,其处理步骤为:建立陀螺简化闭环模型;建立陀螺精确闭环模型;确定陀螺稳定及超调条件;为光纤陀螺选取闭环参数并通过仿真预估陀螺动态性能。本发明设计合理,其通过精确闭环模型可以实现陀螺在振动过程中零偏效应的预估,并为优化陀螺闭环参数提供参考,提高光纤陀螺动态性能。

Description

一种优化光纤陀螺动态性能的方法
技术领域
本发明属于光纤陀螺技术领域,尤其是一种优化光纤陀螺动态性能的方法。
背景技术
光纤陀螺是一种新型的全固态惯性仪表,具有成本低、寿命长、动态范围大和可靠性高等优点,已被广泛的应用于军用及民用领域。当光纤陀螺在振动过程中(无论是随机还是扫频)陀螺的零偏(均值)发生微小偏移,严重降低光纤陀螺的精度。原因是因为光纤陀螺回路中存在延迟、超调等,反馈信号总是滞后Sagnac相位。当由于加速或减速,角速率发生变化时,反馈相位总是试图跟上实际反馈相位,存在一个瞬时残余误差信号。这个残余误差信号在高频振动下变得很大。在偏置相位处于过调制状态下,“残余瞬时信号”很大时,正、负解调周期的采样值将不对称(符号相反,幅值不同),使光纤陀螺的解调环节具有与正弦输入同频的可变增益系数,正弦瞬态误差信号与同频的传函可变增益系数相乘,产生振动零偏效应。
光纤陀螺不仅可以用于测量角速率,还可以用于测量旋转角位移,因此,要求光纤陀螺在整个动态范围内都具有良好的精度。对于同样的标度误差,大输入角速率下,测量误差更大,因此,必须在整个动态范围内有一个精确而稳定的线性标度因数,实现这一标度因数的方法是采用闭环检测方案。在光纤陀螺简化闭环模型中,由于解调模块被简化成了一个增益固定的环节,因此,这一过程无法体现出来。传统的简化模型无法解释陀螺在振动过程中出现的陀螺零偏效应。
发明内容
本发明的目的在于克服现有技术的不足,提供一种设计合理、能够有效提高光纤陀螺的精度和环境适应性的优化光纤陀螺动态性能的方法。
本发明解决现有的技术问题是采取以下技术方案实现的:
一种优化光纤陀螺动态性能的方法,光纤陀螺闭环模型包括依次连接的前向通道环节:sagnac效应、光电转换、滤波放大、AD转换、余弦响应偏置调制与解调、一次积分数字存储、数字截取以及依次连接的后向通道环节:二次积分、D/A转换、功率放大、相位调制和线圈延迟,该方法包括以下步骤:
步骤1、建立陀螺简化闭环模型,该陀螺简化闭环模型计算闭环传递函数在Z域的表达式为:
Figure GDA0002525674590000021
其中,g为前向通道增益;h为后向通道增益;τ为光在环圈中的传播时间;每延迟一个τ的时间,Z域表达为z-1;p为前向通道线路与程序的总的延迟时间除以τ的值;q为后向通道线路与程序的总的延迟时间除以τ的值;
陀螺简化闭环模型中的余弦响应偏置调制与解调环节为:
Figure GDA0002525674590000022
I0为探测器接收到无偏置调制时的光功率;φb为调制相位;
步骤2、建立陀螺精确闭环模型,该精确闭环模型中的余弦响应偏置调制与解调环节为:
Figure GDA0002525674590000023
I0为探测器接收到无偏置调制时的光功率;φ(t)和φ(t-τ)分别为当前陀螺敏感到的相位和一个τ之前敏感到的相位;φb为调制相位;
步骤3、确定陀螺稳定及超调条件,其中确定陀螺稳定的方法为:针对不同的M值求解zM·(z-1)+g·h=0,使得|z|<1的g·h值的范围即为满足系统稳定的增益约束条件,运用MATLAB计算出不同延迟指数M下,保证陀螺闭环系统稳定的回路增益的最大值,画出光纤闭环回路稳定区域与回路总增益g·h和回路延迟指数M的关系曲线;确定超调条件的方法为:运用MATLAB仿真精确闭环模型绘制不同延迟时间下不同增益的幅频特性曲线,画出确保闭环传递函数幅频特性曲线不出尖峰的回路总增益g·h和回路延迟指数M的关系曲线;
步骤4、为光纤陀螺选取闭环参数并通过仿真预估陀螺动态性能,具体方法为:分析不同类型陀螺的相应闭环回路的延迟时间,将具体延迟时间带入到精确闭环模型中,确定其既满足稳定条件又不出现超调现象的闭环增益的范围,根据陀螺实际振动输出要求选择合适的增益。
本发明的优点和积极效果是:
1、本发明针对光纤陀螺通过建立简化的闭环模型,并在简化模型的基础上建立精确闭环模型,根据精确模型判断陀螺稳定工作条件、绘制幅频特性及相频特性曲线;同时,利用精确模型对陀螺进行动态仿真预估陀螺在不同闭环参数下动态性能,最终为陀螺设计参数的选择提供依据,通过调整陀螺闭环参数优化陀螺动态性能。
2、本发明通过精确闭环模型可以实现陀螺在振动过程中零偏效应的预估,并为优化陀螺闭环参数提供参考,提高光纤陀螺动态性能,并达到对陀螺的动态性能进行准确预判。
附图说明
图1是光纤陀螺简化闭环模型示意图;
图2是光纤陀螺简化模型中的解调环节原理图;
图3是光纤陀螺精确模型中的解调环节原理图;
图4是光纤陀螺精确闭环模型示意图;
图5是针对不同延迟时间M下陀螺输出稳定的增益范围示意图;
图6是闭环模型幅频特性曲线出现超调量示意图;
图7是针对不同延迟时间M下陀螺输出稳定条件及超调条件示意图;
图8是陀螺在不同频率的动态输入下输出情况仿真示意图。
具体实施方式
以下结合附图对本发明实施例做进一步详述:
一种优化光纤陀螺动态性能的方法,如图1所示,包括以下步骤:
步骤1、建立陀螺简化闭环模型。
图1给出了闭环模型各环节数学表达式。前向通道用G表示,各环节用G0、G1……G6表示,后向通道用H表示,各环节为H1、H2……H5。其中,前向通道的各个环节依次为:sagnac效应(G0)、光电转换(G1)、滤波放大(G2)、AD转换(G3)、余弦响应偏置调制与解调(G4)、一次积分数字存储(G5)和数字截取(G6);后向通道的各个环节依次为:二次积分(H1)、D/A转换(H2)、功率放大(H3)、相位调制(H4)和线圈延迟(H5)。
计算闭环传递函数为:
Figure GDA0002525674590000041
将公式(1)化简为:
Figure GDA0002525674590000042
其中,
Figure GDA0002525674590000043
g为前向通道增益;
Figure GDA0002525674590000044
h为后向通道增益。如果考虑到前向通道和反馈通道的其它附加延迟,式(2)可以写成:
Figure GDA0002525674590000045
其中(p+2)·τ表示前向通道的时间延迟,q·τ表示反馈通道的时间延迟,τ为光在环圈中的传播时间。
步骤2、建立精确闭环模型。
当输入信号φ为一个快速变化的信号时,闭环系统存在一定的反馈延迟,此时解调环节不能简单的看做
Figure GDA0002525674590000046
如图2所示。为了利用闭环模型在仿真中重现陀螺在振动实验中出现的零偏效应,需要重新考虑闭环模型G4环节即余弦响应偏置调制与解调环节,该环节精确表示为
Figure GDA0002525674590000047
如图3所示。此时精确闭环模型如图4所示。
步骤3、确定陀螺稳定及超调条件。
根据步骤1中公式(1),Z域的特征方程为:
z-1+g·h·z-(p+q+2)=0 (4)
令p+q+2=M,则式(4)变为:
zM·(z-1)+g·h=0 (5)
根据稳定性判据:|z|<1,则系统一定是稳定的。因而针对不同的M值求解式(5),使得|z|<1的g·h值的范围即为满足系统稳定的增益约束条件。运用MATLAB计算出不同延迟指数M下,保证陀螺闭环系统稳定的回路增益的最大值,画出光纤闭环回路稳定区域与回路总增益g·h和回路延迟指数M的关系曲线,如图5所示。
在动态输入条件下,回路跟踪误差大固然是振动过程中陀螺产生零偏的主因,但陀螺回路参数设计不合理,造成闭环传递函数幅频特性存在超调,通常会对振动零偏效应起推波助澜(或放大)作用。因此需要确定陀螺幅频特性曲线(如图6)超调条件,避免超调。运用MATLAB仿真精确闭环模型绘制不同延迟时间下不同增益的幅频特性曲线,画出确保闭环传递函数幅频特性曲线不出尖峰的回路总增益g·h和回路延迟指数M的关系曲线,如图7所示。
步骤4、在限制范围内为光纤陀螺选取合适的闭环参数并通过仿真预估陀螺动态性能。
针对不同类型陀螺,分析其相应的闭环回路的延迟时间,将具体延迟时间带入到精确闭环模型中,确定其既满足稳定条件又不出现超调现象的闭环增益的范围。以某型陀螺参数为例,其闭环增益符合限制条件的最大值为g·h=0.256,分别将该陀螺闭环增益设为g·h,g·h/2,g·h/4,g·h/8带入到精确闭环模型中仿真陀螺在不同频率的动态输入下输出情况,绘制图形如图8所示。根据陀螺实际振动输出要求选择合适的增益。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。

Claims (1)

1.一种优化光纤陀螺动态性能的方法,光纤陀螺闭环模型包括依次连接的前向通道环节:sagnac效应、光电转换、滤波放大、AD转换、余弦响应偏置调制与解调、一次积分数字存储、数字截取以及依次连接的后向通道环节:二次积分、D/A转换、功率放大、相位调制和线圈延迟,其特征在于,该方法包括以下步骤:
步骤1、建立陀螺简化闭环模型,该陀螺简化闭环模型计算闭环传递函数在Z域的表达式为:
Figure FDA0002525674580000011
其中,g为前向通道增益;h为后向通道增益;τ为光在环圈中的传播时间;每延迟一个τ的时间,Z域表达为z-1;p为前向通道线路与程序的总的延迟时间除以τ的值;q为后向通道线路与程序的总的延迟时间除以τ的值;
陀螺简化闭环模型中的余弦响应偏置调制与解调环节为:
Figure FDA0002525674580000012
I0为探测器接收到无偏置调制时的光功率;φb为调制相位;
步骤2、建立陀螺精确闭环模型,该精确闭环模型中的余弦响应偏置调制与解调环节为:
Figure FDA0002525674580000013
I0为探测器接收到无偏置调制时的光功率;φ(t)和φ(t-τ)分别为当前陀螺敏感到的相位和一个τ之前敏感到的相位;φb为调制相位;
步骤3、确定陀螺稳定及超调条件,其中确定陀螺稳定的方法为:针对不同的M值求解zM·(z-1)+g·h=0,使得|z|<1的g·h值的范围即为满足系统稳定的增益约束条件,运用MATLAB计算出不同延迟指数M下,保证陀螺闭环系统稳定的回路增益的最大值,画出光纤闭环回路稳定区域与回路总增益g·h和回路延迟指数M的关系曲线;确定超调条件的方法为:运用MATLAB仿真精确闭环模型绘制不同延迟时间下不同增益的幅频特性曲线,画出确保闭环传递函数幅频特性曲线不出尖峰的回路总增益g·h和回路延迟指数M的关系曲线;
步骤4、为光纤陀螺选取闭环参数并通过仿真预估陀螺动态性能,具体方法为:分析不同类型陀螺的相应闭环回路的延迟时间,将具体延迟时间带入到精确闭环模型中,确定其既满足稳定条件又不出现超调现象的闭环增益的范围,根据陀螺实际振动输出要求选择合适的增益。
CN201710228060.2A 2017-04-10 2017-04-10 一种优化光纤陀螺动态性能的方法 Active CN107014367B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710228060.2A CN107014367B (zh) 2017-04-10 2017-04-10 一种优化光纤陀螺动态性能的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710228060.2A CN107014367B (zh) 2017-04-10 2017-04-10 一种优化光纤陀螺动态性能的方法

Publications (2)

Publication Number Publication Date
CN107014367A CN107014367A (zh) 2017-08-04
CN107014367B true CN107014367B (zh) 2020-09-01

Family

ID=59445200

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710228060.2A Active CN107014367B (zh) 2017-04-10 2017-04-10 一种优化光纤陀螺动态性能的方法

Country Status (1)

Country Link
CN (1) CN107014367B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108120452B (zh) * 2017-12-21 2021-12-24 吉林大学 Mems陀螺仪动态数据的滤波方法
CN110687780B (zh) * 2019-06-18 2022-02-08 蓝箭航天空间科技股份有限公司 运载火箭速率陀螺增益自适应调整方法、调整系统及存储介质
CN111044081B (zh) * 2020-01-03 2022-07-26 中国船舶重工集团公司第七0七研究所 一种自适应小型化光纤环圈测试系统及测试方法
CN112815935A (zh) * 2020-12-31 2021-05-18 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) 光纤陀螺动态特性评估装置及方法
CN113074712B (zh) * 2021-03-23 2023-01-24 安徽华驰动能科技有限公司 一种光纤陀螺环路增益优化方法
CN113074713B (zh) * 2021-03-23 2023-01-24 安徽华驰动能科技有限公司 一种光纤陀螺高动态控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07134039A (ja) * 1993-11-10 1995-05-23 Tokimec Inc 光ファイバジャイロ
CN1510561A (zh) * 2002-12-23 2004-07-07 ƽ 用于设计和使用工作在调制域的模拟电路的系统和方法
CN103984830A (zh) * 2014-05-23 2014-08-13 湖北三江航天红峰控制有限公司 一种光纤陀螺数字闭环控制逻辑的仿真方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7505139B2 (en) * 2007-03-06 2009-03-17 Bergh Ralph A Signal processing for a Sagnac interferometer

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07134039A (ja) * 1993-11-10 1995-05-23 Tokimec Inc 光ファイバジャイロ
CN1510561A (zh) * 2002-12-23 2004-07-07 ƽ 用于设计和使用工作在调制域的模拟电路的系统和方法
CN103984830A (zh) * 2014-05-23 2014-08-13 湖北三江航天红峰控制有限公司 一种光纤陀螺数字闭环控制逻辑的仿真方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"反馈延迟对光纤陀螺振动误差特性的影响";潘雄;《红外与激光工程》;20140831;第43卷(第8期);正文第2608页 *
"采用等效正弦信号输入的光纤陀螺带宽测试方法";陈馨;《惯性技术发展动态发展方向研讨会文集》;20161231;正文第14-17页 *
"闭环模型及参数对光纤陀螺动态误差特性的影响";张书颖;《惯性技术发展动态发展方向研讨会文集》;20161231;正文第97-104页 *

Also Published As

Publication number Publication date
CN107014367A (zh) 2017-08-04

Similar Documents

Publication Publication Date Title
CN107014367B (zh) 一种优化光纤陀螺动态性能的方法
CN109990773B (zh) 一种干涉型光纤陀螺环路增益的检测与闭环控制系统及控制方法
CN111174773B (zh) 一种多相调制解调的光纤陀螺多闭环算法
US9915550B2 (en) Method and apparatus for data fusion of a three-axis magnetometer and three axis accelerometer
CN110375727B (zh) 一种闭环光纤陀螺信号调制方法
CN108168537B (zh) 基于正交解调的谐振式光学陀螺的检测系统和方法
CN104236537A (zh) 基于强度调制器的光源强度噪声抑制数字双闭环方法
CN111220142B (zh) 一种新型超高精度光纤陀螺仪的信号处理方法及装置
CN115143948B (zh) 一种基于光纤陀螺本征频率实时补偿标度因数的方法
US20120312095A1 (en) Gyroscope dynamic motor amplitude compensation for enhanced rate estimation during startup
CN104266663A (zh) 一种运动平台光电跟踪系统二级稳定扰动解耦方法
CN101408425B (zh) 一种光纤陀螺本征频率的跟踪方法
CN108444464B (zh) 一种光纤陀螺仪中的光源相对强度噪声抑制方法
CN112697124B (zh) 闭环谐振式光学陀螺的方波正交解调实现方法及装置
CN105674976B (zh) 光纤陀螺调制解调方法、提高标度因数稳定性方法及装置
CN112815935A (zh) 光纤陀螺动态特性评估装置及方法
CN103471593B (zh) 一种基于gps信息的惯性导航系统测量误差修正方法
CN116184815A (zh) 消除光纤陀螺复位误差的闭环控制方法及系统
CN109459012A (zh) 一种光纤陀螺2π电压参数全数字伺服调整装置及方法
CN116558548A (zh) 一种高q值mems谐振器的稳定幅度控制系统
CN113074713B (zh) 一种光纤陀螺高动态控制方法
Liu et al. A Phase Compensation Method for MEMS Quadruple Mass Gyroscope in Zero Bias Drift
CN115480417A (zh) 一种电光调制器的偏压控制电路系统
CN113063446B (zh) 一种基于测控电路的硅微陀螺标度因数非线性补偿方法
CN104713538A (zh) 一种闭环光纤陀螺最小拍控制方法

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