趙　冰，陳軍波

（1.92941部隊，遼寧葫蘆島125001；2.北京自動化控制設備研究所，北京100074）

光纖陀螺溫度建模與補償方法分析

趙冰1，陳軍波2

（1.92941部隊，遼寧葫蘆島125001；2.北京自動化控制設備研究所，北京100074）

通過光纖陀螺溫度試驗，分析了光纖陀螺的溫度特性；理論上闡述了各項溫度因素對光纖陀螺零偏的影響，并采用逐步回歸分析的方法建立光纖陀螺零偏的溫度數(shù)學模型。通過試驗驗證，采用該模型對光纖陀螺進行溫度漂移的補償，可以有效提高光纖陀螺的測量精度。

光纖陀螺；溫度特性；溫度建模；溫度補償

光纖陀螺是一種新型的固態(tài)慣性測量儀表，與普通機械陀螺和激光陀螺相比，具有結構簡單、體積小、重量輕、耐振動、響應速度快、動態(tài)范圍寬、測量精度高以及抗電磁干擾、無加速度引起的漂移等特點，在航天、航空、航海和兵器領域得到了廣泛的應用[1-4]。光纖陀螺主要由光源、偏振器、PZT相位調制器、光纖耦合器和光纖線圈等組成，其主要組成部件均為溫度敏感器件，環(huán)境溫度對光纖陀螺的性能影響較大；而工程應用要求光纖陀螺應具有較寬的工作溫度范圍及較強的溫度適應性。因此，對光纖陀螺進行溫度特性分析，建立溫度影響的數(shù)學模型，從而對光纖陀螺進行溫度補償，提高陀螺測量精度是實現(xiàn)光纖陀螺工程化應用的必要條件[5]。

本文從溫度試驗入手，分析光纖陀螺的溫度特性。通過理論分析和試驗結果找出影響光纖陀螺精度的各項溫度因素，建立光纖陀螺零偏受溫度影響的數(shù)學模型。利用該模型對光纖陀螺輸出進行溫度補償，可以有效提高光纖陀螺的測量精度。

1　光纖陀螺溫度試驗及溫度特性分析

1.1溫度試驗設計

通過設計溫度試驗，對光纖陀螺的溫度特性進行分析[6-7]。試驗項目包括固定溫度點長時間工作試驗和變溫試驗。

固定溫度點長時間工作試驗：在-40℃～+60℃溫度范圍內，取若干溫度點，測試光纖陀螺在該溫度點環(huán)境下長時間工作（1 h）的輸出特性。

變溫試驗：在全溫度范圍內，選擇不同的溫度變化速率，測試在溫度變化的環(huán)境下，光纖陀螺的輸出特性。

溫度試驗設備組成包括工控機、可調變溫速率溫箱、光纖陀螺和測試臺，如圖1所示。

圖1　光纖陀螺溫度試驗系統(tǒng)組成示意圖Fig.1 Temperature test system schematic diagram of FOG

其中，測試臺實現(xiàn)對光纖陀螺輸出和光纖陀螺內部溫度值的測量和數(shù)據(jù)記錄，溫度傳感器的精度為±0.5℃；溫度控制精度為0.1℃；光纖陀螺常溫精度零偏穩(wěn)定性優(yōu)于0.1(°)/h，零偏重復性不大于0.05(°)/h。

1.2試驗結果分析

對光纖陀螺進行恒溫溫度試驗，得到光纖陀螺零偏隨溫度變化的曲線，2次測試結果如圖2所示。

圖2　恒溫定點試驗結果Fig.2 Test result of constant temperature fixed-point

對光纖陀螺進行-40℃～+60℃的升、降溫測試，得到溫度變化率分別為0.5℃/min、1℃/min、2℃/min的試驗結果，如圖3～5所示。曲線均為500點相鄰值平滑濾波結果。

圖3　溫度變化率0.5℃/min下試驗結果Fig.3 Test result of variable temperature test at0.5℃/min

圖4　溫度變化率1℃/min下試驗結果Fig.4 Test result variable temperature at1℃/min

圖5　溫度變化率2℃/min下試驗結果Fig.5 Test result of variable temperature at2℃/min

從不同溫度變化率對應光纖陀螺輸出曲線上可以看出：

1）在-40℃～+60℃的變化過程中，不同的環(huán)境溫度變化率，同一光纖陀螺零偏變化趨勢相同。

2）升溫過程和降溫過程，光纖陀螺輸出不一致。原因是：首先，光纖陀螺內部溫度相對于環(huán)境溫度有延遲，而升溫和降溫過程相對于環(huán)境溫度（溫箱溫度）的延遲剛好相反，所以造成升降溫陀螺輸出曲線呈環(huán)狀；其次，光纖陀螺的零位漂移不僅僅和溫度有關，還和溫度梯度、溫度變化率等因素有關。

光纖陀螺的輸出零偏主要是由內部器件的測量誤差引起的，而占主要部分的溫度影響是由于光纖陀螺中的溫度敏感元件引起的。光纖陀螺的各個組成部分中，光纖敏感環(huán)圈是關鍵件，光源、Y波導、檢測電路是重要件，耦合器、探測器是一般件。而作為光纖陀螺關鍵件的光纖敏感環(huán)圈是溫度敏感器件，溫度變化將使光纖敏感線圈產(chǎn)生Shupe效應，導致陀螺測量誤差；而重要件的工作特性隨溫度變化也將受影響。

可見，溫度對于光纖陀螺漂移的影響是其對于光纖陀螺內部元件影響的綜合效應，而溫度對陀螺溫度敏感元件的影響是非線性的，各個部件的溫度特性不一致。溫度敏感元件對陀螺漂移影響的不確定性以及目前的工藝條件和元器件的限制使得對內部溫度敏感器件分別進行溫度建模補償較難實現(xiàn)。在現(xiàn)有條件下確定光纖陀螺溫度特性的有效方法是將其視為一個“黑箱”，根據(jù)陀螺的輸出和溫度之間相互關系建立模型。這要求陀螺的輸出和溫度之間存在相關性，同時陀螺漂移要具備一定的重復性。試驗表明，光纖陀螺具備上述2個條件。研究認為，光纖陀螺的溫度漂移主要受溫度、溫度變化率、溫度梯度影響，而溫度梯度對現(xiàn)有光纖陀螺漂移的影響相對較小，所以建立的溫度補償模型暫不考慮溫度梯度項[4]。

2　光纖陀螺溫度補償模型

采用定溫試驗標定模型結合變溫試驗數(shù)據(jù)模型，建立綜合補償模型[8]。

定溫試驗標定模型主要用于補償光纖陀螺在該溫度下的零偏，變溫數(shù)據(jù)模型用于補償變溫條件下的陀螺零偏溫變誤差。模型形式：

2.1定溫標定溫度補償模型

根據(jù)定溫溫度標定試驗得到標定結果，建立定溫標定試驗的溫度多項式補償模型：

對多項式參數(shù)的辨識，辨識算法采用最小二乘算法。確定模型的結構——多項式模型的階次滿足2點要求：由陀螺的零偏溫度曲線可以大致選擇多項式模型的階次；逐次增加模型的階次，并同時計算擬合點的擬合均方根誤差，在擬合誤差允許的情況下，選擇較低階次的多項式模型。

實際計算表明，對于模型參數(shù)求解的正則方程組，當模型階次n較大時，正則方程組往往是病態(tài)的，為了保證計算精度，擬合多項式不宜超過5階，且擬合多項式只是在原始數(shù)據(jù)范圍內估計精度較高，超出這個范圍用擬合多項式做預報，可能得出錯誤的結果。故擬合算法中，待擬合模型的最高階次M設定為5，這對于大多數(shù)情況來說都可以滿足的。

依據(jù)上面的算法，采用數(shù)據(jù)處理軟件Matlab進行擬合算法的程序編寫，建立光纖陀螺的定溫零偏模型，得到多項式模型：

2.2變溫補償模型

2.2.1分析方法

建立變溫數(shù)據(jù)溫度模型時，因影響模型的因素理論上不易分析，所以采用逐步回歸法建立影響變溫漂移的多變量模型。

逐步回歸法是從與y有關的變量中選取對y有顯著性影響的變量來建立回歸方程的一種方法[9-11]。其基本思想是：在所有考慮的因素中，按對y的作用的顯著程度的大小，取最顯著的變量，逐一引入回歸方程，對y作用不顯著的那些變量自始自終都未被引入。另一方面，已被引入回歸方程的變量，在有引入的新變量后若發(fā)現(xiàn)其對y的作用變?yōu)椴伙@著時，則隨時從回歸方程中剔出，直至沒有新的變量能引入方程，且已引入方程的所有變量均不需剔出為止。利用這些變量就可以建立起線性多元回歸方程。

設共有 p-1個自變量：x1,x2,…,xp-1，記y≡xp。對所獲得的n組測量得到的樣本數(shù)據(jù)() xt1,xt2,…,xtp，t=1,2,…,n，計算各個變量的均值-xj及偏差平方和的算術根σj（j=1,2,…,p）：

建立zp關于z1,z2,…,zp-1的線性回歸方程：

由于zj都是標準化后的變量，故：

從而zp關于z1,z2,…,zp-1的線性回歸方程記為：

化為y關于x1,x2,…,xp-1的線性方程：

由變量的關系，可以推導出：

2.2.2動態(tài)溫度補償模型的建立

由前面的分析可知，光纖陀螺的零偏與溫度的關系是復雜的非線性關系，即與溫度、溫度變化率及其交叉項都可能有關系[11]?？紤]三階變量，分別令：作為線性多變量模型的自變量，線性多變量模型：

模型中的常數(shù)項：

式（14）中：Y為陀螺零偏；X為自變量；P為模型的參數(shù)；in為經(jīng)逐步回歸分析后保留在模型內的自變量。

對陀螺進行逐步回歸溫度補償模型建立，得到逐步回歸模型：

3　試驗驗證

利用建立的溫度補償模型對光纖陀螺變溫數(shù)據(jù)進行溫度補償，得到補償效果如圖6～8所示，不同溫度變化率數(shù)據(jù)補償效果見表1。

圖6　補償1℃/min升降溫數(shù)據(jù)效果Fig.6 Compensation1℃/minwarming and cooling effect of data

圖7　補償0.5℃/min升降溫數(shù)據(jù)效果Fig.7 Compensation0.5℃/minwarming and cooling effect of data

圖8　補償2℃/min升降溫數(shù)據(jù)效果Fig.7 Compensation2℃/minwarming and cooling effect of data

表1　不同溫度變化率數(shù)據(jù)補償效果Tab.1 Compensation effect of different temperature rate data

從表1統(tǒng)計的補償效果看，溫度漂移補償模型對于不同溫度變化率的數(shù)據(jù)均有一定的補償效果。因此，采用逐步回歸方法建立的多變量補償模型對光纖陀螺的溫度漂移補償是有效的。

4　結束語

本文通過光纖陀螺的溫度試驗，分析了陀螺溫度特性；利用逐步回歸法建立了光纖陀螺的變溫補償模型。通過試驗驗證可以得出以下結論：

1）根據(jù)逐步回歸分析，光纖陀螺零偏與溫度、溫度變化率以及其交叉項具有相關性，而與溫度變化率的高次項相關性較弱。采用多變量模型對光纖陀螺零偏進行補償，可以有效提高光纖陀螺的測量精度。

2）對于不同的光纖陀螺，其定溫試驗標定模型和變溫數(shù)據(jù)補償模型的參數(shù)和模型的結構可能會有所變化，需要分別對每一個陀螺進行溫度標定。

[1]孫麗，王德釗.光纖陀螺的最新進展[J].航天控制，2013，16（3）：75-80. SUN LI，WANG DEZHAO.Recent development of fiber optic gyro[J].Aerospace Control，2013，16（3）：75-80.（in Chinese）

[2]于明飛.光纖陀螺技術的軍事應用及前景[J].計算機測量技術，2011，27（5）：785-788. YU MINGFEI.The application and prospect of fiber optic gyrocope technology in military[J].Computer Measurement Technology，2011，27（5）：785-788.（in Chinese）

[3]郭秀中.慣性系統(tǒng)陀螺儀理論[M].北京：國防工業(yè)出版社，1996：175-178. GUO XIUZHONG.Gyroscope inertial system theory[M]. Beijing：National Defense Industry Press，1996：175-178.（in Chinese）

[4]劉蘭芳，陳剛，金國良.光纖陀螺儀基本原理與分類[J].現(xiàn)代防御技術，2007，32（4）：36-40. LIU LANFANG，CHEN GANG，JIN GUOLIANG.Principle and classification of fiber optic gyrocope[J].ModernDefense Technology，2007，32（4）：36-40.（in Chinese）

[5]譚顯裕.光纖陀螺的關鍵技術及其軍用研究[J].航空兵器，2003（3）：32-35. TAN XIANYU.The key technology of fiber optic gyro and its military research[J].Aviation Weapon，2003（3）：32-35.（in Chinese）

[6]孟照魁，崔佳濤，章博，等.高精度光纖陀螺溫度實驗研究[J].宇航學報，2007，28（3）：580-583. MENG ZHAOKUI，CUI JIATAO，ZHANG BO，et al.Experiment research of the high precision fiber optic gyrocope tempreature[J].Journal of Astronautics，2007，28（3）：580-583.（in Chinese）

[7]張春喜，李森，陳光見.光纖陀螺慣性測量組合的數(shù)字溫控系統(tǒng)設計[J].中國慣性技術學報，2012，16（5）：541-547. ZHANG CHUNXI，LI SEN，CHEN GUANGJIAN.Digital temperature control system design for FOG inertial measurement unit[J].Journal of Chinese Inerial Technology，2012，16（5）：541-547.（in Chinese）

[8]張鵬飛，龍興武.二頻機抖激光陀螺零偏的溫度特性的逐步回歸分析[J].光學技術，2014，32（5）：738-740. ZHANG PENGFEI，LONG XINGWU.Analysis on temperature characteristic of mechanically dithered RLG’bias with a method of stepwise regression[J].Optical Technique，2014，32（5）：738-740.（in Chinese）

[9]何曉群，劉文卿.應用回歸分析[M].北京：中國人民大學出版社，2001：41-45. HE XIAOQUN，LIU WENQING.Application of regression analysis[M].Beijing：Renmin University of China Press，2001：41-45.（in Chinese）

[10]費業(yè)泰.誤差理論與數(shù)據(jù)處理[M].北京：機械工業(yè)出版社，2005：118-152. FEI YETAI.The error theory and data processing[M].Beijing：Mechanical Industry Press，2005：118-152.（in Chinese）

[11]梅碩基.慣性儀器測試與數(shù)據(jù)分析[M].西安：西北工業(yè)大學出版社，1991：138-186. MEI SHUOJI.Inertial instrument test and data analysis [M].Xi’an：Northwestern Polytechnic University Press，1991：138-186.（in Chinese）

Analysis on Temperature Modeling and Compensation Mothed of FOG

ZHAO Bing1,CHEN Junbo2
（1.The 92941stUnit of PLA,Huludao Liaoning 125001,China; 2.Beijing Institute of Automation and Control Equipment,Beijing 100074,China）

Based on temperature experiment,the temperature characteristic of FOG’s modeling was studied,the factors ef?fected on gyro’s modeling was analyzed.From the aspect of theory the effect on the FOG’s modeling was studied,and some useful methods was proposed in this paper to found the temperature compensating model.The model was used to compensate the FOG’s modeling for reduce gyro’s error.The results of compensation showed that the stability was en?hanced and the compensation model was effective.

FOG;temperature characteristic;temperature modeling;temperature compensation

V241.5；TP212

A

1673-1522（2016）01-0012-05

10.7682/j.issn.1673-1522.2016.01.003

2015-09-16；

2015-12-22

趙冰（1965-），男，高工，大學。