張曉雅，楊志懷，宋麗薇，馬 林

（天津航海儀器研究所，天津 300131）

連續(xù)旋轉(zhuǎn)的光纖陀螺全溫標(biāo)度因數(shù)快速建模補(bǔ)償方法

張曉雅，楊志懷，宋麗薇，馬 林

（天津航海儀器研究所，天津 300131）

在全溫范圍內(nèi)應(yīng)用的光纖陀螺，標(biāo)度因數(shù)誤差是其主要的誤差之一。特別是在大角速率或者高精度應(yīng)用時(shí)，光纖陀螺的標(biāo)度因數(shù)誤差甚至超過(guò)零偏漂移誤差。在實(shí)際使用中，需對(duì)陀螺標(biāo)度因數(shù)在全溫范圍內(nèi)進(jìn)行建模和補(bǔ)償。對(duì)光纖陀螺標(biāo)度因數(shù)誤差機(jī)理進(jìn)行詳細(xì)分析后，提出了一種連續(xù)旋轉(zhuǎn)的光纖陀螺全溫標(biāo)度因數(shù)快速建模補(bǔ)償方法。基于單軸速率轉(zhuǎn)臺(tái)的連續(xù)旋轉(zhuǎn)，可以自動(dòng)快速完成標(biāo)度因數(shù)全溫建模且工程實(shí)現(xiàn)簡(jiǎn)單易行。更重要的是該方法可以有效識(shí)別標(biāo)度因數(shù)在全溫范圍內(nèi)的變化拐點(diǎn)，提高建模和補(bǔ)償?shù)木?。?duì)比試驗(yàn)結(jié)果表明，采用此方法后能精確測(cè)得某型光纖陀螺全溫工作的標(biāo)度因數(shù)真實(shí)拐點(diǎn)為48℃，全溫標(biāo)度因數(shù)補(bǔ)償精度優(yōu)于15×10-6，較按照GJB2426-2004進(jìn)行的多點(diǎn)測(cè)試后補(bǔ)償提高10%左右。

光纖陀螺；標(biāo)度因數(shù)；溫度建模；補(bǔ)償

當(dāng)溫度變化時(shí)光源、Y波導(dǎo)集成光學(xué)調(diào)制器等光學(xué)元器件的平均波長(zhǎng)會(huì)隨之變化，導(dǎo)致光纖陀螺的標(biāo)度因數(shù)變化[1-2]。在全溫范圍內(nèi)應(yīng)用的光纖陀螺，特別是在大角速率或者高精度應(yīng)用時(shí)，光纖陀螺的標(biāo)度因數(shù)誤差會(huì)超過(guò)偏置漂移誤差，嚴(yán)重制約光纖陀螺的環(huán)境適用性。從機(jī)理上消除溫度帶來(lái)的標(biāo)度因數(shù)誤差難度大、成本高，而基于理論機(jī)理的溫度誤差建模和補(bǔ)償是一種簡(jiǎn)單有效的解決方法。特別是在全溫范圍內(nèi)（-40℃～+60℃）對(duì)光纖陀螺的標(biāo)度因數(shù)誤差進(jìn)行建模和補(bǔ)償，對(duì)于提高標(biāo)度因數(shù)穩(wěn)定性指標(biāo)有重要意義。

常用的建模方法有線性回歸模型[3-4]、多種神經(jīng)網(wǎng)絡(luò)模型[5]等。前者結(jié)構(gòu)簡(jiǎn)單，工程實(shí)現(xiàn)容易，缺點(diǎn)是非線性補(bǔ)償能力弱；后者結(jié)構(gòu)復(fù)雜，參數(shù)多，將物理過(guò)程看作一個(gè)黑箱模型，非線性函數(shù)逼近效果能力強(qiáng)，但工程實(shí)現(xiàn)困難。

常用的光纖陀螺標(biāo)度因數(shù)測(cè)試方法是按照GJB2426-2004進(jìn)行一系列的角速率測(cè)試，根據(jù)最小二乘法計(jì)算得到[6]。全溫條件下的光纖陀螺標(biāo)度因數(shù)建模通常是選擇幾個(gè)典型的溫度點(diǎn)，分別進(jìn)行標(biāo)度因數(shù)測(cè)試，得到各溫度點(diǎn)對(duì)應(yīng)的標(biāo)度因數(shù)值作為標(biāo)度因數(shù)補(bǔ)償模型的輸入[3-4]。上述這種光纖陀螺標(biāo)度因數(shù)測(cè)試和建模方法由于頻繁操作速率轉(zhuǎn)臺(tái)效率低，且隨著建模精度需求增加更多的測(cè)試溫度點(diǎn)時(shí)，會(huì)大大增加測(cè)試時(shí)間和成本。本文提出了一種基于連續(xù)旋轉(zhuǎn)的光纖陀螺全溫標(biāo)度因數(shù)快速建模補(bǔ)償方法。基于帶溫箱的單軸速率轉(zhuǎn)臺(tái)，在全溫條件下使其單向勻速旋轉(zhuǎn)即可自動(dòng)快速完成標(biāo)度因數(shù)在全溫條件下的建模，工程實(shí)現(xiàn)簡(jiǎn)單易行。更重要的是該方法可以有效識(shí)別標(biāo)度因數(shù)在全溫范圍內(nèi)的變化拐點(diǎn)，無(wú)盲點(diǎn)地全面反映標(biāo)度因數(shù)隨溫度變化的情況，從而提高建模和補(bǔ)償?shù)木取?/p>

1 光纖陀螺標(biāo)度因數(shù)誤差機(jī)理

目前中高精度光纖陀螺均采用全數(shù)字閉環(huán)光纖陀螺方案，其中光纖陀螺敏感到的輸出信號(hào)與輸入角速率的關(guān)系可以表示為[3]

式中：ΩI為輸入角速度，Ω0為輸出角速度，L為光纖長(zhǎng)度，D為光纖環(huán)直徑，λ為真空中的光波波長(zhǎng)，C為真空中的光速，Kfp為Y波導(dǎo)的調(diào)制系數(shù)，Vpp為階梯波峰值電壓，N為D/A轉(zhuǎn)換器位數(shù)，KSF為光纖陀螺的標(biāo)度因數(shù)。由式(1)可以看出，光纖陀螺的標(biāo)度因數(shù)受到光纖敏感環(huán)圈長(zhǎng)度和直徑、光波波長(zhǎng)、Y波導(dǎo)半波電壓以及電路階梯波峰值電壓的影響，而這些因素又均受到溫度的影響。

全數(shù)字閉環(huán)光纖陀螺方案的優(yōu)點(diǎn)是通過(guò)增加第二反饋回路可以精確跟蹤 Y波導(dǎo)半波電壓隨溫度的變化，即始終保持調(diào)制系數(shù)與Y波導(dǎo)半波電壓的乘積系數(shù)為 2π，從而提高標(biāo)度因數(shù)的穩(wěn)定性[7-8]。因此在理想閉環(huán)反饋狀態(tài)下，標(biāo)度因數(shù)誤差僅由L、D的乘積和光波長(zhǎng)λ決定。

光纖陀螺的標(biāo)度因數(shù)溫度誤差可以表示為[9]

中高精度光纖陀螺通常采用超熒光摻鉺光纖光源(ASE)，其平均波長(zhǎng)漂移有三種與溫度相關(guān)的來(lái)源，可以表示為[10]式中：第一項(xiàng)是有摻鉺光纖的固有溫度系數(shù)導(dǎo)致的固有的平均波長(zhǎng)的變化，第二項(xiàng)和第三項(xiàng)分別表示泵浦波長(zhǎng)和泵浦功率隨溫度變化導(dǎo)致的光源平均波長(zhǎng)變化。通過(guò)優(yōu)化設(shè)計(jì)，可以將ASE光源的溫度系數(shù)做到小于5×10-7/℃。

通過(guò)以上對(duì)標(biāo)度因數(shù)誤差來(lái)源的分析并結(jié)合大量光纖陀螺實(shí)驗(yàn)測(cè)試表明，其標(biāo)度因數(shù)隨溫度變化具有較強(qiáng)的重復(fù)性和規(guī)律性，且通過(guò)溫度的一階最小二乘法進(jìn)行建模和補(bǔ)償即可實(shí)現(xiàn)比較好的全溫標(biāo)度因數(shù)性能。

2 標(biāo)度因數(shù)溫度誤差建模和補(bǔ)償

標(biāo)度因數(shù)的溫度誤差模型可以表示為

式中：K0是與溫度無(wú)關(guān)的標(biāo)度因數(shù)，KT是與溫度相關(guān)的陀螺標(biāo)度因數(shù)誤差系數(shù)，T為陀螺溫度。根據(jù)式(4)模型，光纖陀螺全溫標(biāo)度因數(shù)通常的建模方法是選擇幾個(gè)典型的溫度點(diǎn)，分別按照GJB2426-2004進(jìn)行標(biāo)度因數(shù)測(cè)試，得到的結(jié)果作為標(biāo)度因數(shù)補(bǔ)償模型的輸入，進(jìn)而通過(guò)分段建模的方式進(jìn)行補(bǔ)償[3-4]。

光纖陀螺的輸出可以表示為

式中：0Ω為光纖陀螺的固有零偏，ω為陀螺輸入轉(zhuǎn)速。將光纖陀螺的敏感軸朝天且固定于一勻速單向旋轉(zhuǎn)的單軸速率轉(zhuǎn)臺(tái)，此時(shí)由零偏引起的輸出標(biāo)度因數(shù)測(cè)試誤差可以表示為

圖1 不同轉(zhuǎn)速下陀螺零偏對(duì)標(biāo)度因數(shù)誤差影響Fig.1 Gyro bias’s influence on scale factor error at different rotating speeds

根據(jù)式(6)，圖1給出了不同轉(zhuǎn)速下，陀螺固有零偏對(duì)標(biāo)度因數(shù)測(cè)試的影響曲線。通常中高級(jí)精度光纖陀螺的固有零偏在0.1 (°)/h以下，因此當(dāng)轉(zhuǎn)臺(tái)速率大于50 (°)/s時(shí)，其對(duì)標(biāo)度因數(shù)的誤差影響小于0.6×10-6。

當(dāng)光纖陀螺在以較高速率勻速轉(zhuǎn)動(dòng)時(shí)，陀螺輸出的變化即代表了溫度引起的標(biāo)度因數(shù)變化。設(shè)置溫箱溫度激勵(lì)，讓光纖陀螺遍歷全溫范圍即可完成光纖陀螺的全溫標(biāo)度因數(shù)連續(xù)建模。這種方法可以有效識(shí)別標(biāo)度因數(shù)在全溫范圍內(nèi)的變化拐點(diǎn)，無(wú)盲點(diǎn)地全面反映標(biāo)度因數(shù)隨溫度變化情況，提高建模精度。選定溫度拐點(diǎn)后，采用最小二乘法可以實(shí)現(xiàn)光纖陀螺分段式標(biāo)度因數(shù)建模。

3 試驗(yàn)驗(yàn)證

圖2(a)給出了一只光纖陀螺采用連續(xù)旋轉(zhuǎn)速率法時(shí)的溫度激勵(lì)曲線，圖2(b)給出了在此溫度激勵(lì)下的光纖陀螺輸出曲線，其代表了不同溫度點(diǎn)下的光纖陀螺標(biāo)度因數(shù)變化。圖3給出了同一只光纖陀螺采用常規(guī)國(guó)軍標(biāo)方法進(jìn)行的20℃一個(gè)采樣點(diǎn)的全溫標(biāo)度因數(shù)建模。

從圖2～3中可以看出，連續(xù)旋轉(zhuǎn)法能全面反映各個(gè)溫度下光纖陀螺標(biāo)度情況，有利于標(biāo)度因數(shù)分段補(bǔ)償時(shí)選擇合適的拐點(diǎn)。傳統(tǒng)國(guó)軍標(biāo)方法在測(cè)試點(diǎn)較少情況下不能完全精確反映光纖陀螺的標(biāo)度因數(shù)在全溫條件下的變化情況。按照?qǐng)D3所示的模型建模，當(dāng)補(bǔ)償后的光纖陀螺工作在真實(shí)拐點(diǎn)48℃時(shí)，將發(fā)生較大的標(biāo)度因數(shù)誤差。為避免上述問題，可以在建模時(shí)增加更多的采樣點(diǎn)，與此同時(shí)又會(huì)大大增加測(cè)試時(shí)間和成本，且不能精確反映全溫條件下光纖陀螺標(biāo)度因數(shù)曲線的精確拐點(diǎn)位置。

通過(guò)上述方法對(duì)3只光纖陀螺分別進(jìn)行了全溫標(biāo)度因數(shù)建模，并采用單片F(xiàn)PGA實(shí)現(xiàn)光纖陀螺分段式全溫標(biāo)度因數(shù)補(bǔ)償。最后，采用國(guó)軍標(biāo)測(cè)試方法對(duì)其分別進(jìn)行了多點(diǎn)測(cè)試驗(yàn)證。表1為3支陀螺分別用兩種方法建模后的補(bǔ)償效果。

由于連續(xù)旋轉(zhuǎn)法建模方法更能捕捉精確的各段標(biāo)度因數(shù)拐點(diǎn)，上述實(shí)驗(yàn)中3只陀螺補(bǔ)償后的全溫標(biāo)度因數(shù)重復(fù)性相比多點(diǎn)測(cè)試法建模，補(bǔ)償精度更高，標(biāo)度因數(shù)重復(fù)性指標(biāo)更優(yōu)。

圖2 (a) 溫度激勵(lì)曲線Fig.2(a) Temperature curve

圖2 (b) 連續(xù)旋轉(zhuǎn)法標(biāo)度建模Fig.2(b) Scale factor modeling by continuous rotation method

圖3 常規(guī)多點(diǎn)測(cè)試法標(biāo)度因數(shù)建模Fig.3 Scale factor modeling by conventional mutipoint method

表1 兩種方法建模后全溫標(biāo)度因數(shù)重復(fù)性對(duì)比Tab.1 Compensation results of two modeling methods

4 結(jié) 論

光纖陀螺標(biāo)度因數(shù)溫度誤差主要受到L、D的乘積和光源波長(zhǎng)λ隨溫度變化的影響。陀螺以較大速率勻速轉(zhuǎn)動(dòng)時(shí)，陀螺的輸出變化直接反應(yīng)標(biāo)度因數(shù)隨溫度的變化關(guān)系?；趩屋S速率轉(zhuǎn)臺(tái)，采用連續(xù)旋轉(zhuǎn)法對(duì)光纖陀螺進(jìn)行全溫標(biāo)度因數(shù)，可以在大大提高建模效率的同時(shí)實(shí)現(xiàn)全溫范圍內(nèi)的標(biāo)度因數(shù)變化無(wú)盲點(diǎn)測(cè)試，從而在分段擬合時(shí)可以精確選擇拐點(diǎn)位置，使得建模和補(bǔ)償精度更高。

（References）：

[1] Dagenais D M, Goldberg L, Moeller R P, et al. Wavelength stability characteristics of a high power amplified superfluorescent source[J]. Journal of Lightwave Technology, 1999, 17(8): 1415-1422.

[2] Zhao Dun-hui, Chen Jia-bin, Han Yong-qiang, et al. Temperature compensation of FOG scale factor based on CPSO-BPNN[C]//Control and Decision Conference. 2010: 2898-2901.

[3] Zhang Xiao-feng, Fan Hua-chun, Liang Yu. FOG output drift compensation based-on temperature gradients[J]. Piezoelectrics & Acoustooptics, 2010, 33(1): 38-40.

[4] 汪剛, 吉世濤, 藍(lán)燕, 等. 光纖陀螺標(biāo)度因數(shù)溫度補(bǔ)償?shù)墓こ虘?yīng)用[J]. 導(dǎo)航與控制, 2012, 11(1): 41-45. Wang Gang, Ji Shi-tao, Lan Yan, et al. FOG scale factor temperature compensation of engineering application[J]. Navigation and Control, 2012, 11(1): 41-45.

[5] 那永林, 李立坤, 吳衍記, 等. 光纖陀螺標(biāo)度因數(shù)的測(cè)試誤差分析[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 2012, 20(4): 473-477. Na Yong-lin, Li li-kun, Wu Yan-ji, et al. Measurement error analysis of FOG’s scale factor[J]. Journal of Chinese Inertial Technology, 2012, 20(4): 473-477.

[6] GJB2426A-2004. 光纖陀螺儀測(cè)試方法[S]. 國(guó)防科學(xué)技術(shù)工業(yè)委員會(huì), 2004. GJB2426A-2004. Test methods for fiber optic gyroscope [S]. Commission of Science Technology and Industry for National Defense, 2004.

[7] 楊曄, 陳馨, 張桂才. 基于角速率積分法的光纖陀螺標(biāo)度因數(shù)不對(duì)稱的測(cè)量方法[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 2012 20(3): 343-347. Yang Ye, Chen Xin, Zhang Gui-cai. Measurement of fiberoptic gyroscope scale factor asymmetry based on angular rate integral technology[J]. Journal of Chinese Inertial Technology, 2012, 20(3): 343-347.

[8] Ko J, Fox D. Learning GP-Bayes filters via Gaussian process latent variable models[J]. Autonomous Robots, 2011, 30: 3-23.

[9] Nguyen-Tuong D, Peters J. Incremental online sparsification for model learning in realtime robot control[J]. Neurocomputing, 2011, 74(11): 1859-1867.

[10] 劉元元, 楊功流, 尹洪亮. 基于雙模型的光纖陀螺溫度補(bǔ)償方法[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 2015, 23(1): 131-136. Liu Yuan-yuan, Yang Gong-liu, Yin Hong-Liang. Temperature compensation for fiber optic gyroscope based on dual models[J]. Journal of Chinese Inertial Technology, 2015, 23(1): 131-136.

[11] Wang Xi-jiang, He Zu-yuan, Hotate K, et al. Reduction of polarization-fluctuation induced drift in resonator fiber optic gyro by a resonator with twin 90° polarization-axis rotated splice[J]. Optics Express, 2010, 18(2): 1677-1683.

[12] Yu Xu-hui, Ma Hui-lian, Jin Zhong-he, et al. Improving thermal stability of a resonator fiber optic gyro employing a polarizing resonator[J]. Optics Express, 2013, 21(1): 358-369.

Rapid modeling and compensation method for scale factor of FOG with continuous rotation

ZHANG Xiao-ya, YANG Zhi-huai, SONG Li-wei, Ma Lin

(Tianjin Navigation Instrument Research Institute, Tianjin 300131, China)

The scale-factor error is one of the major errors of FOG, especially when applied in full temperature. Its influences are more than those of bias error in large angular rate or high-precision applications, and therefore should be modeled and compensated in practical application. Based on the analysis of the scale factor’s error mechanism, a modeling and compensation method for FOG scale-factor in full-temperature is presented based on a continuous rotation method, which is easier and more efficient. More importantly, the inflection points can be obtained effectively, and the precisions of the modeling and compensation will be improved. Experiment results show that the real inflection point, which is 48℃, of a certain type of FOG can be accurately measured. The compensation accuracy of scale factor error by this modeling method is better than 15×10-6, which is improved by about 10%.

fiber optic gyro; scale factor; temperature modeling; compensation

U666.1

A

1005-6734(2016)02-0215-03

10.13695/j.cnki.12-1222/o3.2016.02.014

2015-12-10

2016-03-23

國(guó)防預(yù)研項(xiàng)目（51309010102）

張曉雅（1981—），男，高級(jí)工程師，研究方向?yàn)閷?dǎo)航、制導(dǎo)與控制。E-mail: zhangxy1208@163.com