第18卷 第6期 太赫兹科学与电子信息学报Vo1.18，No.6 2020年12月 Journal of Terahertz Science and Electronic Information Technology Dec.，2020 文章编号：2095-4980(2020)06-1003-07
王雪琦，涂刚毅，吴少鹏
(中国船舶重工集团公司第七二四研究所，江苏南京 211106)
摘 要：分布式相参雷达(DCAR)是目前国内外雷达领域的重要研究方向，精确的参数估计是实现其良好相参性能的前提和核心。基于动目标模型，提出一种基于多脉冲积累的相参参数估计
方法。该方法通过对多脉冲信号进行快、慢时间匹配滤波处理，实现多脉冲相参积累；再利用互 相关法进行相参参数估计。仿真分析对比了不同脉冲个数和不同输入信噪比下的参数估计性能和
相参性能，仿真结果表明，该方法具有可行性，且可以有效提高低信噪比情况下的参数估计性能
和相参性能。
关键词：分布式相参；参数估计；动目标；多脉冲积累
中图分类号：TN957.51文献标志码：A doi：10.11805/TKYDA2019182
Coherent parameters estimation method for distributed coherent radar
based on multi-pulse accumulation
WANG Xueqi，TU Gangyi，WU Shaopeng
(No.724 Research Institute of CSIC，Nanjing Jiangshu 211106，China)
Abstract：Distributed Coherent Aperture Radar(DCAR) is an important research direction in the field of radar at home and abroad. Accurate parameter estimation is the premise and core of good coherence
performance. Based on the moving target model, a coherent parameter estimation method based on
multi-pulse accumulation is proposed. The method performs fast-time and slow-time match filtering for
multi-pulse signals, and obtains the results of multi-pulse coherent accumulation. Then the
cross-correlation method is utilized to estimate the coherent parameters. The performance of parameter
estimation and correlation under different numbers of pulses and different input signal-to-noise ratios are
compared by simulation analysis. The simulation results show that the method is feasible and can
effectively improve the performance of parameter estimation and coherence in low signal-to-noise ratio.
Keywords：distributed coherence；parameter estimation；moving target；multi-pulse accumulation
分布式相参雷达(DCAR)因具有较好的探测性能、高角度分辨力、灵活性和机动性等一系列技术优势而成为目前国内外雷达领域研究热点[1-5]。DCAR系统由若干具有机动性的单元雷达和一个中心控制系统组成。这些雷达按照一定的阵列形式进行布局，通过对单元雷达的回波进行信号融合，达到大孔径雷达的性能[5]。 DCAR主要工作于3种模式[2]：独立工作模式、接收相参模式(MIMO模式)和收发相参模式(全相参模式)。独立工作模式下，各单元雷达独立进行目标的搜索和截获；截获目标后转入接收相参模式，各单 元雷达发射正交信号，波束指向同一区域，并根据回波信息来估计相参参数；参数估计达到一定精确度下转入收发相参模式，各单元雷达发射相同的信号波形，并以估计得到的相参参数在收发端进行时延和相位的调整，实现收发全相参。由上述工作过程可知，相参参数的精确估计是实现良好的相参性能的核心和前提。对于相参参数估计问题，目前已经取得了一些成果。文献[6-8]推导了不同结构下参数估计的克拉美罗界，分析了参数估计性能和相参性能；文献[9]提出了2种模式下的相参参数跟踪方法，分析了2种模式的误差影响；文献[10-11]分析了信号的非理想正交性对参数估计的影响，并提出一种基于信号重建的相参参数估计算法。目前，相参参数估计的难点和热点主要 收稿日期：2019-05-20；修回日期：2019-07-21
基金项目：装备发展部装备预先研究项目(41413010201)；船舶重工联合预先研究项目(61411304010101-06)
作者简介：王雪琦(1995-)，男，在读硕士研究生，主要研究方向为分布式相参雷达.email:**************
1004
太赫兹科学与电子信息学报 第18卷
是动目标情况下参数估计和低信噪比情况下的参数估计问题，如何在低信噪比环境下实现动目标的相参合成是DCAR 实现高性能探测的关键技术[12]。
本文首先建立动目标情况下信号模型，然后基于该模型提出一种利用脉冲多普勒处理实现多脉冲积累的相参参数估计方法。通过对多脉冲回波信号进行快、慢时间的匹配滤波，实现多脉冲信号相参积累，估计多普勒频率并利用互相关法进行相参参数估计。同时，本文还分析了多脉冲积累情况下的接收相参模式和收发相参模式的信噪比增益。仿真分析不同脉冲个数和输入信噪比情况下的参数估计精确度和相参性能，结果表明，该方法在低信噪比情况下可以提高相参参数估计精确度和相参合成性能。
1 信号模型
由K 个发射单元和L 个接收单元组成的一般结构的DCAR 系统见图1[6]。假设在接收相参模式下一组脉冲宽度为T p 的窄带正交脉冲信号为()(1,2,,)k s t k K = ，其相关性满足式(1)：
(),
()()
,k
k k A t k k s t s t k k *''
=⎧⊗
-=⎨'
≠⎩ (1) 式中()k A t 为()k s t 自相关结果。对()k s t 以本振信号t
c j2πj e k f t ϕ+进行上变频之后作为发射信号，可以表示为：
t
c j2πj ()()e k
f t k k x t s t ϕ+= (2)
式中：c f 为载频；t k ϕ为第k 个发射单元雷达上变频时引入的本振相位。
假设发射信号的脉冲重复间隔为r T ，相干处理间隔为r NT 。为便于分析，假设在相干处理间隔内目标处于同一个距离单元且目标发射和接收方向的径向速度保持不变。根据DCAR 基线布置准则[13-14]
，对于不同的发射或
接收单元，径向速度认为近似相等。在第(0,1,,1)n n N =- 个脉冲发射时刻，目标到第(1,2,,)k k K = 个发射单
元和第(1,2,,)l l L = 个接收单元的距离可以表示为：
t t t r r r r r ()(0)()(0)k k l l
R n R nT R
R n R nT R ⎧=-⨯⎪⎨=-⨯⎪⎩ (3) 式中：t (0)k R 和r (0)l R 分别为目标到第k 个发射单元和第l 个接收单元的初始距离；t R
和r R 分别为发射方向和接收方向的距离变化率。因此，第n 个周期脉冲从第k 个发射单元经目标散射到达第l 个接收单元的时间可以表示为：
r t r ()()()(0)lk l k lk R
n n n nT c
ττττ=+=-⨯ (4)
式中：c 为光速；t t ()()/k k n R n c τ=和r r
()()/l l n R n c τ=为第n 个周期发射和接收单程时延；(0)[(0)(0)]/lk k l R R c τ=+为
第0个周期脉冲从第k 个发射单元经目标到达第l 个接收单元的时间；t r R
R R =+ 为传播路径的距离变化率。 假设目标只包含单个散射点，且在相干处理间隔内目标对于不同路径电磁波的复散射系数相等，记为ξ。
每个接收单元收到的回波信号都是K 个发射信号的混合回波。在第l 个接收单元处将接收的第n 个周期的混合回波信号以本振信号r j2πj e l
c f t ϕ-进行下变频处理，r l ϕ为第l 个发射单元下变频时引入的本振相位。下变频处理后的信号
可以表示为：
[]t r
j2π()j j 1
(,)()e (,)c lk k l K
f n l k lk l k r t n s t n w t n τϕϕξτ-++==-+∑ (5)
式中：t 为快时间；n 为慢时间(周期脉冲)；(,)l w t n 为第l 个接收单元第n 个周期脉冲回波的噪声信号，是零均值的复高斯过程，且满足
2
(,)(,)()()()l l w E w t n w t n l l t t n n σδδδ*''''''⎡⎤=---⎣⎦ (6)
式中：2
w σ为噪声信号的平均功率；()δ⋅为狄拉克函数。
假设所有的接收信号都以第1个发射单元发射经目标散射到第1个接收单元的传播路径作为参考标准，则令
t r
1111
j2π(0)j j 0e c f τ
ϕϕξξ-++=,t t t 1()()()k k n n n τττ∆=-,r r r 1()()()l l n n n τττ∆=-,t t t 1k k ϕϕϕ∆=-,r r r 1l l ϕϕϕ∆=-,t t t ()2π()k c k k n f n θτϕ∆=-∆+∆,
r r r c ()=2π()l l l n f n θτϕ∆-∆+∆。同时，由于假设目标在相干处理间隔内位于同一个距离单元，对于所有0,1,,1n N =- ，
Fig.1 General architecture for DCAR 图1 一般结构的分布式相参雷达
target
T 1
R L
T 2 … …糖水不等式
T K
R 2
R 1
… …
第6期 王雪琦等：分布式相参雷达多脉冲积累相参参数估计方法 1005
令()lk lk n ττ=，则式(5)可以写为：
r
t
d r j ()j ()j2πr t
0111(,)()()e
(,)l k K
n n f nT l k l k l k r t n s t n n w t n θθξτττ∆+∆+=⎡⎤=--∆-∆+⎣⎦∑ (7) 式中d c /f f R
c = 为目标多普勒频率。由于假设目标在相干处理间隔内接收方向和发射方向的距离变化率保持不变，因此各单元在接收方向和发射方向的距离差保持不变，故令t t ()k k n ττ∆=∆，r r ()l l n ττ∆=∆。由于t ()k n θ∆和r ()l n θ∆与t ()k n τ∆和r ()l n τ∆成线性关系，故令t t ()k k n θθ∆=∆，r r ()l l n θθ∆=∆。t t r
r {,,,}k k l l τθτθ∆∆∆∆统称为相参参数，其中t k τ∆和
t k θ∆称为发射时延差和发射相位差，r l τ∆和r l θ∆称为接收时延差和接收相位差。显然，t r t r 11110ττθθ∆=∆=∆=∆=。根据上述相参参数的定义，单元雷达l 接收的第n 个周期回波信号可以表示为：
r t
d r j j j2πr t 0111
(,)()e (,)l k K
f nT l k l k l k r t n s t w t n θθξτττ∆+∆+==--∆-∆+∑ (8)
对相参参数进行估计并以估计参数对信号进行调整和相参融合就可以实现高的相参性能。
2 相参参数估计
2.1 快时间匹配滤波
对式(8)信号以匹配滤波器()()k k
h t s t *
=-进行快时间匹配滤波处理，匹配滤波输出结果(,)lk y t n 可以表示为： r t
d r j j j2πr t 011(,)(,)()()
e (,)l k
f nT lk l k k l k lk y t n r t n h t A t w t n θθξτττ∆+∆+=⊗=--∆-∆+ (9)
式中(,)(,)()lk l k w t n w t n h t =⊗为噪声的匹配滤波输出。
2.2 慢时间匹配滤波
对快时间匹配滤波之后的信号作慢时间周期采样，得到慢时间采样序列lk Y ，表示如下：
T
T T T (,0),(,1),,(,1)lk lk lk lk y t y t y t N ⎡⎤=⋅⋅⋅-⎣⎦Y (10)
对采样序列做慢时间匹配滤波(离散时间傅里叶变换)，得到时间—多普勒二维域的输出结果为：
r t
d r r
r 1
1
j j j π(1)()j2πj2πr t d r 0110
d r sin[π()](,)(,)e
()e e (,)e sin[π()]l k N N N f f T fnT fnT lk lk k l k
lk n n f f NT Y t f y t n A t w t n f f T θθξτττ--∆+∆-----==-==--∆-∆+-∑∑ (11)
式中f 为慢时间匹配滤波得到的频域变量。由于离散时间傅里叶变换(DTFT)得到的频率变量是连续的，不能直接
计算得出。因此，以离散傅里叶变换(DFT)作为替代。由于DFT 是DTFT 结果在频域的采样，为了尽量减小多普勒频率和相参参数估计误差，应该尽可能增加DFT 点数。
对所有快时间匹配滤波分离的L ×K 路信号都作相同的处理。由式(11)可知，(,)lk Y t f 信号峰值对应的频率为多普勒频率。取每一路最大值所对应的频率作为多普勒频率估计值，为了减小误差，对所有路径求得的多普勒频率估计值取平均：
[]d 11
1ˆarg max (,)L K
lk f l k f Y t f LK ===∑∑ (12) 式中arg max()f
⋅为“⋅”取得最大值所对应频率。
假设多普勒频率估计不存在误差，则相参积累结果可表示为：
r t d r 1
j j j2πr t d 0110
ˆ(,)()e (,)e l k N f nT lk k l k
lk n Y t f N A t w t n θθξτττ-∆+∆-==--∆-∆+∑ (13) 由式(13)可知，N 脉冲信号相参积累，信号
功率增大为单脉冲情况下的N 2倍，而噪声功率根据其分布特性只
能增大为N 倍，因此信噪比提高为原来的N 倍，减小了噪声的影响。
2.3 互相关法参数估计
1d ˆ(,)l Y t f 和1d
ˆ(,)k Y t f 相参积累的结果为： r j r 1d 0111ˆ(,)()e l l l Y t f NA t θξττ∆=--∆ (14) t j t 1d 011ˆ(,)()e k k k k
Y t f NA t θξττ∆=--∆ (15) 对相参积累结果运用互相关法进行参数估计。
1006
太赫兹科学与电子信息学报 第18卷
对d ˆ(,)lk Y t f 和1d
ˆ(,)l Y t f 做互相关运算： t j 2t d 1d 1ˆˆ()(,)(,)||e ()()()k lk lk l k k
X t Y t f Y t f t A t A t θξδτ∆**=⊗-=-∆⊗⊗- (16) 取()lk X t 的峰值时刻作为发射时延差的估计值，峰值相位作为发射相位差的估计值。为了减小误差，降低干扰和噪声的影响，对所有L 组发射参数估计值取平均，得到发射参数估计值：
{}t 1
t 11ˆarg max [()]1ˆphase max[()]L
k lk t
l L
k
lk l X t L X t L τθ==⎧∆=⎪⎪⎨⎪∆=⎪⎩
∑∑ (17) 式中：arg max()t
⋅为“⋅”取得最大值所对应时间变量；max()⋅为“⋅”的最大值；phase()⋅为对“⋅”取相位。
对()lk Y t 和1()k Y t 做互相关运算：
r j 2r d 1d
ˆˆ()(,)(,)||e ()()()l lk lk k l k k Z t Y t f Y t f t A t A t θξδτ∆**=⊗-=-∆⊗⊗- (18) 取()lk Z t 的峰值时刻作为接收时延差的估计值，峰值相位作为接收相位差的估计值。为了减小误差，降低干扰和噪声的影响，对所有K 组接收参数估计值取平均，得到接收参数估计值：
r 1
r 11ˆarg max [()]1ˆphase{max[()]}K
l lk t
k K
l
lk k Z t K Z t K τθ==⎧∆=⎪⎪⎨⎪∆=
⎪⎩
∑∑ (19) 基于多脉冲积累互相关法的DCAR 参数估计流程见图2。
Fig.2 Procedure for estimating coherent parameters of moving targets
图2 动目标相参参数估计流程
3 相参性能分析
在接收相参模式下，利用相参参数估计值t ˆk τ∆,t
ˆk θ∆,r ˆl
τ∆和r
ˆl θ∆对式(13)的信号进行时延和相位调整，并进行累加，得到的累加信号：
r t r t
r t r
ˆˆj j r t d
111
ˆˆj j j j j2πr
t r t 0
11
11
110
ˆˆˆ()(,)e ˆˆ()e
(,)e
e
l k l k l k
L K
lk l k l k L
K
L K N fnT k
l k
lk l
k
l k l k n Y t Y t f N A t w t n θθδθδθθθττξτ
δτδτττ-∆-∆==-+-∆-∆-======+∆+∆=---++∆+∆∑∑∑∑∑∑∑ (20)
式中：t t t ˆk k k δτττ=∆-∆；t t t ˆk k k δθθθ=∆-∆；r r r ˆl l l δτττ=∆-∆；r r r ˆl l l δθθθ=∆-∆。理想情况下，t r t r 0k l k l
δτδτδθδθ====，相参合成信号()Y t 可以表示为：
r t
r ˆˆj j j2πr t 01111
111
ˆˆ()()(,)e e l k L K L K N
fnT k lk l k l k l k n Y t N A t w t n θθξτττ-∆-∆-======-++∆+∆∑∑∑∑∑ (21)
在收发相参模式下，各发射单元发射相同的信号波形()s t ，以估计得到的发射参数t ˆk τ∆,t
ˆk
θ∆对各发射单元的发射信号进行时延和相位调整，得到t
t
ˆj2πj j t ˆ()()e e c k k f t k
k x t s t ϕθτ+-∆'=+∆。在单元雷达l 处，接收信号下变频并进行接收参数的调整，得到的信号可以表示为：
典当管理办法
r t
r j j j2πr t 011r 1
(,)()e (,)l k d K
f nT l l k l k r t n s t nT w t n δθδθξτδτδτ++==----+∑ (22)
对式(22)进行快、慢时间匹配滤波处理之后，相参累加结果可表示为：
r t r 1
j j j2πr t d 0111
ˆ(,)()e
(,)e l k K
N fnT l l k lk k n Y t f N A t w t n δθδθξτδτδτ-+-===---+∑∑ (23)
第6期 王雪琦等：分布式相参雷达多脉冲积累相参参数估计方法 1007
式中()A t 表示()s t 的自相关输出。对所有l 个接收单元的回波信号作相同处理，并进行叠加，得
r t r t r ˆˆj j j j j2πr t r t d 0111
11
11
ˆˆˆ()(,)()e
(,)e e l k l k L
L
K
L
N
fnT l l k l l k l l k l n Y t Y t f N A t w t n δθδθθθξτδτδτττ+-∆-∆-=======---++∆+∆∑∑∑∑∑ (24) 理想情况下，t
r t r 0k l k l
δτδτδθδθ====，累加信号()Y t 可表示为： r t
r ˆˆj j j2πr t 01111
ˆˆ()()(,)e e l k L N
fnT l l k l n Y t NKLA t w t n θθξτττ-∆-∆-===-++∆+∆∑∑ (25)内肋管
假设任意的()k A t 和()A t 峰值功率都相等，记为A P 。因此2种模式下的信号功率相等，都是22()||A NLK P ξ。且
由于时移和相移不会改变高斯过程的分布特性，因此接收相参和收发相参模式下的噪声功率分别为2
w
NLK σ和2
w NL σ，故2种模式下的输出信噪比可表示为：
2MIMO 2
22
FC SN 2SN ||||A
w A w NKL P R NK L P R ξσξσ⎧=⎪⎪
⎨⎪=⎪⎩
(26) 单部雷达工作时的输出信噪比可以表示为：
2single SN
2||A
董振堂w
P R
ξσ=
(27)
因此，信噪比增益为：
MIMO
越野MIMO SN
SN,gain single SN FC
FC 2SN SN,gain single SN =R R NKL
R R R NK L
R ⎧==⎪⎪
⎨⎪=⎪⎩
(28) 由式(28)可以看出，多脉冲信号相参积累在接收相参模式下可以实现NKL 倍的信噪比增益，在收发相参模式
下可以实现2NK L 倍的信噪比增益。
4 仿真分析
以2个接收单元和2个发射单元组成的DCAR 系统进行仿真，仿真参数设置见表1。
表1 仿真参数
Table1 List of simulation parameters
华北水利水电学院学报simulation parameter value simulation parameter parameter value
carrier frequency/GHz 3 distance of transmitting unit 1/km 30
sampling frequency/MHz
200 distance of receiving unit 1/km 30 pulse width/μs
2.56 distance of transmitting unit 2/km 30.015 6 pulse repetition frequency/kHz
10 distance of receiving unit 2/km
30.020 4 scattering coefficient
1 local oscillator phase difference of transmitting unit/(°) 30 radial velocity in the transmitting direction/(m/s) 80 local oscillator phase difference of receiving unit/(°)
45 radial velocity in the receiving direction/(m/s)
120
DFT points
256
仿真采用正交相位编码信号[15]，脉冲个数N =8,16,32，输入信噪比区间[-26,16] dB ，步长为2 dB 。分别进行仿真，每组参数进行100次蒙特卡洛试验，以均方根误差(RMSE)作为参数估计性能的指标，定义如下
RMSE (29) 式中：c M 表示蒙特卡洛试验次数；()x i 和ˆ()x
i 分别表示第i 次试验的真实值和估计值。仿真得到不同输入信噪比情况下的时延差、相位差和多普勒频率估计的均方根误差(RMSE)，见图3。从图3可以看出，相参参数和多普
勒频率的估计精确度随输入信噪比的提高而提高，且在相同输入信噪比条件下，积累脉冲数越多，参数估计精确度越高，表明通过多脉冲积累可以有效提高在低信噪比情况下的参数估计性能。受采样频率的限制，时延差估计误差趋于定值。此外，接收相位差的估计精确度随输入信噪比的提高而提高，
而发射相位趋于定值，这是由于发
射信号非理想正交性影响发射参数的估计[10-11]。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议