摘要
本文主要介紹了機載捷聯慣導系統常用的姿態更新算法——四元素法,并重點介紹了利用四元素法進行姿態更新的一般過程。
關鍵詞:四元素法,連貫導,姿態
1 引言
所謂姿態更新是指將運載體上慣性單元的輸出,實時轉換成運載體的姿態。這里的姿態通常指機體坐標系(b系)相對于導航坐標系(n系)的角位置。
對于一個姿態求解系統而言,其內部的姿態更新算法,無疑是其整個系統的核心,如何根據系統輸入以及應用場合的特點選用一種速度、精度都能充分滿足要求的算法是系統設計人員必須認真考慮的一件事情。目前,常用姿態更新算法主要有歐拉角算法、方向余弦法、四元素法以及旋轉矢量法。
歐拉角算法通過求解歐拉角微分方程直接計算航向角、俯仰角和橫滾角,但由于其微分方程特點,當俯仰角接近90o度時方程出現退化,故其只適用于水平姿態變化不大的情況,而不適用于全姿態的姿態確定。
方向余弦法避免了退化問題,可計算運載體全姿態,但是由于算法計算量大,實時計算困難,工程中很少采用。
四元素法計算量比余弦法小,且算法簡單,易于操作,是比較實用的工程方法。但其對有限轉動引起的不可交換誤差的補償程度不夠,所以只適用于低動態運載體的姿態解算。對于高動態運載體,其算法漂移十分嚴重。
旋轉矢量法根據運載體角速度擬合方式,分為單子樣算法(常數擬合),二子樣算法(直線擬合),三子樣算法(拋物線擬合)。因此可以根據需要采用合適的多子樣算法實現對不可交換誤差做有效補償。旋轉矢量法精度通常優于四元素法,但是通常計算量較四元素法更大,其較四元素法更適合角機動頻繁或者存在嚴重角振動的場合。值得一提的是,當旋轉矢量法采用單子樣計算姿態時,就等同于四元素法。
由于民航飛機以及直升機正常工作狀態常處于低動態,故在此類飛機的捷聯慣導系統中常采用四元素法作為姿態更新算法,本文也將主要針對應用于此類系統的四元素法進行介紹。
2 四元素定義
所謂四元素,即由四個元構成的數:

(2.1)
其中q0、q1、q2、q3是實數,i、j、k既是互相正交的單位向量,又是虛單位,因此四元素既可以看作是四維空間中的一個向量,又可看作一個超復數。
3 四元素性質
在計算運載體姿態時,當只關心機體坐標系相對于導航坐標系的的角位置時,可以認為機體坐標系(b系)是由導航坐標系(n系)經過無中間過程一次性等效旋轉形成的,四元素Q包含了這種等效旋轉的全部信息。其姿態變換公式見下式:

(3.1)
式(3.1)中,為b系中矢量,為b系相對于n系旋轉后,在n系中的矢量投影結果。式中坐標旋轉矩陣具體形式如下:

(3.2)
當已知運載體的航向角Y、俯仰角q、橫滾角g時,機體坐標系相對于導航坐標系的一次性等效旋轉矩陣還可以表示成如下形式:

(3.3)
記

,故當已知姿態變換矩陣時,根據式(3.2)、式(3.3)可以求出當前姿態的四元素以及歐拉角,其中四元素求解方程如下:

(3.4)
式中q0、q1、q2、q3的符號可以按下式確定,其中的值可以任意假設:

(3.5)
由旋轉矩陣式(3.3)得到的歐拉角求解公式如下:

(3.6)
式(3.6)中的航向角和橫滾角的真值按表1和表2確定。
表1
表3.2
4 四元素微分方程的畢卡求解法
捷聯陀螺的輸出通常是角速度,因此,為了計算運載體姿態,我們需要引入四元素微分方程。引入微分方程的好處是,根據上一時刻姿態四元素,通過定時采樣機體坐標系的三軸角增量即可求出新的姿態四元素(即機體坐標系相對于導航坐標系的四元素)。
四元素微分方程見下式:

(4.1)
式(4.1)中:

(4.2)
其中

分別為定時采樣時三個軸的角增量,I為單位矩陣。因此由時刻陀螺儀定時采樣的角增量,結合式(4.1),即可求出

時刻的四元素

。
5 四元素姿態更新流程
利用四元素微分方程進行姿態更新前,需要知道系統的初始四元素,求取初始四元素的方法以捷聯慣導為例,通常是利用捷聯慣導的加速度計在無加速狀態時的輸出結果計算出系統初始的歐拉角,再利用式(3.3)求出初始的姿態旋轉矩陣

,再將

的矩陣元素代入式(3.4)、式(3.5)即可求出初始四元素。
求出初始四元素后,定時采樣系統陀螺儀的角速度即可得到等時間間隔內的系統角增量,代統姿態。
6 結論
本文僅對基于四元素法更新的一般情況進行了討論,實際應用中還必須充分考慮到四元素法由于算法本身問題造成的誤差,例如利用四元素法進行姿態更新時,陀螺儀的角增量為定時采樣的,然而實際陀螺儀的輸出并非完全線性關系,即不為常數,同時四元素法對旋轉軸的旋轉次序是敏感的,即當

的結果相同時,不同的

輸出次序將導致不同的四元素更新結果。而在利用四元素法進行姿態更新時,近似的認為其在等時間內三軸同時產生角增量。這勢必產生前文所述的不可交換誤差。縮短采樣時間間隔可以有效降低近似誤差,然而卻增加了計算量,對系統提出了更高要求,因此實際使用中必須根據應用場合對四元素法進行適當修正,才能得到滿意的結果。