丁宗華 唐志美 代連東 吳健 許正文

(中國(guó)電波傳播研究所 電波環(huán)境特性及模化技術(shù)重點(diǎn)實(shí)驗(yàn)室,青島 266107)

非相干散射雷達(dá)的高信噪比空間物體距離測(cè)量方法

丁宗華 唐志美 代連東 吳健 許正文

(中國(guó)電波傳播研究所 電波環(huán)境特性及?；夹g(shù)重點(diǎn)實(shí)驗(yàn)室,青島 266107)

非相干散射雷達(dá)是目前地基電離層探測(cè)的強(qiáng)大手段,同時(shí)在空間物體探測(cè)方面也具有重要應(yīng)用前景．當(dāng)空間物體回波信噪比很高時(shí),可以利用單個(gè)脈沖進(jìn)行目標(biāo)檢測(cè),而無(wú)需相干積累．非相干散射雷達(dá)常采用相位編碼脈沖,在單個(gè)相位編碼脈沖內(nèi)存在多次相位翻轉(zhuǎn),此相位翻轉(zhuǎn)在接收機(jī)濾波后顯示一定坡度,通過(guò)搜索此坡度曲線上具有最大斜率的坡度點(diǎn)可高精度確定發(fā)射-接收回波脈沖內(nèi)各對(duì)應(yīng)相位翻轉(zhuǎn)的時(shí)間差以及對(duì)應(yīng)的距離,對(duì)這些距離值進(jìn)行加權(quán)最小二乘擬合,即可給出該回波脈沖對(duì)應(yīng)的空間物體距離及誤差估計(jì)值．利用歐洲非相干散射科學(xué)聯(lián)合會(huì)(European Incoherent Scatter Scientific Association,EISCAT)實(shí)測(cè)數(shù)據(jù)分析發(fā)現(xiàn),距離誤差可達(dá)數(shù)十米量級(jí),遠(yuǎn)優(yōu)于以往的基于發(fā)射-接收回波脈沖前沿時(shí)間差方法．

非相干散射雷達(dá);空間物體;高信噪比;距離;相位翻轉(zhuǎn)

引 言

隨著人類(lèi)航天活動(dòng)的日益頻繁,空間碎片等空間物體數(shù)目日益增多,對(duì)人類(lèi)空間活動(dòng)和國(guó)家安全等帶來(lái)了嚴(yán)重影響,同時(shí)為空間和天文科學(xué)研究提出了更多機(jī)遇和挑戰(zhàn),也對(duì)空間物體探測(cè)提出了更高要求．在諸多探測(cè)方法中,地基雷達(dá)具有探測(cè)精度高、探測(cè)覆蓋區(qū)域廣、不受氣象條件限制等優(yōu)點(diǎn),一直是空間物體探測(cè)的主要探測(cè)手段．但是我國(guó)目前專用的地基空間物體探測(cè)手段十分有限,需要不斷加強(qiáng)．

電離層非相干散射雷達(dá)具有發(fā)射功率強(qiáng)(大于1 MW)、天線增益高(高于40 dB)、系統(tǒng)噪聲溫度低(約100 K)等特點(diǎn),是目前地基電離層觀測(cè)的強(qiáng)大手段,同時(shí)在空間碎片等空間物體硬目標(biāo)探測(cè)方面具有重要應(yīng)用價(jià)值．2012年初,中國(guó)電波傳播研究所在云南曲靖建成了我國(guó)首套非相干散射雷達(dá)[1],為有關(guān)實(shí)驗(yàn)研究提供了條件．

由于研究條件所限,國(guó)內(nèi)外的非相干散射雷達(dá)空間物體探測(cè)實(shí)驗(yàn)工作較少,日本、歐洲和俄羅斯等學(xué)者均開(kāi)展過(guò)有關(guān)研究[2-5]．歐洲非相干散射科學(xué)聯(lián)合會(huì)(European Incoherent Scatter Scientific Association,EISCAT)從1990年后期開(kāi)始全面系統(tǒng)地開(kāi)展了非相干散射雷達(dá)的空間碎片探測(cè)實(shí)驗(yàn)研究,研制了專門(mén)的空間碎片信號(hào)處理與數(shù)據(jù)終端機(jī),實(shí)現(xiàn)了對(duì)低地球軌道小尺寸空間碎片的波束駐留探測(cè)能力,近年來(lái)以Markkanen為代表不斷改進(jìn)探測(cè)實(shí)驗(yàn)、數(shù)據(jù)分析和處理方法,推動(dòng)了非相干散射雷達(dá)和空間物體探測(cè)技術(shù)的發(fā)展[6-7]．國(guó)內(nèi)相關(guān)工作甚少,中國(guó)電波傳播研究所參照EISCAT的工作,詳細(xì)介紹了非相干散射雷達(dá)空間物體探測(cè)中的匹配濾波方法,并利用EISCAT非相干散射雷達(dá)原始采樣數(shù)據(jù),提取了空間碎片距離、速度、目標(biāo)散射截面等參數(shù),初步分析了北極地區(qū)空間碎片參數(shù)的統(tǒng)計(jì)分布特征[8-11]．

非相干散射雷達(dá)的空間物體回波被接收后,通過(guò)匹配濾波方法[9-10]測(cè)量接收回波脈沖相對(duì)發(fā)射脈沖的時(shí)延來(lái)計(jì)算空間物體距離,其精度與信號(hào)采樣時(shí)間間隔、信噪比等有關(guān),比如對(duì)1 us采樣間隔,其對(duì)應(yīng)的距離測(cè)量精度約150 m,這對(duì)空間物體的精確定軌來(lái)說(shuō)是不夠的．實(shí)際上發(fā)射和接收脈沖的前后沿常被噪聲和干擾信號(hào)所污染,發(fā)射機(jī)自身相位漂移等都對(duì)脈沖相位和幅度產(chǎn)生影響,這些給基于發(fā)射-接收脈沖前沿/后沿時(shí)延差的距離測(cè)量帶來(lái)了更多誤差．

對(duì)于較大尺寸的空間物體目標(biāo),其回波信噪比較高,此時(shí)不需相干積累即可從單個(gè)回波脈沖實(shí)現(xiàn)空間物體的檢測(cè)和參數(shù)提?。硗夥窍喔缮⑸淅走_(dá)常采用相位編碼脈沖(如交替碼、巴克碼),在單個(gè)相位編碼脈沖內(nèi)存在多次相位翻轉(zhuǎn)(從-180°變?yōu)?80°稱為正相位翻轉(zhuǎn),反之為負(fù)相位翻轉(zhuǎn)),每次發(fā)射-接收相位翻轉(zhuǎn)的時(shí)間差可給出一個(gè)距離測(cè)量,因此從單個(gè)相位編碼脈沖即可得到多個(gè)距離測(cè)量及其誤差估計(jì)．鑒于此,Markkanen[7]提出了基于發(fā)射-接收回波脈沖各相位翻轉(zhuǎn)時(shí)間差的高精度距離測(cè)量方法．本文利用自主的實(shí)驗(yàn)數(shù)據(jù)(2015年10月利用EISCAT 非相干散射雷達(dá)開(kāi)展了自主的空間碎片實(shí)驗(yàn)),國(guó)內(nèi)首次對(duì)該方法進(jìn)行了詳細(xì)分析和驗(yàn)證,同時(shí)探討了該方法的不足．

1 信號(hào)處理與分析方法

1.1 脈沖校準(zhǔn)

為了精確獲得非相干散射雷達(dá)發(fā)射脈沖與接收回波脈沖內(nèi)各相位翻轉(zhuǎn)之間的時(shí)間差,首先需要對(duì)發(fā)射和接收脈沖進(jìn)行校準(zhǔn)．

發(fā)射脈沖幅度通常會(huì)出現(xiàn)輕微變化(EISCAT實(shí)驗(yàn)發(fā)現(xiàn),發(fā)射脈沖幅度是線性下降的,通常約百分之幾),同時(shí)相位漂移對(duì)發(fā)射脈沖相位具有一定影響,這等效為一個(gè)額外的多普勒頻移．為了校準(zhǔn)這些影響,對(duì)發(fā)射脈沖進(jìn)行如下模擬:

z(t)=A·(1-βt)·E(t)·exp(2πivdt)+ε(t) .

(1)

式中:A是復(fù)幅度;β是實(shí)常數(shù);vd為等效多普勒頻移;E(t)為理想交替碼的調(diào)制包絡(luò)(取值為1或-1);ε(t)為噪聲．由非線性最小二乘擬合法可確定以上模型系數(shù)(在擬合時(shí)靠近相位翻轉(zhuǎn)處的數(shù)據(jù)點(diǎn)被忽略),當(dāng)模型參數(shù)確定后,即可得到理想的交替碼發(fā)射脈沖．

對(duì)于接收回波脈沖的校準(zhǔn),首先利用通常的匹配濾波方法[7-8]確定空間物體的距離和速度初值,此處的距離初值用于從原始數(shù)據(jù)中提取空間物體回波脈沖數(shù)據(jù),速度初值用于回波脈沖的多普勒頻移校準(zhǔn)．回波脈沖具有相似的調(diào)制包絡(luò)以及幅度相位變化特征,但是它具有目標(biāo)運(yùn)動(dòng)產(chǎn)生的多普勒頻移．可采用類(lèi)似辦法對(duì)接收回波脈沖進(jìn)行校準(zhǔn)．由于回波信噪比很高,因此可以精確地校準(zhǔn)其幅度和相位變化．

1.2 脈沖相位翻轉(zhuǎn)

發(fā)射或接收回波脈沖的相位發(fā)生翻轉(zhuǎn)表明相鄰采樣數(shù)據(jù)的相位發(fā)生了顯著變化,為此需定義一個(gè)相位變化的門(mén)限值(如1弧度),超過(guò)此門(mén)限表明可能發(fā)生了相位翻轉(zhuǎn),由此初步確定發(fā)射或接收回波脈沖中的相位翻轉(zhuǎn)位置．圖1為EISCAT非相干散射雷達(dá)發(fā)射脈沖和接收脈沖采樣實(shí)部,脈沖寬度為1 920us,采樣時(shí)延間隔為1us,因此數(shù)據(jù)采樣點(diǎn)數(shù)為1 920．由于采用了交替碼相位調(diào)制,在一個(gè)脈沖內(nèi)出現(xiàn)了多次相位翻轉(zhuǎn),圖1中共出現(xiàn)27次相位翻轉(zhuǎn)．

(a) 發(fā)射脈沖實(shí)部采樣

(b) 接收脈沖實(shí)部采樣圖1 非相干散射雷達(dá)發(fā)射和接收脈沖采樣

1.3 脈沖相位翻轉(zhuǎn)時(shí)間的確定

雷達(dá)接收機(jī)實(shí)際上是一個(gè)濾波器,其濾波特性可由脈沖響應(yīng)函數(shù)來(lái)表示,接收機(jī)輸出信號(hào)就是脈沖響應(yīng)和回波信號(hào)的卷積．對(duì)相位翻轉(zhuǎn)來(lái)說(shuō),濾波過(guò)程就相當(dāng)于脈沖響應(yīng)和階梯函數(shù)的卷積(即濾波器的階梯響應(yīng)),從而使得濾波后的相位偏移呈一定斜坡度變化,而不是尖銳變化．?dāng)?shù)據(jù)采樣間隔應(yīng)略小于脈沖響應(yīng)的長(zhǎng)度,這樣保證至少一個(gè)數(shù)據(jù)點(diǎn)位于斜坡曲線內(nèi),稱為斜坡點(diǎn),此斜坡點(diǎn)對(duì)應(yīng)的時(shí)間即為相位翻轉(zhuǎn)時(shí)間．對(duì)理想的接收機(jī)脈沖響應(yīng),正負(fù)相位翻轉(zhuǎn)中的斜度點(diǎn)將關(guān)于零對(duì)稱．圖2為非相干散射雷達(dá)接收機(jī)脈沖響應(yīng)和正負(fù)相位翻轉(zhuǎn)的階梯響應(yīng),可見(jiàn)在脈沖響應(yīng)最大值附近的相位翻轉(zhuǎn)階梯響應(yīng)具有最大斜率．從非相干散射雷達(dá)實(shí)測(cè)數(shù)據(jù)確定相位翻轉(zhuǎn)時(shí)間的具體步驟是:1)首先確定發(fā)射或接收回波脈沖中相位變化(兩個(gè)相鄰采樣數(shù)據(jù)相位之差的絕對(duì)值)超過(guò)一定門(mén)限的數(shù)據(jù)點(diǎn),搜索這些數(shù)據(jù)點(diǎn)附近幅度變化最大的數(shù)據(jù)點(diǎn),視為初值斜坡點(diǎn)．2)在階梯響應(yīng)函數(shù)曲線上搜索與此初值斜坡點(diǎn)幅度對(duì)應(yīng)的時(shí)間,視為修正量．3)初值斜坡點(diǎn)與修正量一起確定了最終的斜坡點(diǎn)時(shí)間,也就是最終的相位翻轉(zhuǎn)時(shí)間．

(a) 接收機(jī)脈沖響應(yīng)

(b) 相位翻轉(zhuǎn)對(duì)應(yīng)的階梯脈沖響應(yīng)圖2 非相干散射雷達(dá)接收機(jī)脈沖響應(yīng)和 相位翻轉(zhuǎn)階梯響應(yīng)

下面分析相位翻轉(zhuǎn)時(shí)間的誤差,它包括系統(tǒng)誤差(階梯響應(yīng)函數(shù)模型誤差)和隨機(jī)噪聲誤差．假設(shè)接收機(jī)濾波器的階梯響應(yīng)為p(t),相位翻轉(zhuǎn)時(shí)間τs對(duì)應(yīng)的信號(hào)值為s,設(shè)時(shí)延τ和信號(hào)值s為隨機(jī)變量,ε為零均值高斯噪聲,RSN為信噪比,σε為噪聲標(biāo)準(zhǔn)偏差,則

s=p(τ)+ε.

(2)

(3)

當(dāng)s已知時(shí),τ滿足條件概率密度分布,可表示為

(4)

當(dāng)信噪比很高時(shí)(通常對(duì)于大尺寸空間物體)σε很小,D(τ|s)具有窄的尖峰．這里στ為相位翻轉(zhuǎn)時(shí)間的標(biāo)準(zhǔn)偏差,可表示為

(5)

式(5)表明已知信噪比時(shí),若斜坡點(diǎn)時(shí)間對(duì)應(yīng)于脈沖響應(yīng)最大值,則可高精度地確定相位翻轉(zhuǎn)時(shí)間．這也對(duì)應(yīng)于階梯響應(yīng)的最大斜率處,此時(shí)階梯響應(yīng)變化引起的相位翻轉(zhuǎn)時(shí)間變化量最?。?/p>

1.4 距離測(cè)量

(6)

假定在單個(gè)脈沖反射期間(脈沖寬度為1 920us)的空間物體視線速度不變,則單個(gè)脈沖反射期的空間物體距離隨時(shí)間滿足線性變化,對(duì)單個(gè)脈沖內(nèi)多次相位翻轉(zhuǎn)對(duì)應(yīng)的多個(gè)距離測(cè)量值進(jìn)行加權(quán)最小二乘擬合,通過(guò)擬合可得到一個(gè)距離及誤差估計(jì)值,此距離估計(jì)值相對(duì)于以往的接收-發(fā)射脈沖時(shí)延法來(lái)說(shuō)精度具有很大提高．

2 實(shí)測(cè)數(shù)據(jù)分析與討論

2015年10月我們利用EISCAT500MHz非相干散射雷達(dá)(海拔445m,地理經(jīng)度16°01'E,地理緯度78°09'N)開(kāi)展了首次北極空間物體凝視觀測(cè)實(shí)驗(yàn),實(shí)驗(yàn)中采用64位交替碼,碼元寬度30us,脈沖寬度1 920us,脈沖重復(fù)周期為20ms,數(shù)據(jù)采樣率為1MHz,雷達(dá)波束方位角和仰角分別固定為90°和75°,天線增益約42.5dB,半功率波束寬度為1.1°,系統(tǒng)噪聲溫度約70K,每個(gè)脈沖周期內(nèi)所有數(shù)據(jù)均被連續(xù)采樣,這樣單個(gè)重復(fù)周期內(nèi)的數(shù)據(jù)為20 000個(gè)復(fù)采樣點(diǎn),接收機(jī)濾波器脈沖響應(yīng)為高斯型,其脈沖長(zhǎng)度為2.2us．

圖3為某個(gè)脈沖內(nèi)的原始采樣數(shù)據(jù),從中可見(jiàn)在第4 300個(gè)采樣位置處(通過(guò)匹配濾波檢測(cè),計(jì)算其距離約644.40km)出現(xiàn)強(qiáng)散射回波,其信噪比約130,可認(rèn)為此處存在高信噪比的空間物體目標(biāo)．

圖3 單個(gè)脈沖周期內(nèi)的采樣數(shù)據(jù)

圖4為單個(gè)發(fā)射和接收脈沖內(nèi)的27個(gè)斜坡點(diǎn)的位置修正量,其中圓圈代表正相位翻轉(zhuǎn),星號(hào)為負(fù)相位翻轉(zhuǎn)．

(a) 發(fā)射脈沖斜坡點(diǎn)

(b) 接收脈沖斜坡點(diǎn)圖4 發(fā)射和接收脈沖的斜坡點(diǎn)位置

從圖4中27個(gè)相位翻轉(zhuǎn)可計(jì)算27個(gè)距離測(cè)量值,對(duì)這些測(cè)量值采用加權(quán)最小二乘擬合,見(jiàn)圖5,可得到此脈沖內(nèi)的空間物體距離擬合值及其誤差估計(jì)．橫坐標(biāo)零時(shí)刻表示具有最小擬合誤差的某相位翻轉(zhuǎn)對(duì)應(yīng)的空間物體反射時(shí)間(它對(duì)應(yīng)最終距離),縱軸表示測(cè)量值相對(duì)最終距離的偏差,不帶誤差棒表示該測(cè)量值對(duì)應(yīng)的斜坡點(diǎn)離階梯響應(yīng)最大斜率處較遠(yuǎn)(即誤差較大),擬合時(shí)舍棄不帶誤差棒的數(shù)據(jù),由此可得該脈沖對(duì)應(yīng)的空間物體最終距離(644 337.4m)和誤差估計(jì)(18.50m)．

圖5 由單個(gè)脈沖內(nèi)各相位翻轉(zhuǎn)計(jì)算的距離及其誤差

圖6為采用以往的匹配濾波方法[9-10]分析上文相同數(shù)據(jù)得到每個(gè)脈沖內(nèi)的信噪比和空間物體距離,脈沖寬度1 920us,脈沖重復(fù)周期為20ms．空間物體穿越雷達(dá)波束期間的回波信噪比和距離一般滿足平穩(wěn)、連續(xù)變化特征,我們?nèi)】臻g物體穿越整個(gè)波束期間的最大信噪比時(shí)刻(圖6(a)中圓圈表示)對(duì)應(yīng)的距離643.35km作為空間物體距離估計(jì)值．由于單個(gè)脈沖重復(fù)周期內(nèi)僅得到一個(gè)距離測(cè)量值,因此無(wú)法給出其誤差估計(jì)．此處計(jì)算信噪比時(shí),我們利用每個(gè)脈沖重復(fù)周期內(nèi)遠(yuǎn)離空間物體回波的采樣數(shù)據(jù)(數(shù)據(jù)長(zhǎng)度為1 920點(diǎn))計(jì)算背景噪聲．由于原始數(shù)據(jù)中存在干擾信號(hào)等原因,影響了對(duì)空間物體的檢測(cè)、信噪比和距離計(jì)算,甚至出現(xiàn)目標(biāo)的錯(cuò)誤檢測(cè)(虛警)和參數(shù)的錯(cuò)誤計(jì)算．

(a) 不同脈沖重復(fù)周期對(duì)應(yīng)的信噪比

(b) 不同脈沖重復(fù)周期對(duì)應(yīng)的空間物體距離圖6 采用以往匹配濾波方法提取的空間物體參數(shù)

從圖4～6可知,以往匹配濾波方法利用發(fā)射-接收脈沖的時(shí)間差作為空間物體距離測(cè)量值,不能對(duì)其誤差進(jìn)行定量估計(jì),其測(cè)量誤差由脈沖波形、信噪比和采樣間隔等決定,一般為數(shù)百米以上．而本文方法首先利用單個(gè)脈沖內(nèi)各發(fā)射-接收相位翻轉(zhuǎn)時(shí)間差計(jì)算其對(duì)應(yīng)的空間物體距離及其誤差,對(duì)高的信噪比其均方根誤差為十米量級(jí),因此本文方法不僅測(cè)量精度高,而且可給出定量的誤差估計(jì)．

以上結(jié)果表明本方法得到的距離精度可達(dá)數(shù)十米量級(jí),遠(yuǎn)優(yōu)于以往的基于發(fā)射-接收脈沖前沿時(shí)間差的距離測(cè)量方法．但是若將本方法用于空間物體精密定軌和碰撞預(yù)警等,尚需進(jìn)一步考慮和修正以下誤差．首先，由于非相干散射雷達(dá)主要用于電離層觀測(cè),對(duì)雷達(dá)站地理坐標(biāo)、天線三軸中心等精度要求不高,這直接影響了對(duì)空間物體距離、位置等測(cè)量精度．其次,空間物體散射信號(hào)在電離層和接收機(jī)系統(tǒng)里的傳播時(shí)延等可能導(dǎo)致數(shù)十米以上誤差，特別是在電離層劇烈擾動(dòng)和較低仰角探測(cè)時(shí),雷達(dá)至空間物體的信號(hào)傳播路徑不再是直線,由電離層折射產(chǎn)生的距離誤差更大,可到數(shù)百米．再次,有一些因素會(huì)影響到距離擬合誤差，比如,假定信號(hào)模型中幅度為常數(shù)或規(guī)則變化,這僅適用于無(wú)尺寸點(diǎn)目標(biāo)．較大尺寸的不規(guī)則目標(biāo)具有多個(gè)強(qiáng)的散射中心,干涉效應(yīng)會(huì)對(duì)回波脈沖幅度和相位產(chǎn)生調(diào)制．此外,本文方法的實(shí)際距離測(cè)量精度還需結(jié)合空間物體精確軌道數(shù)據(jù)進(jìn)行分析驗(yàn)證．

3 結(jié) 論

非相干散射雷達(dá)是目前地基電離層監(jiān)測(cè)的強(qiáng)大手段,在空間碎片等探測(cè)方面也具有重要應(yīng)用潛力．以往的非相干散射雷達(dá)空間碎片探測(cè)研究中,主要利用發(fā)射和接收脈沖前后沿時(shí)間差來(lái)測(cè)量距離,其測(cè)量精度受脈沖前/后沿波形、采樣間隔等影響,測(cè)量精度一般為數(shù)百米量級(jí)．

本文參照EISCAT的最新工作,結(jié)合自主實(shí)測(cè)數(shù)據(jù)詳細(xì)介紹了基于相位編碼脈沖內(nèi)多次相位翻轉(zhuǎn)的高信噪比空間物體距離測(cè)量方法．該方法基于兩點(diǎn)考慮:1)相位編碼脈沖內(nèi)存在多次相位翻轉(zhuǎn),因此基于發(fā)射-接收脈沖之間的多次相位翻轉(zhuǎn)時(shí)間差測(cè)量可實(shí)現(xiàn)多次距離測(cè)量．2)假定單個(gè)相位編碼脈沖(本文脈沖寬度為1 920us)內(nèi)空間物體速度為常數(shù),其距離隨時(shí)間線性變化．因此對(duì)單個(gè)脈沖內(nèi)的多次距離測(cè)量進(jìn)行最小二乘擬合,可得到單個(gè)脈沖內(nèi)的最佳距離估計(jì)值和誤差估計(jì)值．

本方法可一定程度克服以往方法中由于噪聲和干擾所致脈沖波形畸變引起的較大誤差,距離誤差可達(dá)數(shù)十米量級(jí),這對(duì)利用我國(guó)曲靖非相干散射雷達(dá)開(kāi)展空間物體的高精度探測(cè),提升我國(guó)的空間物體監(jiān)測(cè)預(yù)警能力具有重要現(xiàn)實(shí)指導(dǎo)意義．下一步,需結(jié)合空間物體精確軌道數(shù)據(jù)進(jìn)一步分析非相干散射雷達(dá)的空間物體測(cè)量誤差以及修正方法,以期逐步提高非相干散射雷達(dá)的空間物體參數(shù)實(shí)際測(cè)量精度．

[1] 丁宗華, 代連東, 董明玉,等.非相干散射雷達(dá)進(jìn)展:從傳統(tǒng)體制到EISCAT3D[J].地球物理學(xué)進(jìn)展,2014,29(5): 2376-2381.

DINGZH,DAILD,DONGMY,etal.Progressoftheincoherentscatteringradar:fromthetraditionalradartothelatestEISCAT3D[J].Progressingeophysics, 2014, 29(5): 2376-2381. (inChinese)

[2]SATOT,IKEDAK,KIMURAI,etal.Shapeofspacedebrisestimatedfromradarcrosssectionvariation[J].Journalofspacecraftandrocket, 1994, 31(4): 665-670.

[3]MARKKANENJ,LEHTINENM,HUUSKONENA,etal.Measurementsofsmall-sizedebriswithbackscatterofradiowaves[R].Kiruna:EISCAT, 2002.

[4]MARKKANENJ,LEHTINENM,LANDGRAFM.RealtimespacedebrismonitoringwithEISCAT[J].Advancesinspaceresearch, 2005, 35: 1197-1209.

[5]MARKKANENJ,POSTILAM,VANEYKENA.SmallsizespacedebrisdatacollectionwithEISCATradarfacilities[R].Kiruna:EISCAT, 2006.

[6]NYGRENT,MARKKANENJ,AIKIOA,etal.High-precisionmeasurementofsatellitevelocityusingtheEISCATradar[J].Annalesgeophysicae, 2012, 30: 1555-1565.

[7]MARKKANENJ,NYGRENT,MARKKANENM,etal.High-precisionmeasurementofsatelliterangeandvelocityusingtheEISCATradar[J].Annalesgeophysicae, 2013, 31: 859-870.

[8] 金旺, 吳振森, 吳健, 等.非相干散射雷達(dá)探測(cè)空間碎片實(shí)驗(yàn)研究[J].電波科學(xué)學(xué)報(bào), 2011, 26(5): 956-960.

JINW,WUZS,WUJ,etal.Spacedebrisexperimentstudyusingtheincoherentscatterradar[J].Chinesejournalofradioscience, 2011, 26(5): 956-960. (inChinese)

[9] 丁宗華, 許正文, 趙振維. 非相干散射雷達(dá)的空間碎片參數(shù)提取[J]. 現(xiàn)代雷達(dá),2014, 36(2): 1-5.

DINGZH,XUZW,ZHAOZW.Theparameterestimationofspacedebrisfromtheincoherentscatterradar[J].Modernradar, 2014, 36(2): 1-5. (inChinese)

[10]丁宗華, 代連東, 趙振維, 等.非相干散射雷達(dá)的空間碎片信號(hào)處理與參數(shù)估算方法研究[J]. 電波科學(xué)學(xué)報(bào), 2014, 29(3): 417-423.

DINGZH,DAILD,ZHAOZW,etal.Theinvestigationofspacedebris'ssignalprocessingandparameterestimationbasedontheincoherentscatteringradar[J].Chinesejournalofradioscience, 2014, 29(3): 417-423.(inChinese)

[11]丁宗華, 代連東, 付佗, 等. 曲靖非相干散射雷達(dá)空間碎片初步探測(cè)結(jié)果[C]//第八屆全國(guó)空間碎片學(xué)術(shù)交流會(huì)論文集， 北京, 2015.

DINGZH,DAILD,FUT,etal.ThepreliminaryresultsofthespacedebrisobservationbytheQujingincoherentscatterradar[C]//Proceedingsofthe8thspacedebris,Beijing, 2015.(inChinese)

The range measurement of space object with high signal to noise ratio based on the incoherent scatter radar

DING Zonghua TANG Zhimei DAI Liandong WU Jian XU Zhengwen

(NationalKeyLaboratoryofElectromagneticEnvironment,ChinaResearchInstituteofRadiowavePropagation,Qingdao266107,China)

The incoherent scatter radar(ISR) is a main tool for monitoring the ionosphere on the ground, and can also be used to detect the space objects. When the signal to noise ratio is a bit high, only one pulse is needed to detect the space object and the coherent integration is unnecessary. There usually exist many phase flips in the radar phase-modulated pulses. When the echo pulse with phase flips is filtered by the receiver, the resulting phase flips are not sharp but consist of many slopes. By investigating the echo amplitude variation within this slope, the time difference between the transmission and reception sub-pulse and then the range measurements are known. Then all these range measurements during one pulse are applied to a weighted least-square fitting, which gives the range and its error estimation. By analyzing the EISCAT raw data, it is found that the range error is about several tens of meters,which is far better than that from the frequently-used method based on time difference between the transmission and reception pulse.

incoherent scatter radar; space object; high signal to noise ratio; range; phase flip

10.13443/j.cjors.2016092702

2016-09-27

十二五民用航天預(yù)研基金

TN958

A

1005-0388(2016)06-1081-06

丁宗華 (1978-),男,湖北人,博士,高級(jí)工程師,研究方向?yàn)殡婋x層非相干散射探測(cè)與電波傳播．

唐志美 (1993-),女,貴州人， 碩士研究生,研究方向?yàn)榭臻g碎片信號(hào)處理與數(shù)據(jù)分析．

代連東 (1988-),男,河北人，碩士,工程師,研究方向?yàn)榉窍喔缮⑸涮綔y(cè)與實(shí)驗(yàn)．

丁宗華, 唐志美, 代連東,等. 非相干散射雷達(dá)的高信噪比空間物體距離測(cè)量方法 [J].電波科學(xué)學(xué)報(bào),2016,31(6):1081-1086.

DING Z H, TANG Z M, DAI L D,et al. The range measurement of space object with high signal to noise ratio based on the Incoherent Scatter Radar [J].Chinese journal of radio science,2016,31(6):1081-1086.(in Chinese).DOI:10.13443/j.cjors.2016092702

聯(lián)系人： 丁宗華 E-mail：zhdingmou@163.com

DOI 10.13443/j.cjors.2016092702