CN111174773B - 一种多相调制解调的光纤陀螺多闭环算法 - Google Patents

一种多相调制解调的光纤陀螺多闭环算法 Download PDF

Info

Publication number
CN111174773B
CN111174773B CN201911405551.5A CN201911405551A CN111174773B CN 111174773 B CN111174773 B CN 111174773B CN 201911405551 A CN201911405551 A CN 201911405551A CN 111174773 B CN111174773 B CN 111174773B
Authority
CN
China
Prior art keywords
modulation
light intensity
fiber
closed
loop
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
CN201911405551.5A
Other languages
English (en)
Other versions
CN111174773A (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 CN201911405551.5A priority Critical patent/CN111174773B/zh
Publication of CN111174773A publication Critical patent/CN111174773A/zh
Priority to PCT/CN2020/103238 priority patent/WO2021135171A1/zh
Application granted granted Critical
Publication of CN111174773B publication Critical patent/CN111174773B/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

Abstract

本发明公开了一种多相调制解调的光纤陀螺多闭环算法。通过光纤陀螺的调制波形的合理设计,在调制解调周期内分为多个相位,在不同的相位上实现不同调制幅度,并对时间段的陀螺输出信号进行采样,结合不同的闭环处理方法分别实现角速率、半波电压参数、信号强度等参数的解调,实现多参数多闭环。本发明光纤陀螺多闭环算法提高光纤陀螺的线性度、偏置稳定性等,无需增加或改变现有硬件,通过调制解调方法以实时监控光纤陀螺的工作状态,提高可靠性。

Description

一种多相调制解调的光纤陀螺多闭环算法
技术领域
本发明涉及光纤陀螺信号处理方法中如何同时实现多个参数的调制解调的问题,尤其是涉及到一种多相调制解调的光纤陀螺多闭环算法。
背景技术
光纤陀螺是一种新型角速率测量仪,由于其具有全固态、带宽大及具有多种协议数字输出的优点,被广泛用于导航和姿态控制系统中。光纤陀螺的工作原理是基于光学萨格奈克效应的光纤干涉仪,即当环形干涉仪旋转时,产生一个正比于旋转角速率大小的萨格奈克相位差(也常称为萨格奈克相移或萨格奈克相位),通过检测该相位差,结合预先测量得到的标度因数可以推算得到环形干涉仪所在系统的角速度,而标度因数一般在光纤陀螺制作完毕之后通过转台标定而得,是一个常数,所以只需要检测得到萨格奈克相位差,则结合标度因数的大小则立刻可得到相应的角速度的大小。
萨格奈克相位差φSag和系统角速率Ω的关系如下:
Figure BDA0002348529180000011
其中,L为光纤陀螺光纤环光纤长度,D为光纤环直径,λ光纤陀螺所用光源波长,c为真空中光速。
光相位无法直接测量,需要通过干涉转换为光强测量,常通过对萨格奈克干涉仪添加一个±π/2的方波相位调制,以同时检测角速度的方向和大小,并达到最大的灵敏度,此时的光纤陀螺输出强度信号I应该表示为:
Figure BDA0002348529180000012
其中,I0为光源光信号强度,φsag为角速度对应的光相位。
高精度光纤陀螺中,为了实行超高精度的测量,采用的为闭环控制和测量的方式,并且随着精度的提高,不仅仅要求系统中实现一个闭环。在高精度光纤陀螺中,可以通过多个闭环回路,实现多闭环算法,达到不同的的测量目标。
在超高精度的闭环光纤陀螺中,需要对系统的半波电压参数、角速度、光强进行实时测量。为了满足高精度光纤陀螺的应用,需要对一系列的参数进行闭环调制解调或者进行解算,从而使得光纤陀螺工作在相位偏置点,使得采用的电光相位调制系数和真实波导参数一致,并工作在合理的光强稳定点上,这些算法需要在一个软件里面实现,并通过调制解调的方式融合在一起。现有技术缺少了这样的多闭环算法进行实时控制。
发明内容
鉴于现有高精度光纤陀螺对于多个参数具有闭环需求,同时要实现多个闭环算法和参数解算算法的现状,本发明的目的在于提供一种多相位调制的光纤陀螺多闭环算法,通过调制波形的多相位设计和对应调制解调算法的设计,同时完成多个光纤陀螺工作参数的解算,从而实现多闭环控制,从而最终为光纤陀螺的工作提供一个稳定的工作点,满足稳定性和高精度的需求。
本发明所采用的技术方案的步骤如下:
将光纤陀螺的一个调制解调周期分为多个相位,实现多相位调制,并对实时的陀螺输出信号进行采集;根据采样得到的陀螺输出信号按照预先设计的特定解调方式进行解调,从而得到被测的角速度、光强和半波电压参数等信息;然后通过这些信息,反馈调节光纤陀螺相位调制器调制的多相调制的幅度,从而实现光纤陀螺的角速度、光强和半波电压参数的多闭环控制。
所述的光纤陀螺中包含有相位调制器和光源。
所述多相位调制的一个调制解调周期包括六个调制相位,六个调制相位按照时间顺序顺次分别记为A、B、C、D、E、F:
A、D调制电压相同均为V1,V1=0.7Vπ,Vπ表示光纤陀螺所采用的相位调制器的半波电压参数,调制深度相同均为M1,调制持续时间相同,符号相反;
B、E调制电压相同均为V2,V2=0.75Vπ,调制深度相同均为M2,调制持续时间相同,符号相反;
C、F调制电压相同均为V3,V3=1.25Vπ,调制深度相同均为M3,调制持续时间相同,符号相反;
A、B、C三个调制相位持续的总时间和D、E、F三个调制相位持续的总时间相同,均等于光纤陀螺的渡越时间。
所述的光纤陀螺的渡越时间通过预先测量得到。调制电压和调制相位之间的比例系数由光纤陀螺所采用的相位调制器的半波电压参数Vπ所决定。
上述V1、V2、V3选择大小不同,均不等于π相位的整数倍,特别的选取 V1=0.7Vπ,V2=0.75Vπ,V3=1.25Vπ。其中半波电压参数Vπ初始值由光纤陀螺所采用的相位调制器的标称参数获得,或者通过预先测量获得,并通过以下半波电压参数Vπ闭环算法实时更新。
在多相位调制的一个调制解调周期内,A、B、C、D、E、F六个调制相位下的陀螺输出信号,分别记为SA、SB、SC、SD、SE、SF,对应的数字采样结果分别记为SA(n)、SB(n)、SC(n)、SD(n)、SE(n)、SF(n),其中n为第n个调制解调周期,n的起始值为1,按照半波电压参数闭环算法实时更新半波电压参数 Vπ提高光纤陀螺的测量性能,按照角速度闭环算法实时更新光纤陀螺的角速度提高输出准确性能,按照光强解算算法输出光强变量调制光源的光强输出进行光纤陀螺的内部补偿,这样共同能实现帮助光纤陀螺提高线性度、偏置稳定性。
采用半波电压参数闭环算法实时更新半波电压参数Vπ具体为:
(1)第n个调制解调周期内的采样结果分别为SA(n)、SB(n)、SC(n)、SD(n)、 SE(n)、SF(n),半波电压参数记为Vπ(n),计算半波电压参数中间变量 D(n)=[SB(n)+SE(n)-SC(n)-SF(n)];
(2)根据半波电压参数中间变量D(n)的大小,按照半波电压参数闭环幅度 dVπ更新下一个第n+1个调制解调周期的半波电压参数Vπ(n+1):
如果D(n)=0,则Vπ(n+1)=Vπ(n);
如果D(n)>0,则Vπ(n+1)=Vπ(n)+dVπ;
如果D(n)<0,则Vπ(n+1)=Vπ(n)-dVπ;
其中,半波电压参数闭环幅度dVπ为预设的参数,为正数,dVπ取值满足: dVπ<0.1,典型值可为dVπ=0.001;
(3)重复以上步骤,得到各个调制解调周期下的半波电压参数Vπ(n),作为半波电压参数的测量值,以半波电压参数的测量值为光纤陀螺内部解算的实际参数,进而提高光纤陀螺的测量性能。
按照角速闭环算法实时更新光纤陀螺的角速度,具体为:
(1)第n个调制解调周期内的采样结果分别为SA(n)、SB(n)、SC(n)、SD(n)、 SE(n)、SF(n),角速度记为R(n),R(n)的初始值R(1)取值为0,计算角速度中间变量RD(n):
RD(n)=[SA(n)+SB(n)+SF(n)-SC(n)-SD(n)-SE(n)];
(2)根据角速度中间变量RD(n)的大小,在n>1时,按照角速度闭环幅度 dR(n)更新角速度R(n):
如果RD(n)=0,则R(n)=R(n-1);
如果RD(n)>0,则R(n)=R(n-1)–dR(n);
如果RD(n)<0,则R(n)=R(n-1)+dR(n);
其中,角速度闭环幅度dR(n)为预设的参数,为正数,有两种取值方式:
A.固定取值:dR(n)=Const,Const为正值常数,典型值可取Const=10-6
B.固定比例:dR(n)=k×|RD(n)|,其中k为角速度比例系数,和闭环光纤陀螺的光强和电路放大倍数相关,通过光纤陀螺的光电参数预先标定和实验测试得到;
(3)在每一个调制解调周期中重复以上闭环过程,将实时的角速度R(n)进行输出,通过光纤陀螺的通信输出接口发送给外部应用系统,进而提高光纤陀螺的精度。
具体实施的应用系统例如为惯性导航系统。
按照光强解算算法输出光强变量调制光源的光强输出,如下:
(1)第n个调制解调周期内的采样结果分别为SA(n)、SB(n)、SC(n)、SD(n)、 SE(n)、SF(n),计算光强变量ID(n):
ID(n)=[SA(n)-SB(n)+SD(n)-SE(n)];
这样处理获得的光强变量ID(n)正比于于闭环光纤陀螺的光强、比例系数和电路放大倍数相关,在确定的一个光纤陀螺系统中,电路放大倍数不变的情况下,光强变量ID(n)的变化情况代表了光强的变化情况。
(2)将光强变量ID(n)作为光纤陀螺的光强监控值进行光强反馈,通过增加光强度调制器,将光强度调制器串联接在光纤陀螺内部光源的输出端,利用光强变量ID(n)控制光强度调制器的强度调制:
若光强变量ID(n)大于预先设置的光强信号阈值ID0,则通过控制光强度调制器增加光强损耗,降低光强;
若光强变量ID(n)小于预先设置的光强信号阈值ID0,则通过控制光强度调制器减少光强损耗,增大光强;
若光强变量ID(n)等于预先设置的光强信号阈值ID0,则通过控制光强度调制器保持光强不变;
从而实现光强的闭环反馈。
所述的半波电压参数闭环算法、角速度闭环算法和光强解算算法构成了多闭环算法,均采用一个调制解调周期内的六个相位对应的采样值,三者相互不干扰,在实现中采用FPGA作为数据处理器,实现多闭环算法的并行计算。
本发明原理情况如下:
作为导航系统的核心,光纤陀螺的精度决定着导航系统的导航精度,高精度陀螺的应用能够大幅度提高导航系统的性能,而同样的高精度陀螺的稳定性和可靠性也对应地影响着导航系统的稳定性和可靠性。通过闭环技术的引进,中低精度的光纤陀螺以及被广泛应用于各类系统中,而随着对陀螺精度和性能的提高,需要进一步对光纤陀螺的工作状态进行闭环控制,使其不仅仅在角速度对应的光学相位实现闭环,还需要对其他光强、电光相位调制系数等参数进行控制和测量,从而进一步提高其性能和稳定性。
本发明是通过设计合适的调制解调方式,具体的为多相调制解调的测量,将一个调制解调周期按时间顺序分为6个相位,并在不同的相位上采用不同的具有特定规律的调制电压,同时根据陀螺角速度的影响机理、电光相位调制系数的影响机理以及光强的影响机理的不同,通过在解调算法上的设计,同时实现这三者的分开独立解算,从而实现多闭环算法。所采用的采集不需增加硬件,具有可靠性好扩展性灵活的特点。
对于角速度产生的光相位,通过两个偏置相位大小一致,方向相反的相位调制,对应于这两个相位的采用信号差值反应了相位的大小;在开环陀螺中根据这个大小可以实现角速度的解算,在闭环陀螺中根据这个差值的大小可以实现相位残差的解算,从而调整闭环相位,在闭环良好的情况下根据闭环相位得到角速度的大小,这是传统光纤陀螺中方波调制解调的原理。
而根据两个相位和为2π的相位调制,如0.75π和1.25π,对应的调制电压分别为0.75Vπ和1.25Vπ,理论上在相位闭环下起对应的信号大小一致,满足 cos0.75π=cos1.25π。实际调制过程中,相位是通过电压提供的,实际调制的电压分别为0.75Vπ和1.25Vπ,如果此时Vπ和实际的半波电压(记为Vπ0)不准,则 0.75Vπ和1.25Vπ调制电压对应的调制相位偏离0.75π和1.25π,此时两个相位的对应的输出信号不等,根据这个不准可以实现半波电压偏离的估计,进而可以调整估计的半波电压值Vπ,直到此时两个相位对应输出的信号一致;在本发明中,在0.75π和1.25π之外增加了-0.75π和-1.25π的调制,结合四个信号的采样值实现半波电压估计值的结算,可以进一步降低对于转动相位等的敏感性,提高测量性能,改进测量稳定性。
同样的,根据两个相位具有一定偏差的相位调制,其信号采样的差值和光强相关,根据这个采样值可以实现光强的结算。例如通过0.75Vπ和0.70Vπ的相位偏置,其光强差值dI为:
dI=I0[1+cos(0.70π)]-I0[1+cos(0.75π)]=I0[cos(0.70π)-cos(0.75π)]
其中,dI表示光强差值,I0表示光源光信号强度,也简称为光强。
假设系统中的电路系统增益不变,则对应两个采样值正比于光强I0,根据采样值差可以反应光强的大小情况。通过不同调制解调周期下的测量,可以实现光强的实时监控和跟踪,在需要光强闭环的应用情况下,可以增加光强调制器,利用该信号实现强度的闭环反馈。
在本发明中,在0.75π和0.70π之外增加了-0.75π和-0.70π的调制,结合四个信号的采样值实现光强信号的结算,可以进一步降低对于转动相位等的敏感性,提高测量性能,改进测量稳定性。
本发明具有的有益效果是:
首次提出一种多相调制解调的光纤陀螺多闭环算法,该方法通过对调制解调周期的多相分解和设计设计,通过具有不同特征的六个调制相位的对陀螺信号进行相位调制,并针对不同的闭环和解算量设计对应的解调和闭环算法,同时实现了三个参数的闭环控制和解算,从而能够为高精度光纤陀螺的研制提供一个多闭环算法,提供一个稳定的工作点。
该方法不需增加额外器件,在相同的成本下提高了应用系统的可靠性,提高了光纤陀螺的稳定性,对于高精度导航应用系统具有重要的应用价值和推广价值。
附图说明
图1:不同调制幅度点对应的光强幅度图;
图2:多相调制解调波形图;
图3:半波电压闭环算法实现的原理图;
图4:角速度闭环算法实现的原理图;
图5:光强测量和解调的示意图图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
在光纤陀螺干涉仪的相位和光强对应曲线如图1,在光纤陀螺中,通过相位偏置的方式使得在零转速输入的情况下,干涉仪两束光之间依旧存在一定的相位差;在方波相位偏置(方波相位调制解调或方波调制解调)数字闭环光纤陀螺仪中,偏置点成对称选取,如图中的A和D,B和E,C和F,这三对偏置点相对 0相对对称,相位幅度大小相同,符号相反。
通过一对偏置点的光强的检测,可以实现角速度的解算(在闭环光纤陀螺中,实现角速度闭环残差的结算);在现有陀螺中,典型的方波调制只采用一对偏置点,也即只有一个相位调制相位。
在本发明中,将通过三对偏置点的选取,实现数字闭环光纤陀螺的角速度闭环、电光相位调制系数闭环和光强的结算等。
在本发明基础上,后续可以进一步扩展,通过其他偏置点的选取,实现更多和陀螺相关的参数的结算。
(1)、
半波电压参数闭环算法情况如图3所示。在数字闭环光纤陀螺中,我们通过调制电压,实现光相位的调制,实际输出控制的量为调制电压值。调制电压和光相位之间的比例关系极为电光相位调制系数,有时候又以半波电压Vπ来表示这个系数,半波电压Vπ指的是产生相位所需要施加的调制电压值。Vπ的值决定于所采用的电光相位调制器,其典型值为2V~5V,一般出厂时候均有对应的标称值。而实际该值受环境的温度,湿度等参数的影响,在实际过程中需要对其进行测试,或者实时标定。
在本发明中,通过图2中的调制波形,可以实现对于Vπ的实时测量和闭环,使得Vπ的测量值(Vπt)最终等于真实的Vπ值(记为Vπ0)。在Vπt=Vπ0的时候,如图1所示,此时:
SB=I0[1+cos(0.75Vπt*π/Vπ0)]=I0[1+cos(0.75Vπ0*π/Vπ0)]=I0[1+cos(0.75π)]
SC=I0[1+cos(1.25Vπt*π/Vπ0)]=I0[1+cos(1.25Vπ0*π/Vπ0)]=I0[1+cos(0.75π)]
SE=I0[1+cos(-0.75Vπt*π/Vπ0)]=I0[1+cos(0.75Vπ0*π/Vπ0)]=I0[1+cos(0.75π)]
SF=I0[1+cos(1.75Vπt*π/Vπ0)]=I0[1+cos(1.75Vπ0*π/Vπ0)]=I0[1+cos(0.75π)]
这四个采样值相同,满足SB=SC=SE=SF,而如果此时当前的Vπ的测量值 (Vπt)和真实Vπ值不等,存在误差k,也即Vπt=(1+k)Vπ0,则此时的采样值分别为:
在干涉相应图上表现为图4 ,表现为工作点此时均向0点运动。
SB=I0[1+cos(0.75Vπt*π/Vπ0)]=I0[1+cos(0.75π+0.75kπ)]
SC=I0[1+cos(1.25Vπt*π/Vπ0)]=I0[1+cos(0.75π-0.125kπ)]
SE=I0[1+cos(-0.75Vπt*π/Vπ0)]=I0[1+cos(0.75π+0.75kπ)]
SF=I0[1+cos(1.75Vπt*π/Vπ0)]=I0[1+cos(0.75π-0.125kπ)]
此时SB=SE≠SF=SC,记中间变量为:
D=(SB+SE)-(SF+SC)
D中包括了k的信息,从而可以根据k的大小,实时调整Vπt的大小,在 D=0的情况下,此时k=0,Vπt=Vπ0,从而实现半波电压的闭环。
在数字闭环光纤陀螺中,如果令n表示调制解调周期的序号,那么上面式子也可写为离散的表示方式,如下:
SB(n)=I0[1+cos(0.75π+0.75kπ)]
SC(n)=I0[1+cos(0.75π-0.125kπ)]
SE(n)=I0[1+cos(0.75π+0.75kπ)]
SF(n)=I0[1+cos(0.75π-0.125kπ)]
D(n)=(SB(n)+SE(n))-(SF(n)+SC(n));
根据D(n)的大小,按照等步长或者等比例的方式调整下一个周期的半波电压测量值Vπt(n)。
(2)、
多相调制解调波形如图2,图中画出了两个调制解调周期的波形,一个调制解调周期的周期长度为2τ,τ为光纤陀螺的渡越时间,可以预先测试得到,或者根据光纤陀螺的光纤长度和折射率计算得到。一个调制解脱周期内部顺序包括了6个调制相位A、B、C、D、E、F,并且这六个相位调制幅度有所不同,其中A和D,B和E,C和F两两之间分别构成对称调制相位对,每对两两之间持续的时间相同,并且调制电压幅度大小一样,分别记为V1、V2、V3,符号相反,虽然在应用中常让这三对对称调制相位持续时间相同,均为τ/3,但这不是必须要求,可以根据具体系统需求进行配置。同样A、B、C之间的顺序也非必须要求,可具体根据需求配置。
理论上,三个调制相位对的幅度可以任选,但实际中由于实现和稳定性等原因,V1、V2、V3根据所采用的集成电光相位调制器(简称电光相位调制器,相位调制器或者Y波导等)的电光相位调制系数的电压选取,一般选取确定比例的值。作为一个实用的例子,假设光纤陀螺的电光相位调制系数标称值已知,其半波电压记为Vπ,可选比例系数分别为0.7,0.75,1.25,也即V1=0.7Vπ, V2=0.75Vπ,V2=1.25Vπ。
对于每一对对称调制相位,均可按照典型方波调制的光纤陀螺的信息处理方法,通过对相反调制相位下对应的采样值实现角速度的结算(开环陀螺)或者角速度残差的结算(闭环陀螺)。在闭环陀螺中,通过使得每一对对称调制相位两个相反调制相位下的采样差为零实现闭环控制,在实现闭环控制时候,对应的反馈相位调制则正比于角速度,此即为典型的方波调制解调算法。
角速度闭环算法情况如图4所示,对于闭环光纤陀螺,如果此时反馈相位和实际角速度对应的相位之间不等,产生了相位残差,则此时可以根据A、B、 C、E、D、F的采样值获得。如图表示的是在有相位残差的情况下,此时A、B、 C、E、D、F六个工作点的净相位(偏置相位+反馈相位+角速度信号)均向右偏离,此时A、B、C、E、D、F六个点3对对称调制相位的采样结果SA、SB、SC、 SE、SD、SF不同,根据他们的采样值可以得到相位残差的大小,进而可以调制下一个周期的闭环反馈相位,从而最终重新获得闭环的结果,在重新获得闭环的情况下,角速度对应的光相位等于反馈相位,同时三对对称调制的采样结果相同,也即SA=SD,SB=SE,SC=SF。
在数字化光纤陀螺中,1个调制解调周期为一个计算周期,在数字化下可以将第n个调制解调周期的采样值具体记为SA(n)、SB(n)、SC(n)、SE(n)、SD(n)、 SF(n),对应的中间变量和计算结果也均可增加后缀(n)。
(2)、
光强测量解调如图5所示,通过具有不同调制深度的相位A、B、D和E的采样结果,实现光强的闭环测量。这四个调制相位上对应的光强度信号大小分别记为SA、SB、SD、SE,则:
SA=I0[1+cos(0.70π)]
SB=I0[1+cos(0.75π)]
SD=I0[1+cos(0.70π)]
SE=I0[1+cos(0.75π)]
记中间变量为:
ID=(SA+SD)-(SB+SE)
ID中包括了I0的信息,和I0成正比,从而可以根据ID的大小,实时获得光强I0的大小。
在数字闭环光纤陀螺中,如果令n表示调制解调周期的序号,那么上面式子也可写为离散的表示方式,如下:
SA(n)=I0(n)[1+cos(0.70π)]
SB(n)=I0(n)[1+cos(0.75π)]
SD(n)=I0(n)[1+cos(0.70π)]
SE(n)=I0(n)[1+cos(0.75π)]
中间变量为:
ID(n)=[SA(n)+SD(n)]-[SB(n)+SE(n)]=2[cos(0.70π)-cos(0.75π)]I0(n)
从而光强的实时解调值为:
I0(n)=ID(n)*2[cos(0.70π)-cos(0.75π)]
可见光强的实时值正比于ID(n),考虑到实际光纤陀螺系统中的光电转换过程也为线性过程,实际采样值为光强ID(n)乘以一个和电路成正比的系数,比例系数和电路放大倍数相关,在确定的一个光纤陀螺系统中,电路放大倍数不变的情况下,ID(n)的变化情况代表了光强的变化情况,ID(n)解算结果,作为光纤陀螺光强监控值;根据闭环光纤陀螺的设计功能需求,在需要进行光强反馈的系统中,通过增加光强度调制器,利用ID(n)控制光强度调制器的强度调制,从而实现光强的闭环反馈;在不需要光强反馈的系统,可以通过ID(n)实时反应系统光强信号的大小。
由此实施可见,本发明能够为光纤陀螺提供一个稳定的工作点,提高光纤陀螺的线性度、偏置稳定性等;不需增加或改变现有硬件,通过调制解调方法以实时监控光纤陀螺的工作状态,从而提高其应用系统的可靠性,具有很强的工业应用价值。

Claims (6)

1.一种多相调制解调的光纤陀螺多闭环方法,所述的光纤陀螺中包含有相位调制器和光源,其特征在于:将光纤陀螺的一个调制解调周期分为多个相位,实现多相位调制,并对实时的陀螺输出信号进行采集;根据采样得到的陀螺输出信号按照预先设计的解调方式进行解调,从而得到被测的角速度、光强和半波电压参数信息;然后通过这些信息,反馈调节光纤陀螺相位调制器调制的多相调制的幅度,从而实现光纤陀螺的角速度、光强和半波电压参数的多闭环控制;
所述多相位调制的一个调制解调周期包括六个调制相位,六个调制相位按照时间顺序顺次分别记为A、B、C、D、E、F:
A、D调制电压相同均为V1,V1=0.7Vπ,Vπ表示光纤陀螺所采用的相位调制器的半波电压参数,调制深度相同均为M1,调制持续时间相同,符号相反;
B、E调制电压相同均为V2,V2=0.75Vπ,调制深度相同均为M2,调制持续时间相同,符号相反;
C、F调制电压相同均为V3,V3=1.25Vπ,调制深度相同均为M3,调制持续时间相同,符号相反;
A、B、C三个调制相位持续的总时间和D、E、F三个调制相位持续的总时间相同,均等于光纤陀螺的渡越时间。
2.根据权利要求1所述的一种多相调制解调的光纤陀螺多闭环方法,其特征在于:在多相位调制的一个调制解调周期内,A、B、C、D、E、F六个调制相位下的陀螺输出信号,分别记为SA、SB、SC、SD、SE、SF,对应的数字采样结果分别记为SA(n)、SB(n)、SC(n)、SD(n)、SE(n)、SF(n),其中n为第n个调制解调周期,按照半波电压参数闭环算法实时更新半波电压参数Vπ提高光纤陀螺的测量性能,按照角速度闭环算法实时更新光纤陀螺的角速度提高输出准确性能,按照光强解算算法输出光强变量调制光源的光强输出进行光纤陀螺的内部补偿,这样共同能实现帮助光纤陀螺提高线性度、偏置稳定性。
3.根据权利要求2所述的一种多相调制解调的光纤陀螺多闭环方法,其特征在于:采用半波电压参数闭环算法实时更新半波电压参数Vπ具体为:
(1)第n个调制解调周期内的采样结果分别为SA(n)、SB(n)、SC(n)、SD(n)、SE(n)、SF(n),半波电压参数记为Vπ(n),计算半波电压参数中间变量D(n)=[SB(n)+SE(n)-SC(n)-SF(n)];
(2)根据半波电压参数中间变量D(n)的大小,按照半波电压参数闭环幅度dVπ更新下一个第n+1个调制解调周期的半波电压参数Vπ(n+1):
如果D(n)=0,则Vπ(n+1)=Vπ(n);
如果D(n)>0,则Vπ(n+1)=Vπ(n)+dVπ;
如果D(n)<0,则Vπ(n+1)=Vπ(n)-dVπ;
其中,半波电压参数闭环幅度dVπ为预设的参数,为正数,dVπ取值满足:dVπ<0.1;
(3)重复以上步骤,得到各个调制解调周期下的半波电压参数Vπ(n),作为半波电压参数的测量值,以半波电压参数的测量值为光纤陀螺内部解算的实际参数,进而提高光纤陀螺的测量性能。
4.根据权利要求2所述的一种多相调制解调的光纤陀螺多闭环方法,其特征在于:按照角速闭环算法实时更新光纤陀螺的角速度,具体为:
(1)第n个调制解调周期内的采样结果分别为SA(n)、SB(n)、SC(n)、SD(n)、SE(n)、SF(n),角速度记为R(n),R(n)的初始值R(1)取值为0,计算角速度中间变量RD(n):
RD(n)=[SA(n)+SB(n)+SF(n)-SC(n)-SD(n)-SE(n)];
(2)根据角速度中间变量RD(n)的大小,在n>1时,按照角速度闭环幅度dR(n)更新角速度R(n):
如果RD(n)=0,则R(n)=R(n-1);
如果RD(n)>0,则R(n)=R(n-1)–dR(n);
如果RD(n)<0,则R(n)=R(n-1)+dR(n);
其中,角速度闭环幅度dR(n)为预设的参数,为正数,有两种取值方式:
A.固定取值:dR(n)=Const,Const为正值常数;
B.固定比例:dR(n)=k×|RD(n)|,其中k为角速度比例系数;
(3)在每一个调制解调周期中重复以上闭环过程,将实时的角速度R(n)进行输出,通过光纤陀螺的通信输出接口发送给外部应用系统,进而提高光纤陀螺的精度。
5.根据权利要求2所述的一种多相调制解调的光纤陀螺多闭环方法,其特征在于:按照光强解算算法输出光强变量调制光源的光强输出,如下:
(1)第n个调制解调周期内的采样结果分别为SA(n)、SB(n)、SC(n)、SD(n)、SE(n)、SF(n),计算光强变量ID(n):
ID(n)=[SA(n)-SB(n)+SD(n)-SE(n)];
(2)将光强变量ID(n)作为光纤陀螺的光强监控值进行光强反馈,通过增加光强度调制器,将光强度调制器串联接在光纤陀螺内部光源的输出端,利用光强变量ID(n)控制光强度调制器的强度调制:
若光强变量ID(n)大于预先设置的光强信号阈值ID0,则通过控制光强度调制器增加光强损耗,降低光强;
若光强变量ID(n)小于预先设置的光强信号阈值ID0,则通过控制光强度调制器减少光强损耗,增大光强;
若光强变量ID(n)等于预先设置的光强信号阈值ID0,则通过控制光强度调制器保持光强不变;从而实现光强的闭环反馈。
6.根据权利要求2所述的一种多相调制解调的光纤陀螺多闭环方法,其特征在于:所述的半波电压参数闭环算法、角速度闭环算法和光强解算算法均采用一个调制解调周期内的六个相位对应的采样值,三者相互不干扰,在实现中采用FPGA作为数据处理器,实现多闭环方法的并行计算。
CN201911405551.5A 2019-12-30 2019-12-30 一种多相调制解调的光纤陀螺多闭环算法 Active CN111174773B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201911405551.5A CN111174773B (zh) 2019-12-30 2019-12-30 一种多相调制解调的光纤陀螺多闭环算法
PCT/CN2020/103238 WO2021135171A1 (zh) 2019-12-30 2020-07-21 一种多相调制解调的光纤陀螺多闭环方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911405551.5A CN111174773B (zh) 2019-12-30 2019-12-30 一种多相调制解调的光纤陀螺多闭环算法

Publications (2)

Publication Number Publication Date
CN111174773A CN111174773A (zh) 2020-05-19
CN111174773B true CN111174773B (zh) 2022-04-05

Family

ID=70655881

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911405551.5A Active CN111174773B (zh) 2019-12-30 2019-12-30 一种多相调制解调的光纤陀螺多闭环算法

Country Status (2)

Country Link
CN (1) CN111174773B (zh)
WO (1) WO2021135171A1 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111174773B (zh) * 2019-12-30 2022-04-05 浙江大学 一种多相调制解调的光纤陀螺多闭环算法
CN111998843B (zh) * 2020-08-10 2022-10-21 中国科学院光电技术研究所 一种低噪声数字双闭环光纤陀螺信号调制解调方法
CN112291172B (zh) * 2020-12-24 2021-04-09 武汉慧联无限科技有限公司 一种网关及信息解调方法
CN113720357A (zh) * 2021-09-16 2021-11-30 北京控制工程研究所 一种3s光纤imu真空全温条件下陀螺标度因数标定及补偿方法
CN114295115B (zh) * 2022-01-07 2023-03-31 北京思卓博瑞科技有限公司 一种提高光纤陀螺动态范围的方法及装置
CN114719837B (zh) * 2022-05-27 2022-09-16 北京邮电大学 用于数字闭环光纤陀螺仪的六状态调制解调方法和系统
CN115077567B (zh) * 2022-08-18 2022-11-04 中国船舶重工集团公司第七0七研究所 一种基于波导复位误差的标度因数补偿系统及方法
CN116026303A (zh) * 2023-03-29 2023-04-28 航天星云(北京)科技有限公司 光纤陀螺电路板、调试方法、光纤陀螺及惯性导航系统
CN116045946B (zh) * 2023-03-31 2023-06-09 中国船舶集团有限公司第七〇七研究所 光纤陀螺及降低其标度因数非线性度的方法
CN116045957B (zh) * 2023-03-31 2023-06-16 中国船舶集团有限公司第七〇七研究所 一种基于光纤陀螺扩频采样的误差消除方法

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60157012A (ja) * 1984-01-27 1985-08-17 Fujitsu Ltd 閉ル−プ光フアイバジヤイロシステム
GB2179140B (en) * 1985-08-14 1989-08-16 Stc Plc Sagnac effect device
FR2654827B1 (fr) * 1989-11-17 1992-03-20 Photonetics Dispositif de mesure a fibre optique, gyrometre, centrale de navigation et de stabilisation.
US5052808A (en) * 1990-02-15 1991-10-01 Litton Systems, Inc. Method and apparatus for interferometric rotation sensor phase modulation, intensity demodulation, and control
CN101701819B (zh) * 2009-11-23 2011-05-11 浙江大学 一种双轴复用光纤陀螺及其信号调制解调方法
CN101706278B (zh) * 2009-11-23 2011-05-11 浙江大学 一种可监控光功率大小的调制解调方法
CN101881669B (zh) * 2010-06-22 2011-11-09 浙江大学 集成电光相位调制器半波电压的实时闭环测量和追踪方法
WO2012008955A1 (en) * 2010-07-14 2012-01-19 Ralph Bergh Occasional calibration phase-difference modulation for sagnac interferometer
CN104089614A (zh) * 2014-05-21 2014-10-08 浙江大学 一种光纤陀螺中光强和相位调制项的实时同步测量方法
CN105973219B (zh) * 2016-04-29 2018-08-17 重庆华渝电气集团有限公司 光纤陀螺调制解调方法
IT201700091341A1 (it) * 2017-08-07 2017-11-07 Civitanavi Systems S R L Schema di modulazione della fase ottica di un MIOC di un giroscopio in fibra ottica di tipo interferometrico.
CN108132062B (zh) * 2017-12-05 2020-09-18 北京控制工程研究所 基于四分频四态调制的光纤陀螺光路功率检测方法及装置
CN108519079B (zh) * 2018-03-28 2021-11-26 北京航空航天大学 一种双闭环光纤陀螺六态调制降低串扰的方法
CN208350099U (zh) * 2018-04-19 2019-01-08 中国航空工业集团公司西安飞行自动控制研究所 基于六态调制的闭环光纤陀螺光功率自检测装置
CN110006417B (zh) * 2019-03-15 2021-04-20 湖北三江航天红峰控制有限公司 一种数字闭环光纤陀螺随机四态调制方法
CN109990773B (zh) * 2019-04-11 2020-01-07 北京航空航天大学 一种干涉型光纤陀螺环路增益的检测与闭环控制系统及控制方法
CN111174773B (zh) * 2019-12-30 2022-04-05 浙江大学 一种多相调制解调的光纤陀螺多闭环算法

Also Published As

Publication number Publication date
WO2021135171A1 (zh) 2021-07-08
WO2021135171A9 (zh) 2021-09-02
CN111174773A (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
CN111174773B (zh) 一种多相调制解调的光纤陀螺多闭环算法
CN109990773B (zh) 一种干涉型光纤陀螺环路增益的检测与闭环控制系统及控制方法
JPH03210417A (ja) ファイバー光測定装置、ジャイロメータ、セントラルナビゲーション、及び安定化システム
CN103411601B (zh) 一种实现光路差分的双干涉式光纤陀螺仪调制解调方法
CN1118684C (zh) 开环光纤陀螺信号检测方法及装置
US20130044328A1 (en) Systems and methods for environmentally insensitive high-performance fiber-optic gyroscopes
CN100541127C (zh) 采用不对称方波调制方法测量干涉式光纤陀螺本征频率的方法
CN103105177B (zh) 低成本光纤陀螺本征频率测量方法
CN106871931B (zh) 一种闭环光纤陀螺温度补偿方法
CN115143948A (zh) 一种基于光纤陀螺本征频率实时补偿标度因数的方法
CN115077567A (zh) 一种基于波导复位误差的标度因数补偿系统及方法
CN114993281A (zh) 基于低相干光的全数字化闭环角速度检测系统及方法
CN113720321A (zh) 一种光纤陀螺温度补偿系统
CN115752416A (zh) 用于检测角速度信号的谐振式光纤陀螺仪
CN101101213B (zh) 一种补偿数字闭环光纤陀螺标度因数非线性度的方法
CN103940593B (zh) 一种基于数字锯齿波的光纤陀螺渡越时间快速测量方法
CN100561128C (zh) 光纤陀螺渡越时间在线精密测量方法
CN114111754B (zh) 消除复位误差的光纤陀螺闭环控制方法、系统及存储介质
CN105674976A (zh) 光纤陀螺调制解调方法、提高标度因数稳定性方法及装置
CN104296740A (zh) 一种基于反熔丝fpga的光纤陀螺主控板晶振选取方法
CN113074713B (zh) 一种光纤陀螺高动态控制方法
CN114705176A (zh) 一种方波调制光纤陀螺全温启动2π电压参数调整方法
CN104613955A (zh) 一种兼顾大角速率和高精度的光纤陀螺
CN104713538A (zh) 一种闭环光纤陀螺最小拍控制方法
CN102944232B (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