一種衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法
【專(zhuān)利摘要】本發(fā)明公開(kāi)了一種衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法，步驟如下：(1)運(yùn)用拉格朗日多項(xiàng)式插值法，采用切比雪夫節(jié)點(diǎn)，即取不等間距的點(diǎn)將衛(wèi)星一天的軌道預(yù)報(bào)結(jié)果進(jìn)行插值，再由衛(wèi)星根據(jù)拉格朗日多項(xiàng)式計(jì)算衛(wèi)星軌道；(2)采用先找到俯仰角曲線(xiàn)下降段與坐標(biāo)軸交點(diǎn)再判斷時(shí)間窗口的方法，通過(guò)大步長(zhǎng)搜索俯仰角曲線(xiàn)與坐標(biāo)軸交點(diǎn)，采用變步長(zhǎng)搜索計(jì)算可見(jiàn)時(shí)間窗口。該方法針對(duì)未來(lái)星上搭載的自主任務(wù)規(guī)劃平臺(tái)，適用于低主頻小內(nèi)存的衛(wèi)星平臺(tái)，采用拉格朗日差值的方式實(shí)現(xiàn)了星上快速軌道計(jì)算，并設(shè)計(jì)了變步長(zhǎng)搜索與線(xiàn)性擬合相結(jié)合的衛(wèi)星對(duì)目標(biāo)可見(jiàn)時(shí)間窗口的快速計(jì)算方法，提供一種更加高效的衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法。
【專(zhuān)利說(shuō)明】
一種衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)領(lǐng)域，具體涉及一種衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法，適用 于低主頻小內(nèi)存平臺(tái)上的衛(wèi)星。
【背景技術(shù)】
[0002] 星上任務(wù)規(guī)劃平臺(tái)計(jì)算任務(wù)的目的與地面平臺(tái)的一致，但與地面任務(wù)規(guī)劃平臺(tái)可 以使用高性能計(jì)算機(jī)或集群計(jì)算不同，受星上能源系統(tǒng)、熱控系統(tǒng)和處理器性能等條件的 制約，目前星上任務(wù)規(guī)劃平臺(tái)的存儲(chǔ)器和處理器主頻與地面使用的計(jì)算機(jī)相比差距很大。 任務(wù)規(guī)劃過(guò)程包含許多復(fù)雜的迭代計(jì)算求解過(guò)程，其中訪(fǎng)問(wèn)預(yù)報(bào)所需要的計(jì)算資源最多， 因此同樣的計(jì)算過(guò)程在星上任務(wù)規(guī)劃平臺(tái)需要很長(zhǎng)的時(shí)間；另一方面衛(wèi)星在軌飛行速度快 (太陽(yáng)同步軌道衛(wèi)星約為7km/s)，而應(yīng)急任務(wù)出現(xiàn)的非常突然，從獲得目標(biāo)信息到執(zhí)行拍攝 任務(wù)之間給衛(wèi)星留出的規(guī)劃時(shí)間非常有限，在這種計(jì)算資源與時(shí)間矛盾下，迫切需要更高 效的衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)計(jì)算方法。
[0003] 衛(wèi)星對(duì)目標(biāo)可見(jiàn)時(shí)間窗口，是指衛(wèi)星飛行過(guò)程中，在其機(jī)動(dòng)能力限制范圍內(nèi)，衛(wèi)星 成像載荷對(duì)目標(biāo)點(diǎn)可見(jiàn)弧段所對(duì)應(yīng)的開(kāi)始和結(jié)束時(shí)間范圍。衛(wèi)星繞地球飛行過(guò)程中，由于 地球遮擋和衛(wèi)星機(jī)動(dòng)能力限制，衛(wèi)星對(duì)目標(biāo)不是時(shí)時(shí)可見(jiàn)的，而是在某個(gè)圈次某個(gè)弧段對(duì) 目標(biāo)可見(jiàn)，超出這個(gè)弧段衛(wèi)星對(duì)目標(biāo)不可見(jiàn)，也就無(wú)法對(duì)其進(jìn)行成像。
[0004] 衛(wèi)星對(duì)目標(biāo)的是否可見(jiàn)需要判斷衛(wèi)星與目標(biāo)之間是否有遮擋、衛(wèi)星能否通過(guò)姿態(tài) 機(jī)動(dòng)使成像載荷指向目標(biāo)，判斷方法為在t時(shí)刻建立衛(wèi)星和目標(biāo)兩點(diǎn)的連線(xiàn)，通過(guò)求連線(xiàn)與 地球表面的交點(diǎn)判斷是否被地球遮擋，再求衛(wèi)星對(duì)目標(biāo)的指向姿態(tài)(r，p，y)，判斷是否超出 衛(wèi)星的姿態(tài)機(jī)動(dòng)能力范圍，如果都滿(mǎn)足條件
[0005] 1" G [ ~Tmax , Tmax] /^\ P ^ [ -Pmax , Pmax] /^\ Y ^ [ -Ymax , Ymax]
[0006] 則說(shuō)明該時(shí)刻衛(wèi)星對(duì)目標(biāo)可見(jiàn)，其中rmax、pmax、y max分別表示衛(wèi)星最大滾動(dòng)角、最大 滾動(dòng)角和最大偏航角。
[0007] 衛(wèi)星對(duì)目標(biāo)的可見(jiàn)時(shí)間窗口的計(jì)算是用At的步長(zhǎng)在一定時(shí)間范圍內(nèi)[Tbegin，Tend] 遍歷判斷衛(wèi)星對(duì)目標(biāo)是否可見(jiàn)，搜索到衛(wèi)星進(jìn)出可見(jiàn)弧段的時(shí)刻點(diǎn)，就得出了衛(wèi)星在該時(shí) 間范圍內(nèi)對(duì)目標(biāo)的可見(jiàn)時(shí)間窗口 [TWbegin，TWend]，具體流程如附圖1所示。
[0008] 該計(jì)算流程需要瞬時(shí)衛(wèi)星的位置信息，在地面平臺(tái)中星歷數(shù)據(jù)可以預(yù)先通過(guò)高精 度星歷預(yù)報(bào)模塊推算出來(lái)，實(shí)際使用中通過(guò)查表的方式得到，花費(fèi)時(shí)間可以忽略;對(duì)于大時(shí) 間粒度星歷數(shù)據(jù)可以采用切比雪夫差值法，差值參數(shù)實(shí)時(shí)計(jì)算得到，在目前常見(jiàn)的計(jì)算機(jī) 上計(jì)算時(shí)間也可忽略不計(jì)。
[0009] 以存儲(chǔ)和主頻都極其有限的航天級(jí)DSP產(chǎn)品為例，一天24小時(shí)的秒級(jí)星歷數(shù)據(jù)至 少為24 X 60 X 60 X 7 X 8B~5MB(每秒數(shù)據(jù)為7個(gè)雙精度浮點(diǎn)數(shù)），無(wú)法在星上空間存儲(chǔ)，而且 大量數(shù)據(jù)的給上注也帶來(lái)很大挑戰(zhàn);采用切比雪夫差值的方法，由于其參數(shù)計(jì)算涉及高階 迭代計(jì)算，而且算法中調(diào)用次數(shù)極多，在多次實(shí)驗(yàn)中都無(wú)法在可接受的時(shí)間內(nèi)計(jì)算出時(shí)間 窗口（大于20分鐘），故該方法也不可行；另外還可以采用簡(jiǎn)化軌道外推模型，實(shí)現(xiàn)星上實(shí)時(shí) 軌道計(jì)算，但是這種方法在較遠(yuǎn)的時(shí)間計(jì)算精度過(guò)低。因此，開(kāi)發(fā)一種適合星上DSP計(jì)算的 小數(shù)據(jù)量、低誤差的快速星歷預(yù)報(bào)方法勢(shì)在必行。
[0010] 可見(jiàn)時(shí)間窗口計(jì)算極慢的另一個(gè)原因是搜索算法效率極低，采用固定At的步長(zhǎng) 在時(shí)間窗口搜索過(guò)程中，如果選取的At過(guò)小，則搜索次數(shù)多，搜索時(shí)間長(zhǎng);如果選取的At 過(guò)大，則容易遺漏時(shí)間較短的窗口。因此考慮采用更高效的搜索算法以縮短計(jì)算時(shí)間。

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

[0011] 為解決現(xiàn)有衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)計(jì)算領(lǐng)域的上述一個(gè)或多個(gè)問(wèn)題，本發(fā)明的目的在于提 供一種應(yīng)用于低主頻小內(nèi)存平臺(tái)上的，更加高效的衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法。
[0012] 為實(shí)現(xiàn)上述目的，本發(fā)明提供以下技術(shù)方案：
[0013] 一種衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法，步驟如下：
[0014] (1)運(yùn)用拉格朗日多項(xiàng)式插值法，采用切比雪夫節(jié)點(diǎn)，即取不等間距的點(diǎn)將衛(wèi)星一 天的軌道預(yù)報(bào)結(jié)果進(jìn)行插值，再由衛(wèi)星根據(jù)拉格朗日多項(xiàng)式計(jì)算衛(wèi)星軌道；
[0015] (2)采用先找到俯仰角曲線(xiàn)下降段與坐標(biāo)軸交點(diǎn)再判斷時(shí)間窗口的方法，通過(guò)大 步長(zhǎng)搜索俯仰角曲線(xiàn)與坐標(biāo)軸交點(diǎn)，采用變步長(zhǎng)搜索計(jì)算可見(jiàn)時(shí)間窗口。
[0016] 更進(jìn)一步地，步驟（1)的拉格朗日多項(xiàng)式插值法的具體操作方法為：將衛(wèi)星在 J2000慣性坐標(biāo)系下的位置速度參數(shù)(X,:/, 隨時(shí)間t的變化視為函數(shù)，其中 時(shí)間為自變量，利用拉格朗日多項(xiàng)式對(duì)該函數(shù)進(jìn)行插值擬合，設(shè)在連續(xù)區(qū)間[a，b]上有n+1 個(gè)不同點(diǎn)xq，xi, . . .，xn，分別對(duì)應(yīng)的函數(shù)值為yo，yi,. . .，yn，引入記號(hào)
[0017] COn+i(x) = (x-xo) (x-Xl)-- - (x-Xn)
[0018] 貝lj: ? 'n+i(xk) = (xk-xo)."(xk_xk-i)(xk-xk+i)."(xk_xn)
[0019] 拉格朗日多項(xiàng)式為
[0021 ]取不等間距的點(diǎn)進(jìn)行插值為
[0022]更進(jìn)一步地，步驟(2)的具體操作方法為:采用大步長(zhǎng)判斷連續(xù)兩點(diǎn)是否下降且跨 零點(diǎn)，再通過(guò)二分法快速找到交點(diǎn)時(shí)刻，判斷該時(shí)刻是否有時(shí)間窗口，如果沒(méi)有則繼續(xù)搜 索，如果有則開(kāi)始搜索時(shí)間窗口的起止時(shí)刻，由于出現(xiàn)時(shí)間窗口的位置曲線(xiàn)斜率大時(shí)間窗 口短，因此采用直線(xiàn)擬合的方式快速跳到窗口起止點(diǎn)附近，再通過(guò)簡(jiǎn)單搜索即可找到精確 起止時(shí)刻點(diǎn)，即得出目標(biāo)可見(jiàn)時(shí)間窗口。
[0023]本發(fā)明具有以下有益效果：
[0024] (1)由于完整星歷數(shù)據(jù)過(guò)大、簡(jiǎn)單模型外推方法精度低，所以仍然考慮差值擬合的 方法，由于通過(guò)大粒度星歷數(shù)據(jù)實(shí)時(shí)計(jì)算差值參數(shù)的復(fù)雜度高，故考慮采用固定差值參數(shù)， 這樣可以在地面系統(tǒng)中分段計(jì)算出差值參數(shù)，通過(guò)指令上注的方式將差值參數(shù)上注到衛(wèi) 星，可以顯著降低星上任務(wù)規(guī)劃平臺(tái)計(jì)算軌道數(shù)據(jù)的復(fù)雜度。但是采用拉格朗日插值法會(huì) 出現(xiàn)龍格現(xiàn)象。因此為了使整體誤差控制在一定范圍，不能采用等距節(jié)點(diǎn)，采用切比雪夫節(jié) 點(diǎn)，即取不等間距的點(diǎn)進(jìn)行插值，使得接近邊界區(qū)域的點(diǎn)比中間區(qū)域的點(diǎn)更密集。利用該方 法將衛(wèi)星一天的軌道預(yù)報(bào)結(jié)果進(jìn)行插值，可以?xún)H僅上傳多項(xiàng)式參數(shù)，再由衛(wèi)星根據(jù)拉格朗 日多項(xiàng)式計(jì)算衛(wèi)星軌道，能夠有效減少上傳的數(shù)據(jù)量。
[0025] (2)衛(wèi)星對(duì)目標(biāo)的時(shí)間窗口的基本計(jì)算方法，其計(jì)算所需時(shí)間與算法循環(huán)次數(shù)成 正比，與搜索步長(zhǎng)成反比，因此減少搜索次數(shù)可以有效減少計(jì)算時(shí)間。但是采取加大搜索步 長(zhǎng)的方式可能遺漏長(zhǎng)度較短的時(shí)間窗口，因此可以考慮通過(guò)數(shù)據(jù)分析得出規(guī)律，再基于該 規(guī)律采用變步長(zhǎng)搜索計(jì)算可見(jiàn)時(shí)間窗口。
[0026] (3)將帶有快速時(shí)間窗口計(jì)算的任務(wù)規(guī)劃預(yù)處理算法移植到DSP上可以很快地計(jì) 算出時(shí)間窗口，完成一個(gè)四條帶拼接任務(wù)的處理從原來(lái)的大于20分鐘降低到20秒鐘，而且 精度也可以滿(mǎn)足要求。
【附圖說(shuō)明】
[0027] 為了更清楚地說(shuō)明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖做簡(jiǎn)單的介紹，顯而易見(jiàn)地，下面描述中的附圖是本發(fā)明 的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù) 這些附圖獲得其他的附圖。
[0028] 附圖1為本發(fā)明【背景技術(shù)】的可見(jiàn)時(shí)間窗口計(jì)算流程圖。
[0029] 附圖2為本發(fā)明實(shí)施例的兩種節(jié)點(diǎn)選取方式對(duì)比，其中a為等距節(jié)點(diǎn)，b為切比雪夫 節(jié)點(diǎn)。
[0030] 附圖3為本發(fā)明實(shí)施例的衛(wèi)星對(duì)目標(biāo)指向姿態(tài)曲線(xiàn)。
[0031] 附圖4為本發(fā)明實(shí)施例的拉格朗日差值的誤差曲線(xiàn)。
【具體實(shí)施方式】
[0032] 下面通過(guò)具體實(shí)施例對(duì)本發(fā)明作進(jìn)一步的說(shuō)明。
[0033]本發(fā)明公開(kāi)了一種衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法，步驟如下：
[0034] (1)運(yùn)用拉格朗日多項(xiàng)式插值法，采用切比雪夫節(jié)點(diǎn)，即取不等間距的點(diǎn)將衛(wèi)星一 天的軌道預(yù)報(bào)結(jié)果進(jìn)行插值，再由衛(wèi)星根據(jù)拉格朗日多項(xiàng)式計(jì)算衛(wèi)星軌道；
[0035] (2)采用先找到俯仰角曲線(xiàn)下降段與坐標(biāo)軸交點(diǎn)再判斷時(shí)間窗口的方法，通過(guò)大 步長(zhǎng)搜索俯仰角曲線(xiàn)與坐標(biāo)軸交點(diǎn)，采用變步長(zhǎng)搜索計(jì)算可見(jiàn)時(shí)間窗口。
[0036] 其中，步驟（1)為基于拉格朗日插值的衛(wèi)星軌道快速預(yù)報(bào)方法。由于完整星歷數(shù)據(jù) 過(guò)大、簡(jiǎn)單模型外推方法精度低，所以仍然考慮插值擬合的方法，由于通過(guò)大粒度星歷數(shù)據(jù) 實(shí)時(shí)計(jì)算插值參數(shù)的復(fù)雜度高，故考慮采用固定插值參數(shù)，這樣可以在地面系統(tǒng)中分段計(jì) 算出插值參數(shù)，通過(guò)指令上注的方式將插值參數(shù)上注到衛(wèi)星，可以顯著降低星上任務(wù)規(guī)劃 平臺(tái)計(jì)算軌道數(shù)據(jù)的復(fù)雜度。
[0037] 拉格朗日多項(xiàng)式插值法是一種利用多項(xiàng)式對(duì)函數(shù)進(jìn)行近似的常用方法，將衛(wèi)星在 J2000慣性坐標(biāo)系下的位置速度參數(shù)(mi,>',i)隨時(shí)間t的變化視為函數(shù)，其中時(shí)間為自 變量，利用拉格朗日多項(xiàng)式對(duì)該函數(shù)進(jìn)行插值擬合，設(shè)在連續(xù)區(qū)間[a，b]上有n+1個(gè)不同點(diǎn) xo，xi,…，x n，分別對(duì)應(yīng)的函數(shù)值為yo，yi,…，yn，引入記號(hào)
[0038] ? n+l(x) = (x-xo) (x-Xl)-- - (x-Xn)
[0039] 貝 1J:
[0040] 〇 7 n+1 ( Xk ) = ( Xk-X0 ) ? ? ? ( Xk-Xk-1) ( Xk-Xk+1) ? ? ? ( Xk-Xn )
[0041] 拉格朗日多項(xiàng)式為
[0043]但是采用拉格朗日插值法會(huì)出現(xiàn)龍格現(xiàn)象。因此為了使整體誤差控制在一定范 圍，不能采用等距節(jié)點(diǎn)，采用切比雪夫節(jié)點(diǎn)，即取不等間距的點(diǎn)進(jìn)行插值，使得接近邊界區(qū) 域的點(diǎn)比中間區(qū)域的點(diǎn)更密集：
[0045]利用等距節(jié)點(diǎn)和切比雪夫節(jié)點(diǎn)分別對(duì)函數(shù)f (x) zQ+x2)^1進(jìn)行插值，取n=10,結(jié) 果對(duì)比如附圖2所示，其中，a為等距節(jié)點(diǎn)，b為切比雪夫節(jié)點(diǎn)，可以看出，由于采用了切比雪 夫節(jié)點(diǎn)，邊界處取點(diǎn)更密集，龍格現(xiàn)象不再明顯，因此本文亦采用切比雪夫節(jié)點(diǎn)。
[0046]利用該方法將衛(wèi)星一天的軌道預(yù)報(bào)結(jié)果進(jìn)行插值，可以?xún)H僅上傳多項(xiàng)式參數(shù)，再 由衛(wèi)星根據(jù)拉格朗日多項(xiàng)式計(jì)算衛(wèi)星軌道，能夠有效減少上傳的數(shù)據(jù)量。
[0047] 例如一天24小時(shí)總共86400秒，將一天時(shí)間共分為16段，其中每段5400秒，對(duì)每一 個(gè)參數(shù)，在每一段進(jìn)行次數(shù)為10的拉格朗日插值，產(chǎn)生一組由11個(gè)雙精度浮點(diǎn)數(shù)構(gòu)成的系 數(shù)，這樣原本5M數(shù)據(jù)壓縮到16乂11乂6乂88~80數(shù)據(jù)，上傳和占用內(nèi)存空間的數(shù)據(jù)量大大 減少。將拉格朗日多項(xiàng)式完全展開(kāi)，求出多項(xiàng)式系數(shù)。
[0048]其中，步驟(2)為快速時(shí)間窗口搜索方法。衛(wèi)星對(duì)目標(biāo)的時(shí)間窗口的基本計(jì)算方法 中，其計(jì)算所需時(shí)間與算法循環(huán)次數(shù)成正比，與搜索步長(zhǎng)成反比，因此減少搜索次數(shù)可以有 效減少計(jì)算時(shí)間。但是采取加大搜索步長(zhǎng)的方式可能遺漏長(zhǎng)度較短的時(shí)間窗口，因此可以 考慮通過(guò)數(shù)據(jù)分析得出規(guī)律，再基于該規(guī)律采用變步長(zhǎng)搜索計(jì)算可見(jiàn)時(shí)間窗口。
[0049]計(jì)算一天內(nèi)衛(wèi)星對(duì)目標(biāo)的指向姿態(tài)（忽略是否可見(jiàn)），取北京時(shí)間2013年2月2日上 午八點(diǎn)至次日上午八點(diǎn)，目標(biāo)點(diǎn)煒經(jīng)高度坐標(biāo)為(60°，120°，0° )，假設(shè)衛(wèi)星側(cè)擺俯仰能力均 為45°，將衛(wèi)星對(duì)目標(biāo)指向的側(cè)擺角、俯仰角、衛(wèi)星的Z軸坐標(biāo)和時(shí)間窗口判定隨時(shí)間的變化 曲線(xiàn)繪制在一個(gè)圖內(nèi)，如附圖3所示(為方便觀(guān)察已將每個(gè)縱坐標(biāo)進(jìn)行歸一化處理），從圖中 可以看出在一天內(nèi)時(shí)間窗口出現(xiàn)的時(shí)間很短（約600秒），所以采用固定小步長(zhǎng)搜索方式效 率極低，而采用大步長(zhǎng)搜索容易跨過(guò)整個(gè)時(shí)間窗口。另外，從圖中可以看出兩個(gè)姿態(tài)角的變 化周期性不嚴(yán)格(與Z軸坐標(biāo)為基準(zhǔn)），所以也很難采用周期性的預(yù)測(cè)或擬合方法。
[0050] 但是從圖中可以總結(jié)出來(lái)一些規(guī)律并給出合理的解釋?zhuān)?br>[0051] (1)時(shí)間窗口總是出現(xiàn)在俯仰角曲線(xiàn)下降段與坐標(biāo)軸交點(diǎn)處
[0052] (2)時(shí)間窗口總是出現(xiàn)在俯仰角曲線(xiàn)振幅較大段；
[0053] (3)時(shí)間窗口是否出現(xiàn)與側(cè)擺曲線(xiàn)關(guān)系不大；
[0054] (4)俯仰角曲線(xiàn)的周期不會(huì)與Z軸坐標(biāo)曲線(xiàn)周期相差很多。
[0055]因此可以考慮采用先找到俯仰角曲線(xiàn)下降段與坐標(biāo)軸交點(diǎn)再判斷時(shí)間窗口的方 法，通過(guò)大步長(zhǎng)搜索俯仰角曲線(xiàn)與坐標(biāo)軸交點(diǎn)，實(shí)際中采用大步長(zhǎng)判斷連續(xù)兩點(diǎn)是否下降 且跨零點(diǎn)，再通過(guò)二分法快速找到交點(diǎn)時(shí)刻，判斷該時(shí)刻是否有時(shí)間窗口，如果沒(méi)有則繼續(xù) 搜索，如果有則開(kāi)始搜索時(shí)間窗口的起止時(shí)刻，由于出現(xiàn)時(shí)間窗口的位置曲線(xiàn)斜率大時(shí)間 窗口短，因此可以采用直線(xiàn)擬合的方式快速跳到窗口起止點(diǎn)附近，再通過(guò)簡(jiǎn)單搜索即可找 到精確起止時(shí)刻點(diǎn)，由于衛(wèi)星執(zhí)行觀(guān)測(cè)動(dòng)作的時(shí)間以秒為單位，故算法只精確到秒。
[0056]以下為一個(gè)具體的應(yīng)用實(shí)施例。
[0057]采用拉格朗日差值方法方法對(duì)一天24小時(shí)的軌道預(yù)報(bào)結(jié)果進(jìn)行插值，得到的6個(gè) 參數(shù)的誤差曲線(xiàn)如附圖4所示，從圖中可以看出采用該方法計(jì)算的星歷數(shù)據(jù)誤差較低(24小 時(shí)內(nèi)位置誤差在200米以?xún)?nèi)），優(yōu)于上海小衛(wèi)星中心的吳會(huì)英等人發(fā)明星上軌道外推算法 (18小時(shí)誤差在10公里以?xún)?nèi)），而且采用多項(xiàng)式計(jì)算更為簡(jiǎn)單快速，適合衛(wèi)星上計(jì)算能力較 低的環(huán)境，是一種較好的衛(wèi)星星上軌道快速推導(dǎo)方法。
[0058]對(duì)于快速時(shí)間窗口搜索方法，選取維度為0°、30°、60°和80°的四個(gè)目標(biāo)點(diǎn)，利用標(biāo) 準(zhǔn)算法、基于精確星歷和差值星歷的快速搜索算法計(jì)算目標(biāo)點(diǎn)的時(shí)間窗口，時(shí)間范圍為北 京時(shí)間2013-2-2 08:00:00~2013-2-2 17:00:00，統(tǒng)計(jì)其計(jì)算過(guò)程中搜索次數(shù)和時(shí)間窗口 的誤差，為避免計(jì)算平臺(tái)不同造成的時(shí)間統(tǒng)計(jì)不準(zhǔn)確問(wèn)題，試驗(yàn)中用搜索次數(shù)代替執(zhí)行時(shí) 間比較算法效率，結(jié)果表1如所示。
[0059] 表1時(shí)間窗口計(jì)算結(jié)果對(duì)比
[0062] 從表中可以看出，搜索次數(shù)與存在時(shí)間窗口的數(shù)量正相關(guān)，快速搜索算法搜索次 數(shù)遠(yuǎn)遠(yuǎn)小于標(biāo)準(zhǔn)時(shí)間窗口計(jì)算方法;星上任務(wù)規(guī)劃平臺(tái)時(shí)間精度為1秒，因此精確算法、基 于差值星歷與基于精確星歷的快速搜索算法結(jié)果幾乎一致，說(shuō)明星歷誤差對(duì)時(shí)間窗口計(jì)算 的影響小于系統(tǒng)的時(shí)間分辨率，該時(shí)間窗快速搜索算法正確有效。
[0063] 將帶有快速時(shí)間窗口計(jì)算的任務(wù)規(guī)劃預(yù)處理算法移植到DSP上可以很快地計(jì)算出 時(shí)間窗口，完成一個(gè)四條帶拼接任務(wù)的處理從原來(lái)的大于20分鐘降低到20秒鐘，而且精度 也可以滿(mǎn)足要求。
[0064] 本文中應(yīng)用了具體實(shí)施例對(duì)本發(fā)明的原理及實(shí)施方式進(jìn)行了闡述，以上實(shí)施例的 說(shuō)明只是用于幫助理解本發(fā)明的方法及其核心思想。以上所述僅是本發(fā)明的優(yōu)選實(shí)施方 式，應(yīng)當(dāng)指出，由于文字表達(dá)的有限性，而客觀(guān)上存在無(wú)限的具體結(jié)構(gòu)，對(duì)于本技術(shù)領(lǐng)域的 普通技術(shù)人員來(lái)說(shuō)，在不脫離本發(fā)明原理的前提下，還可以做出若干改進(jìn)、潤(rùn)飾或變化，也 可以將上述技術(shù)特征以適當(dāng)?shù)姆绞竭M(jìn)行組合;這些改進(jìn)潤(rùn)飾、變化或組合，或未經(jīng)改進(jìn)將發(fā) 明的構(gòu)思和技術(shù)方案直接應(yīng)用于其它場(chǎng)合的，均應(yīng)視為本發(fā)明的保護(hù)范圍。
【主權(quán)項(xiàng)】
1. 一種衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法，其特征在于步驟如下： (1) 運(yùn)用拉格朗日多項(xiàng)式插值法，采用切比雪夫節(jié)點(diǎn)，即取不等間距的點(diǎn)將衛(wèi)星一天的 軌道預(yù)報(bào)結(jié)果進(jìn)行插值，再由衛(wèi)星根據(jù)拉格朗日多項(xiàng)式計(jì)算衛(wèi)星軌道； (2) 采用先找到俯仰角曲線(xiàn)下降段與坐標(biāo)軸交點(diǎn)再判斷時(shí)間窗口的方法，通過(guò)大步長(zhǎng) 搜索俯仰角曲線(xiàn)與坐標(biāo)軸交點(diǎn)，采用變步長(zhǎng)搜索計(jì)算可見(jiàn)時(shí)間窗口。2. 根據(jù)權(quán)利要求1所述的衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法，其特征在于，步驟（1)的拉格朗 日多項(xiàng)式插值法的具體操作方法為：將衛(wèi)星在J2000慣性坐標(biāo)系下的位置速度參數(shù)隨時(shí)間t的變化視為函數(shù)，其中時(shí)間為自變量，利用拉格朗日多項(xiàng)式 對(duì)該函數(shù)進(jìn)行插值擬合，設(shè)在連續(xù)區(qū)間[a，b]上有n+1個(gè)不同點(diǎn)χο,χι, . . .，χη，分別對(duì)應(yīng)的函 數(shù)值為y。，yi，…，yn，引入記號(hào) ω n+i (X ) = ( χ-χο)(X-Xi)…（χ-χη) 則：ω ' n+l ( Xk ) = ( Xk-XO )…（Xk_Xk-l)( Xk_Xk+l)…（Xk-Xn )拉格朗日多項(xiàng)式為 取不等間距的點(diǎn)進(jìn)行插43. 根據(jù)權(quán)利要求1或2所述的衛(wèi)星訪(fǎng)問(wèn)預(yù)報(bào)快速計(jì)算方法，其特征在于，步驟(2)的具體 操作方法為:采用大步長(zhǎng)判斷連續(xù)兩點(diǎn)是否下降且跨零點(diǎn)，再通過(guò)二分法快速找到交點(diǎn)時(shí) 亥IJ，判斷該時(shí)刻是否有時(shí)間窗口，如果沒(méi)有則繼續(xù)搜索，如果有則開(kāi)始搜索時(shí)間窗口的起止 時(shí)刻，由于出現(xiàn)時(shí)間窗口的位置曲線(xiàn)斜率大時(shí)間窗口短，因此采用直線(xiàn)擬合的方式快速跳 到窗口起止點(diǎn)附近，再通過(guò)簡(jiǎn)單搜索即可找到精確起止時(shí)刻點(diǎn)，即得出目標(biāo)可見(jiàn)時(shí)間窗口。
【文檔編號(hào)】G06F17/50GK105893659SQ201610186066
【公開(kāi)日】2016年8月24日
【申請(qǐng)日】2016年6月2日
【發(fā)明人】陳英武, 袁駔, 何磊, 楊文沅, 賀仁杰, 姚鋒, 邢立寧, 劉曉路, 王沛, 陳成, 陳盈果, 王濤, 劉慧慧
【申請(qǐng)人】中國(guó)人民解放軍國(guó)防科學(xué)技術(shù)大學(xué)