星載脈沖線性調(diào)頻波雷達(dá)測(cè)高方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于星球探測(cè)領(lǐng)域，特別是針對(duì)星載探地雷達(dá)對(duì)地高度信息的需要，在不 增加其他載荷的前提下提供準(zhǔn)確的高度信息。
【背景技術(shù)】
[0002] 近年來(lái)，隨著人類對(duì)宇宙空間探索的逐漸深入，衛(wèi)星被越來(lái)越多地用于對(duì)外星球 地形、地貌、地質(zhì)構(gòu)造的探測(cè)。為實(shí)現(xiàn)探測(cè)，衛(wèi)星一般搭載星載探地雷達(dá)，使用雷達(dá)回波實(shí)現(xiàn) 對(duì)探測(cè)目標(biāo)的檢測(cè)和分析。
[0003] 星載探地雷達(dá)在使用過(guò)程中需要有衛(wèi)星高度信息的支持，如果高度不夠準(zhǔn)確，會(huì) 導(dǎo)致開(kāi)窗時(shí)間錯(cuò)誤，不能準(zhǔn)確采集回波信號(hào)，難以實(shí)現(xiàn)探測(cè)目標(biāo)，因此需要有比較精確的衛(wèi) 星高度信息。在目前的火星探測(cè)衛(wèi)星上，有些星載雷達(dá)裝載有激光測(cè)高儀或者Ku、S、C、波 段的測(cè)高儀，這需要額外的載荷。

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

[0004] 為了解決衛(wèi)星軌道高度不精確，本發(fā)明在不適用其他載荷的前提下提供了一種星 載脈沖線性調(diào)頻波雷達(dá)測(cè)高的方法來(lái)得到衛(wèi)星飛行的高度。
[0005] (二）技術(shù)方案
[0006] 本發(fā)明提供一種星載脈沖線性調(diào)頻波雷達(dá)測(cè)高方法，用于繞星球飛行的衛(wèi)星，該 方法包括步驟：
[0007] 步驟A，根據(jù)地面高程數(shù)據(jù)選擇一段平坦地區(qū)作為測(cè)高起始地區(qū)；
[0008] 步驟B，根據(jù)衛(wèi)星姿態(tài)調(diào)節(jié)雷達(dá)天線方向，使天線波束垂直于星球表面；
[0009] 步驟C，發(fā)射脈沖線性調(diào)頻雷達(dá)信號(hào)，信號(hào)的載頻為F，帶寬為B，脈沖寬度為T(mén)111，脈 沖重復(fù)周期（PRT)為T(mén)。；
[0010] 步驟D，根據(jù)衛(wèi)星飛行軌道提供的高度信息R。，估計(jì)雷達(dá)信號(hào)回波返回衛(wèi)星的時(shí)間
-，其中c為光速；
[0011] 步驟E，根據(jù)步驟D計(jì)算的回波延時(shí)時(shí)間開(kāi)窗采集回波信號(hào)，采樣時(shí)間窗多出 20us，前后各留IOus作余量長(zhǎng)為1"+20 · 10 6,開(kāi)始采樣時(shí)間為τ。-1〇 · 10 6,采樣率為Fs ;
[0012] 步驟F，對(duì)步驟E采集的回波信號(hào)進(jìn)行匹配濾波，累加匹配濾波的輸出結(jié)果，并分 析累加后信號(hào)的信噪比SNR ;
[0013] 步驟G，判斷步驟F中輸出信號(hào)的信噪比是否滿足要求，若滿足則繼續(xù)進(jìn)行下一 步，若不滿足則回到步驟B繼續(xù)進(jìn)行發(fā)射、采集、匹配濾波和累加；以及
[0014] 步驟Η，根據(jù)匹配濾波的結(jié)果得到一個(gè)脈沖信號(hào)，脈沖信號(hào)距離采樣時(shí)間窗起點(diǎn)的 時(shí)間為τ i，根據(jù)脈沖信號(hào)在時(shí)間軸上的位置計(jì)算從信號(hào)采集開(kāi)始到出現(xiàn)脈沖的時(shí)間，結(jié)合 步驟E中設(shè)定的開(kāi)窗時(shí)間推導(dǎo)出衛(wèi)星距離星球表面的高度信息〗
[0015] 根據(jù)本發(fā)明的實(shí)施例，該方法還可以包括：步驟I，衛(wèi)星慣性導(dǎo)航系統(tǒng)可以提供星 球運(yùn)行的速度 ν'(?)，根據(jù)地面高程數(shù)據(jù)可以得到星下點(diǎn)位置的地形高度變化曲線d(l)，再 結(jié)合步驟H計(jì)算的高度信息H。，可以實(shí)時(shí)推導(dǎo)出衛(wèi)星高度信息h (t)。
[0016] 根據(jù)本發(fā)明的實(shí)施例，該方法還可以包括：步驟J，下一次出現(xiàn)星球表面為平坦地 區(qū)時(shí)，重復(fù)上述步驟得到新的高度信息，校正實(shí)時(shí)跟蹤高度誤差。
[0017] 根據(jù)本發(fā)明的實(shí)施例，該方法還可以包括：步驟K，根據(jù)記錄下來(lái)的數(shù)據(jù)可以做出 衛(wèi)星飛行的高度軌跡圖和地面高程變化圖。
[0018] (三）有益效果
[0019] 本發(fā)明可以在不增加其他載荷的前提下，利用星載探地雷達(dá)實(shí)現(xiàn)對(duì)衛(wèi)星高度的計(jì) 算，并實(shí)時(shí)跟蹤衛(wèi)星的高度。相比于現(xiàn)有技術(shù)，節(jié)省了衛(wèi)星的負(fù)重和成本。
【附圖說(shuō)明】
[0020] 圖1示意性示出了根據(jù)本發(fā)明實(shí)施例的星載脈沖線性調(diào)頻波雷達(dá)測(cè)高方法的流 程圖；
[0021] 圖2示意性示出了根據(jù)本發(fā)明實(shí)施例的星載脈沖線性調(diào)頻波雷達(dá)測(cè)高方法的匹 配濾波輸出結(jié)果圖。
【具體實(shí)施方式】
[0022] 為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，以下結(jié)合針對(duì)火星的具體實(shí) 施例，并參照附圖，對(duì)本發(fā)明作進(jìn)一步的詳細(xì)說(shuō)明。
[0023] 根據(jù)本發(fā)明的實(shí)施例，通過(guò)繞火星飛行的衛(wèi)星上的星載探地雷達(dá)發(fā)射天線向火星 表面周期地發(fā)射脈沖線性調(diào)頻信號(hào)，該信號(hào)在傳到火星表面時(shí)發(fā)生反射。接收天線通過(guò)預(yù) 估發(fā)射信號(hào)的傳輸時(shí)間，選擇適合的開(kāi)窗時(shí)間對(duì)回波信號(hào)進(jìn)行采樣接收，并通過(guò)對(duì)雷達(dá)回 波信號(hào)的延時(shí)進(jìn)行分析，得到衛(wèi)星的高度信息。
[0024] 圖1示意性示出了根據(jù)本發(fā)明實(shí)施例的火星星載脈沖線性調(diào)頻波雷達(dá)測(cè)高方法 的流程圖。圖2示意性示出了該方法中回波信號(hào)匹配濾波的輸出結(jié)果圖。參照?qǐng)D1，根據(jù)本 發(fā)明實(shí)施例的星載脈沖線性調(diào)頻波雷達(dá)測(cè)高的方法的步驟包括：
[0025] 步驟A，根據(jù)火星地面高程數(shù)據(jù)選擇一段平坦地區(qū)作為測(cè)高起始地區(qū)。盡量選取地 形起伏抖動(dòng)不超過(guò)50米，且變化平滑的區(qū)域。
[0026] 步驟B，根據(jù)衛(wèi)星姿態(tài)調(diào)節(jié)雷達(dá)天線方向，使天線波束垂直于火星表面。天線波束 垂直于火星表面才可以得到信號(hào)強(qiáng)度比較強(qiáng)的回波信號(hào)，而且所測(cè)的延時(shí)才是衛(wèi)星到火星 表面星下點(diǎn)的延時(shí)。
[0027] 步驟C，發(fā)射脈沖線性調(diào)頻雷達(dá)信號(hào)，信號(hào)的載頻為F，帶寬為B，脈沖寬度為T(mén)111，脈 沖重復(fù)周期（PRT)為T(mén)。。因此發(fā)送的脈沖線性調(diào)頻雷達(dá)信號(hào)x(t)為：
[0029] 其中，t表示時(shí)間變量，η為自然數(shù)，少為初始相位，信號(hào)的載頻為F，帶寬為B，脈沖 寬度為T(mén)ni,脈沖重復(fù)周期（PRT)為T(mén)。。相位函數(shù)為：
[0030] 對(duì)相位求導(dǎo)得到頻率函數(shù)：
[0031]
[0032] 步驟D，根據(jù)衛(wèi)星飛行軌道提供的高度信息R。，估計(jì)雷達(dá)信號(hào)回波返回衛(wèi)星的時(shí)間
其中c為光速。衛(wèi)星在繞火星飛行的過(guò)程中會(huì)提供衛(wèi)星的飛行高度信息，但是所 提供的飛行高度信息不夠精確，有上千米的誤差。因此τ。只是雷達(dá)信號(hào)發(fā)射到接受大概 所需的延時(shí)，不夠準(zhǔn)確，需要做進(jìn)一步的處理，得到更準(zhǔn)確的延時(shí)。
[0033] 步驟Ε，根據(jù)步驟D計(jì)算的回波延時(shí)時(shí)間開(kāi)窗采集回波信號(hào)，由于衛(wèi)星飛行軌道高 度信息不夠準(zhǔn)確，因此采樣時(shí)間窗需要比脈沖信號(hào)寬度T ni多出20us，前后各留IOus作余 量。因此時(shí)間窗長(zhǎng)為!"+20 · 10 6,開(kāi)始采樣時(shí)間為τ。-1〇 · 10 6。采樣率為Fs。
[0034] 設(shè)雷達(dá)信號(hào)經(jīng)過(guò)的準(zhǔn)確延時(shí)為τ，其中τ <T。。則回波信號(hào)的表達(dá)式為：
選擇的米樣時(shí)間為：[τ。_10 · 10 6， τ。+1"+10 · 10 6]，采樣率為Fs。采樣點(diǎn)數(shù)為：N = Fs · (Tm+20 · 10 6)。得到的離散采樣信 號(hào)為：
[0035]
[0036] 其中，k = 0,1 …N-1。
[0037] 步驟F，對(duì)步驟E采集的回波信號(hào)進(jìn)行匹配濾波，并累加匹配濾波的輸出結(jié)果，并 分析累加后信號(hào)的信噪比SNR。
[0038] 匹配濾波器為Μ? = -〇:，其中上標(biāo)*表示取共輒。匹配濾波器的輸出函數(shù)為
[0040] 如圖2所示，匹配濾波器的輸出函數(shù)可以近似為一脈沖信號(hào) [0041] sout (t) ^Tn- Sinc(Ba-T1))
[0042] 將得到的脈沖信號(hào)結(jié)果存儲(chǔ)起來(lái)，并且與之前得到的結(jié)果進(jìn)行累加得到累加后的 信號(hào)。脈沖的位置s n (t)出現(xiàn)在開(kāi)