沈愛(ài)國(guó),姜秋喜
(1.電子工程學(xué)院雷達(dá)對(duì)抗系, 合肥230037;2.安徽省電磁制約重點(diǎn)實(shí)驗(yàn)室, 合肥230037)
合成孔徑雷達(dá)(Synthetic Aperture Radar,SAR)在軍事領(lǐng)域應(yīng)用范圍不斷擴(kuò)大,對(duì)SAR干擾技術(shù)的研究逐漸成為電子對(duì)抗領(lǐng)域的研究熱點(diǎn)[1-8]。為了滿足不斷提高的分辨率要求,SAR雷達(dá)工作的信號(hào)帶寬也越來(lái)越寬,已經(jīng)達(dá)到幾百M(fèi)Hz甚至GHz量級(jí)。傳統(tǒng)的雷達(dá)干擾技術(shù)在干擾信號(hào)產(chǎn)生過(guò)程中需要利用所截獲雷達(dá)脈沖的全部信息,干擾信號(hào)帶寬至少與雷達(dá)發(fā)射信號(hào)帶寬相等,也就是說(shuō)干擾信號(hào)的瞬時(shí)帶寬不小于雷達(dá)信號(hào)帶寬。若利用傳統(tǒng)的干擾技術(shù)對(duì)寬帶SAR實(shí)施干擾,不僅對(duì)干擾機(jī)發(fā)射設(shè)備提出了高要求,同時(shí)也有可能對(duì)己方的電子設(shè)備造成影響,帶來(lái)電磁兼容等方面的問(wèn)題??梢?jiàn),要對(duì)寬帶SAR實(shí)施有效干擾,常規(guī)的干擾技術(shù)已經(jīng)難以奏效,必須尋求新的干擾方法。
文獻(xiàn)[9-10]指出帶寬大于SAR工作帶寬10%的噪聲干擾在處理階段就很難將其去除而不對(duì)原有圖像產(chǎn)生顯著影響,帶寬大于SAR工作帶寬50%的欺騙干擾信號(hào)即可產(chǎn)生接近于SAR分辨率的欺騙干擾圖像,在不對(duì)圖像進(jìn)行精確測(cè)量的圖像判讀階段,能夠達(dá)到一定的欺騙干擾效果。因此,對(duì)于大帶寬的SAR,可以考慮降低帶寬的干擾。文獻(xiàn)[11]指出,在干擾雷達(dá)的波形中,噪聲干擾是效果最差的干擾波形。與信號(hào)完全不相關(guān)的噪聲得不到任何相關(guān)增益,而全相關(guān)的目標(biāo)能得到全部處理增益。即便有一點(diǎn)相關(guān)性的干擾波形也比隨機(jī)噪聲要有效。根據(jù)這一結(jié)論,文獻(xiàn)[12]提出了子帶脈沖循環(huán)轉(zhuǎn)發(fā)干擾技術(shù)并分析了該干擾技術(shù)對(duì)SAR的壓制干擾效果。子帶脈沖循環(huán)轉(zhuǎn)發(fā)干擾在一定程度上提高了干擾脈沖的密集程度,改善了子帶脈沖干擾的干擾效果。但循環(huán)轉(zhuǎn)發(fā)的最大問(wèn)題在于其轉(zhuǎn)發(fā)的子帶脈沖之間最小間隔應(yīng)不小于子帶脈沖的寬度,限制了提高子帶脈沖干擾信號(hào)密集程度的能力。本文研究了利用子帶脈沖對(duì)梳狀譜信號(hào)進(jìn)行調(diào)制的干擾技術(shù)。
雷達(dá)發(fā)射的線性調(diào)頻(LFM)脈沖信號(hào)可以表示為
式中:fc為載頻;Tp為脈沖寬度。則雷達(dá)的基帶信號(hào)可以表示為
設(shè)子帶脈沖寬度為τ,采樣起始時(shí)刻為t1,則子帶脈沖信號(hào)可以表示為
可以選擇兩種極端的情況來(lái)討論子帶脈沖信號(hào)的基本特征,截取LFM信號(hào)兩端,可以得到的子帶脈沖信號(hào)分別為
圖1a)~圖1c)分別給出了式(2)的雷達(dá)發(fā)射脈沖基帶信號(hào)和式(4)定義的兩個(gè)子帶脈沖信號(hào)的時(shí)域波形,其中雷達(dá)信號(hào)脈沖寬度T0=20 μs,子帶脈沖寬度τ=Tp/4,雷達(dá)信號(hào)帶寬為200 MHz,仿真時(shí)間長(zhǎng)度為40 μs。從圖中可以明顯看出子帶脈沖只是雷達(dá)發(fā)射脈沖的一部分,保持了雷達(dá)發(fā)射脈沖的幅度特征和變化規(guī)律。為了便于對(duì)子帶脈沖信號(hào)的特性加深理解,圖1d)~圖1f)中同時(shí)給出了上述時(shí)域波形對(duì)應(yīng)信號(hào)的時(shí)頻分布特性,可以明顯看出子帶脈沖信號(hào)的時(shí)頻分布特性保持了原始雷達(dá)脈沖信號(hào)的時(shí)頻分布規(guī)律,其差別僅僅是子帶脈沖信號(hào)所占據(jù)的時(shí)間寬度和頻率寬度都只是雷達(dá)發(fā)射脈沖信號(hào)的一部分。從子帶脈沖的信號(hào)表達(dá)式以及其時(shí)域波形和時(shí)頻分布特性均可以得出結(jié)論:子帶脈沖信號(hào)仍是一個(gè)線性調(diào)頻信號(hào),其調(diào)頻斜率與雷達(dá)發(fā)射的脈沖信號(hào)一致,信號(hào)的脈沖寬度和信號(hào)帶寬決定于參數(shù),中心頻率決定于子帶脈沖截取的起始位置和脈沖寬度。本文只是對(duì)子帶脈沖信號(hào)的基本特性進(jìn)行簡(jiǎn)單闡釋,子帶脈沖與干擾效果緊密相關(guān)的更多特性與其匹配濾波輸出有關(guān)。
圖1 子帶寬信號(hào)波形及時(shí)頻分布
文獻(xiàn)[13]通過(guò)嚴(yán)格的數(shù)學(xué)推導(dǎo)得到式(2)所示的具有矩形包絡(luò)的線性調(diào)頻信號(hào)頻譜函數(shù)的精確表達(dá)式,根據(jù)匹配濾波器的定義得到對(duì)應(yīng)的匹配濾波器傳遞函數(shù)。但是嚴(yán)格意義上的傳遞函數(shù)不僅要求通帶內(nèi)幅頻特性存在菲涅耳起伏,而且要求具有復(fù)雜的相頻特性,很難綜合出來(lái)。因此,文獻(xiàn)[13]提出了利用全通的相位校正網(wǎng)絡(luò)代替匹配濾波器進(jìn)行分析,由于矩形包絡(luò)的線性調(diào)頻信號(hào)具有近似矩形的頻譜和平方率的相頻特性,滿足要求的全通相位校正網(wǎng)絡(luò)的傳遞函數(shù)的幅頻特性和相頻特性為
由于全通相位校正網(wǎng)絡(luò)只對(duì)信號(hào)相位有校正作用,故可以看成是一相位均衡器。根據(jù)式(5)可以得到該相位校正網(wǎng)絡(luò)的沖擊響應(yīng)為
將式(6)與式(2)的信號(hào)求卷積,可以得到輸出信號(hào)為
得到輸出信號(hào)的幅度系數(shù)為
比較式(7)和式(9)可以看出,二者之間具有一定的差別。首先是兩者的包絡(luò)有差別,如果線性調(diào)頻性的時(shí)帶積較大,則壓縮后時(shí)寬只占Tp的很小部分。如果只考慮峰值附近很小t值范圍內(nèi)的形狀,滿足t?Tp,則二者就近似相等。其次,用相位校正網(wǎng)絡(luò)壓縮線性調(diào)頻信號(hào),將帶來(lái)殘留線性調(diào)頻項(xiàng)e-jπkt2,該殘留調(diào)頻的斜率和原先信號(hào)大小相等,符號(hào)相反。這個(gè)效應(yīng)只有在時(shí)帶積D較小時(shí),才對(duì)壓縮信號(hào)中的載波有可觀測(cè)的影響。在時(shí)帶積較大時(shí),由于壓縮信號(hào)持續(xù)時(shí)間較窄,殘留項(xiàng)可以忽略不計(jì)。
上述分析表明,全通相位均衡網(wǎng)絡(luò)可以作為對(duì)線性調(diào)頻信號(hào)匹配濾波器的良好近似,特別是線性調(diào)頻信號(hào)的時(shí)帶積較大時(shí)更是如此。寬帶SAR具有很大的時(shí)間帶寬積,適合于用相位均衡網(wǎng)絡(luò)對(duì)信號(hào)的匹配濾波輸出進(jìn)行分析。另外,由于式(3)所示的子帶脈沖信號(hào)起始時(shí)刻可以隨意選取,子帶脈沖信號(hào)具有不對(duì)稱特性,用全通的相位校正網(wǎng)絡(luò)進(jìn)行近似分析更為合適。將式(3)與式(6)卷積,可以得到子帶脈沖信號(hào)的匹配濾波輸出為
對(duì)式(10)進(jìn)行變換,可以得到
根據(jù)式(11)可以得到子帶脈沖信號(hào)匹配濾波輸出信號(hào)的幅度函數(shù)為
比較式(8)、式(12)的幅度系數(shù)可以看出,子帶脈沖信號(hào)在經(jīng)過(guò)匹配濾波器之后其包絡(luò)仍具有sinc函數(shù)形式,且峰值位置與雷達(dá)信號(hào)匹配濾波輸出的峰值位置一致,在t=0時(shí)達(dá)到最大值;子帶脈沖信號(hào)的最大值為原始雷達(dá)信號(hào)的τ/Tp,即子帶脈沖信號(hào)的匹配濾波輸出信號(hào)幅度與子帶信號(hào)的寬度成正比。進(jìn)一步研究子帶脈沖輸出信號(hào)的寬度,由于雷達(dá)信號(hào)匹配濾波輸出的第一零點(diǎn)出現(xiàn)在t=1/kTp時(shí)刻,根據(jù)式(10)可以得到子帶脈沖信號(hào)匹配濾波輸出的第一零點(diǎn)出現(xiàn)在t=1/kτ時(shí)刻,通常定義該值為輸出脈沖的主瓣寬度,即雷達(dá)信號(hào)和子帶脈沖信號(hào)匹配濾波輸出信號(hào)的主瓣寬度分別為1/kTp和1/kτ。
圖2給出了圖1中雷達(dá)發(fā)射脈沖和兩個(gè)子帶脈沖信號(hào)匹配濾波輸出波形,該處的輸出脈沖信號(hào)幅度以雷達(dá)信號(hào)輸出為基準(zhǔn)進(jìn)行了歸一化處理。從仿真結(jié)果可以看出:脈沖壓縮輸出峰值位置保持不變,但幅度降低了,對(duì)信號(hào)輸出峰值歸一化后為0.25;子帶脈沖匹配濾波輸出的主瓣寬度也展寬了,雷達(dá)發(fā)射脈沖濾波輸出信號(hào)的第一零點(diǎn)位于t=0.005 μs,對(duì)應(yīng)的主瓣寬度等于信號(hào)帶寬的倒數(shù);而子帶脈沖輸出信號(hào)的第一零點(diǎn)位于t=0.002 μs,正好等于子帶寬脈沖所占頻帶寬度的倒數(shù)。與理論分析結(jié)果完全一致。
圖2 子帶脈沖信號(hào)匹配濾波輸出
最簡(jiǎn)單的梳狀譜信號(hào)可以理解為一定頻帶范圍內(nèi)多個(gè)單頻信號(hào)的線性組合,表示為
式中:N為單頻信號(hào)的數(shù)目;[f1,f2]為信號(hào)的頻率范圍。利用式(3)的子帶脈沖信號(hào)對(duì)該梳狀譜信號(hào)進(jìn)行調(diào)制,得到
子帶脈沖調(diào)制梳狀譜信號(hào)的實(shí)際效果相當(dāng)于將子帶脈沖信號(hào)根據(jù)梳狀譜分布進(jìn)行頻率搬移后再線性組合,因此在頻域討論更容易理解。令h(t)的頻譜為H(f),可以得到
令s1(t)的頻譜為S1(f),則可以得到調(diào)制后的信號(hào)頻譜為
即干擾信號(hào)為子帶脈沖信號(hào)頻移后的線性組合。
考慮子帶脈沖調(diào)制梳狀譜干擾信號(hào)的匹配濾波輸出特性,式(3)的子帶脈沖信號(hào)頻移fd后得到
與式(6)的沖擊響應(yīng)函數(shù)進(jìn)行卷積運(yùn)算,得到頻移后的子帶脈沖信號(hào)匹配濾波輸出為
根據(jù)上式可以得到頻移子帶脈沖信號(hào)匹配濾波輸出信號(hào)的幅度函數(shù)為
比較式(19)與式(12)可以看出,對(duì)子帶脈沖延遲fd后,只要頻移后的子帶脈沖信號(hào)在雷達(dá)信號(hào)帶寬范圍之內(nèi),其匹配濾波輸出信號(hào)幅度仍保持原子帶脈沖匹配濾波輸出特性,只是在時(shí)間上位移fd/k。為了保證頻移后的子帶脈沖信號(hào)不再發(fā)生脈壓損失,需要對(duì)移頻范圍進(jìn)行限制。假定雷達(dá)發(fā)射信號(hào)的頻率范圍為f1~f2,截取的子帶脈沖信號(hào)頻率范圍為f1≤fi1-fi2≤f2,那么頻移fd的范圍為f1-fi1~f2-fi2,對(duì)應(yīng)頻移fd的變化范圍為
根據(jù)頻移與延遲之間的關(guān)系,由上式得到移頻后的子帶脈沖壓縮峰值時(shí)間延遲為
式中:td1對(duì)應(yīng)于將子帶脈沖向低頻端調(diào)制,壓縮峰值滯后于原始子帶脈沖壓縮峰值;td2對(duì)應(yīng)于將子帶脈沖向高頻端調(diào)制,壓縮峰值超前于原始子帶脈沖壓縮峰值。對(duì)于一個(gè)選定的子帶脈沖信號(hào),其頻移范圍還將受到進(jìn)一步的限制,對(duì)應(yīng)的延時(shí)范圍也受到限制。根據(jù)式(21)可以知道,確定子帶脈沖起始時(shí)刻t1后,移頻導(dǎo)致的延時(shí)覆蓋寬度為
也就是說(shuō),對(duì)于帶寬一定的子帶脈沖信號(hào),在不影響其脈壓輸出信號(hào)幅度的情況下,對(duì)其進(jìn)行頻移產(chǎn)生的延時(shí)范圍是確定。利用該結(jié)論分析式(14)中的信號(hào)可以知道,子帶脈沖調(diào)制的梳狀譜干擾脈沖壓縮后的峰值分布于td2~td1之間,最大覆蓋寬度為Δtd。那么,假定梳狀譜信號(hào)的頻率在取值范圍內(nèi)均勻分布,則調(diào)制后的信號(hào)脈壓峰值也將在對(duì)應(yīng)范圍內(nèi)均勻分布,此時(shí)各個(gè)脈壓峰值之間的時(shí)間間隔為
為了使子帶調(diào)制梳狀譜信號(hào)脈壓后各個(gè)脈壓峰值相互之間不受影響,或者說(shuō)可以忽略相互間副瓣的影響,需要各峰值的間隔要遠(yuǎn)大于1/kτ,即要求滿足
因此,在設(shè)置梳狀譜信號(hào)時(shí),需要對(duì)頻率的點(diǎn)數(shù)和分布特征進(jìn)行考慮。當(dāng)Δti?1/kτ時(shí),可以忽略副瓣的影響,此時(shí)根據(jù)線性系統(tǒng)疊加原理可以得到式(14)中子帶調(diào)制梳狀譜干擾信號(hào)匹配濾波輸出信號(hào)幅度為
上述分析表明,子帶調(diào)制的梳狀譜干擾信號(hào),經(jīng)距離向脈沖壓縮后可以形成多個(gè)脈沖壓縮峰值。但隨著載機(jī)的移動(dòng),干擾信號(hào)的壓縮峰值分布與回波信號(hào)的彎曲特性不一致,在進(jìn)行距離遷移校正時(shí)不能得到完全校正,因此方位向的干擾信號(hào)之間不能對(duì)正,表現(xiàn)為輸出圖像的展寬。為了驗(yàn)證子帶調(diào)制梳狀譜干擾技術(shù)對(duì)寬帶SAR干擾的有效性,我們通過(guò)仿真試驗(yàn)對(duì)該干擾技術(shù)的效果進(jìn)行了定量分析。
設(shè)雷達(dá)信號(hào)的中心頻率為500 MHz,信號(hào)帶寬為200 MHz,脈沖寬度為 1 μs,脈沖重復(fù)頻率為 200 Hz。設(shè)成像場(chǎng)景為200 m×200 m的正方形區(qū)域,雷達(dá)距離場(chǎng)景的中心斜距為10 km,載機(jī)速度為100 m/s,方位向分辨率為1 m。干擾機(jī)位于場(chǎng)景的中心位置,在干擾機(jī)位置有一個(gè)理想的點(diǎn)目標(biāo)。在回波信中分別加入隨機(jī)噪聲干擾、子帶脈沖循環(huán)轉(zhuǎn)發(fā)干擾、子帶調(diào)制梳狀譜干擾信號(hào)及其循環(huán)轉(zhuǎn)發(fā)干擾信號(hào),計(jì)算點(diǎn)目標(biāo)圖像指標(biāo)在不同干信比條件下的變化情況。圖3給出了該點(diǎn)目標(biāo)圖像指標(biāo)的變化曲線,其中圖3a)為點(diǎn)目標(biāo)圖像的距離向峰值旁瓣比,圖3b)為方位向峰值旁瓣比,圖3c)為距離向積分旁瓣比,圖3d)為方位向積分旁瓣比。圖中的干擾1為隨機(jī)噪聲干擾,干擾2為子帶脈沖重復(fù)轉(zhuǎn)發(fā)干擾,干擾3為子帶調(diào)制梳狀譜干擾,干擾4為子帶調(diào)制梳狀譜重復(fù)轉(zhuǎn)發(fā)干擾。從圖中可以看出,對(duì)所有的干擾信號(hào),目標(biāo)的圖像指標(biāo)都隨著干擾信號(hào)強(qiáng)度的增加而下降,最終都達(dá)到一定的穩(wěn)定狀態(tài)。圖像指標(biāo)進(jìn)入穩(wěn)定狀態(tài)說(shuō)明圖像特征已經(jīng)完全由干擾信號(hào)決定,繼續(xù)增加干擾功率,圖像指標(biāo)不再變化。
圖3 不同干擾信號(hào)對(duì)點(diǎn)目標(biāo)圖像性能影響對(duì)比曲線
可以根據(jù)圖像指標(biāo)進(jìn)入穩(wěn)態(tài)的干信比門(mén)限值來(lái)判定不同干擾之間性能的優(yōu)劣。分析圖中所示的仿真結(jié)果可以發(fā)現(xiàn),要完全壓制點(diǎn)目標(biāo)圖像,在該仿真參數(shù)設(shè)置條件下,隨機(jī)噪聲干擾需要的干信比約為44 dB,而子帶脈沖重復(fù)轉(zhuǎn)發(fā)干擾和子帶調(diào)制梳狀譜干擾所需要的干信比幾乎相當(dāng),約為32 dB,比射頻噪聲降低12 dB左右,與理論分析的結(jié)論基本一致。將子帶調(diào)制梳狀譜干擾進(jìn)行5次重復(fù)轉(zhuǎn)發(fā),得到了圖中干擾4所示的圖像指標(biāo)變化曲線,可以看出此時(shí)干擾效果得到了進(jìn)一步加強(qiáng),這是由于重復(fù)轉(zhuǎn)發(fā)提高了相關(guān)干擾信號(hào)的密集程度。進(jìn)一步分析可以發(fā)現(xiàn),由于仿真中重復(fù)轉(zhuǎn)發(fā)與子帶調(diào)制梳狀譜干擾所采用的子帶脈沖是相同的,因此這兩種干擾信號(hào)對(duì)目標(biāo)圖像性能的影響基本一致,這也說(shuō)明其性能指標(biāo)的改善主要來(lái)自于子帶脈沖信號(hào)所獲得的相干處理增益,驗(yàn)證了理論分析的正確性。
本文研究了對(duì)寬帶SAR的子帶調(diào)制梳狀譜干擾技術(shù),詳細(xì)分析了子帶脈沖信號(hào)形式及其匹配濾波輸出特性,由于子帶脈沖信號(hào)保持了原始雷達(dá)脈沖信號(hào)的調(diào)頻斜率,與雷達(dá)信號(hào)具有一定的相干性,因此獲得的部分處理增益,從而節(jié)省了干擾功率。子帶脈沖信號(hào)的梳狀譜調(diào)制相當(dāng)于在頻域?qū)ψ訋}沖進(jìn)行復(fù)制和搬移,提高了干擾信號(hào)的密集程度,進(jìn)而改善了干擾效果。理論分析和仿真實(shí)驗(yàn)均表明,子帶脈沖調(diào)制梳狀譜干擾相比于噪聲干擾來(lái)說(shuō),降低了對(duì)干擾功率的需求,具有良好的干擾效果。
[1]Goj W W.Synthetic aperture radar and electronic warfare[M].[S.l.]:Artech House,1993.
[2]Dumper K,Cooper P S,Wons A F,et al.Spaceborne synthetic aperture radar and noise jamming[C]//IEEE International Radar Conference.Edinburgh:IEEE Press,1997:411-414.
[3]Huang Hongxu,Zhou Yiyu,Wu Jing,et al.A frequencybased inter/intra partly coherent jamming style to SAR[C]//Internation Conference on Signal Processing System.Dalian.China:IEEE Press,2010:434-437.
[4]Ye Wei,Ruan Hang,Zhang Shuxian,et al.Study of noise jamming based on convolution modulation to SAR[C]//Internation Conference on CMCE.Changchun.China:IEEE Press,2010:169-172.
[5]吳曉芳,王雪松,梁景修.SAR-GMTI高逼真勻速運(yùn)動(dòng)假目標(biāo)調(diào)制干擾[J].宇航學(xué)報(bào),2012,33(10):1472-1479.Wu Xiaofang,Wang Xuesong,Liang Jingxiu.Modulation jamming method for high-vivid false uniformly-moving targets against SAR-GMTI[J].Journal of Astronautics,2012,33(10):1472-1479.
[6]張立峰,吳建新.基于總體最小二乘的多通道SAR-GMTI方法[J].現(xiàn)代雷達(dá),2013,35(4):26-29.Zhang Lifeng,Wu Jianxin.A method of multi-channel SARGMTI based on total least squares[J].Modern Radar,2013,35(4):26-29.
[7]趙 博,楊軍,孫光才,等.一種虛假大場(chǎng)景SAR快速轉(zhuǎn)發(fā)欺騙干擾方法研究[J].電子與信息學(xué)報(bào),2012,34(4):963-968.Zhao Bo,Yang Jun,Sun Guangcai,et al.A method of SAR fast repeater deception jamming for large false scene[J].Journal of Electronics& Information Technology,2012,34(4):963-968.
[8]韓國(guó)強(qiáng),李永禎,邢世其,等.對(duì)新型SAR欺騙干擾效果的評(píng)估方法[J].宇航學(xué)報(bào),2011,32(9):1994-2002.Han Guoqiang,Li Yongzhen,Xing Shiqi,et al.Research on an evaluation method for new deceptive jamming effect on SAR[J].Journal of Astronautics,2011,32(9):1994-2002.
[9]楊 光,李紹濱,姜義成.輻射定標(biāo)SAR欺騙干擾的影響改進(jìn)方法[J].現(xiàn)代雷達(dá),2010,32(1):26-29.Yang guang,Li Shaobin,Jiang Yicheng.Impact of radiometric calibration to SAR deception jamming and new improved jamming method[J].Modern Radar,2010,32(1):26-29.
[10]吳彥鴻,王宏艷,賈 鑫.對(duì)SAR電子干擾的初步分析[J].現(xiàn)代雷達(dá),2007,29(9):101-104.Wu Yanhong,Wang Hongyan,Jia Xin.Elementary problem study of ECM against SAR[J].Modern Radar,2007,29(9):101-104.
[11]Schleher D C.Electronic warfare in the information age[M].Norwood,Massachusetts:Artech House,1999.
[12]陳立明,沈愛(ài)國(guó),姜秋喜.寬帶SAR的子帶脈沖循環(huán)轉(zhuǎn)發(fā)干擾[J].現(xiàn)代防御技術(shù),2010,38(4):100-104.Chen Liming,Shen Aiguo,Jiang Qiuxi.Jamming wideband LFM signal with sub-band pulse[J].Modern Defence Technology,2010,38(4):100-104.
[13]張澄波.綜合孔徑雷達(dá)[M].北京:科學(xué)出版社,1989.Zhang Chenbo.Synthetic aperture radar[M].Beijing:Publishing House of Science,1989.