慣組中光纖陀螺溫度漂移的二維插值建模研究

白子杰，蔡春龍，牟玉濤

(北京航天時(shí)代光電科技有限公司，北京 100094)

0 引言

光纖陀螺是基于薩格納克(Sagnac)效應(yīng)進(jìn)行角速度測(cè)量的光學(xué)慣性儀表。因其具備可靠性高、體積功耗小、生產(chǎn)工藝相對(duì)簡(jiǎn)單等優(yōu)點(diǎn)，被越來越廣泛地應(yīng)用于航空航天及武器裝備領(lǐng)域中。以光纖陀螺為基礎(chǔ)的捷聯(lián)慣性導(dǎo)航系統(tǒng)也正逐步向高精度發(fā)展。

光纖陀螺工作于慣性測(cè)量組合中，輸出慣組自身三軸方向的角速度，用于對(duì)應(yīng)慣導(dǎo)系統(tǒng)的對(duì)準(zhǔn)與導(dǎo)航計(jì)算。光纖陀螺因本身儀表的精度限制和所處不同工作環(huán)境的影響，輸出與真實(shí)值會(huì)存在誤差，其中零偏是主要誤差之一。光纖陀螺的工作原理和器件構(gòu)成使其對(duì)溫度較為敏感，溫度相關(guān)變量的變化會(huì)使其輸出產(chǎn)生很大波動(dòng)，從而使零偏產(chǎn)生溫度漂移。零偏的溫度漂移會(huì)對(duì)慣組所在慣導(dǎo)系統(tǒng)的導(dǎo)航精度產(chǎn)生影響，因此該誤差是光纖慣導(dǎo)向高精度發(fā)展過程中必須解決的問題之一。

20世紀(jì)80年代初期D. M. Shupe發(fā)現(xiàn)在干涉型光纖陀螺中，光纖環(huán)溫度變化時(shí)會(huì)產(chǎn)生熱致非互易性相位延遲，并命名為Shupe效應(yīng)[1]。Shupe效應(yīng)從理論上分析了溫度漂移產(chǎn)生的原因。在此之后F.Mohr研究證實(shí)了使用四極對(duì)稱繞法繞制光纖環(huán)可以極大地降低光纖陀螺中因溫度變化產(chǎn)生的相位誤差，但是這種方法繞制的光纖環(huán)過于理想，在實(shí)際應(yīng)用中并不能完美實(shí)現(xiàn)。

目前，針對(duì)溫度漂移誤差可以使用溫度控制或溫度建模補(bǔ)償來消除。當(dāng)光纖慣組的工作環(huán)境溫度范圍較大時(shí)，溫度控制的升溫幅度較大、啟動(dòng)時(shí)間過長(zhǎng)，升溫相對(duì)較難實(shí)現(xiàn)且不利于系統(tǒng)的快速響應(yīng)。溫度建模補(bǔ)償則能更好地解決此問題。目前，溫度建模多采用多項(xiàng)式模型的擬合方法，如最小二乘法，因其原理簡(jiǎn)單、用法成熟被廣泛應(yīng)用于各種產(chǎn)品中。同時(shí)隨著計(jì)算機(jī)技術(shù)的發(fā)展，神經(jīng)網(wǎng)絡(luò)也得到了人們的重視。單獨(dú)使用各種不同的神經(jīng)網(wǎng)絡(luò)或與其他方法組合建立非線性模型的方法也得到了研究。張瀟等在研究分類建模的過程中使用了最小二乘法建模[2]；周琪等研究了利用自適應(yīng)神經(jīng)模糊推理系統(tǒng)(Adaptive Network-based Fuzzy Inference System，ANFIS)對(duì)光纖陀螺溫度漂移誤差進(jìn)行建模補(bǔ)償[3]；胡國(guó)欣等在研究基于跟蹤微分器的建模時(shí)，使用了徑向基函數(shù)(Radial Basis Function，RBF)神經(jīng)網(wǎng)絡(luò)的方法進(jìn)行建模[4]；李健等研究了基于小波神經(jīng)網(wǎng)絡(luò)的光纖陀螺溫度建模補(bǔ)償[5]。

但是多項(xiàng)式模型的擬合方法在較復(fù)雜的輸出情況下，擬合效果有限；而神經(jīng)網(wǎng)絡(luò)的方法計(jì)算量較大，所需時(shí)間較長(zhǎng)，且無法直觀得到所建立模型，不利于在工程上的應(yīng)用。本文經(jīng)研究發(fā)現(xiàn)，可以利用二維插值模型對(duì)溫度漂移誤差進(jìn)行補(bǔ)償[6]，并對(duì)比多項(xiàng)式模型進(jìn)行分析。

本文首先通過研究光纖陀螺的工作原理及溫度漂移誤差產(chǎn)生的原因，分析了相關(guān)變量對(duì)溫度漂移誤差的影響。在對(duì)比使用多項(xiàng)式模型中最小二乘法建模補(bǔ)償?shù)慕Y(jié)果后，基于二維插值的方法對(duì)光纖陀螺的輸出進(jìn)行建模補(bǔ)償，并通過試驗(yàn)對(duì)其結(jié)果進(jìn)行了驗(yàn)證。

1 溫度漂移誤差機(jī)理分析

光纖陀螺是基于薩格納克效應(yīng)進(jìn)行工作的，即沿閉合光路相向傳播的兩束光波之間的相位差正比于閉合光路法向的輸出角速度。在光纖陀螺相對(duì)慣性空間存在轉(zhuǎn)動(dòng)時(shí)，兩束相向傳播的光傳播路程發(fā)生變化，從而產(chǎn)生時(shí)間差Δt，最終引起相應(yīng)的相位差，其相位差如式(1)

(1)

其中，λ0為光的波長(zhǎng)，Ω為轉(zhuǎn)速，Δφ為相位差，l為光纖環(huán)長(zhǎng)度，r為光纖環(huán)半徑。

光纖陀螺所在工作環(huán)境溫度的不同以及在工作中自身溫度的不斷上升，使其輸出發(fā)生變化，即產(chǎn)生溫度漂移誤差。光纖線圈是光纖陀螺的重要部件之一，當(dāng)光纖線圈上某段光纖的溫度隨著時(shí)間發(fā)生變化時(shí)，自身的折射率也會(huì)隨之變化。相向傳播的光在光纖線圈中經(jīng)過該段光纖的時(shí)間不同，因其折射率的不同，它們所經(jīng)歷的光程也不同，從而產(chǎn)生相位誤差，被稱為非互易性相位延遲，即Shupe效應(yīng)。這種相位誤差無法與薩格納克效應(yīng)相區(qū)分，從而產(chǎn)生了零偏的溫度漂移誤差[7-9]。

假設(shè)光纖線圈的總長(zhǎng)度為L(zhǎng)，其中因溫度變化產(chǎn)生的相位延遲為

(2)

其中，k為真空中的波數(shù)，即波長(zhǎng)λ的倒數(shù)，理論物理中定義為k=2π/λ，n為折射率，α為光纖線性膨脹系數(shù)，L為光纖總長(zhǎng)，ΔT為光纖溫度的變化。

當(dāng)兩束光波分別沿順時(shí)針(Clock Wise，CW)與逆時(shí)針(Counter Clock Wise，CCW)傳播時(shí)，達(dá)到初始端時(shí)刻為t，則經(jīng)過l處的時(shí)刻分別為

(3)

其中，c為真空中光速，c0為光在光纖中的傳播速度c0=c/n。

代入(2)式可求出兩束光分別產(chǎn)生的相位延遲ΔφCW與ΔφCCW。相減后，再經(jīng)過整理得到相位差公式

(4)

為了抑制Shupe效應(yīng)，減少相位誤差，在實(shí)際中可以采用對(duì)稱法來繞制光纖線圈，這樣可以很大程度上保證光纖環(huán)對(duì)稱部分的一小段光纖具有相同的溫度，從而具有相同的折射率。但是這種方法要求光纖線圈繞制處于絕對(duì)理想的狀態(tài)中，在實(shí)際中無法完美實(shí)現(xiàn)，并不能完全解決Shupe效應(yīng)帶來的誤差。

同時(shí)，本文使用的某型號(hào)光纖慣組要求工作環(huán)境溫度范圍較大。在溫度不同的環(huán)境下，光纖陀螺內(nèi)部的光源、光纖環(huán)、探測(cè)器等光學(xué)器件以及慣組中其他結(jié)構(gòu)零件，由于熱脹冷縮造成的形變和擠壓等也會(huì)產(chǎn)生漂移誤差。慣組在通電啟動(dòng)后，光纖陀螺等各部件會(huì)產(chǎn)生熱量使自身溫度升高并向外傳播，此過程光纖陀螺的溫度持續(xù)升高，直到工作一定時(shí)間后，溫度才會(huì)達(dá)到穩(wěn)態(tài)。在溫度由低升高到最終穩(wěn)態(tài)的過程中，溫度變化率在從大向小發(fā)生變化，相應(yīng)也造成了溫度漂移誤差。最后慣組在工作時(shí)，內(nèi)部會(huì)產(chǎn)生一個(gè)變化的溫度場(chǎng)，若當(dāng)前環(huán)境溫度變化較快時(shí)，內(nèi)部器件因受熱不均也會(huì)產(chǎn)生溫度延遲，造成漂移誤差[13]。

本文使用的光纖慣組陀螺在某溫度下工作一段時(shí)間，其溫度、溫度變化率、角速度原始輸出如圖1所示，可見在慣組正常工作時(shí)陀螺因溫度的變化使輸出產(chǎn)生了較大的漂移誤差，遠(yuǎn)遠(yuǎn)高于陀螺本身的精度，因此必須針對(duì)溫度相關(guān)信息對(duì)陀螺進(jìn)行建模補(bǔ)償，抑制溫度漂移誤差。

圖1 光纖陀螺溫度、溫變率與原始輸出Fig.1 Temperature, temperature change rate and original output of fiber optic gyroscope

2 誤差模型的分析與建立

光纖陀螺安裝于慣組中，陀螺轉(zhuǎn)動(dòng)將兩束光傳播產(chǎn)生的相位差轉(zhuǎn)化為電信號(hào)，然后轉(zhuǎn)換為脈沖信號(hào)輸出，信號(hào)輸出在DSP中經(jīng)過標(biāo)度因數(shù)等參數(shù)的計(jì)算后得到慣組的角速度。光纖陀螺在工作時(shí)其輸出與輸入角速度存在如下簡(jiǎn)化關(guān)系

D=Kω+D0+ΔDT+ε

(5)

其中，K為陀螺標(biāo)度因數(shù)，ω為輸入角速度，D0為零偏，ΔDT為零偏溫度漂移誤差，ε為隨機(jī)誤差。其中ΔDT主要由Shupe效應(yīng)產(chǎn)生的非互易性相位延遲所造成。

光纖陀螺在工作時(shí)環(huán)境較為復(fù)雜,其輸入輸出存在某種特定的關(guān)系，同樣溫度、溫度變化率等溫度相關(guān)變量與溫度漂移誤差也存在某種特定的關(guān)系。通過第1節(jié)的理論分析得出，溫度與溫度變化率能夠直接造成光纖陀螺漂移誤差的產(chǎn)生，但因?yàn)樵趹T組中溫度對(duì)其他部件也造成影響，多種效應(yīng)共同作用產(chǎn)生最終的漂移誤差，所以很難理論分析出具體的數(shù)學(xué)關(guān)系。因此需要使用系統(tǒng)辨識(shí)的方法，把慣組視為一個(gè)黑箱，通過對(duì)慣組試驗(yàn)所得數(shù)據(jù)進(jìn)行分析，建立數(shù)學(xué)模型，使用模型擬合出這種特定的關(guān)系，從而用數(shù)學(xué)的方法補(bǔ)償?shù)艄饫w陀螺的溫度漂移誤差，使其輸出更加逼近真實(shí)情況。

2.1 多項(xiàng)式模型

傳統(tǒng)建模方法多使用多項(xiàng)式模型進(jìn)行擬合，最小二乘法是其中一種在工程中被廣泛應(yīng)用的多項(xiàng)式擬合方法。所謂最小二乘，是一種數(shù)學(xué)上的優(yōu)化方法，使擬合出的模型與實(shí)際數(shù)據(jù)的偏差平方和最小。但在實(shí)際試驗(yàn)中發(fā)現(xiàn)，僅使用單一的溫度T與溫度變化率ΔT作為輸入時(shí)，往往擬合效果并不理想。分析發(fā)現(xiàn)其他的溫度相關(guān)變量，如溫度與溫度變化率的不同冪次項(xiàng)、溫度與溫度變化率的乘積，以及各自不同冪次項(xiàng)與對(duì)方的乘積都會(huì)對(duì)溫度漂移誤差產(chǎn)生影響，因此應(yīng)在自變量的選擇中添加上述變量。

在MATLAB中對(duì)光纖陀螺輸出與溫度、溫度的不同冪次項(xiàng)、溫度變化率、溫度變化率的不同冪次項(xiàng)、溫度與溫度變化率的乘積等變量進(jìn)行逐步回歸分析[14]。在分析中將變量逐個(gè)引入模型，每次引入一個(gè)新變量進(jìn)行F檢驗(yàn)，同時(shí)對(duì)已經(jīng)引入的變量進(jìn)行t檢驗(yàn)，在檢驗(yàn)的過程中刪除不顯著的變量，使模型中任意一個(gè)變量都是顯著的。當(dāng)最終引入所有變量后，不再有新的變量引入時(shí)，模型中剩下的所有變量均為顯著。

根據(jù)顯著性分析所得結(jié)果判斷，選取下列各項(xiàng)作為自變量

建立陀螺輸出與上述變量相關(guān)的模型為

Y(T,ΔT)=Y0+f1(T)+f2(T2)+f3(T3)+

f4(ΔT)+f5(ΔT2)+f6(ΔT3)+

f7(T×ΔT)+f8(T2×ΔT)+

f9(T×ΔT2)+f10(T3×ΔT)+

f11(T×ΔT3)

(6)

對(duì)上述模型進(jìn)行最小二乘擬合：設(shè)X是某一確定性常值向量，無法直接觀測(cè)到X，只能觀測(cè)到X各個(gè)分量的線性組合Z。其中關(guān)系如下

Z=HX+V

(7)

(8)

經(jīng)過推導(dǎo)得出使式(8)成立的條件為

(9)

(10)

因此，構(gòu)造矩陣H包含所有溫度相關(guān)變量自變量，同時(shí)構(gòu)造參數(shù)矩陣X包含模型中溫度相關(guān)變量的系數(shù)，y為陀螺原始輸出，可得

y=H×X+ε

(11)

利用最小二乘法得

(12)

此時(shí)，求得了此多項(xiàng)式模型中的所有參數(shù)，代入其他數(shù)據(jù)即可完成模型的建立。

2.2 二維插值模型

在以往的試驗(yàn)分析中，使用多項(xiàng)式模型往往不能很好地適應(yīng)不同位置下的輸出，尤其在試驗(yàn)數(shù)據(jù)較為復(fù)雜的情況下，會(huì)產(chǎn)生很大的誤差。因?yàn)榘褢T組在不同溫度下工作的數(shù)據(jù)放在一起時(shí)波動(dòng)較大，在不同溫度下升溫速率也不相同，許多關(guān)鍵的數(shù)據(jù)點(diǎn)在多項(xiàng)式擬合中會(huì)被當(dāng)作誤差所忽略。所以在最小二乘建模中，僅在選取的自變量上優(yōu)化也不能很好地?cái)M合出整個(gè)試驗(yàn)范圍內(nèi)的輸出結(jié)果。針對(duì)這一規(guī)律，本文發(fā)現(xiàn)使用插值法更適合在該情況下的模型建立。插值法會(huì)使模型一定通過所選中的數(shù)據(jù)點(diǎn)，保留數(shù)據(jù)中這些具有特點(diǎn)的真實(shí)結(jié)果，同時(shí)也能使模型更好地符合真實(shí)情況。

根據(jù)之前理論分析，使用插值法建模時(shí)需要使用溫度T與溫度變化率ΔT這2個(gè)變量進(jìn)行二維插值。一般的二維插值使用的數(shù)據(jù)為等間距且單調(diào)變化的網(wǎng)格形式，但是慣組實(shí)際工作時(shí)的數(shù)據(jù)并不能按照溫度與溫變率這2個(gè)方向排列成標(biāo)準(zhǔn)的網(wǎng)格形式，無論是溫度還是溫變率均難以在試驗(yàn)中控制為不同的標(biāo)準(zhǔn)大小。因此在建模時(shí)只能使用基于不規(guī)則的散點(diǎn)數(shù)據(jù)進(jìn)行插值的方法。

在插值法建模中，第一要考慮已知數(shù)據(jù)點(diǎn)的選擇，需要讓選出的數(shù)據(jù)點(diǎn)大致均勻地分布在需要估計(jì)出的點(diǎn)的周圍，在數(shù)據(jù)密集處選擇較多的數(shù)據(jù)點(diǎn)，在數(shù)據(jù)稀疏處選擇較少的數(shù)據(jù)點(diǎn)，分布不均的數(shù)據(jù)往往會(huì)造成較大的插值誤差。第二要注意數(shù)據(jù)點(diǎn)的篩選，過多的數(shù)據(jù)點(diǎn)并不會(huì)使插值的準(zhǔn)確度提高，反而會(huì)因?yàn)閿?shù)據(jù)過多造成冗余，掩蓋住真實(shí)的結(jié)果，對(duì)某一次的試驗(yàn)數(shù)據(jù)過度擬合后，對(duì)其他位置試驗(yàn)數(shù)據(jù)的重復(fù)適應(yīng)性會(huì)相應(yīng)變差。過多的數(shù)據(jù)點(diǎn)也會(huì)造成數(shù)據(jù)計(jì)算量過大，增加計(jì)算時(shí)間。慣組在啟動(dòng)初期溫度變化較大，相同時(shí)間內(nèi)數(shù)據(jù)較少且變化較為規(guī)律，百秒平滑后每個(gè)數(shù)據(jù)均可選擇；隨著溫度變化幅度變小后，溫度在較小范圍內(nèi)波動(dòng)，相同時(shí)間內(nèi)數(shù)據(jù)變多，為了更好地反映數(shù)據(jù)變化規(guī)律，每隔2～3個(gè)數(shù)據(jù)擇一個(gè)[15]。

經(jīng)過多次測(cè)試，在每個(gè)啟動(dòng)溫度點(diǎn)下根據(jù)溫度變化的快慢且參考變化趨勢(shì)的規(guī)律性，選擇適量的原始數(shù)據(jù)點(diǎn)，畫出以溫度與溫變率為自變量，陀螺輸出為因變量的三維曲面圖(其中x、y軸為溫度與溫變率，z軸為陀螺原始輸出)，圖2所示即為二維插值模型的基礎(chǔ)。

圖2 陀螺原始輸出曲面Fig.2 Gyro original output surface

在圖2的基礎(chǔ)模型上進(jìn)行插值，得出自變量范圍內(nèi)所需新的溫度與溫變率下陀螺的輸出值。相較于一維插值，二維插值更為復(fù)雜，但許多一維插值的方法無法適用。對(duì)于曲面上的插值，首先把數(shù)據(jù)點(diǎn)投影至自變量所在平面，需要根據(jù)已知的數(shù)據(jù)點(diǎn)表示出自變量平面內(nèi)任意一點(diǎn)，當(dāng)使用的數(shù)據(jù)點(diǎn)是標(biāo)準(zhǔn)的網(wǎng)格形式時(shí)較容易實(shí)現(xiàn)。但是因?yàn)閼T組的數(shù)據(jù)為不規(guī)則的散點(diǎn)，不能根據(jù)橫縱坐標(biāo)分割成標(biāo)準(zhǔn)的矩形，因此要重新對(duì)x、y軸所在的二維平面進(jìn)行分割。Delaunay三角剖分法可以根據(jù)一定的算法，建立空間中散點(diǎn)的聯(lián)系，最終使平面形成三角網(wǎng)格的形式，然后可以依據(jù)此三角網(wǎng)格進(jìn)行插值。

首先，把溫度與溫變率所在平面組成的平面點(diǎn)集按Delaunay方法進(jìn)行三角割分，然后進(jìn)行數(shù)據(jù)的插值計(jì)算。假設(shè)經(jīng)過分割后平面內(nèi)某3個(gè)點(diǎn)組成三角形，3個(gè)點(diǎn)分別為A、B、C，其中每個(gè)點(diǎn)對(duì)應(yīng)的值為f(A)、f(B)、f(C),此三角形內(nèi)任意一點(diǎn)X可以用相對(duì)此三角形的重心坐標(biāo)來表示

X=lA+mB+nC

(13)

其中，l、m、n為坐標(biāo)系數(shù)，當(dāng)0

X點(diǎn)輸出線性插值表示為

f(X)=l·f(A)+m·f(B)+n·f(C)

(14)

因此，想要求得任意點(diǎn)的線性插值，需要求出式(14)中的系數(shù)l、m、n。假設(shè)三角形3個(gè)頂點(diǎn)為3個(gè)已知數(shù)據(jù)點(diǎn)，其中每個(gè)頂點(diǎn)溫度表示為T1、T2、T3，溫變率表示為ΔT1、ΔT2、ΔT3，陀螺輸出表示為y1、y2、y3。所求插值點(diǎn)X的溫度為T、溫變率為ΔT、陀螺輸出為y。

陀螺輸出可使用自變量線性表示為

yi=aTi+bΔTi+c

(15)

以此建立三元一次方程組

(16)

求解可得

(17)

對(duì)于三角形任意一點(diǎn)X的輸出，可表示為

(18)

將式(16)代入式(17)中，可求得系數(shù)l、m、n

(19)

將式(19)代入式(13)中，即可求得任意點(diǎn)X在二維線性插值下的陀螺輸出。依照此方法可以根據(jù)圖2中的原始模型，插值出任意定義域范圍內(nèi)陀螺的輸出值，以此進(jìn)行建模補(bǔ)償。

3 試驗(yàn)設(shè)計(jì)與驗(yàn)證

根據(jù)模型需求搭建試驗(yàn)平臺(tái)，試驗(yàn)用光纖陀螺采用了光電公司某型號(hào)光纖慣組產(chǎn)品中的光纖陀螺，其精度達(dá)到5‰。試驗(yàn)以此慣組為基礎(chǔ)，同時(shí)還包含溫箱、直流穩(wěn)壓電源、通信電纜、數(shù)據(jù)串口盒、測(cè)試電腦等設(shè)備。試驗(yàn)具體設(shè)計(jì)為：使慣組在不同溫度下啟動(dòng)并工作2h，直到內(nèi)部陀螺溫度趨于穩(wěn)定并記錄采集數(shù)據(jù)。

首先進(jìn)行第1組試驗(yàn)，將慣組放置于溫箱中，在-40℃～60℃范圍內(nèi)設(shè)置以下溫度點(diǎn)：-40℃，-25℃，-10℃，5℃，20℃，35℃，45℃，60℃。在第1個(gè)溫度點(diǎn)充分保溫，在確保慣組中的陀螺處于設(shè)定溫度點(diǎn)后，慣組上電進(jìn)行測(cè)試。2h后慣組斷電，在第2個(gè)溫度點(diǎn)再次充分保溫后上電測(cè)試。依次完成所有溫度點(diǎn)的測(cè)試并保存好數(shù)據(jù)后，改變慣組位置進(jìn)行第2組試驗(yàn)，同樣按照設(shè)定的溫度重復(fù)進(jìn)行上述步驟，用于驗(yàn)證模型效果。實(shí)驗(yàn)平臺(tái)如圖3所示。

圖3 試驗(yàn)平臺(tái)Fig.3 Test platform

本慣組采樣頻率為1000Hz，所得數(shù)據(jù)在疊加了白噪聲后，其真實(shí)的輸出結(jié)果被淹沒在噪聲中無法看出變化，同時(shí)過多的原始數(shù)據(jù)會(huì)造成冗余對(duì)于插值建模反而不利。因此對(duì)于原始數(shù)據(jù)進(jìn)行百秒平滑，減少噪聲對(duì)慣性器件輸出的影響，消除了不穩(wěn)定的噪聲對(duì)輸出數(shù)據(jù)的影響，之后對(duì)經(jīng)過平滑的溫度數(shù)據(jù)使用差分法求得對(duì)應(yīng)的溫度變化率。經(jīng)過平滑后的陀螺溫度、溫度變化率、輸出數(shù)據(jù)如圖4所示。

圖4 經(jīng)平滑處理后的試驗(yàn)數(shù)據(jù)Fig.4 Test data after smoothing

使用平滑處理后的第1組試驗(yàn)數(shù)據(jù)建立多項(xiàng)式模型，之后把第2組試驗(yàn)中相同陀螺的輸出數(shù)據(jù)輸入該模型。其模型擬合效果及補(bǔ)償后的結(jié)果如圖5所示。

圖5 最小二乘建模結(jié)果Fig.5 Modeling results of the least squares model

通過補(bǔ)償結(jié)果可見，以此方法擬合的模型，對(duì)相同陀螺的不同試驗(yàn)數(shù)據(jù)補(bǔ)償?shù)慕Y(jié)果并不理想，尤其是在-25℃下的補(bǔ)償結(jié)果。這種在不同溫度下啟動(dòng)的工作數(shù)據(jù)組合在一起波動(dòng)較大，而傳統(tǒng)多項(xiàng)式模型在此情況下擬合能力有限，無法滿足更高精度的要求。

之后使用二維插值方法建模，按2.2節(jié)中圖2所選擇的數(shù)據(jù)點(diǎn)作為插值的基礎(chǔ)模型。把另一組驗(yàn)證數(shù)據(jù)的溫度和溫變率作為插值點(diǎn)的自變量輸入模型，得到插值點(diǎn)的輸出值。此模型補(bǔ)償后結(jié)果如圖6所示。

圖6 二維插值建模結(jié)果Fig.6 Modeling results of the two-dimensional interpolation model

通過圖6與圖5對(duì)比可以看出，二維插值模型的補(bǔ)償效果明顯優(yōu)于最小二乘法的多項(xiàng)式模型擬合。慣組在-25℃啟動(dòng)時(shí)的數(shù)據(jù)也得到了較好的補(bǔ)償。

表1列出了陀螺輸出數(shù)據(jù)補(bǔ)償前與多項(xiàng)式模型和二維插值模型進(jìn)行溫度建模補(bǔ)償后的結(jié)果對(duì)比。

表1 補(bǔ)償前后陀螺零偏穩(wěn)定性

從表1中可以看出，通過二維插值建模對(duì)陀螺輸出數(shù)據(jù)進(jìn)行補(bǔ)償，很大程度上減小了因溫度造成的零偏漂移誤差，使陀螺的零偏穩(wěn)定性由補(bǔ)償前的0.0153(°)/h減小到0.0051(°)/h；使慣組在剛剛啟動(dòng)工作溫度并未達(dá)到穩(wěn)態(tài)勢(shì)時(shí)，輸出仍較為穩(wěn)定，繼而使慣導(dǎo)系統(tǒng)在上電短時(shí)間內(nèi)對(duì)準(zhǔn)導(dǎo)航精度得到提升。

4 結(jié)論

目前，光纖慣組中陀螺輸出因受溫度影響而產(chǎn)生溫度漂移誤差，很大程度上限制了其向高精度的進(jìn)一步發(fā)展,本文經(jīng)研究分析:

1)從光纖陀螺溫度漂移誤差產(chǎn)生的機(jī)理入手，分析了Shupe效應(yīng)的影響以及產(chǎn)生溫度漂移誤差的事實(shí)。

2)建立光纖陀螺的溫度誤差模型，對(duì)比了傳統(tǒng)的多項(xiàng)式擬合方法，得出的補(bǔ)償結(jié)果并不理想。

3)提出了一種使用二維插值法對(duì)慣組中光纖陀螺溫度漂移誤差進(jìn)行建模補(bǔ)償?shù)姆椒?。?jīng)過試驗(yàn)的設(shè)計(jì)與驗(yàn)證，有效地補(bǔ)償了慣組輸出因溫度相關(guān)信息而產(chǎn)生的誤差項(xiàng)，使陀螺的零偏穩(wěn)定性大大減小，且優(yōu)于傳統(tǒng)多項(xiàng)式模型，在工程中具有一定的應(yīng)用價(jià)值。