星敏感器是以恒星為參照系,以星空為工作對象的高精度空間姿態測量裝置,通過探測天球上不同位置的恒星并進行解算,為衛星、洲際戰略導彈、宇航飛船等航空航天飛行器提供準確的空間方位和基準,并且與慣性陀螺一樣都具有自主導航能力,具有重要的應用價值。
星敏感器的研究發展與應用已歷經半個多世紀,隨著新材料,新器件的出現和工藝技術的進步,精度提高,功耗減小,成本降低,應用領域日益廣泛的新型星敏感器不斷推出。因此,及時收集整理分析比較國外星敏感器的信息,有利于國內有關姿態測量控制技術的發展。
根據具體計算方法的不同,星敏感器中常用的姿態結算算法可以分為靜態確定性姿態解算算法和動態濾波估計算法。
確定性姿態解算算法
確定性姿態解算算法是指根據一組矢量觀測值,求出衛星本體坐標系相對于慣性坐標系的方向余弦矩陣。由確定性算法求解出來的結果具有明確的幾何和物理意義,而且只需通過一次測量就可以得到衛星的瞬時姿態。所以靜態確定性算法具有穩定性高、計算速度快、占用內存少等優點,也是目前星敏感器中主要使用的姿態確定算法。
但是直接求解Wahba問題比較困難,而且很難獲得最優解。1968年,Davenport提出了q-方法,利用四元數參數化姿態矩陣,將Wahba問題轉化為K矩陣的最大特征值求解問題,極大的推動了靜態確定性姿態解算算法的發展。后來研究者又提出了TRIAD算法、Euler-q算法、QUEST算法和FORM算法等。Shuster指出利用TRIAD法求解時,在兩個觀測矢量測量精度不對等時將無法得到最優的結果[32]。Baritzhack提出利用兩次TRIAD法進行加權處理得到更為準確的姿態矩陣的方法[33]。Markley從FOAM方法入手推導了目前形式最為簡潔的雙矢量觀測情形下的閉合解形式,并且分析了該算法的方差。
QUEST算法是最小二乘意義下的最優四元數估計,該算法最早應用于1979年的MAGSAT任務,也是目前解決Wahba問題的最常用算法。Shuster提出了QUEST測量模型并證明了其對于小視場敏感器是比較精確的,并利用QUEST測量模型推導了TRIAD法和QUEST法的方差陣,從理論上證明了QUEST法優于TRIAD法。
動態濾波估計算法
在實際應用中,衛星軌道參數的測量誤差和星敏感器的安裝誤差都會給觀測矢量的測量帶來不確定性的誤差,而且這些誤差很難克服。為了滿足高精度姿態控制的需要,可以利用動態濾波估計方法來結算航天器的姿態信息。動態濾波估計方法是利用航天器根據姿態動力學方程,建立狀態方程和觀測方程,根據觀測信息得到一定準則下的最優估計方法得到航天器的真實姿態。相對于靜態確定性算法,動態濾波估計算法利用了更多的觀測信息,能提供統計意義下的最優解,可以避免不確定因素對航天器姿態的影響,提高姿態確定的精度。
擴展卡爾曼濾波(EKF)算法是航天器最常用的實時姿態確定算法[36]。根據姿態參數的選取不同和觀測量的不同形式,常見的有乘性擴展卡爾曼濾波(MEKF)和加性擴展卡爾曼濾波(AEKF),其中MEKF被廣泛應用于各種航天器姿態確定任務并且發展最為成熟。但是EKF魯棒性不強,易于發散,對于非線性特性較強的估計問題時常不能得到最優解。Julier和Uhlmann利用UT變換取代了局部線性化[37],提出Unscented卡爾曼濾波器(UKF),在初始誤差較大的情況下,依然有良好的收斂性,得到比較優的結果。不管是EKF還是UKF算法,都是基于系統的隨機部分服從高斯分布的假設條件上提出來的,在星敏感器的姿態動力學模型中存在有不確定性的力矩模型誤差情況下,其結果的有效性難以保證。
基于QUEST算法,Shuster提出了濾波QUEST算法[38],該算法利用姿態分布矩陣B的傳播實現了卡爾曼濾波的遞推處理功能。Bar-Itzhack也對QUEST算法進行了擴展,提出RE-QUEST算法[39],通過K矩陣的傳播來實現遞推功能。從本質上來說,濾波QUEST算法和RE-QUEST算法在數學上是等效的。但是由于這兩種算法的精度相對于EKF算法來說比較差,并沒有在工程中廣泛使用。近些年來,隨著新的濾波估計方法的出現,越來越多的算法用在了星敏感器動態姿態估計中,如粒子濾波算法[40]、高斯濾波算法[41]和多模自適應估計算法[42]等。動態濾波估計方法的優點是可以利用先驗的知識來逼近統計意義下的最優解,但是由于采用了非線性的方法,算法復雜性較高,在實際使用中尚有一定的困難。
無論是靜態確定性姿態解算方法還是動態濾波算法都在星敏感器產品中得到了實踐和應用。表 5中對比分析了星敏感器姿態解算算法的特點。對于有陀螺的姿態確定系統,目前最實用也最常用的是乘性擴展卡爾曼濾波(MEKF)的方法,對于無陀螺的姿態確定系統,可以采用QUEST或者預測濾波估計的方法。