CN113074712A - 一种光纤陀螺环路增益优化方法 - Google Patents
一种光纤陀螺环路增益优化方法 Download PDFInfo
- Publication number
- CN113074712A CN113074712A CN202110307328.8A CN202110307328A CN113074712A CN 113074712 A CN113074712 A CN 113074712A CN 202110307328 A CN202110307328 A CN 202110307328A CN 113074712 A CN113074712 A CN 113074712A
- Authority
- CN
- China
- Prior art keywords
- fiber
- loop
- optic gyroscope
- closed
- gain
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000005457 optimization Methods 0.000 title claims abstract description 29
- 239000013307 optical fiber Substances 0.000 title description 6
- 238000001514 detection method Methods 0.000 claims abstract description 66
- 238000013461 design Methods 0.000 claims abstract description 45
- 239000000835 fiber Substances 0.000 claims abstract description 38
- 230000007613 environmental effect Effects 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 25
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000001629 suppression Effects 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000003111 delayed effect Effects 0.000 claims description 2
- 241000764238 Isis Species 0.000 claims 1
- 230000008030 elimination Effects 0.000 claims 1
- 238000003379 elimination reaction Methods 0.000 claims 1
- 239000000284 extract Substances 0.000 claims 1
- 238000012360 testing method Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 7
- 230000006872 improvement Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 241000238876 Acari Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C19/00—Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
- G01C19/58—Turn-sensitive devices without moving masses
- G01C19/64—Gyrometers using the Sagnac effect, i.e. rotation-induced shifts between counter-rotating electromagnetic beams
- G01C19/72—Gyrometers 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
本发明公开了一种光纤陀螺环路增益优化方法,包括如下步骤:步骤一、设计检测方法从噪声背景中提取闭环误差信号;步骤二、基于闭环误差模型推导光纤陀螺的动态方程;步骤三、针对非线性条件,进行环路增益优化设计;步骤四、考虑噪声以及环境扰动,进一步优化设计光纤陀螺系统闭环参数;步骤五、基于上述优化结果进行检测电路与参数设计。本发明的方法,提高了闭环光纤陀螺系统的动态性能与检测精度,并且进一步提出了在外界扰动与噪声条件下的环路增益优化设计方法。
Description
技术领域
本发明属于光纤陀螺环路增益优化技术领域,涉及基于光纤陀螺闭环误差特点的环路增益分布原理及相关控制器硬件电路与控制参数设计方法。
背景技术
光纤陀螺是一种高精度的惯性传感器件,在制导、导航等军用与民用领域起着非常重要的作用。光纤陀螺基于Sagnac效应,Sagnac效应描述了光纤环的旋转使光相位发生变化的现象,当光纤环发生旋转时,其相位会产生一个正比于旋转角速度的变化值。本发明可以通过检测光相位的变化情况获得系统的角速度。
在光纤陀螺系统中,信号检测系统是重要组成部分。由于开环检测存在着精度较低、动态范围有限的缺陷,本发明采用光纤陀螺闭环检测方案,选用频率对应于光纤环本征频率的方波信号进行调制,其中τ为光的渡越时间。数字闭环检测实现了良好的线性与较大的动态范围,受到广泛应用。
如图1所示,光纤陀螺采用最小互易结构,主要包括光源、耦合器、Y波导、光纤环、探测器、前置放大器、A/D转换器、数字控制器与D/A转换器等重要组件。光源发出光传输到Y波导并在此起偏,之后以顺时针、逆时针方向进入光纤线圈。当系统转动时,会随之产生同旋转角速度成正比的Sagnac相位差,之后转换为光强的变化,再被探测器捕获。为提高系统灵敏度,采用方波调制的闭环方案。前置放大器放大微弱信号,由A/D转换器转换为离散数字量,数字信号处理电路进行解调,由D/A转换器形成阶梯波,作用在Y波导上最终形成反馈相位实现闭环工作。
光纤陀螺环路增益的优化对于高频角速度的高精度的测量意义十分重大,通过优化光纤陀螺环路增益,能满足光纤陀螺系统中的高角速度测量要求。通过优化控制算法,对光纤陀螺环路增益进行优化能提高光纤陀螺系统的动态性能并提高光纤陀螺系统的检测精度,因此环路增益的设计对于提高光纤陀螺的动态性能和检测精度具有重要意义。此外,检测电路环路增益的分布直接影响闭环误差信号的检测,这在相关研究中被忽略。
发明内容
为了解决上述技术问题,本发明提出一种光纤陀螺环路增益优化方法,对于满足光纤陀螺系统的高角速度测量要求,对于实现光纤陀螺的动态性能和检测精度有重大意义。本发明提出用于光纤陀螺系统的一种环路增益的优化设计方法,以及一种信号检测方法,以同时提高光纤陀螺闭环系统的快速追踪能力和检测精度。首先,本发明从理论上分析了光纤陀螺环路增益的分布对光纤陀螺性能的影响并建立了一种环路增益设计方法的理论模型。其次,根据环路增益的理论模型分析建立了闭环光纤系统的动态方程。第三,本发明给出了光纤陀螺环路增益闭环控制参数的设计优化方法,提高了光纤陀螺的动态性能。第四,针对外部的扰动和噪声,提出了光纤陀螺环路增益闭环参数的优化设计方法。最后,通过实验证明了该方法的有效性。
本发明的技术方案为:一种光纤陀螺环路增益优化方法,包括如下步骤:
步骤一、设计检测方法从噪声背景中提取闭环误差信号;
步骤二、基于闭环误差模型推导光纤陀螺的动态方程;
步骤三、针对非线性条件,进行环路增益优化设计;
步骤四、考虑噪声以及环境扰动,进一步优化设计光纤陀螺系统闭环参数;
步骤五、基于上述优化结果进行检测电路与参数设计。
有益效果:
本发明根据光纤陀螺闭环误差的特点,从理论上研究了环路增益的分布原理;基于获得的光纤陀螺闭环误差模型,设计了光纤陀螺的动态方程,从而分别提高了光纤陀螺闭环误差的高动态性能与外界干扰条件下的高检测精度。本发明实现了在较高的检测精度下同时具有较高的动态性能。采用本发明提供的光纤陀螺环路增益设计方法,光纤陀螺系统在动态性能与检测精度性能上都有较大改善;本发明所给出的光纤陀螺环路增益闭环控制参数的设计优化方法,能有效提高光纤陀螺的动态性能,在此基础上,针对外部的扰动和噪声,仅一步提出了光纤陀螺环路增益闭环参数的优化设计方法,使得系统在存在外部扰动与噪声的情况下依然能保持较好的性能。
附图说明
图1为基于Sagnac效应的闭环光纤陀螺原理方案;
图2为本发明闭环系统信号检测原理示意图;
图3为环路增益优化前光纤陀螺系统静态检测精度;
图4为环路增益优化后光纤陀螺系统静态检测精度。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅为本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域的普通技术人员在不付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
首先,如图1,明确闭环光纤陀螺的信号检测原理。如图2是闭环光纤陀螺系统信号检测原理示意图。光电检测器将光信号转换为电信号以进行进一步处理。闭环误差信号被调制为高频载波,以抑制检测电路中的1/f噪声。本发明设计了一种数字解调方案用于高精度地提取光纤陀螺闭环误差。闭环光纤陀螺的前向通道由光电探测器,前置放大器,A/D转换器和解调器组成。从闭环控制器输出的反馈相位被转换为模拟信号,并通过D/A转换器施加在Y波导上。
光电探测器探测到的干涉强度P可以表示为:
其中,P0是干扰强度的直流分量,是由Sagnac效应引起的相位延迟,φm(t)是调制相位,φf(t)是阶梯波反馈相位,τ是光通过到返回Y波导所用的传输时间。当光线通过并返回到Y波导时,反馈电压分别为φf(t-τ)和φf(t)。因此,本发明获得了反馈相位本发明定义为光纤陀螺的闭环误差,有
方波调制可以将余弦函数转换为正弦函数,以利用正弦函数在零点附近的近似线性来提高信号检测的灵敏度。传输时间τ决定了方波的调制周期T=2τ,调制信号施加到Y波导上的调制相位在2kτ≤t<(2k+1)τ时为而在(2k+1)τ≤t<2(k+1)τ时为0,其中k=0,1,…,n。
假设前向通道的增益为k1,本发明可以得到调制后的信号SA(t)为:
其中fm(t)表示幅度为±1的调制方波。本发明可以看到,调制方波将闭环误差调制为周期信号。然后,可以通过以下方式获得实际信号:
其中w(t),主要是白噪声,本发明认为w(t)是平均值为零的宽带噪声。
闭环检测方案将闭环误差转换为包含大量噪声的微弱的闭环误差信号。实际信号中包含的噪声给高精度地提取闭环误差信号带来了很大的困难,这进一步增加了设计信号检测方案时的复杂性。因此,本发明必须通过相应的控制算法解决改善动态性能与提高检测精度之间的矛盾。本发明设计的控制算法是通过线性矩阵不等式(LMI)方法获得的,本发明通过Lyapunov函数从理论上分析了光纤陀螺闭环增益的分配原则。根据闭环增益的分配原则以及光纤陀螺闭环误差的特性,优化了光纤陀螺环路增益的设计方法。
进一步,分析光纤陀螺闭环误差的相关检测方法,总结出了增加前向通道增益对于提高光纤陀螺的检测精度具有重要意义。考虑到调制的闭环误差是周期信号,并且检测到的闭环误差信号与噪声互不相关,因此可以设计相关检测方法来从强噪声背景中高精度地提取出微弱的闭环误差信号。由于解调周期仅为几微秒,因此本发明可以假设闭环误差在一个解调周期内保持不变。设计一个解调方波fd(t),该方波与fm(t)具有相同的幅度,频率和相位,并作为fm(t)的参考信号。进一步本发明可以得到相关检测的输出:
其中T表示方波的调制周期。
将调制方波和参考波信号可以展开为傅立叶级数:
因此,可以将相关检测方法应用于消除干涉强度的直流分量。
此外,由于等式(3)中假设w(t)是白噪声,并且有fm(t)fd(t)等于1。根据统计定律,由于w(t)是平均值为零的随机信号,可以得出为零。本发明可以推导出解调部分的总输出为从以上分析可以得出结论,通过增加前向通道的增益k1值,可以在抑制白噪声的同时提高闭环误差的检测精度。因此,闭环误差的信噪比的提升可以为随后的控制算法设计提供更多的便利。
之后,为了实现光纤陀螺的良好追踪性能和快速动态响应,本发明分析了系统参数的设计,以确保光纤陀螺系统满足快速稳定性。本发明将信噪比提升后的解调信号发送到控制器。可以得出闭环光纤陀螺的动态方程为:
本发明将控制器设计为二阶,并假设Kc=[kc1 kc2],kc1>0以及kc2>0。然后,考虑受控系统,光纤陀螺的动态方程可写为:
以下定理用于设计闭环检测参数k1,k2和闭环控制参数kc1,kc2,使上述系统动态方程(5)满足动态快速稳定条件。
证明:本发明为系统(5)选择以下Lyapunov函数,
易得:
本发明知道矩阵Ac的多项式特征值满足方程:
然后,本发明得出结论:系统(5)是指数稳定的。
定理1为闭环光纤陀螺系统参数的设计优化提供了充分的条件,从而保证了具有时滞的光纤陀螺的状态变量快速指数性收敛于零。由定理1,能发现 不仅同反馈控制变量kc1和kc2相关,还与环路增益,比如前向通道增益k1以及信号检测电路反馈通道增益k2相关。然后,本发明可以看到,系统的闭环性能取决于闭环检测的设计参数以及反馈控制参数的设计与优化。只要满足β=0,则系统(5)渐近稳定。因此,定理1为光纤陀螺的环路增益设计提供了理论指导。
进一步,考虑噪声及环境扰动影响下的环路增益方法。首先将光纤陀螺闭环检测系统离散化,并考虑没有外界环境扰动条件下的环路增益优化设计方法。
设计鲁棒控制算法,最关键的就是要求出其反馈增益矩阵Kc,用来保证状态变量x(k+1)的指数稳定性。x(k+1)经相应变换可表示成如下的差分方程:
非线性响应会影响到光纤陀螺的稳定性及动态性能,而且光路噪声的干扰会影响光纤陀螺的检测精度。基于此,针对光纤陀螺,本专利采用鲁棒控制算法。一方面对系统状态方程进行稳定性控制,通过反馈相位来近似抵消Sagnac相位,使系统始终工作在正弦函数的零点附近,提高系统的稳定性;另一方面对光纤陀螺系统的输出信号进行稳定性控制,以获取系统的指数稳定性并达到规定的噪声抑制水平。
x(k+1)=Ax(k)+Bk1sin(-k2Kcx(k))+Bk1sin(-k2Kcx(k-1))+Bw(k)
(8)
下面详细说明如何设计反馈增益矩阵Kc,首先引入以下定义:
定义1给定一常数0<α<1,若存在一常数κ≥0,如果控制器的噪声干扰w(k)=0时的状态变量满足:
定义2给定一常数γ>0,如果系统是指数稳定的,并且在零初始条件下对所有的非零噪声干扰w(k)∈L2[0,∞),满足则称光纤陀螺系统既满足指数稳定性的动态性能且具有γ抑制度的H∞噪声抑制性能。其中,w(s)和x(s)分别为系统噪声和系统状态。
首先,忽略噪声干扰w(k)的影响,通过设计控制器反馈增益矩阵保证光纤陀螺系统具有良好的快速动态跟踪性能。
定理2如果存在正定矩阵P,Q,R∈Rn×n,反馈增益矩阵Kc∈R1×n,及正常数α,ε1,ε2,使得下式成立:
证明:针对上述系统,有如下的分段李雅普诺夫方程:
V(k)=xT(k)Px(k)+xT(k-1)Qx(k-1)+ηT(k)Rη(k)
上式中P,Q的定义如定理1中所规定,η(k)=x(k)-x(k-1)。
计算ΔV(k)=V(k+1)-V(k),有如下方程:
ΔV(k)+αV(k)=xT(k+1)Px(k+1)-(1-α)xT(k)Px(k)+xT(k)Qx(k)
-(1-α)xT(k-1)Qx(k-1)+ηT(k+1)Rη(k+1)-(1-α)ηT(k)Rη(k)
可得到,对任意的常数ε1>0,ε2>0,有式:
此外,有下式:
式中λ1=minλ(P)是矩阵P的最小特征值,λ2=maxλ(P)+maxλ(Q)+4maxλ(R),maxλ(·)是最大特征值。
在实际工程应用中,光纤陀螺系统不可避免的会受到外部环境的扰动影响,设计满足H∞噪声抑制的控制器对提高系统的抗干扰能力具有重要的意义,并且能够提高光纤陀螺在实际应用中的鲁棒性。利用定理2的结果,讨论存在扰动的情况下光纤陀螺系统的性能。下述定理3给出了光纤陀螺系统获得指数稳定性及期望H∞性能的充分条件。
定理3如果存在正定矩阵P,Q,R∈Rn×n,反馈增益矩阵Kc∈R1×n,及正常数α,ε1,ε2,使得下式成立:
证明:证明方法类似于定理2有:
由φ+δ1 TPδ1+δ2 TRδ2<0,可得M<0。于是有:
ΔV(k)<-αV(k)-(xT(k)x(k)-γ2wT(k)w(k))
令Г(t)=xT(k)x(k)-γ2wT(k)w(k),则易得:
在零初始条件下有:
于是可以得出关于光纤陀螺系统具有噪声抑制度为γ的期望H∞性能。
接下来,研究用于实现闭环检测参数的硬件电路,并设计闭环控制参数。实现的闭环信号检测方法的硬件由光电检测器,开关电路,前置放大器,A/D转换器,控制器,积分器,D/A转换器和D/A驱动电路组成。由于数字方波对噪声不敏感,因此也可以在以数字方波为参考信号的数字场中实现调制和解调波。这对于改善所检测的闭环误差信号的信噪比是有利的。因此,本发明通过使用数字信号处理单元FPGA来实现解调和控制器。
当调制方波扫过余弦函数的零点时,调制信号在和之间会产生尖峰。为了满足检测光强度和转换效率一定乘积,光电检测器输出的峰值可以达到1V以上。当光电检测器的输出被直接放大K2次时,尖峰将超过A/D转换器的最大输入电压。因此,本发明需要设计一个模拟开关电路以消除检测器输出中的尖峰,以实现正向链增益k1的设计值。
在控制回路中,数字反馈由控制器计算并累加以形成反馈相位斜坡。然后在D/A转换器和驱动电路之后,将反馈相位斜坡施加到Y波导以生成反馈电压为了抵消本发明获得光纤陀螺输出的确定值由于D/A转换器的最大量程受限,数字相位斜波会在数字量为时定期重置,对应于光相位为2π时,其中nDA=16是D/A转换器的分辨率。这样,可通过来计算反馈链的增益。
根据本发明的实施例,基于上面的推导分析过程,提出一种用于光纤陀螺的环路增益的方法,提高了闭环光纤陀螺系统的动态性能与检测精度,并且进一步提出了在外界扰动与噪声条件下的环路增益优化设计方法。具体包括如下步骤:
第一步,设计相关检测方法来从强噪声背景中提取微弱的闭环误差信号。设计一个解调方波fd(t),该方波与fm(t)具有相同的幅度,频率和相位,并作为fm(t)的参考信号。获得相关检测的输出V0:
其中T表示方波的调制周期。
将调制方波和参考波信号展开为傅立叶级数:
因此,可以将相关检测方法应用于消除干涉强度的直流分量。
此外,由于等式(3)中假设w(t)是白噪声,并且知道fm(t)fd(t)等于1。然后,根据统计定律,由于w(t)是平均值为零的随机信号,可以得出 为零。进一步可以推导出解调部分的总输出为从以上分析可以得出结论,通过增加前向通道的k1值,可以在抑制白噪声的情况下提高闭环误差的检测精度。因此,闭环误差的信噪比的提升可以为随后的控制算法设计提供更多的便利。
基于闭环误差模型,光纤陀螺的动态方程可写为:
第三步,为了实现高动态性能,采用本发明提出的环路增益优化方法。
惯导系统要求光纤陀螺应该能够准确地测量高频角速度。因为高频角速度零点附近的信号等效于光纤陀螺的斜坡信号。当系统阶数大于两阶时,高频输入的稳态误差e(∞)为零。换句话说,当控制器被设计为具有两阶或两阶以上时,可以在光纤陀螺的稳态下精确地追踪输入角速度的变化。因此,本发明将控制器设计为二阶,并假设Kc=[kc1kc2],kc1>0以及kc2>0。然后,为了找到满足动态快速稳定条件(5)的闭环检测参数k1,k2和闭环控制参数kc1,kc2的设计原理,由定理1:
定理1为闭环光纤陀螺系统参数的设计原理提供了充分的条件,从而保证了具有时滞的光纤陀螺的状态变量快速指数性收敛于零。由定理1,能发现 不仅同反馈控制变量kc1和kc2相关,还与环路增益,比如前向通道增益k1以及信号检测电路反馈通道增益k2相关。然后,本发明可以看到,系统的闭环性能取决于闭环检测的设计参数以及反馈控制参数。只要满足β=0,则系统(5)渐近稳定。
第四步,考虑噪声以及环境扰动,进一步优化设计闭环参数。
在实际工程应用中,光纤陀螺系统不可避免的会受到外部环境的扰动影响,设计满足H∞噪声抑制的控制器对提高系统的抗干扰能力具有重要的意义,并且能够提高光纤陀螺在实际应用中的鲁棒性。由定理3:
定理3如果存在正定矩阵P,Q,R∈Rn×n,反馈增益矩阵Kc∈R1×n,及正常数α,ε1,ε2,使得下式成立:
其中,δ1=[A 0 Bk1 Bk1 B],δ2=[A-I 0 Bk1 Bk1 B],
利用定理3的结果,讨论存在扰动的情况下光纤陀螺系统的环路增益参数设计,可以得到光纤陀螺系统获得指数稳定性及期望H∞性能的充分条件。
第五步实现检测电路与参数设计。
实现检测电路并仿真波形。实现的闭环信号检测方法的硬件由光电检测器,开关电路,前置放大器,A/D转换器,控制器,积分器,D/A转换器和D/A驱动电路组成。40MHz振荡器用于生成时钟和参考信号。开关电路时钟,A/D采样时钟和D/A时钟分别为开关电路,A/D转换器和D/A转换器提供时钟信号。并在调制解调过程中应用了调制方波和解调方波。开关时钟是通过峰值消除开关电路中的分频计数器获得的。由开关电路时钟控制的开关电路的通断状态可以产生周期性的序列信号,以消除光电探测器输出的尖峰。为了对前置放大器后的闭环误差信号进行采样,设计了一种特殊的A/D采样时钟,可以避免开关电路的截止状态。由于A/D转换的时间花费了若干个个A/D时钟节拍,因此A/D转换后的调制闭环误差出现了延迟。因此,解调方波fd(t)被设计成具有对应A/D转换的延迟,以匹配调制的闭环误差的时序。A/D采样时钟通过计数器和逻辑乘积单元获得。通过分频器产生调制方波,并设计一个延迟单元以获得解调方波。数字反馈由积分器累加,并将结果加到调制方波中。
对给定的正常数α=0.1,可以应用LMI工具箱中的求解器feasp来判断该线性矩阵不等式的可行性,并在可行的情况下,利用得到的可行解得系统指数稳定的状态反馈控制律u(k)=-k2Kcx(k)。
已知光纤陀螺的参数为:干涉强度I0=2×10-2mW,探测器光电转换效率γ为0.91A/W,跨阻R为200K,探测器增益为γR=1.82×105,AD位数为12位,AD的基准电压为2伏,半周期采点个数Ndem为15,控制算法增益为1/1024,则由此可得光纤陀螺前向通道增益k1=4197.648。假设系统前向通道增益变化范围为±20%,所以不确定参数的变化范围||Δk1||≤839.5296。由于DA位数为16位,来计算反馈链的增益。
其中,γ和R分别是光电检测器的转换效率和跨阻,K2是前置放大器的增益,Nad和Vref分别是A/D转换器的分辨率和参考电压,Ndem是每个解调周期的采样数。
根据定理1环路增益优化方法,本发明优化前向通道增益的唯一可调变量是前置放大器的增益K2。通过理论分析,本发明可以知道前向通道增益k1是影响闭环误差检测精度的重要因素。幅频特性说明了输出信号和输入信号之间的误差。在实际应用中,要求光纤陀螺应实时反映角速度幅度而不会出现失真。因为高精度光纤陀螺比例因子的可重复性误差应在电源频率内低于0.1%。根据定理1的理论分析,还确定了前置放大器的增益K2=38.44。
进一步根据定理3,优化设计控制器的反馈矩阵Kc。
Kc=[-6.67 -1.22]。
第六步,实验验证效果。
静态精度测试主要测试常温时零偏的稳定性(也称为零漂)。零偏的稳定性是衡量光纤陀螺测试精度的主要指标。为了比较改进的方案对测试精度的提高,相同的测试环境下,将0.01deg/h光纤陀螺的原系统和优化系统分别进行了静态精度测试,如图3-4所示。由0.01deg/h光纤陀螺测试数据可见,经过优化以后系统的零偏稳定性的测试结果有了很大改善。
由以上试验数据可知,合理分配环路增益,可提高系统对噪声的抑制能力。试验数据表明,优化系统提高了系统的检测精度。以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明实施例揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (8)
1.一种光纤陀螺环路增益优化方法,其特征在于:包括如下步骤:
步骤一、设计检测方法从噪声背景中提取闭环误差信号;
步骤二、基于闭环误差模型推导光纤陀螺的动态方程;
步骤三、针对非线性条件,进行环路增益优化设计;
步骤四、考虑噪声以及环境扰动,进一步优化设计光纤陀螺系统闭环参数;
步骤五、基于上述优化结果进行检测电路与参数设计。
4.根据权利要求1所述的一种光纤陀螺环路增益优化方法,其特征在于:所述步骤三,在考虑光学效应非线性条件下,通过下述定理实现光纤陀螺系统指数稳定的高动态性能的环路增益优化设计:
5.根据权利要求4所述的一种光纤陀螺环路增益优化方法,其特征在于:所述步骤四,在考虑外界扰动与噪声的条件下,所建立的闭环光纤陀螺的动态方程为:
x(k+1)=Ax(k)+Bk1sin(-k2Kcx(k))+Bk1sin(-k2Kcx(k-1))+Bw(k)
其中,w(k)是平均值为零的宽带噪声。
7.根据权利要求1所述的一种光纤陀螺环路增益优化方法,其特征在于:设计的控制算法通过Lyapunov函数分析光纤陀螺闭环增益的分配原则,并以此为基础对光纤陀螺环路增益进行优化。
8.根据权利要求1所述的一种光纤陀螺环路增益优化方法,其特征在于:所述步骤五,实现检测电路与参数设计具体包括如下步骤:
设计检测电路并仿真波形进行闭环信号检测的硬件,包括:光电检测器,开关电路,前置放大器,A/D转换器,控制器,积分器,D/A转换器和D/A驱动电路;
振荡器用于生成时钟和参考信号,开关电路时钟,A/D采样时钟和D/A时钟分别为开关电路,为A/D转换器和D/A转换器提供时钟信号,并在调制解调过程中应用了调制方波和解调方波,开关时钟是通过峰值消除开关电路中的分频计数器获得的,由开关电路时钟控制的开关电路的通断状态产生周期性的序列信号,以消除光电探测器输出的尖峰;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110307328.8A CN113074712B (zh) | 2021-03-23 | 2021-03-23 | 一种光纤陀螺环路增益优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110307328.8A CN113074712B (zh) | 2021-03-23 | 2021-03-23 | 一种光纤陀螺环路增益优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113074712A true CN113074712A (zh) | 2021-07-06 |
CN113074712B CN113074712B (zh) | 2023-01-24 |
Family
ID=76613316
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110307328.8A Active CN113074712B (zh) | 2021-03-23 | 2021-03-23 | 一种光纤陀螺环路增益优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113074712B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115790566A (zh) * | 2023-01-31 | 2023-03-14 | 中国船舶集团有限公司第七〇七研究所 | 一种防止光纤陀螺闭环失效的调节方法及光纤陀螺 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6473182B1 (en) * | 2000-07-27 | 2002-10-29 | Northrup Grumman Corporation | Automatic gain control of a closed loop fiber optic gyroscope using non-linear control laws |
CN101458095A (zh) * | 2009-01-06 | 2009-06-17 | 北京航空航天大学 | 一种适用于全数字闭环光纤陀螺的前向增益估计系统 |
CN103217562A (zh) * | 2013-03-04 | 2013-07-24 | 北京航空航天大学 | 抑制光学电压互感器非线性特性的交直流电压检测方法及装置 |
CN107014367A (zh) * | 2017-04-10 | 2017-08-04 | 中国船舶重工集团公司第七0七研究所 | 一种优化光纤陀螺动态性能的方法 |
CN107796300A (zh) * | 2016-09-02 | 2018-03-13 | Abb瑞士股份有限公司 | 使用环路增益来确定干涉对比度的闭环干涉型传感器 |
CN109990773A (zh) * | 2019-04-11 | 2019-07-09 | 北京航空航天大学 | 一种干涉型光纤陀螺环路增益的检测与闭环控制系统及控制方法 |
-
2021
- 2021-03-23 CN CN202110307328.8A patent/CN113074712B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6473182B1 (en) * | 2000-07-27 | 2002-10-29 | Northrup Grumman Corporation | Automatic gain control of a closed loop fiber optic gyroscope using non-linear control laws |
CN101458095A (zh) * | 2009-01-06 | 2009-06-17 | 北京航空航天大学 | 一种适用于全数字闭环光纤陀螺的前向增益估计系统 |
CN103217562A (zh) * | 2013-03-04 | 2013-07-24 | 北京航空航天大学 | 抑制光学电压互感器非线性特性的交直流电压检测方法及装置 |
CN107796300A (zh) * | 2016-09-02 | 2018-03-13 | Abb瑞士股份有限公司 | 使用环路增益来确定干涉对比度的闭环干涉型传感器 |
CN107014367A (zh) * | 2017-04-10 | 2017-08-04 | 中国船舶重工集团公司第七0七研究所 | 一种优化光纤陀螺动态性能的方法 |
CN109990773A (zh) * | 2019-04-11 | 2019-07-09 | 北京航空航天大学 | 一种干涉型光纤陀螺环路增益的检测与闭环控制系统及控制方法 |
Non-Patent Citations (2)
Title |
---|
张勇等: "高精度数字闭环光纤陀螺的调制增益控制", 《哈尔滨工程大学学报》 * |
李敏等: "基于伪随机辨识的光纤陀螺前向增益估计", 《光电工程》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115790566A (zh) * | 2023-01-31 | 2023-03-14 | 中国船舶集团有限公司第七〇七研究所 | 一种防止光纤陀螺闭环失效的调节方法及光纤陀螺 |
Also Published As
Publication number | Publication date |
---|---|
CN113074712B (zh) | 2023-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111174773B (zh) | 一种多相调制解调的光纤陀螺多闭环算法 | |
CN107389097B (zh) | 光纤陀螺Sagnac光纤环本征频率跟踪测量方法 | |
CN109990773A (zh) | 一种干涉型光纤陀螺环路增益的检测与闭环控制系统及控制方法 | |
CN102520374B (zh) | 一种具有双调制特征的光纤磁场传感器数字闭环检测装置 | |
CN111999739A (zh) | 相位调制测距测速相干激光雷达方法及装置 | |
CN111220142B (zh) | 一种新型超高精度光纤陀螺仪的信号处理方法及装置 | |
US7990299B2 (en) | Optical sampling method for analog to digital conversion | |
WO2019116549A1 (ja) | 測距装置及び制御方法 | |
CN113074712B (zh) | 一种光纤陀螺环路增益优化方法 | |
CN104677398A (zh) | 一种基于频漂动态补偿的φ-otdr增敏方法及装置 | |
CN110375727A (zh) | 一种闭环光纤陀螺信号调制方法 | |
CN1851403A (zh) | 谐振式微光机电陀螺的闭环光频检测装置 | |
CN113720321A (zh) | 一种光纤陀螺温度补偿系统 | |
CN101975584B (zh) | 一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法 | |
CN100533063C (zh) | 光纤陀螺热噪声、散粒噪声及强度噪声分离的方法 | |
CN109029412B (zh) | 一种光纤陀螺闭环反馈回路非线性度测试方法及系统 | |
CN110986912B (zh) | 信号处理电路信号干扰检测方法、存储介质及计算机设备 | |
CN102901515A (zh) | 一种光纤陀螺渡越时间的在线快速测量方法 | |
JP2951547B2 (ja) | 光波距離計および距離測定方法 | |
CN116952212A (zh) | 一种大量程光纤陀螺仪及其测量方法 | |
CN115077510B (zh) | 一种基于交替振幅方波的干涉式光纤陀螺调制方法 | |
CN114325030B (zh) | 一种基于光学电流互感器的双方波调制开环解调方法 | |
CN110987014B (zh) | 光纤陀螺信号处理电路信号干扰检测方法、存储介质及计算机设备 | |
EP0919795B1 (en) | Method for determining the phase difference of light waves propagated over two paths | |
CN1945209A (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 |