修正的星載PD雷達雜波仿真模型

王之立，于 澤

(北京航空航天大學 電子信息工程學院，北京 100191)

0 引言

雜波仿真在雷達仿真中占有重要地位，尤其對于一些目標檢測的算法驗證、雷達性能指標的驗證都需要有較為精確的雜波仿真。早在20世紀60年代，國外就出現了對于機載脈沖多普勒(PD)雷達雜波仿真的相關研究[1]。隨著技術和國際形勢的發展，機載雷達探測范圍小、留空時間短等弊端逐漸暴露。星載PD雷達成為世界各國研究的新方向，文獻[2-4]闡述了各國的研究現狀。

文獻[5-7]介紹了PD雷達雜波仿真的相關方法，由于PD雷達仿真需要模擬雜波在距離—多普勒域對目標的遮蓋，往往需要考慮距離向和方位向的混疊，因此仿真覆蓋區域很大，數據量往往很大。如文獻[8-9]中所述，為了盡可能地降低運算量，通常雜波仿真模型在機載條件下做了如下簡化：忽略載機與雜波單元在同一個相干處理時間中的空間相對運動，只保留相對運動導致的多普勒相位。雜波劃分方法通常有矩形網格映像法和等距離—多普勒網格映像法，將地面劃分成網格，每個網格中的雜波等效為一個雜波單元。機載雷達平臺運動速度慢，在一個相干處理時間內載機與雜波單元的相對運動往往不會超過一個最小分辨單元，因此簡化的機載PD雜波模型廣泛適用。而對于星載條件，由于衛星的飛行速度遠大于飛機，一個相干處理時間內衛星與雜波單元的相對運動往往超過幾個甚至幾十個最小分辨單元，雜波會在距離時域及多普勒頻域出現二維展寬。此時，簡化的機載PD雜波模型已經無法準確地對展寬進行模擬。

目前已有的星載PD雷達仿真方法往往基于簡化的機載PD雜波模型，考慮衛星高度導致的地球曲率，但對雜波在距離時域—多普勒頻域產生的二維展寬問題沒有進行分析[10]，文獻[11]分析了星載條件下地球曲率和地球自轉對雜波數學模型的影響，但沒有對仿真模型進行分析。本文在未經簡化的原始雜波模型的基礎上，對現有的雜波模型進行了修正，提出了一種修正的星載PD雷達雜波仿真模型。該模型相對未經簡化的原始雜波模型在運算量上有大幅的進步，使得其能適用于星載PD雷達仿真的大數據量；相對于簡化的機載PD雜波模型，該模型修正了衛星高速運動時導致的一些誤差，使得其能適用于星載條件。

1 現有雜波模型

1.1 發射信號模型

Chirp信號是當前使用頻繁的發射信號形式[12]，

每個發射脈沖為：

p(t)=a(t)cos[ωt-φ(t)]，

式中，a(t)為矩形窗函數，滿足

τp為脈沖寬度；Tprt為脈沖重復周期。考慮到通常采用線性調頻信號，有

φ(t)=-πbt2，

b為調頻率。

1.2 原始雜波回波模型

對于某個斜距為R(t)的地面雜波單元，其回波信號經由接收機解調后可以表示為：

式中，A(t)為隨時間變化的幅度系數，滿足[13]

P為發射功率；η為天線效率；λ為載波波長；σ為電磁波反射物的RCS；Gr和Gs分別為雷達的收發天線方向圖。

考慮到Gs(t)，Gr(t)，R(t)隨時間t變化緩慢，由于信號收發時間間隔短，基于停走模型，s(t)可做如下近似：

如文獻[12]所述，利用變量置換，將以時間t為變量的一元函數s(t)轉換成以快時間τ和慢時間η為變量的二元函數。令

t=η+τ，η=kTprt，

則有

該模型基于物理實際過程，結果相對準確，但對于PD雷達仿真而言，由于雜波單元數目太大，其運算量往往難以實際運用。因此，實際機載PD雜波仿真時往往對此模型進行了簡化，以降低運算量。

1.3 簡化的機載PD雜波模型

式中，η1為某個相干處理時間的開始方位向慢時刻；ηNa為某個相干處理時間的終止方位向慢時刻；A1為η1時刻目標回波對應的幅度；vd1為η1時刻平臺與目標的徑向相向速度；R1為η1時刻平臺與目標的回波距離。

機載模型中的近似條件是相干處理時間內飛機與雜波單元的相對運動不超出最小網格，可以忽略。當模型運用到星載條件下時，由于衛星速度可以達到飛機速度的幾十倍，這一假設不再成立。相干處理時間內雜波單元與衛星的相對運動往往能跨越幾個甚至幾十個最小網格，從而導致雜波回波在距離—多普勒域的展寬。因而，本節所述的簡化雜波模型不再適用于星載條件下的雜波仿真。

2 修正的星載PD雜波模型

2.1 基本原理

常見的地面雜波劃分有2種方法：矩形網格映像法[15]和等距離—多普勒網格映像法[16]。本文提出的修正模型可運用于2種網格映像法中，為了方便說明，本文采用更為直觀的矩形網格映像法，且網格在地表劃分方向分別為沿飛行方向和垂直于飛行方向在地面的投影，每個網格單元的大小需要能夠保證相鄰網格單元的間隔不超過雷達距離和多普勒分辨率。地表連續分布的雜波散射體可以離散地等效成一個個的雜波單元點。在相干處理時間內，衛星與地面發生相對運動，矩形網格與衛星保持相對靜止也與地面發生相對運動。因此，雜波單元與衛星的相對運動可等效為地面的雜波單元在地面網格中發生相對運動。地表網格的俯視圖如圖1所示，其中實心點為某個雜波單元，其在地面網格中的等效相對運動可用虛箭頭線表示。

圖1 等效雜波單元點在網格中相對運動示意

與傳統的網格映像法不同，本文中的網格不僅僅用于初始時刻劃分各個雜波單元，還用于標識各個雜波單元在相干處理時間內與衛星的相對位置。文中修正模型的思路便是將雜波單元在網格中的連續運動離散成多個網格間的運動，各網格內的仿真模型仍然沿用機載PD雜波模型，從而可以解決雜波在距離—多普勒域的展寬問題。

2.2 回波模型

在相干處理時間內，對于某個雜波單元，假設其跨越了k個網格，則將相干處理時間劃分成對應的k段，使得雜波單元在每段時間內都不跨越網格。在每段時間內，按照機載仿真模型，假設相對靜止只考慮相對運動帶來的多普勒頻率；在各段之間考慮相對運動，更新衛星與各個雜波單元的相對位置。

假設在某段相干處理時間內，方位向慢時間為η1≤η≤ηNa。期間目標相對于平臺在k個網格間運動，對應網格的斜距分別為R1，R2，R3，...，Rk，對應網格的平臺徑向速度分別為vd1，vd2，vd3，...，vdk。

在方位向慢時間η1≤η≤ηNa內，雜波單元對應的回波為：

式中，si(η，τ)為雜波單元在第i段時間ηi_start≤η≤ηi_end內對應的回波，在ηi_start≤η≤ηi_end時間內雜波單元在第i個網格內，si(η，τ)滿足

式中，Ai為第i個網格處目標回波對應的幅度；vdi為第i個網格處平臺與目標之間的徑向相向速度；Ri為第i個網格處平臺與目標間的斜距；ηi為雜波單元進入Ri對應網格時的慢時間；Ei為第i個網格的相位補償，用于補償進入第i個網格時的初始相位，滿足

式中，Rηi_start為第i段時間初始慢時刻ηi_start對應的目標斜距。

如上式所示，當雜波單元在各個網格間移動時，除了移動雜波單元到所對應的網格，為了后續方法與機載雷達網格映像法一致，還要對該雜波單元的回波添加Ei·exp{-j·4π/λ·vdi·(ηi-η1)}的相位。

3 仿真及結果分析

脈沖重復頻率fPRF=5 kHz，載波波長λ=0.03 m，衛星對地速度V=7 000 m/s，衛星軌道高度H=550 km，相干處理時間內方位向采樣點數Na=200，采樣率Fs=15 MHz，信號的調頻率b=1×1012Hz/s，脈沖寬度τ=2.1 μs。

分別用簡化的機載PD雜波模型、原始雜波模型和修正的星載PD雜波模型對下視角α=45°，斜視角β=45°的一處雜波單元進行仿真，其結果如圖2所示。

圖2 雜波回波距離時域-方位時域圖

從圖2中可以看出，雜波單元在相干處理時間內與衛星發生相對移動，導致其回波偏移了多個距離門。原始雜波模型在相干處理時間內描繪了斜距的連續變化，機載PD雜波模型則沒有表現出斜距的變化，修正的星載PD雜波模型離散地將平臺與雜波單元的斜距變化表現了出來。

PD雷達的雜波仿真一個主要的目的是模擬雜波在距離—多普勒域對目標回波的遮蓋。在星載條件下，由于雜波單元與衛星的相對運動，導致了雜波回波在距離—多普勒域的二維展寬，如圖3所示。機載PD雜波模型無法準確地在星載條件下進行仿真，而修正的星載PD雜波模型則解決了這個問題，精確度更高。

圖3 雜波回波距離時域-方位頻域圖

4 結束語

在星載條件下，對現有機載PD雷達雜波模型的簡化條件進行了分析，通過仿真對其存在的誤差做了分析論證。針對誤差的成因，細化了簡化條件，提出了一個修正的星載PD雷達雜波模型。修正模型細化了相干處理時間內雷達與雜波單元的相對位置變化，相較于傳統模型更接近理論情況。仿真結果表明，修正模型的誤差在距離向和方位向均低于分辨率，因此不會對后期處理造成影響，保證了雜波仿真的精確度。