孟 喆,劉光炎
(南京電子技術(shù)研究所, 南京210039)
由兩顆或兩顆以上衛(wèi)星構(gòu)成的主星收發(fā)信號,其他輔星被動接收回波信號的多基站星載SAR系統(tǒng)是將衛(wèi)星編隊和SAR技術(shù)有機結(jié)合的新體制天基遙感系統(tǒng),是全天時、全天候,高效獲取全球無縫高精度地面數(shù)字高程模型(DEM)的優(yōu)選手段。收發(fā)分置單航過多基站SAR系統(tǒng)對地成像已在航天飛機雷達地形測繪使命(Shuttle Radar Topography Mission,SRTM)上得到驗證,該系統(tǒng)得到2 000個30 m×30 m,絕對高度分辨率15 m的數(shù)字高程模型,完成對地表近80%信息獲取;基于兩個TerraSAR-X雷達衛(wèi)星組成的Tan-DEM-X多基站星載系統(tǒng)首次實現(xiàn)了單航過干涉成像[1-2]。
單個星載SAR系統(tǒng)的頻率源對信號進行上下變頻、調(diào)制、解調(diào)。頻率源不存在頻率偏差問題,并且頻率源噪聲中的低頻分量被消除。與單個星載SAR不同,相位同步是多基站SAR系統(tǒng)的特殊問題,非相干的頻率源被置于不同的衛(wèi)星上,頻率源噪聲的低頻分量無法消除,導致引入的相位誤差不可忽視。為了保證多基站SAR系統(tǒng)的正常運行,實現(xiàn)主輔星載SAR系統(tǒng)的相位同步,需要在主輔星SAR系統(tǒng)之間建立相應的相位同步鏈路。
工程上容易實現(xiàn)的衛(wèi)星同步鏈路有直達波方案和雙向?qū)鞣桨福?-4]。與前者相比,后者可以通過差分抵消空間鏈路及部分雷達設(shè)備和同步天線等引起的相位誤差。本文基于主星發(fā)射,主輔星接收的多基站SAR系統(tǒng)模型,通過對主星與輔星雷達載波相位的雙向?qū)鳎o以相位補償處理,在主星與輔星之間實現(xiàn)相位同步。
相位同步原則是在兩星或多星間建立的參考相位與主輔星載波相位相近,補償由不同衛(wèi)星上雷達的頻率源之間相位噪聲和頻率差引起的多基干涉相位誤差,下面的分析以兩星為例展開。在主星與輔星之間建立一個相位同步鏈路,可以將包含雷達載波頻率偏差和相位噪聲的信號在主星與輔星的雷達設(shè)備之間進行交換。
如圖1所示,一次雙向?qū)鞯倪^程為主星在t時刻發(fā)射脈寬為Tp的同步信號,t+τ12時刻由輔星接收。在t+τsys時刻輔星發(fā)射同步脈沖,經(jīng)過τ21后主星接收。圖2為星上雙向?qū)魍芥溌肥疽鈭D,以主星發(fā)射信號為例,以晶振為基準的頻率源產(chǎn)生的本振信號與信號產(chǎn)生的基帶信號經(jīng)過調(diào)制放大后變?yōu)橐詅syn為周期的發(fā)射信號,經(jīng)過空間傳輸?shù)竭_輔星并由其接收;間隔一定時間后,輔星返回一個同步信號,主星按類似過程接收。忽略大氣對信號傳輸?shù)挠绊?,?jīng)過主輔星解調(diào)的同步脈沖數(shù)據(jù)被分別下傳至地面,利用短時間τsys內(nèi)正、反向傳輸?shù)目臻g鏈路的路徑基本相同,對消去空間鏈路引入的相位,就得到一個周期為τsyn、重復頻率為fsyn的主、被動雷達載波相位差的采樣序列,通過對該補償信號進行N(N=PRF/fsyn)倍插值處理變?yōu)榕cSAR信號PRF相同的相位差,對輔星進行相位補償,即實現(xiàn)多基站雷達間的相位同步,示意圖如圖3所示。
圖1 雙向?qū)魍芥溌肥疽鈭D
圖2 星上雙向?qū)魍芥溌肥疽鈭D
圖3 地面處理流程
由于雷達的頻率源是以晶振作為基準源,晶振的穩(wěn)定度是評定同步鏈路性能好壞的重要因素之一,在此假設(shè)超穩(wěn)晶振模型采用單邊帶相位密度譜的冪律模型形式[7-9]
為了分析方便,記φji代表衛(wèi)星i發(fā)射、衛(wèi)星j接收到的解調(diào)相位。
假設(shè)主輔星雷達載波頻率為
式中:f0為標稱頻率;Δfi為第i個頻率源的頻率偏差。則t時刻衛(wèi)星i的雷達載波相位為
式中:t0為數(shù)據(jù)獲取的起始時刻。為了文中討論方便,令t0=0,其中φ0i為載波初相,nφi(t)為載波相位噪聲。
依據(jù)圖1b)和圖2所示,主輔星在一次對傳中,收/發(fā)相應時刻相位見表1。
表1 一次對傳主輔星發(fā)射/接收時刻相位
由表1可知,t+τ12時刻輔星的解調(diào)相位為
假定同步信號雙向?qū)鬟^程中,忽略雙向信號傳輸時延差異,即 τ12=τ21=τ。
對于雙向傳輸?shù)耐芥溌范?,相當于是對每一個衛(wèi)星解調(diào)的相位以fsys進行采樣。
采樣時刻 tk=k/fsyn,其中 k=1,2,…,kmax,kmax=[Tdata·fsys],Tdata為數(shù)據(jù)獲取時間。則補償?shù)南辔粸?/p>
從式(6)中也可以看出,采用解調(diào)相位差作為輔星的補償相位的優(yōu)點是只要在一次對傳過程中涉及到的天線、鏈路路徑、所有收發(fā)硬件等共有部分將會相互抵消[5]。
同步鏈路的建立就是為了獲取相位噪聲與頻率偏差,但是從式(6)可以看出,經(jīng)過補償后的SAR信號還存在剩余的相位誤差。除了接收機噪聲誤差之外,由于用有限的fsyn對相位噪聲譜Sφ(f)采樣,還會引入插值誤差和混疊誤差[6]。
插值誤差是由于在|f|<fsyn/2以外的頻率分量由于采樣無法重構(gòu)造成的誤差,混疊誤差來自于采樣頻率是fsyn整數(shù)倍的周期譜在|f|<fsyn/2帶內(nèi)疊加形成。接收機噪聲引入的誤差是由于在接收機噪聲帶寬內(nèi)的噪聲譜以fsyn采樣,頻譜折疊到|f|<fsyn/2內(nèi)引起的。
考慮到影響同步鏈路的一些因素,可以將這一過程用數(shù)學模型表示如圖4所示。
圖4 同步鏈路剩余誤差模型
模型中的各傳遞函數(shù)含義如下:
(1)HLP(f):脈沖交替同步傳遞函數(shù),相當于對每一個頻率源在不同時刻采樣,考慮到τ≤τsys,忽略τ的影響。
(2)Hsyn(f)專門用于濾出補償相位,消除濾波器的失配誤差。
(3)Haz(f)為方位脈壓傳遞函數(shù),采用方位壓縮時可以降低剩余誤差[4]。
分析同步鏈路的剩余誤差,需要知道每一個誤差源的功率譜密度表達式[6],見式(10)~式(12)。
(1)接收機噪聲功率譜密度為
(2)插值誤差的功率譜密度函數(shù)為
(3)混疊誤差的功率譜密度函數(shù)為
根據(jù)已有的研究和信號與系統(tǒng)原理可知,未經(jīng)方位壓縮的各個誤差信號的方差,表達式中不含方位脈壓傳遞函數(shù)[4,7-8],詳見式(13)~ 式(15)。
(1)接收機噪聲方差為
因此總剩余誤差表達式為
依據(jù)上面的理論分析,下面通過仿真來模擬分析同步鏈路誤差對相位同步的影響。
假設(shè)仿真模擬的試驗參數(shù)見表2。
表2 仿真主要參數(shù)
圖5為方位壓縮前后,不同信噪比(SNR)情況下接收機噪聲誤差隨同步頻率的變化情況,由圖可見,當SNR固定時,方位壓縮前接收機噪聲誤差隨著同步頻率的增加不變;方位壓縮以后,隨著同步頻率的增加,接收機噪聲誤差變小。圖6為方位壓縮前后,不同頻段下,插值誤差隨同步頻率的變化情況,同理可見,插值誤差隨同步頻率的增加而減小。圖7為方位壓縮前后,不同頻段下,混疊誤差隨同步頻率的變化情況,由圖可見,混疊誤差也隨同步頻率的增加而減小。由圖還可以看到,方位壓縮有利于改善接收機的噪聲誤差、插值誤差和混疊誤差。
圖5 接收機噪聲誤差
圖6 插值誤差
圖7 混疊誤差
圖8~圖10為不同SNR下,方位壓縮前后不同頻段的同步鏈路總誤差隨同步頻率的變化情況。仿真結(jié)果表明,無論是否經(jīng)過方位壓縮,信噪比的大小對于總剩余誤差的貢獻非常明顯,不同頻段,隨著同步頻率的增加,信噪比越大,總剩余誤差抑制性越好。這是因為,接收機的噪聲誤差在總誤差中所占比重隨著同步頻率的增加而變大。同時,由仿真分析可以看到,在信噪比較大的情況下,不同頻段,同步鏈路剩余誤差的抑制不同,相對而言,L頻段對同步鏈路的剩余誤差抑制性最好。
從仿真分析可以得到:在相同同步頻率下,不同頻段的系統(tǒng),方位壓縮后可以較好地抑制同步鏈路的總剩余誤差。
圖8 SNR=10 dB時3個頻段總誤差
圖9 SNR=30 dB時3個頻段總誤差
圖10 SNR=50 dB時3個頻段總誤差
綜合上面的分析可知,相位同步鏈路系統(tǒng)中,接收機的噪聲誤差在總誤差中所占比重隨著同步頻率的增加而變大,主要原因是插值誤差和混疊誤差隨著同步頻率的增加受到的抑制更加明顯,而接收機噪聲受到的抑制相對較小。因此,提高接收通道的信噪比可以有效抑制同步鏈路的剩余誤差。
對多基站星載SAR系統(tǒng)而言,雙向?qū)鞣桨笇崿F(xiàn)相位同步是一個有益的探索,方案除了考慮基準頻率的高一致性外,還需要考慮同步鏈路的可行性、鏈路系統(tǒng)的高信噪比和相位誤差補償?shù)囊籽a償性,以期實現(xiàn)相位同步的高性能??梢姡p向同步對傳方案是實現(xiàn)多星SAR系統(tǒng)相位同步的一個選擇。
[1]Fritz T,Breit H,Balss U et al.Interferonetric processing of tan DEM-X data[C]//IEEE International Symposium on Geo Science and Remote Sensing.Vancouner BC,Canada:IEEE Press,2011.
[2]Martinez N Y,Eineder M,Brcic R,et al.TanDEM-X mission:SAR image coregistration aspects[C]//8th European Conference on Synthetic Aperture Radar.Eurogress Aachen,Germany:[s.n.],2010:576-579.
[3]Wang W Q.Multi-antenna synthetic aperture radar[M].[S.l.]:CRC Press,2013.
[4]Younis M,Metzig R,Krieger G.Performance prediction of a phase synchronization link for bistatic SAR[J].IEEE Geoscience and Remote Sensing Letters,2006,3(3):429-433.
[5]Younis M,Metzig R.Bachmann G K M,et al.Performance prediction and verification for the synchronization link of Tan-DEM-X[C]//IEEE Geoscience and Remote Sensing Symposium.Barcelone,Spain:IEEE Press,2007:5206-5209.
[6]Krieger G,Moreira A,F(xiàn)iedler H,et al.TanDEM-X:a satellite formation for high-resolution SAR interferometry[J].IEEE Transactions on Geoscience and Remote Sensing,2007,45(11):3317-3341.
[7]Krieger G,Younis M.Impact of oscillator noise in bistatic and multistatic SAR[J].IEEE Geoscience and Remote Sensing Letters,2006,3(3):424-428.
[8]張升康,楊汝良.振蕩器相位噪聲對雙站SAR成像影響分析[J]. 測試技術(shù)學報,2008,1(22):7-12.Zhang Shengkang,Yang Ruliang.Analysis of oscillator phase noise effects on bistatic SAR imaging[J].Journal of Test and Measurement Technology,2008,22(1):7-12.