白丝美女被狂躁免费视频网站,500av导航大全精品,yw.193.cnc爆乳尤物未满,97se亚洲综合色区,аⅴ天堂中文在线网官网

首頁(yè) / 專(zhuān)利庫(kù) / 天文學(xué) / 星系 / 銀河系 / 太陽(yáng)系 / 一種巡航段自主導(dǎo)航位置偏差修正方法

一種巡航段自主導(dǎo)航位置偏差修正方法

閱讀:318發(fā)布:2020-09-20

專(zhuān)利匯可以提供一種巡航段自主導(dǎo)航位置偏差修正方法專(zhuān)利檢索,專(zhuān)利查詢(xún),專(zhuān)利分析的服務(wù)。并且本 發(fā)明 公開(kāi)的一種巡航段自主導(dǎo)航 位置 偏差修正方法,屬于深空探測(cè)技術(shù)領(lǐng)域。本發(fā)明建立探測(cè)器的動(dòng) 力 學(xué)模型并對(duì)其進(jìn)行線性化處理;建立自主導(dǎo)航測(cè)量模型,光學(xué)導(dǎo)航相機(jī)獲取可見(jiàn)行星的灰度圖像,通過(guò)星載 圖像處理 完成可見(jiàn)行星質(zhì)心的識(shí)別和 跟蹤 ;確定探測(cè)器的初始軌道;通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置;通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正可見(jiàn)行星質(zhì)心 角 位置觀測(cè)信息,獲得更精確的觀測(cè)模型;根據(jù)動(dòng)力學(xué)模型及測(cè)量模型,基于非線性系統(tǒng)濾波 算法 解算探測(cè)器實(shí)時(shí)狀態(tài),實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。本發(fā)明能夠提高巡航段自主導(dǎo)航方法估計(jì) 精度 、濾波收斂速度,實(shí)現(xiàn)快速精確估計(jì),為深空探測(cè)巡航段導(dǎo)航方案設(shè)計(jì)提供支持。,下面是一種巡航段自主導(dǎo)航位置偏差修正方法專(zhuān)利的具體信息內(nèi)容。

1.一種巡航段自主導(dǎo)航位置偏差修正方法,其特征在于:包括如下步驟,
步驟1:建立巡航段探測(cè)器的動(dòng)學(xué)模型;
探測(cè)器在太陽(yáng)系中運(yùn)動(dòng),主要受到太陽(yáng)引力場(chǎng)的作用,其動(dòng)力學(xué)滿(mǎn)足開(kāi)普勒二體方程;
為了盡可能精確地建立探測(cè)器在巡航段的動(dòng)力學(xué)模型,模型攝動(dòng)項(xiàng)考慮了大行星的引力影響,太陽(yáng)光壓攝動(dòng)影響和探測(cè)器推力因素,故探測(cè)器在日心慣性坐標(biāo)系下的動(dòng)力學(xué)方程表達(dá)式如公式(1)所示:
式中r和v分別表示探測(cè)器的位置和速度矢量;為第i個(gè)攝動(dòng)行星的位置矢量;為第i個(gè)攝動(dòng)行星相對(duì)探測(cè)器的位置矢量,即 μs為太陽(yáng)引力常數(shù),μi為第i個(gè)攝動(dòng)行星的引力常數(shù);nt為攝動(dòng)行星的個(gè)數(shù);CR為探測(cè)器表面的反射系數(shù),Ssrp為太陽(yáng)輻射光壓因子,m為探測(cè)器的質(zhì)量,k為推力系數(shù),最后一項(xiàng)a表示其他未建模加速度;
選取探測(cè)器位置和速度作為狀態(tài)變量,則 則探測(cè)
器的狀態(tài)方程如公式(2)所示:
步驟2:建立巡航段自主導(dǎo)航測(cè)量模型;
相機(jī)成像的過(guò)程采用小孔成像的模型,定義可見(jiàn)行星的質(zhì)心f1在相機(jī)坐標(biāo)系下的位置坐標(biāo)為rp=[xc?yc?zc]T,則其在相機(jī)像平面內(nèi)的像原像素坐標(biāo)如公式(3)所示:
其中:f為相機(jī)焦距,zc為目標(biāo)點(diǎn)沿著相機(jī)基準(zhǔn)線到相機(jī)成像平面的距離;
相機(jī)系與探測(cè)器本體系重合,考慮到相機(jī)測(cè)量過(guò)程中的姿態(tài)偏差,即定義xc,yc方向的姿態(tài)偏差分別為θ1,θ2,則相機(jī)測(cè)量過(guò)程中的隨機(jī)旋轉(zhuǎn)將對(duì)特征點(diǎn)的位置測(cè)量產(chǎn)生影響,故實(shí)際的位置坐標(biāo)如公式(4)所示:
因?yàn)樽藨B(tài)偏差較小,故簡(jiǎn)化為公式(5)
忽略分母影響,則簡(jiǎn)化為公式(6):
選擇可見(jiàn)行星質(zhì)心的角位置作為相機(jī)的觀測(cè)量,即測(cè)量矢量Z包括方位角 和俯仰角如公式(7)所示:
其中: 和υφ分別代表俯仰角和方位角的測(cè)量噪聲,可見(jiàn)行星的質(zhì)心位置矢量為R,則俯仰角和方位角計(jì)算如公式(8)所示:
步驟3:確定探測(cè)器初始軌道;
Nobs次觀測(cè)的初始軌道確定問(wèn)題寫(xiě)成公式(9)所示形式:
其中:rk,Rk分別代表第k次航天器和可見(jiàn)行星在日心慣性系下的位置矢量;未知的距離ρk是待求量,表示探測(cè)器到可見(jiàn)行星的單位視線矢量;將待求量ρk表示成高階多項(xiàng)式的形式,通過(guò)擬合方法求解滿(mǎn)足公式(9)約束的待求量ρk的近似解,即確定探測(cè)器初始軌道;
步驟4:光行差效應(yīng)校正;
通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置矢量;通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正步驟2中如公式(7)所示的角位置觀測(cè)信息,從而獲得校正后更精確的真實(shí)觀測(cè)模型,即完成光行差效應(yīng)校正;
步驟5:根據(jù)步驟1中的探測(cè)器的動(dòng)力學(xué)模型、步驟3中確定的探測(cè)器初始軌道及步驟4修正后真實(shí)的測(cè)量模型,基于非線性系統(tǒng)濾波算法解算探測(cè)器實(shí)時(shí)狀態(tài)信息,實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航;
根據(jù)步驟1得到的巡航段探測(cè)器動(dòng)力學(xué)模型、步驟3中確定的探測(cè)器初始軌道及步驟4修正后真實(shí)的測(cè)量模型,通過(guò)導(dǎo)航濾波對(duì)探測(cè)器的狀態(tài)進(jìn)行估計(jì);由于動(dòng)力學(xué)模型及測(cè)量模型均呈現(xiàn)非線性,故選用非線性濾波器,解算探測(cè)器實(shí)時(shí)狀態(tài)信息,實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。
2.如權(quán)利要求1所述的一種巡航段自主導(dǎo)航位置偏差修正方法,其特征在于:步驟4具體實(shí)現(xiàn)方法如下,
步驟4.1:通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置;
光行差校正指由于光到達(dá)探測(cè)器的過(guò)程中天體的運(yùn)動(dòng),導(dǎo)致天體真實(shí)位置和觀測(cè)位置之間產(chǎn)生的位移,對(duì)其進(jìn)行校正;針對(duì)導(dǎo)致天體真實(shí)位置和觀測(cè)位置之間產(chǎn)生的位移進(jìn)行校正問(wèn)題,需要對(duì)可見(jiàn)行星位置矢量進(jìn)行進(jìn)一步修正;定義觀測(cè)時(shí)間為tk,k=1,2...,Nobs,各個(gè)測(cè)量方向的單位矢量為 因?yàn)樘綔y(cè)器不能確定自己相對(duì)于太陽(yáng)和地
球的距離,而僅已知觀測(cè)可見(jiàn)行星的位置矢量Rk,和相對(duì)可見(jiàn)行星的觀測(cè)方向的單位矢量因此基于時(shí)間tk-δtk對(duì)可見(jiàn)行星真實(shí)位置矢量進(jìn)行求解,其中δtk是光從可見(jiàn)行星到達(dá)探測(cè)器所需的時(shí)間;首先確定δtk,如公式(18)所示:
其中c代表光速;對(duì)行星的位置矢量進(jìn)行校正,如公式(19)所示:
Rupdated=R(tk-δtk)???????????(19)
可見(jiàn)行星位置矢量的更新由式(18)和(19)確定,即完成可見(jiàn)行星位置矢量的校正;
步驟4.2:通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正步驟2中如公式(7)所示的角位置觀測(cè)信息,得到更精確的測(cè)量模型;
定義探測(cè)器速度方向和觀測(cè)方向 之間的夾角是θobs,真實(shí)方向和觀測(cè)方向之間的畸變角為ε,則修正后滿(mǎn)足精度要求的真實(shí)天體觀測(cè)方向 如公式(20)所示:
其中:是探測(cè)器速度的單位矢量,θtrue是探測(cè)器速度和真實(shí)觀測(cè)方向之間的夾角,且θtrue=θobs+ε;不同時(shí)刻的畸變角ε表示成公式(21)的形式,
其中:c和v分別代表光速大小和探測(cè)器速度大??;
觀測(cè)方向和真實(shí)方向之間的關(guān)系如公式(22)所示,
可見(jiàn)行星和探測(cè)器之間視線方向的夾角和探測(cè)器的速度矢量是已知的,因此真實(shí)觀測(cè)方向可以通過(guò)式(22)求解,結(jié)合公式(7)從而得到修正后真實(shí)的測(cè)量模型,至此完成光行差效應(yīng)校正。
3.如權(quán)利要求1或2所述的一種巡航段自主導(dǎo)航位置偏差修正方法,其特征在于:
為進(jìn)一步提高步驟3中探測(cè)器初始軌道確定的求解效率和求解精度,步驟3中將待求量ρk表示成三階多項(xiàng)式的形式,通過(guò)最小二乘擬合方法求解滿(mǎn)足公式(9)約束的待求量ρk的近似解,選如下方法實(shí)現(xiàn):
待求量ρk在日心慣性系下三軸分量如公式(10)所示:
考慮三階多項(xiàng)式的形式,將公式(10)表示成如公式(11)所示的形式:
公式(11)也應(yīng)該滿(mǎn)足公式(9)中的幾何約束,所以,
tk表示第k次的測(cè)量時(shí)間;將公式(12)中等式左邊移到右邊,定義殘差如公式(13)所示:
進(jìn)一步將公式(13)表示成公式(14)所示形式:
為了得到最好的擬合,需要調(diào)整多項(xiàng)式擬合的系數(shù)ai,bi,和ci,使得殘差φx,φy,和φz最小,所以有,
因此,通過(guò)擬合方法確定待求多項(xiàng)式中的系數(shù),為進(jìn)一步提高求解效率,選最小二乘法確定待求多項(xiàng)式中的系數(shù),如公式(16):
{a}={a1?a2?a3?a4}T,={b1?b2?b3?b4}T,{c}={c1?c2?c3?c4}T,3×3的矩陣T如公式(17):
至此,如公式(11)所示的三階多項(xiàng)式的系數(shù)已經(jīng)全部求解,根據(jù)系數(shù)能夠得到滿(mǎn)足精度要求的待求量ρk的近似解,即完成探測(cè)器初始軌道確定。
4.如權(quán)利要求1或2所述的一種巡航段自主導(dǎo)航位置偏差修正方法,其特征在于:步驟5所述的非線性濾波器擴(kuò)展卡爾曼濾波(EKF)、無(wú)跡卡爾曼濾波(UKF)、模型預(yù)測(cè)濾波器(MPF)或粒子濾波(PF)。
5.如權(quán)利要求4所述的一種巡航段自主導(dǎo)航位置偏差修正方法,其特征在于:為進(jìn)一步提高導(dǎo)航濾波精度及收斂速度,步驟5所述的非線性濾波器選擴(kuò)展卡爾曼濾波(EKF)。
6.一種巡航段自主導(dǎo)航位置偏差修正方法,其特征在于:建立巡航段探測(cè)器的動(dòng)力學(xué)模型并對(duì)動(dòng)力學(xué)模型進(jìn)行線性化處理;建立自主導(dǎo)航測(cè)量模型,光學(xué)導(dǎo)航相機(jī)獲取可見(jiàn)行星的灰度圖像,通過(guò)星載圖像處理軟件完成可見(jiàn)行星質(zhì)心的識(shí)別和跟蹤;通過(guò)基于空間的初始軌道確定算法確定探測(cè)器的初始軌道,所述的探測(cè)器的初始軌道包括探測(cè)器的位置矢量和速度矢量;通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置;通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正可見(jiàn)行星質(zhì)心角位置觀測(cè)信息,從而獲得更精確的觀測(cè)模型,即完成位置偏差修正,即光行差效應(yīng)校正;根據(jù)探測(cè)器巡航段的動(dòng)力學(xué)模型及測(cè)量模型,基于非線性系統(tǒng)濾波算法解算探測(cè)器實(shí)時(shí)狀態(tài),實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。

說(shuō)明書(shū)全文

一種巡航段自主導(dǎo)航位置偏差修正方法

技術(shù)領(lǐng)域

[0001] 本發(fā)明涉及一種巡航段自主導(dǎo)航位置偏差修正方法,屬于深空探測(cè)技術(shù)領(lǐng)域。

背景技術(shù)

[0002] 深空探測(cè)是人類(lèi)了解宇宙太陽(yáng)系的形成和演化、探索生命起源的主要途徑,巡航段是深空探測(cè)最長(zhǎng)的軌道段,該階段具有動(dòng)學(xué)模型和星歷數(shù)據(jù)不精確,飛行時(shí)間長(zhǎng)等特點(diǎn)。巡航段飛行任務(wù)中的姿態(tài)信息由星敏感器確定,位置估計(jì)仍然依賴(lài)于地面通信站。然而在此飛行過(guò)程中,如果完全依靠地面測(cè)控站進(jìn)行探測(cè)器導(dǎo)航,將嚴(yán)重增加地面站的負(fù)擔(dān);另外,當(dāng)通信系統(tǒng)發(fā)生臨時(shí)性或永久性故障失效時(shí),自主導(dǎo)航性能的優(yōu)劣對(duì)探測(cè)任務(wù)的成敗起著至關(guān)重要的作用。此外,深空探測(cè)巡航段的導(dǎo)航精度直接關(guān)系到后續(xù)目標(biāo)天體的捕獲、接近、繞飛等探測(cè)任務(wù)的成敗,因此巡航段高精度的自主導(dǎo)航方法研究有重要意義。
[0003] 基于以上特點(diǎn)和要求,參考深空1號(hào)(Deep?Space?1)任務(wù),星塵(Stardust)任務(wù),智能1號(hào)(SMART-1)任務(wù)等眾多深空探測(cè)任務(wù),光學(xué)敏感器成為巡航段探測(cè)器自主確定位置和速度的主要敏感器。探測(cè)器在星際巡航段一般采用光學(xué)敏感器觀測(cè)天體質(zhì)心的自主導(dǎo)航方法,該方法能夠?qū)崿F(xiàn)探測(cè)器位置和速度的自主估計(jì)。然而星際巡航過(guò)程中,由于探測(cè)器和被觀測(cè)天體在軌運(yùn)行過(guò)程中相距很遠(yuǎn)且存在相對(duì)運(yùn)動(dòng),光從被觀測(cè)天體到達(dá)探測(cè)器的過(guò)程中,會(huì)導(dǎo)致探測(cè)器觀測(cè)天體的方向(視方向)與其真實(shí)方向之間產(chǎn)生視向量偏差,即天體真實(shí)位置和觀測(cè)位置之間產(chǎn)生位置矢量偏差,故導(dǎo)航誤差隨著探測(cè)器和被觀測(cè)天體之間距離的增加而增大,因此需要對(duì)此偏差進(jìn)行校正。

發(fā)明內(nèi)容

[0004] 針對(duì)現(xiàn)有技術(shù)中利用小行星測(cè)量的巡航段自主導(dǎo)航方法中,小行星數(shù)量有限、受篩選準(zhǔn)則影響、星歷誤差大、以及探測(cè)器和被觀測(cè)天體之間存在相對(duì)運(yùn)動(dòng)的問(wèn)題,本發(fā)明公開(kāi)的一種巡航段自主導(dǎo)航位置偏差修正方法,要解決的技術(shù)問(wèn)題是提高巡航段自主導(dǎo)航方法的估計(jì)精度、濾波收斂速度,實(shí)現(xiàn)探測(cè)器狀態(tài)的快速精確估計(jì),為深空探測(cè)巡航段導(dǎo)航方案設(shè)計(jì)提供技術(shù)支持。
[0005] 本發(fā)明的目的是通過(guò)下述技術(shù)方案實(shí)現(xiàn)的。
[0006] 本發(fā)明公開(kāi)的一種巡航段自主導(dǎo)航位置偏差修正方法,建立巡航段探測(cè)器的動(dòng)力學(xué)模型并對(duì)動(dòng)力學(xué)模型進(jìn)行線性化處理。建立自主導(dǎo)航測(cè)量模型,光學(xué)導(dǎo)航相機(jī)獲取可見(jiàn)行星的灰度圖像,通過(guò)星載圖像處理軟件完成可見(jiàn)行星質(zhì)心的識(shí)別和跟蹤。通過(guò)基于空間的初始軌道確定算法確定探測(cè)器的初始軌道,所述的探測(cè)器的初始軌道包括探測(cè)器的位置矢量和速度矢量。通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置。通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正可見(jiàn)行星質(zhì)心位置觀測(cè)信息,從而獲得更精確的觀測(cè)模型,即完成位置偏差修正,即光行差效應(yīng)校正。根據(jù)探測(cè)器巡航段的動(dòng)力學(xué)模型及測(cè)量模型,基于非線性系統(tǒng)濾波算法解算探測(cè)器實(shí)時(shí)狀態(tài),實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。
[0007] 本發(fā)明公開(kāi)的一種巡航段自主導(dǎo)航位置偏差修正方法,包括如下步驟:
[0008] 步驟1:建立巡航段探測(cè)器的動(dòng)力學(xué)模型。
[0009] 探測(cè)器在太陽(yáng)系中運(yùn)動(dòng),主要受到太陽(yáng)引力場(chǎng)的作用,其動(dòng)力學(xué)滿(mǎn)足開(kāi)普勒二體方程。為了盡可能精確地建立探測(cè)器在巡航段的動(dòng)力學(xué)模型,模型攝動(dòng)項(xiàng)考慮了大行星的引力影響,太陽(yáng)光壓攝動(dòng)影響和探測(cè)器推力因素,故探測(cè)器在日心慣性坐標(biāo)系下的動(dòng)力學(xué)方程表達(dá)式如公式(1)所示:
[0010]
[0011] 式中r和v分別表示探測(cè)器的位置和速度矢量;為第i個(gè)攝動(dòng)行星的位置矢量;為第i個(gè)攝動(dòng)行星相對(duì)探測(cè)器的位置矢量,即 μs為太陽(yáng)引力常數(shù),μi為第i個(gè)攝動(dòng)行星的引力常數(shù);nt為攝動(dòng)行星的個(gè)數(shù)。CR為探測(cè)器表面的反射系數(shù),Ssrp為太陽(yáng)輻射光壓因子,m為探測(cè)器的質(zhì)量,k為推力系數(shù),最后一項(xiàng)a表示其他未建模加速度。
[0012] 選取探測(cè)器位置和速度作為狀態(tài)變量,則則探測(cè)器的狀態(tài)方程如公式(2)所示:
[0013]
[0014] 步驟2:建立巡航段自主導(dǎo)航測(cè)量模型。
[0015] 相機(jī)成像的過(guò)程采用小孔成像的模型,定義可見(jiàn)行星的質(zhì)心f1在相機(jī)坐標(biāo)系下的位置坐標(biāo)為rp=[xc?yc?zc]T,則其在相機(jī)像平面內(nèi)的像原像素坐標(biāo)如公式(3)所示:
[0016]
[0017] 其中:f為相機(jī)焦距,zc為目標(biāo)點(diǎn)沿著相機(jī)基準(zhǔn)線到相機(jī)成像平面的距離。
[0018] 相機(jī)系與探測(cè)器本體系重合,考慮到相機(jī)測(cè)量過(guò)程中的姿態(tài)偏差,即定義xc,yc方向的姿態(tài)偏差分別為θ1,θ2,則相機(jī)測(cè)量過(guò)程中的隨機(jī)旋轉(zhuǎn)將對(duì)特征點(diǎn)的位置測(cè)量產(chǎn)生影響,故實(shí)際的位置坐標(biāo)如公式(4)所示:
[0019]
[0020] 因?yàn)樽藨B(tài)角偏差較小,故簡(jiǎn)化為公式(5)
[0021]
[0022] 忽略分母影響,則簡(jiǎn)化為公式(6):
[0023]
[0024] 選擇可見(jiàn)行星質(zhì)心的角位置作為相機(jī)的觀測(cè)量,即測(cè)量矢量Z包括方位角 和俯仰角 如公式(7)所示:
[0025]
[0026] 其中: 和υφ分別代表俯仰角和方位角的測(cè)量噪聲,可見(jiàn)行星的質(zhì)心位置矢量為R,則俯仰角和方位角計(jì)算如公式(8)所示:
[0027]
[0028] 步驟3:確定探測(cè)器初始軌道。
[0029] Nobs次觀測(cè)的初始軌道確定問(wèn)題寫(xiě)成公式(9)所示形式:
[0030]
[0031] 其中:rk,Rk分別代表第k次航天器和可見(jiàn)行星在日心慣性系下的位置矢量。未知的距離ρk是待求量, 表示探測(cè)器到可見(jiàn)行星的單位視線矢量。將待求量ρk表示成高階多項(xiàng)式的形式,通過(guò)擬合方法求解滿(mǎn)足公式(9)約束的待求量ρk的近似解,即確定探測(cè)器初始軌道。
[0032] 為進(jìn)一步提高步驟3中探測(cè)器初始軌道確定的求解效率和求解精度,步驟3中將待求量ρk表示成三階多項(xiàng)式的形式,通過(guò)最小二乘擬合方法求解滿(mǎn)足公式(9)約束的待求量ρk的近似解,優(yōu)選如下方法實(shí)現(xiàn):
[0033] 待求量ρk在日心慣性系下三軸分量如公式(10)所示:
[0034]
[0035] 考慮三階多項(xiàng)式的形式,將公式(10)表示成如公式(11)所示的形式:
[0036]
[0037] 公式(11)也應(yīng)該滿(mǎn)足公式(9)中的幾何約束,所以,
[0038]
[0039] tk表示第k次的測(cè)量時(shí)間。將公式(12)中等式左邊移到右邊,定義殘差如公式(13)所示:
[0040]
[0041] 進(jìn)一步將公式(13)表示成公式(14)所示形式:
[0042]
[0043] 為了得到最好的擬合,需要調(diào)整多項(xiàng)式擬合的系數(shù)ai,bi,和ci,使得殘差φx,φy,和φz最小,所以有,
[0044]
[0045] 因此,通過(guò)擬合方法確定待求多項(xiàng)式中的系數(shù),為進(jìn)一步提高求解效率,選最小二乘法確定待求多項(xiàng)式中的系數(shù),如公式(16):
[0046]
[0047] {a}={a1?a2?a3?a4}T,={b1?b2?b3?b4}T,{c}={c1?c2?c3?c4}T,3×3的矩陣T如公式(17):
[0048]
[0049] 至此,如公式(11)所示的三階多項(xiàng)式的系數(shù)已經(jīng)全部求解,根據(jù)系數(shù)能夠得到滿(mǎn)足精度要求的待求量ρk的近似解,即完成探測(cè)器初始軌道確定。
[0050] 步驟4:光行差效應(yīng)校正。
[0051] 通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置矢量。通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正步驟2中如公式(7)所示的角位置觀測(cè)信息,從而獲得更精確的觀測(cè)模型,即完成光行差效應(yīng)校正。
[0052] 步驟4具體實(shí)現(xiàn)方法如下:
[0053] 步驟4.1:通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置。
[0054] 光行差校正指由于光到達(dá)探測(cè)器的過(guò)程中天體的運(yùn)動(dòng),導(dǎo)致天體真實(shí)位置和觀測(cè)位置之間產(chǎn)生的位移,對(duì)其進(jìn)行校正。針對(duì)導(dǎo)致天體真實(shí)位置和觀測(cè)位置之間產(chǎn)生的位移進(jìn)行校正問(wèn)題,需要對(duì)可見(jiàn)行星位置矢量進(jìn)行進(jìn)一步修正。定義觀測(cè)時(shí)間為tk,k=1,2...,Nobs,各個(gè)測(cè)量方向的單位矢量為 k=1,2...,Nobs。因?yàn)樘綔y(cè)器不能確定自己相對(duì)于太陽(yáng)和地球的距離,而僅已知觀測(cè)可見(jiàn)行星的位置矢量Rk,和相對(duì)可見(jiàn)行星的觀測(cè)方向的單位矢量 因此基于時(shí)間tk-δtk對(duì)可見(jiàn)行星真實(shí)位置矢量進(jìn)行求解,其中δtk是光從可見(jiàn)行星到達(dá)探測(cè)器所需的時(shí)間。首先確定δtk,如公式(18)所示:
[0055]
[0056] 其中c代表光速。對(duì)行星的位置矢量進(jìn)行校正,如公式(19)所示:
[0057] Rupdated=R(tk-δtk)???(19)
[0058] 可見(jiàn)行星位置矢量的更新由公式(18)和(19)確定,即完成可見(jiàn)行星位置矢量的校正。
[0059] 步驟4.2:通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正步驟2中如公式(7)所示的角位置觀測(cè)信息,得到更精確的測(cè)量模型。
[0060] 定義探測(cè)器速度方向和觀測(cè)方向 之間的夾角是θobs,真實(shí)方向和觀測(cè)方向之間的畸變角為ε,則修正后滿(mǎn)足精度要求的真實(shí)天體觀測(cè)方向 如公式(20)所示:
[0061]
[0062] 其中:是探測(cè)器速度的單位矢量,θtrue是探測(cè)器速度和真實(shí)觀測(cè)方向之間的夾角,且θtrue=θobs+ε。不同時(shí)刻的畸變角ε表示成公式(21)的形式,
[0063]
[0064] 其中:c和v分別代表光速大小和探測(cè)器速度大小。
[0065] 觀測(cè)方向和真實(shí)方向之間的關(guān)系如公式(22)所示,
[0066]
[0067] 可見(jiàn)行星和探測(cè)器之間視線方向的夾角和探測(cè)器的速度矢量是已知的,因此真實(shí)觀測(cè)方向可以通過(guò)公式(22)求解,結(jié)合公式(7)從而得到修正后真實(shí)的測(cè)量模型,至此完成光行差效應(yīng)校正。
[0068] 步驟5:根據(jù)步驟1中的探測(cè)器的動(dòng)力學(xué)模型、步驟3中確定的探測(cè)器初始軌道及步驟4修正后真實(shí)的測(cè)量模型,基于非線性系統(tǒng)濾波算法解算探測(cè)器實(shí)時(shí)狀態(tài)信息,實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。
[0069] 根據(jù)步驟1得到的巡航段探測(cè)器動(dòng)力學(xué)模型、步驟3中確定的探測(cè)器初始軌道及步驟4修正后真實(shí)的測(cè)量模型,通過(guò)導(dǎo)航濾波對(duì)探測(cè)器的狀態(tài)進(jìn)行估計(jì)。由于動(dòng)力學(xué)模型及測(cè)量模型均呈現(xiàn)非線性,故選用非線性濾波器,解算探測(cè)器實(shí)時(shí)狀態(tài)信息,實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。
[0070] 所述的非線性濾波器擴(kuò)展卡爾曼濾波(EKF)、無(wú)跡卡爾曼濾波(UKF)、模型預(yù)測(cè)濾波器(MPF)或粒子濾波(PF)等。
[0071] 為進(jìn)一步提高導(dǎo)航濾波精度及收斂速度,所述的非線性濾波器優(yōu)選擴(kuò)展卡爾曼濾波(EKF)。
[0072] 有益效果:
[0073] 1、現(xiàn)有技術(shù)中利用光學(xué)敏感器測(cè)量小行星的導(dǎo)航方法中,近地小行星數(shù)量有限且受到嚴(yán)格篩選準(zhǔn)則的約束,可用導(dǎo)航小行星數(shù)量難以滿(mǎn)足任務(wù)要求,加之地面站觀測(cè)得到的小行星星歷信息有較大誤差,導(dǎo)致該方法導(dǎo)航精度較低。本發(fā)明公開(kāi)的一種巡航段自主導(dǎo)航位置偏差修正方法,通過(guò)選取星歷信息更精確的可見(jiàn)行星為觀測(cè)天體,提高導(dǎo)航算法的估計(jì)精度,滿(mǎn)足巡航段自主導(dǎo)航的精度需求。
[0074] 2、針對(duì)巡航段,探測(cè)器和被觀測(cè)天體二者在軌運(yùn)行過(guò)程中相距很遠(yuǎn),再加上探測(cè)器和被觀測(cè)天體速度的存在,導(dǎo)致光從被觀測(cè)天體到達(dá)探測(cè)器的過(guò)程中,天體真實(shí)位置和觀測(cè)位置之間產(chǎn)生位置偏差(光行差效應(yīng)),本發(fā)明公開(kāi)的一種巡航段自主導(dǎo)航位置偏差修正方法針對(duì)該問(wèn)題對(duì)光行差效應(yīng)進(jìn)行校正,即通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置;通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正角位置觀測(cè)信息,從而獲得行星質(zhì)心更精確的測(cè)量模型,提高導(dǎo)航方法的估計(jì)精度。
[0075] 3、本發(fā)明公開(kāi)的一種巡航段自主導(dǎo)航位置偏差修正方法,采用非線性濾波器對(duì)探測(cè)器的狀態(tài)進(jìn)行估計(jì),能夠提高自主導(dǎo)航算法的精度及濾波收斂速度。附圖說(shuō)明
[0076] 圖1為巡航段自主導(dǎo)航位置偏差修正方法的流程圖;
[0077] 圖2為具體實(shí)施例中采用該自主導(dǎo)航方法時(shí),探測(cè)器在日心慣性坐標(biāo)系下的導(dǎo)航誤差曲線。(圖2a為探測(cè)器x方向位置導(dǎo)航誤差曲線、圖2b為探測(cè)器y方向位置導(dǎo)航誤差曲線、圖2c為探測(cè)器z方向位置導(dǎo)航誤差曲線、圖2d為探測(cè)器x方向速度導(dǎo)航誤差曲線、圖2e為探測(cè)器y方向速度導(dǎo)航誤差曲線、圖2f為探測(cè)器z方向速度導(dǎo)航誤差曲線、)具體實(shí)施方式
[0078] 為了更好的說(shuō)明本發(fā)明的目的和優(yōu)點(diǎn),下面結(jié)合附圖和實(shí)例對(duì)發(fā)明內(nèi)容做進(jìn)一步說(shuō)明。
[0079] 實(shí)施例1:
[0080] 本實(shí)例針對(duì)深空探測(cè)巡航段進(jìn)行仿真驗(yàn)證。利用光學(xué)相機(jī)測(cè)量探測(cè)器至可見(jiàn)行星質(zhì)心的視線信息,同時(shí)對(duì)光行差效應(yīng)進(jìn)行校正,然后利用擴(kuò)展卡爾曼濾波器(EKF),對(duì)探測(cè)器的位置、速度狀態(tài)進(jìn)行聯(lián)合估計(jì),實(shí)現(xiàn)高精度實(shí)時(shí)自主導(dǎo)航。
[0081] 本實(shí)例公開(kāi)的一種巡航段自主導(dǎo)航位置偏差修正方法,包括如下步驟:
[0082] 步驟1:建立巡航段探測(cè)器的動(dòng)力學(xué)模型。
[0083] 探測(cè)器在太陽(yáng)系中運(yùn)動(dòng),主要受到太陽(yáng)引力場(chǎng)的作用,其動(dòng)力學(xué)滿(mǎn)足開(kāi)普勒二體方程。為了盡可能精確地建立探測(cè)器在巡航段的動(dòng)力學(xué)模型,模型攝動(dòng)項(xiàng)考慮大行星的引力影響,太陽(yáng)光壓攝動(dòng)影響和探測(cè)器推力因素,故探測(cè)器在日心慣性坐標(biāo)系下的動(dòng)力學(xué)方程表達(dá)式如公式(1)所示:
[0084]
[0085] 式中r和v分別表示探測(cè)器的位置和速度矢量; 為第i個(gè)攝動(dòng)行星的位置矢量;為第i個(gè)攝動(dòng)行星相對(duì)探測(cè)器的位置矢量,即 μs為太陽(yáng)引力常數(shù),μi為第i個(gè)攝動(dòng)行星的引力常數(shù);nt為攝動(dòng)行星的個(gè)數(shù)。CR為探測(cè)器表面的反射系數(shù),Ssrp為太陽(yáng)輻射光壓因子,m為探測(cè)器的質(zhì)量,k為推力系數(shù),最后一項(xiàng)a表示其他未建模加速度。
[0086] 選取探測(cè)器位置和速度作為狀態(tài)變量,則則探測(cè)器的狀態(tài)方程如公式(2)所示:
[0087]
[0088] 步驟2:建立巡航段自主導(dǎo)航測(cè)量模型。
[0089] 相機(jī)成像的過(guò)程采用小孔成像的模型,定義可見(jiàn)行星的質(zhì)心f1在相機(jī)坐標(biāo)系下的位置坐標(biāo)為rp=[xc?yc?zc]T,則其在相機(jī)像平面內(nèi)的像原像素坐標(biāo)如公式(3)所示:
[0090]
[0091] 其中:f為相機(jī)焦距,zc為目標(biāo)點(diǎn)沿著相機(jī)基準(zhǔn)線到相機(jī)成像平面的距離。
[0092] 相機(jī)系與探測(cè)器本體系重合,考慮到相機(jī)測(cè)量過(guò)程中的姿態(tài)偏差,即定義xc,yc方向的姿態(tài)偏差分別為θ1,θ2,則相機(jī)測(cè)量過(guò)程中的隨機(jī)旋轉(zhuǎn)將對(duì)特征點(diǎn)的位置測(cè)量產(chǎn)生影響,故實(shí)際的位置坐標(biāo)如公式(4)所示:
[0093]
[0094] 因?yàn)樽藨B(tài)角偏差較小,故簡(jiǎn)化為公式(5):
[0095]
[0096] 忽略分母影響,則簡(jiǎn)化為公式(6):
[0097]
[0098] 選擇可見(jiàn)行星質(zhì)心的角位置作為相機(jī)的觀測(cè)量,即測(cè)量矢量Z包括方位角 和俯仰角 如公式(7)所示:
[0099]
[0100] 其中: 和υφ分別代表俯仰角和方位角的測(cè)量噪聲,可見(jiàn)行星的質(zhì)心位置矢量為R,則俯仰角和方位角計(jì)算如公式(8)所示:
[0101]
[0102] 步驟3:確定探測(cè)器初始軌道
[0103] Nobs次觀測(cè)的初始軌道確定問(wèn)題寫(xiě)成公式(9)所示形式,
[0104]
[0105] 其中:rk,Rk分別代表第k次航天器和可見(jiàn)行星在日心慣性系下的位置矢量。未知的距離ρk是待求量, 表示探測(cè)器到可見(jiàn)行星的單位視線矢量。為進(jìn)一步提高步驟3中探測(cè)器初始軌道確定的求解效率和求解精度,步驟3中將待求量ρk表示成三階多項(xiàng)式的形式,通過(guò)最小二乘擬合方法求解滿(mǎn)足公式(9)約束的待求量ρk的近似解,選如下方法實(shí)現(xiàn):
[0106] 待求量ρk在日心慣性系下三軸分量如公式(10)所示:
[0107]
[0108] 考慮三階多項(xiàng)式的形式,將公式(10)表示成如公式(11)所示的形式:
[0109]
[0110] 公式(11)也應(yīng)該滿(mǎn)足公式(9)中的幾何約束,所以,
[0111]
[0112] tk表示第k次的測(cè)量時(shí)間。將公式(12)中等式左邊移到右邊,定義殘差如公式(13)所示:
[0113]
[0114] 進(jìn)一步將公式(13)表示成公式(14)所示形式:
[0115]
[0116] 為了得到最好的擬合,需要調(diào)整多項(xiàng)式擬合的系數(shù)ai,bi,和ci,使得殘差φx,φy,和φz最小,所以有,
[0117]
[0118] 因此,通過(guò)擬合方法確定待求多項(xiàng)式中的系數(shù),為進(jìn)一步提高求解效率,選最小二乘法確定待求多項(xiàng)式中的系數(shù),如公式(16):
[0119]
[0120] {a}={a1?a2?a3?a4}T,={b1?b2?b3?b4}T,{c}={c1?c2?c3?c4}T,3×3的矩陣T如公式(17):
[0121]
[0122] 至此,如公式(11)所示的三階多項(xiàng)式的系數(shù)已經(jīng)全部求解,根據(jù)系數(shù)能夠得到滿(mǎn)足精度要求的待求量ρk的近似解,即完成探測(cè)器初始軌道確定。
[0123] 步驟4:光行差效應(yīng)校正。
[0124] 通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置矢量。通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正步驟2中如公式(7)所示的角位置觀測(cè)信息,從而獲得更精確的觀測(cè)模型,即完成光行差效應(yīng)校正。
[0125] 步驟4具體實(shí)現(xiàn)方法如下:
[0126] 步驟4.1:通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)位置信息校正可見(jiàn)行星的位置。
[0127] 光行差校正指由于光到達(dá)探測(cè)器的過(guò)程中天體的運(yùn)動(dòng),導(dǎo)致天體真實(shí)位置和觀測(cè)位置之間產(chǎn)生的位移,對(duì)其進(jìn)行校正。針對(duì)這個(gè)問(wèn)題,需要對(duì)可見(jiàn)行星位置矢量進(jìn)行進(jìn)一步修正。假設(shè)觀測(cè)時(shí)間為tk,k=1,2...,Nobs,各個(gè)測(cè)量方向的單位矢量為 k=1,2...,Nobs。因?yàn)樘綔y(cè)器不能確定自己相對(duì)于太陽(yáng)和地球的距離,而僅已知觀測(cè)可見(jiàn)行星的位置矢量Rk,和相對(duì)可見(jiàn)行星的觀測(cè)方向的單位矢量 因此基于時(shí)間tk-δtk對(duì)可見(jiàn)行星真實(shí)位置矢量進(jìn)行求解,其中δtk是光從可見(jiàn)行星到達(dá)探測(cè)器所需的時(shí)間。首先確定δtk,如公式(18)所示:
[0128]
[0129] 其中c代表光速。對(duì)行星的位置矢量進(jìn)行校正,如公式(19)所示:
[0130] Rupdated=R(tk-δtk)???(19)
[0131] 可見(jiàn)行星位置矢量的更新由公式(18)和(19)確定,即完成可見(jiàn)行星位置矢量的校正。
[0132] 步驟4.2:通過(guò)可見(jiàn)行星和探測(cè)器之間的相對(duì)速度信息校正步驟2中如公式(7)所示的角位置觀測(cè)信息,得到更精確的測(cè)量模型。
[0133] 定義探測(cè)器速度方向和觀測(cè)方向 之間的夾角是θobs,真實(shí)方向和觀測(cè)方向之間的畸變角為ε,則修正后滿(mǎn)足精度要求的真實(shí)天體觀測(cè)方向 如公式(20)所示:
[0134]
[0135] 其中:是探測(cè)器速度的單位矢量,θtrue是探測(cè)器速度和真實(shí)觀測(cè)方向之間的夾角,且θtrue=θobs+ε。不同時(shí)刻的畸變角ε表示成公式(21)的形式:
[0136]
[0137] 其中:c和v分別代表光速大小和探測(cè)器速度大小。
[0138] 觀測(cè)方向和真實(shí)方向之間的關(guān)系如公式(22)所示,
[0139]
[0140] 可見(jiàn)行星和探測(cè)器之間視線方向的夾角和探測(cè)器的速度矢量是已知的,因此真實(shí)觀測(cè)方向可以通過(guò)公式(22)求解,結(jié)合公式(7)從而得到修正后真實(shí)的測(cè)量模型,至此完成光行差效應(yīng)校正。
[0141] 步驟5:根據(jù)步驟1中的探測(cè)器的動(dòng)力學(xué)模型、步驟3中確定的探測(cè)器初始軌道及步驟4修正后真實(shí)的測(cè)量模型,基于非線性系統(tǒng)濾波算法解算探測(cè)器實(shí)時(shí)狀態(tài)信息,實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。
[0142] 根據(jù)步驟1得到的巡航段探測(cè)器動(dòng)力學(xué)模型、步驟3中確定的探測(cè)器初始軌道及步驟4修正后真實(shí)的測(cè)量模型,通過(guò)導(dǎo)航濾波對(duì)探測(cè)器的狀態(tài)進(jìn)行估計(jì)。由于動(dòng)力學(xué)模型及測(cè)量模型均呈現(xiàn)非線性,故選用非線性濾波器,為進(jìn)一步提高導(dǎo)航濾波精度及收斂速度,所述的非線性濾波器優(yōu)選擴(kuò)展卡爾曼濾波(EKF)。解算探測(cè)器實(shí)時(shí)狀態(tài)信息,實(shí)現(xiàn)實(shí)時(shí)導(dǎo)航。
[0143] 地球是繞太陽(yáng)運(yùn)行的可見(jiàn)天體,選取地球?yàn)楸挥^測(cè)的可見(jiàn)天體,探測(cè)器通過(guò)觀測(cè)4
地球進(jìn)行導(dǎo)航。對(duì)于初始軌道確定部分,用于觀測(cè)的次數(shù)Nobs=8,時(shí)間間隔T=3×10s。導(dǎo)航相機(jī)的焦距為200mm,視場(chǎng)角為3°,分辨率為1024×1024,圖像處理精度為0.1個(gè)像素,相機(jī)測(cè)量誤差為σc=10-4rad,姿態(tài)偏差θ1,2=10-4rad??紤]探測(cè)器規(guī)劃調(diào)度和測(cè)量信息處理的時(shí)間,取濾波周期為ΔT=1500s,仿真時(shí)間設(shè)為T(mén)=1.5×105s。探測(cè)器初始位置誤差為
105km,速度誤差為10km/s,探測(cè)器的初始位置和速度矢量(r和v),地球的初始位置和速度矢量(R和V)如表1所示,
[0144] 表1?初始狀態(tài)仿真參數(shù)
[0145]
[0146] 基于位置偏差修正的巡航段自主導(dǎo)航方法的數(shù)值仿真結(jié)果分別如圖2所示,圖中分別給出了探測(cè)器位置和速度的導(dǎo)航估計(jì)誤差曲線。從仿真結(jié)果可以看出,修正后的導(dǎo)航方法有較高的導(dǎo)航精度和較快的濾波收斂速度,能夠?qū)崿F(xiàn)對(duì)探測(cè)器的位置及速度進(jìn)行實(shí)時(shí)估計(jì),最終得到高精度的狀態(tài)估計(jì)信息。
[0147] 本發(fā)明保護(hù)范圍不僅局限于實(shí)施例,實(shí)施例僅用于解釋本發(fā)明,凡與本發(fā)明在相同原理和構(gòu)思條件下的變更或修改均在本發(fā)明公開(kāi)的保護(hù)范圍之內(nèi)。
高效檢索全球?qū)@?/div>

專(zhuān)利匯是專(zhuān)利免費(fèi)檢索,專(zhuān)利查詢(xún),專(zhuān)利分析-國(guó)家發(fā)明專(zhuān)利查詢(xún)檢索分析平臺(tái),是提供專(zhuān)利分析,專(zhuān)利查詢(xún),專(zhuān)利檢索等數(shù)據(jù)服務(wù)功能的知識(shí)產(chǎn)權(quán)數(shù)據(jù)服務(wù)商。

我們的產(chǎn)品包含105個(gè)國(guó)家的1.26億組數(shù)據(jù),免費(fèi)查、免費(fèi)專(zhuān)利分析。

申請(qǐng)?jiān)囉?/a>

QQ群二維碼
意見(jiàn)反饋