王喜龍，程遠(yuǎn)，吳汪洋

（1．海軍駐431廠軍代表室，遼寧葫蘆島125004;2．中國艦船研究院，北京 100192）

0 引言

計(jì)算機(jī)技術(shù)的飛速發(fā)展使得捷聯(lián)系統(tǒng)的設(shè)計(jì)與分析工作（特別是系統(tǒng)的誤差分析工作）可以首先在計(jì)算機(jī)進(jìn)行，在此基礎(chǔ)上再進(jìn)行系統(tǒng)的硬件（包括陀螺、加速度計(jì)與計(jì)算機(jī)等）及系統(tǒng)軟件（包括各種計(jì)算機(jī)算法及不同迭代周期的選擇等）的設(shè)計(jì)和選擇。但由于條件所限，在不能提供系統(tǒng)硬件的情況下為了進(jìn)行捷聯(lián)慣導(dǎo)系統(tǒng)的研究，開發(fā)慣導(dǎo)仿真器是必不可少的。

1 捷聯(lián)慣導(dǎo)系統(tǒng)仿真器的實(shí)現(xiàn)

本文針對飛行器用捷聯(lián)慣性導(dǎo)航系統(tǒng)的工作環(huán)境進(jìn)行仿真，其主要目的是仿真出符合實(shí)際慣導(dǎo)系統(tǒng)的誤差，捷聯(lián)慣導(dǎo)系統(tǒng)仿真器包括4個模塊:飛行軌跡發(fā)生器，慣性敏感元件仿真器（陀螺仿真器和加速度計(jì)仿真器），導(dǎo)航計(jì)算機(jī)仿真器和誤差處理器。仿真結(jié)構(gòu)如圖1所示［1］。

1.1 飛行軌跡發(fā)生器［2－4］

飛行軌跡數(shù)據(jù)生成是為了測試驗(yàn)證導(dǎo)航軟硬件系統(tǒng)實(shí)時計(jì)算性能而必需的輸入信號數(shù)據(jù)。對于捷聯(lián)慣導(dǎo)仿真系統(tǒng)，需要2類飛行軌跡數(shù)據(jù):飛行姿態(tài)數(shù)據(jù)（姿態(tài)角、姿態(tài)角轉(zhuǎn)動速率）和飛行軌跡數(shù)據(jù)（位置、速度、加速度）。實(shí)際飛行時，這2類數(shù)據(jù)的變化是有關(guān)聯(lián)的。如果根據(jù)飛機(jī)模型和飛行力學(xué)原理，用數(shù)學(xué)的方法來產(chǎn)生一個與實(shí)際完全相符的飛行軌跡數(shù)據(jù)文件，將會因?yàn)楦黠w行軌跡量之間的耦合關(guān)系，使數(shù)學(xué)模型的建立非常復(fù)雜，甚至無法建立起精確的數(shù)學(xué)模型。

圖1 數(shù)字仿真器原理方框圖Fig．1Principle block diagram of digital simulation

為了簡化建模，認(rèn)為2類數(shù)據(jù)按照各自獨(dú)立的規(guī)律變化，二者之間的變化沒有關(guān)系，這樣就可以很方便地用數(shù)學(xué)公式產(chǎn)生出任意時刻精確的理想飛行軌跡數(shù)據(jù)。測試用飛行軌跡數(shù)學(xué)模型的基礎(chǔ)是相互獨(dú)立的3個姿態(tài)角時間函數(shù)、3個位置時間函數(shù)。這樣計(jì)算生成的飛行軌跡雖然不是實(shí)際飛行軌跡，但可以很方便地得到任意時刻的精確測試用飛行軌跡數(shù)據(jù)，用來檢驗(yàn)捷聯(lián)慣導(dǎo)系統(tǒng)仿真器的正確性和精度。

1.1.1 飛行姿態(tài)數(shù)據(jù)生成

設(shè)定飛行姿態(tài)角是隨時間變化的余弦函數(shù)（也可以根據(jù)需要設(shè)定為其他函數(shù)）:

機(jī)體系相對于地理系的轉(zhuǎn)動角速率?在機(jī)體系中的投影為:

1.1.2 飛行軌跡數(shù)據(jù)生成

設(shè)定位置（緯度、經(jīng)度、高度）的變化規(guī)律為時間的正弦函數(shù)（也可以根據(jù)需要設(shè)定為其他函數(shù)），并計(jì)算出位置的一階和二階導(dǎo)數(shù):

計(jì)算出地球直角坐標(biāo)系中飛行載體位置所對應(yīng)的直角坐標(biāo)系R（x，y，z），并求直角坐標(biāo)的一階和二階導(dǎo)數(shù)，

其中，R′為導(dǎo)航系相對于地球的速度在地球坐標(biāo)系中的投影。

在地球系中的相對于地球加速度在導(dǎo)航系中的投影和在導(dǎo)航系中的相對于地球加速度在導(dǎo)航系中的投影為:

其中，R″為相對于地球加速度在地球系中的投影。

1.2 慣性敏感元件仿真器

慣性敏感元件仿真器給出陀螺與加速度計(jì)的仿真模型。當(dāng)只研究計(jì)算機(jī)的算法誤差時，則不考慮敏感元件的誤差;當(dāng)研究敏感元件的誤差時，其誤差也通過仿真器給出。

1.2.1 陀螺仿真器的數(shù)學(xué)模型［5－6］

職業(yè)教育主要培養(yǎng)面向企業(yè)一線具有較強(qiáng)實(shí)踐性和職業(yè)性的技能型人才，要求教師不僅具備專業(yè)理論知識，還需要有一定的實(shí)踐能力。近年來，隨著我國社會經(jīng)濟(jì)發(fā)展，產(chǎn)業(yè)結(jié)構(gòu)調(diào)整等新情況的出現(xiàn)，職業(yè)教育師資隊(duì)伍總量不足、培養(yǎng)機(jī)制不完善等，因此，有必要分析職院?！半p師型”教師隊(duì)伍建設(shè)存在的問題以及提出解決路徑。

理想角速率陀螺儀測量的是機(jī)體坐標(biāo)系（b系）相對與慣性坐標(biāo)系（i系）的轉(zhuǎn)動角速率在機(jī)體坐標(biāo)系中的投影?。從飛行軌跡數(shù)據(jù)中，可以得到機(jī)體坐標(biāo)系相對于導(dǎo)航坐標(biāo)系（n系）的轉(zhuǎn)動角速率在機(jī)體坐標(biāo)系中的投影。通過飛行軌跡數(shù)據(jù)中的水平速度、緯度、高度可以計(jì)算出導(dǎo)航坐標(biāo)系相對于慣性坐標(biāo)系的轉(zhuǎn)動角速率在導(dǎo)航坐標(biāo)系中的投影，通過姿態(tài)角算出從導(dǎo)航坐標(biāo)系到機(jī)體坐標(biāo)系之間的轉(zhuǎn)換矩陣乘上轉(zhuǎn)換矩陣得到相加，得到陀螺儀模型的理想輸出。

完整的陀螺模型包括理想輸出量、隨機(jī)誤差量（隨機(jī)常值誤差、白噪聲誤差、一階馬爾可夫過程誤差）、刻度因子誤差，用數(shù)學(xué)公式表示為:

式中:UG為完整陀螺模型的輸出;KG為刻度因子; δKG為刻度因子誤差;為陀螺模型理想輸出;εGb為陀螺隨機(jī)常數(shù)誤差;εGr為陀螺一階馬爾可夫過程誤差;wGg為陀螺白噪聲誤差。

1.2.2 加速度計(jì)仿真器的數(shù)學(xué)模型［7－8］

完整的加速度計(jì)模型包括理想輸出量、隨機(jī)誤差量（隨機(jī)常值誤差、白噪聲誤差、一階馬爾可夫過程誤差）、刻度因子誤差，可用數(shù)學(xué)公式表示為:

1.3 導(dǎo)航計(jì)算機(jī)仿真器

導(dǎo)航計(jì)算機(jī)的仿真可以根據(jù)系統(tǒng)對計(jì)算機(jī)的要求及現(xiàn)有的機(jī)載計(jì)算機(jī)的性能進(jìn)行。在仿真中可以對計(jì)算精度、計(jì)算時間、計(jì)算字長及單精度與雙精度計(jì)算等方面取得必要的數(shù)據(jù)，以便選擇合適的機(jī)載計(jì)算機(jī)及計(jì)算機(jī)算法。

1.4 誤差處理器

由導(dǎo)航計(jì)算機(jī)仿真器計(jì)算出的帶誤差的導(dǎo)航參數(shù)與飛行軌跡發(fā)生器產(chǎn)生的精確導(dǎo)航參數(shù)進(jìn)行比較，得出計(jì)算的位置、姿態(tài)、地速、高度等參數(shù)的導(dǎo)航誤差。可以通過誤差曲線圖分析仿真結(jié)果是否正確，完成全部的捷聯(lián)系統(tǒng)數(shù)字仿真器工作。

利用C語言設(shè)計(jì)的捷聯(lián)慣導(dǎo)系統(tǒng)仿真器已在微機(jī)上實(shí)現(xiàn)。仿真環(huán)境是在靜態(tài)條件下，選擇地理坐標(biāo)系為數(shù)字平臺。在對飛行軌跡數(shù)據(jù)仿真的算法公式中:θ，λ和h分別為緯度、經(jīng)度和高度;φ，γ和ψ分別為俯仰角、橫滾角和航向角。仿真過程中，式（1）和式（3）所涉及參數(shù)設(shè)定的缺省值為［9］:

飛機(jī)位置、姿態(tài)角的誤差曲線圖分別如圖2～圖4所示。

仿真結(jié)果表明，捷聯(lián)慣導(dǎo)仿真器工作正常，慣導(dǎo)的位置誤差、姿態(tài)誤差隨時間積累，長時間工作后慣導(dǎo)精度無法滿足導(dǎo)航或制導(dǎo)的精度要求。以上仿真誤差曲線是在陀螺儀與加速度計(jì)沒有誤差，也沒有初始對準(zhǔn)誤差的理想條件下生成的捷聯(lián)慣導(dǎo)仿真誤差曲線，產(chǎn)生誤差的主要因素是計(jì)算誤差?？梢钥闯觯S著時間的增長，誤差也在呈發(fā)散趨勢增長。因此，必須采取有效手段，對慣導(dǎo)系統(tǒng)進(jìn)行修正，以提高導(dǎo)航精度。

2 結(jié)語

本文所設(shè)計(jì)的捷聯(lián)慣導(dǎo)系統(tǒng)仿真器利用面向?qū)ο蟮乃枷脒M(jìn)行模塊分解，具有很強(qiáng)的擴(kuò)展性。該仿真器大大縮短了從理論方案到實(shí)際應(yīng)用之間的距離，為工程應(yīng)用提供了一種全面觀察捷聯(lián)系統(tǒng)在不同應(yīng)用條件下精度、性能的強(qiáng)有力的工具。結(jié)合不同的使用環(huán)境，可以對其中具體的模塊進(jìn)行相應(yīng)的重新設(shè)計(jì)。比如軌跡發(fā)生器仿真模塊、慣性敏感元件誤差模塊都可以依據(jù)系統(tǒng)具體的工作環(huán)境，對其中的各種參數(shù)或公式進(jìn)行重新選擇，從而達(dá)到更加真實(shí)的仿真效果。

［1］陳哲．捷聯(lián)慣導(dǎo)系統(tǒng)原理［M］．北京:宇航出版社，1986．206－211．

［2］李麗娟，王琪，羅超．VC環(huán)境下捷聯(lián)慣導(dǎo)模擬器的開發(fā)［J］．應(yīng)用科技，2006，33（3）:46－49．

［3］劉柱，高偉，曹潔．捷聯(lián)慣導(dǎo)系統(tǒng)仿真器的設(shè)計(jì)［J］．應(yīng)用科技，2004，31（7）:40－42．

［4］LEE J G，PACK H W．Multiposition alignment of strapdown inertial navigation system［J］．IEEE Transactions on Aerospace and Electronic Systems，1993，29（4）:1323－1328．

［5］GAI E．The century of inertial navigation technology［J］．Aerospace Conference Proceedings，2000，11（3）:59－60．

［6］LOVREN N，PIEPER J K A．Strapdown inertial navigation system for flat-earch model［J］．IEEE Transaction on Aerospace and Electronic Systems，1997，33（1）:214－224．

［7］吳小蘭．易于實(shí)現(xiàn)的捷聯(lián)式慣性導(dǎo)航系統(tǒng)仿真［J］．彈箭與制導(dǎo)學(xué)報，2006，26（4）:89－91．

［8］黃金明，王立文，鄭淑濤．飛行模擬器綜合導(dǎo)航系統(tǒng)建模與仿真［J］．系統(tǒng)仿真學(xué)報，2006，18（s2）:665－668．

［9］趙玉霞．捷聯(lián)慣導(dǎo)系統(tǒng)仿真算法的研究及其實(shí)現(xiàn)［D］．大連:大連理工大學(xué)，2005．