一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计

苏海; 张晓娟; 叶盛波; 张群英; 方广有

doi:10.12000/JR16076

一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计

DOI: 10.12000/JR16076

苏海^①,②,,
张晓娟^①, ,,
叶盛波^①,,
张群英^①,,
方广有^①,

①.
中国科学院电磁辐射与探测技术重点实验室北京 100190
②.
中国科学院大学北京 100190

基金项目: 国家自然科学基金(61172017)，中国科学院国防科技创新基金

详细信息

作者简介:
苏海：苏海(1991–)，男，广西北流人，现为中国科学院电子学研究所在读硕士研究生，主要研究方向为伪随机编码体制超宽带雷达技术。Email: suhai_214@163.com

张晓娟(1964–)，女，现为中国科学院电子学研究所研究员，博士生导师，主要研究方向为微波成像新方法及应用，微波遥感、天线、计算电磁学、电磁散射与逆散射。E-mail: xjzhang@mail.ie.ac.cn

叶盛波(1983–)，男，湖北武汉人，现为中国科学院电子学研究所副研究员，研究方向为超宽带雷达系统、探地雷达研制、雷达信号处理方法研究。E-mail: sbye@mail.ie.ac.cn

张群英(1970–)，女，陕西西安人，现为中国科学院电子学研究所研究员，博士生导师，研究方向为微波探测新方法研究、超宽带雷达系统、雷达信号处理。E-mail: qyzhang@mail.ie.ac.cn

方广有(1963–)，男，现为中国科学院电子学研究所研究员，博士生导师，中国科学院电子学研究所副所长，研究方向为超宽带雷达系统、测月雷达、太赫兹成像、地球物理探测方法、电磁场理论。E-mail: gyfang@mail.ie.ac.cn

通讯作者:
张晓娟 xjzhang@mail.ie.ac.cn

中图分类号: TN958
计量
- 文章访问数: 3064
- HTML全文浏览量: 621
- PDF下载量: 836
- 被引次数: 10
出版历程
- 收稿日期: 2016-05-30
- 修回日期: 2017-02-12
- 网络出版日期: 2017-02-28

Design of an Ultra-wideband Pseudo Random Coded MIMO Radar Based on Radio Frequency Switches

Su Hai^{①,②
,},
Zhang Xiaojuan^{①
, ,},
Ye Shengbo^{①
,},
Zhang Qunying^{①
,},
Fang Guangyou^{①
,}

①.
The Key Laboratory of Electromagnetic Radiation and Sensing Technology,Chinese Academy of Sciences,Beijing 100190,China
②.
University of Chinese Academy of Sciences,Beijing 100190,China

Funds: The National Natural Science Foundation of China (61172017), The National Defense Science and Technology Innovation Fundation of Chinese Academy of Sciences

摘要

摘要: 多发多收超宽带雷达能实时探测到目标的距离向和方位向信息，在地质勘测、生命救援、穿墙目标跟踪等领域具有重要的应用价值。该文提出一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计方法，研制出一套伪随机编码超宽带MIMO雷达系统。为了降低系统的成本和复杂性，采用射频开关切换的方法来实现多发多收。由于射频开关的耐压值有限，采用峰值功率为18 dBm的伪随机编码信号作为雷达发射信号。最后通过混合采样方法实现超宽带信号的采样接收，并且调用FPGA内部的DSP硬核来进行实时脉冲压缩，以缓解后续数据处理运算量大的问题。实验结果表明，该雷达系统能实时探测到目标的方位向和距离向信息。
- 超宽带雷达 /
- 伪随机编码 /
- 开关切换 /
- 多发多收 /
- FPGA
Abstract: A Multiple-Input Multiple-Output (MIMO) ultra-wideband radar can detect the range and azimuth information of targets in real time. It is widely used for geological surveys, life rescue, through-wall tracking, and other military or civil fields. This paper presents the design of an ultra-wideband pseudo random coded MIMO radar that is based on Radio Frequency (RF) switches and implements a MIMO radar system. RF switches are employed to reduce cost and complexity of the system. As the switch pressure value is limited, the peak power of the transmitting signal is 18 dBm. The ultra-wideband radar echo is obtained by hybrid sampling, and pulse compression is computed by Digital Signal Processors (DSPs) embedded in an Field-Programmable Gate Array (FPGA) to simplify the signal process. The experiment illustrates that the radar system can detect the range and azimuth information of targets in real time.
- Ultra-wideband radar /
- Pseudo random code /
- Radio Frequency (RF) switch /
- Multiple-Input Multiple-Output (MIMO) /
- Field-Programmable Gate Array (FPGA)

HTML全文

1. 引言

非合作无源双基地雷达(Non-cooperative Passive Bistatic Radar, NPBR)因具有造价低、抗截获、反隐身等诸多优势^[1–6]而成为国内外研究热点。目前，NPBR研究的外辐射源信号主要包括通信电台信号、电视广播信号、手机信号以及雷达信号等。无论针对哪种外辐射源信号，在NPBR中，3大同步(包括空间同步、时间同步以及相位同步)问题始终是制约目标有效探测的关键问题。其中，空间同步是指接收和发射天线同时照射相同空域，空间同步关系到接收到的回波是否持续含有目标信息，即接收回波的有效性；时间同步是指接收天线必须准确获知发射天线发射波形时刻，时间同步关系到获取目标运动参数的准确性；相位同步是指发射和接收天线接收到的信号能够在长时间内严格相参，相位同步关系到目标相参积累性能，进而影响目标的探测性能。

本文基于雷达外辐射源信号模型，重点讨论与时间同步及相位同步相关的几个问题，具体包括，发射波形参数估计、直达波到达时间估计以及基于波形参数估计的长时间相参积累等问题。其中，发射波形参数估计主要包括脉冲宽度、脉冲重复间隔、载频、带宽等脉间捷变的参数估计。基于估计的波形参数值，进一步给出了捷变波形的GRFT (Generalized Radon Fourier Transform)长时间相参积累弱目标探测算法，最后通过数值实验验证了本文算法的有效性。

本文内容安排如下：第1节简单介绍了NPBR研究现状及存在的关键问题，明确了本文研究范畴；第2节建立了雷达外辐射源信号模型，给出了基于直达波的NPBR参数估计方法，并提出了一种鲁棒性高的脉冲提取方法；第3节给出了基于GRFT的捷变波形长时间积累算法；第4节利用实测数据对本文算法进行了验证；最后对本文内容进行了总结。

2. NPBR直达波参数估计

图1给出了NPBR工作示意图。NPBR接收的信号主要包括两个部分：(1)直接接收到的辐射源旁瓣信号，称为直达波信号；(2)目标前向散射的回波信号。NPBR合理布设的情况下，直达波信号信噪比将远大于目标回波信噪比，十分有利于发射波形参数的估计。因此，本文将基于直达波信号对辐射源发射波形及波达时间等参数进行估计。

图 1 NPBR示意图

Figure 1. Sketch map of NPBR

下载: 全尺寸图片幻灯片

2.1 信号模型

设外辐射源发射如式(1)所示的LFM脉冲串信号，脉冲串信号脉冲重复间隔、脉宽、载频及带宽均为脉间捷变。

$\begin{aligned}&\! s(\tau ,n) \!=\!\! A{\rm{rect}}\!\!\left( {\frac{\tau }{{{T_{\rm{p}}}\left( n \right)}}} \!\right)\!\exp \!\left\{ {{\rm j}{{π}}\! \left( {2{f_{\rm{c}}}\left( n \right)\tau \!+\! \gamma \left( n \right){\tau ^2}} \right)} \!\right\}\!,\\& \quad\quad\quad \ \ n = 0,1, ·\!·\!·,N - 1\end{aligned}$

(1)

式中， $\tau$ 为快时间，A为幅度，N为脉冲数， ${T_{\rm{p}}}(n)$ 为脉冲宽度， ${f_{\rm{c}}}\left( n \right)$ 为载频， $\gamma (n) = {B_{\rm{s}}}(n)/{T_{\rm{p}}}(n)$ 为调频率， ${B_{\rm{s}}}\left( n \right)$ 为信号带宽，n为脉冲序号。则直达波信号模型可写为：

$\begin{aligned}{s_{\rm{r}}}(\tau ,n) = & {A_1}{\rm{rect}}\left( {\frac{{\tau - {\tau _0}}}{{{T_{\rm{p}}}\left( n \right)}}} \right)\exp \left\{ { - {\rm j}2{{π}} {f_{\rm{c}}}\left( n \right){\tau _0}} \right\}\\& \cdot \!\exp \left\{ {{\rm j}{{π}} \left[ {2{f_{\rm{I}}}\left( n \right)\tau \!+\! \gamma \left( n \right){{\left( {\tau \!-\! {\tau _0}} \right)}^2}} \right]} \right\}\end{aligned}$

(2)

式中， ${\tau _0}{\rm{ = }}2{d_0}/{\rm c}$ 为直达波波达时间，d₀为外辐射源与接收天线之间的距离， ${f_{\rm{I}}}\left( n \right){\rm{ = }}{f_{\rm{c}}}\left( n \right) - {f_{{\rm{down}}}}$ 为中频频率， ${f_{{\rm{down}}}}$ 为下变频参考频率。

显然，式(2)中 ${\tau _0}$ , ${T_{\rm{p}}}\left( n \right)$ , $\gamma (n)$ 以及 ${f_{\rm{c}}}\left( n \right)$ 为待估计的未知参数。

2.2 直达波参数估计

根据待估计参数，可将直达波参数估计过程分为如图2所示的两个部分：脉冲提取和脉冲参数估计。

图 2 直达波参数估计流程图

Figure 2. Flow chart of parameter estimation of the direct waveform

下载: 全尺寸图片幻灯片

2.2.1 脉冲提取脉冲提取通常分为3个步骤，包括时域直达波初提取、自适应带通滤波以及时域直达波脉冲精提取。直达波信号中，通常可能包含大量同频段电台、通信等干扰成分，因此原始直达波信号在时域进行脉冲提取误差较大，需要在初提取的基础上进行带通滤波。而由于直达波信号中心频率、带宽未知，因此无法直接设计带通滤波器。考虑到LFM信号频率响应近似为矩形，本文采用通过脉冲提取方法提取信号的频率响应，从而获得信号中心频率与带宽，进而设计带通滤波器滤除直达波干扰信号。因此，脉冲提取方法的性能决定了直达波脉冲提取的效果。

脉冲串波形信息主要包含脉冲的上升沿、下降沿以及脉内调制信息。直达波脉冲提取可以等效于直达波脉冲的上升沿以及下降沿的提取。为此，本小节给出一种高鲁棒性的脉冲提取方法，该方法基本主要步骤包括：

步骤1 自适应计算噪声阈值门限，对超过噪声阈值门限的信号样本索引号进行差分处理；

步骤2 提取差分值大于最小脉冲间距的索引号作为预选上升沿；

步骤3 将信号序列翻转，用相同的差分处理方法以及信号序列翻转前后的对应关系获得预选下降沿。上升沿与下降沿一一对应，组成预选脉冲；

步骤4 从预选脉冲中剔除不符合预设脉冲条件的脉冲，将剩余脉冲作为最终提取结果。

通常，噪声采样幅值无法连续超过噪声门限，因此无法形成与直达波脉冲宽度和脉冲重复间隔等特征相近的脉冲。本质上，该方法正是利用了噪声与直达波脉冲的这种特征区别，因此具有较好的鲁棒性。

2.2.2 脉冲参数估计图2中，脉冲参数估计部分包括以下几个步骤：调频率和中心频率的估计、直达波脉冲对齐、直达波波达时间及载频估计。

在2.2.1小节脉冲提取的基础上，可以获得直达波信号的脉冲宽度及脉冲重复间隔。本小节针对提取的任意脉冲，给出脉内调制信息的提取方法。

(a) 调频率和中心频率的估计

每个提取的脉冲信号均为一个LFM信号，对于LFM信号的调频率及中心频率的估计方法较多。例如LvD^[7,8], GRFT^[9–11]等。考虑到GRFT参数估计的最优性^[12,13]，这里采用GRFT对LFM信号进行参数估计。设定中心频率范围为 ${\tilde f_k} \in \left[ {{f_{\min }},{f_{\max }}} \right]$ 及调频率范围为 ${\tilde \gamma _l} = \left[ {{\gamma _{\min }},{\gamma _{\max }}} \right]$ ，则中心频率及调频率估计值可由式(3)得到。式(3)中，L和K分别为中心频率和调频率搜索点数， ${s_{{\rm{pulse}}}}(\tau ,n)$ 为提取的第n个脉冲信号。

$\begin{aligned}& {\left[ {{{\tilde f}_{\rm{I}}}\left( n \right),\tilde \gamma \left( n \right)} \right]}\\& \ \ = \mathop {\rm argmax}\limits_{{{\tilde f}_k},{{\tilde \gamma }_l}} \left\{ {\sum\limits_{l = 1}^L {\sum\limits_{k = 1}^K {{s_{{\rm{pulse}}}}(\tau ,n)} } } \right.\\& \ \ \quad \cdot {\exp \left\{ { - {\rm j}{π} \left[ {2{{\tilde f}_k}\tau + {{\tilde \gamma }_l}{\tau ^2}} \right]} \right\}} \Biggr\}\\& \ \ = \mathop {\rm argmax}\limits_{{{\tilde f}_k},{{\tilde \gamma }_l}} \left\{ {\sum\limits_{l = 1}^L {\sum\limits_{k = 1}^K {\tilde A\exp } } \left\{{\rm j}{π} \left[ {2\left( {{f_{\rm{I}}}\left( n \right) - {{\tilde f}_k}} \right)\tau } \right. \right.} \right.\\& \ \ \quad{ { { + \left( {\gamma \left( n \right) - {{\tilde \gamma }_l}} \right){\tau ^2}} \Bigr]} \Big\}} \Biggr\}\end{aligned}$

(3)

式(3)中， $\tilde A$ 为与 $\tau$ 无关的复幅度。在实际中，由于提取脉冲的上升沿和下降沿可能存在较大误差，导致中心载频的估计值 ${\tilde f_{\rm{I}}}\left( n \right)$ 与真实估计值存在一个不可忽略的偏差，记为 $\Delta f\left( n \right)$ ，即 ${\tilde f_{\rm{I}}}\left( n \right) = {f_{\rm{I}}}\left( n \right) -$ $\Delta f\left( n \right)$ 。中心频率偏差将在后续脉压中产生峰值位置的偏移。

(b) 直达波脉冲对齐

利用获取的脉冲宽度、中心频率以及调频率，可构建如式(4)所示的参考信号，并对式(2)所示的直达波脉冲信号进行脉冲压缩。脉压结果如式(5)所示。

$\begin{aligned}\!\!\!\!\!\!\!\!\!\!\!\!\! {s_{{\rm{ref}}}}(\tau ,n) =& {\rm{rect}}\left( {\frac{\tau }{{{{\tilde T}_{\rm{p}}}\left( n \right)}}} \right)\\ & \cdot\exp \left\{ {{\rm j}{{π}} \left[ {2{{\tilde f}_{\rm{I}}}\left( n \right)\tau + \tilde \gamma \left( n \right){\tau ^2}} \right]} \right\}\end{aligned}$

(4)

$\begin{aligned}&{s_{{\rm{PC}}}}(\tau ,n)\\ &\ \ = {s_{\rm{r}}}(\tau ,n) \otimes s_{{\rm{ref}}}^*(\tau ,n)\\&\ \ = {A_1}\exp \left\{ { - {\rm j}2{{π}} {\tau _0}{f_{\rm{c}}}\left( n \right)} \right\} \\ & \ \ \quad\cdot\int_{ - \infty }^\infty {{\rm{rect}}\left( {\frac{{u - {\tau _0}}}{{{T_{\rm{p}}}\left( n \right)}}} \right){\rm{rect}}\left( {\frac{{\tau - u}}{{{{\tilde T}_{\rm{p}}}\left( n \right)}}} \right)} \\& \ \ \quad\cdot \exp \left\{ {{\rm j}{{π}} \left[ {2{f_{\rm{I}}}\left( n \right)u + \gamma \left( n \right){{\left( {u - {\tau _0}} \right)}^2}} \right]} \right\}\\&\ \ \quad \cdot \exp \left\{ { - {\rm j}{{π}} \left[ {2{{\tilde f}_{\rm{I}}}\left( n \right)\left( {\tau - u} \right) + \tilde \gamma \left( n \right){{\left( {\tau - u} \right)}^2}} \right]} \right\}\!\!{\rm d}u\\& \ \ \approx {\rm{rect}}\left( {\frac{{\tau + {T_{\rm{p}}}\left( n \right)/2 - {\tau _0}}}{{{{\tilde T}_{\rm{p}}}\left( n \right)}}} \right) \\ & \ \ \quad \cdot \exp \left\{ {{\rm j}{{π}} \left[ {\left( {\tau + {\tau _0}} \right)\Delta f\left( n \right) - 2{\tau _0}{f_{\rm{c}}}\left( n \right)} \right]} \right\}\\& \ \ \quad\cdot{\rm{sinc}}\left\{ {\left( {\tau \!-\! {\tau _0} \!+\! {{\tilde T}_{\rm{p}}}\left( n \right)} \right){{π}} \tilde \gamma \left( n \right)\left[ {\tau \!-\! {\tau _0} \!+\! \Delta {\tau _0}\left( n \right)} \right]} \right\}\end{aligned}$

(5)

式(5)中，

$\Delta {\tau _0}\left( n \right) = \Delta f\left( n \right)/\tilde \gamma \left( n \right)$

(6)

显然，由于中心频率估计误差 $\Delta f\left( n \right)$ 的存在，脉压后，峰值位置存在一个与脉冲号相关偏移量 $\Delta {\tau _0}\left( n \right)$ ，即不同脉冲的峰值位置处于不同距离单元。利用式(6)中不同脉冲峰值的位置偏移与频率估计误差之间的关系，可实现脉冲间相对频率估计误差的补偿，并实现脉冲对齐。记各脉冲与第1个脉冲的峰值偏差为：

$\Delta {\tau _{\rm{r}}}\left( n \right){\rm{ = }}\Delta {\tau _0}\left( n \right) - \Delta {\tau _0}\left( 0 \right)$

(7)

由式(7)，定义相对频率偏差为：

$\Delta {f_{\rm{r}}}\left( n \right) = \Delta {\tau _{\rm{r}}}\left( n \right)\tilde \gamma \left( n \right)$

(8)

重新构建脉压匹配函数为：

$\begin{aligned}{\hat s_{{\rm{ref}}}}(\tau ,n) =& {\rm{rect}}\left( {\frac{\tau }{{{{\tilde T}_{\rm{p}}}\left( n \right)}}} \right)\\&\! \cdot\exp\! \left\{ {{\rm j}\!{{π}}\! \left\{ {2\!\left[ {{{\tilde f}_{\rm{I}}}\left( n \right) \!-\! \Delta {f_{\rm{r}}}\left( n \right)} \right]\!\!\tau \!+\! \tilde \gamma \left( n \right){\tau ^2}} \right\}} \right\}\end{aligned}$

(9)

则脉压结果更新为：

$\begin{aligned}&s'\!\!_{{\rm{PC}}}(\tau ,n)\\& \ \ = {s_{\rm{r}}}(\tau ,n) \otimes \hat s_{{\rm{ref}}}^*(\tau ,n)\\ & \ \ \approx{A_0}{\rm{rect}}\left( {\frac{{\tau + {T_{\rm{p}}}\left( n \right)/2 - {\tau _0}}}{{{{\tilde T}_{\rm{p}}}\left( n \right)}}} \right) \\ & \ \ \quad \cdot\exp \left\{ {{\rm j}{{π}} \left[ {\left( {\tau + {\tau _0}} \right)\Delta f\left( 0 \right) - 2{\tau _0}{f_{\rm{c}}}\left( n \right)} \right]} \right\}\\ & \ \ \quad \cdot {\mathop{\rm sinc}\nolimits} \left\{ {\left( {\tau \!-\! {\tau _0} + {{\tilde T}_{\rm{p}}}\left( n \right)} \right){{π}} \tilde \gamma \left( n \right)\left[ {\tau \!-\! {\tau _0} \!+\! \Delta {\tau _0}\left( 0 \right)} \right]}\right\}\end{aligned}$

(10)

由式(10)可知，脉冲峰值位置均为：

$\tau = {\tau _{{\rm{peak}}}} = {\tau _0} - \Delta {\tau _0}\left( 0 \right)$

(11)

显然，与脉冲号无关，实现了脉冲对齐。

由式(7)和式(8)可知

$\begin{aligned}\Delta {f_{\rm{r}}}\left( n \right) &= \left( {\Delta {\tau _0}\left( n \right) - \Delta {\tau _0}\left( 0 \right)} \right)\tilde \gamma \left( n \right)\\ &{\rm{ = }}\Delta f\left( n \right) - \Delta {\tau _0}\left( 0 \right)\tilde \gamma \left( n \right)\\ & = \Delta f\left( n \right) - \Delta f\left( 0 \right)\tilde \gamma \left( n \right)/\tilde \gamma \left( 0 \right)\end{aligned}$

(12)

将式(12)代入式(10)中相位项，得到

$\begin{aligned}&\hat \phi \left( {n,{\tau _{{\rm{peak}}}}} \right)\\ &\ \ = \exp \left\{ {{\rm j}{{π}} \left[ {\left( {{\tau _{{\rm{peak}}}} + {\tau _0}} \right)\Delta f\left( 0 \right) - 2{\tau _0}{f_{\rm{c}}}\left( n \right)} \right]} \right\}\\ &\ \ = \exp \Big\{ {{\rm j}{{π}}\Big\{ {\left( {2{\tau _{{\rm{peak}}}} + \Delta f\left( 0 \right)/\tilde \gamma \left( 0 \right)} \right)\Delta f\left( 0 \right)} } \\ &\ \ \quad - 2\left( {{\tau _{{\rm{peak}}}} + \Delta f\left( 0 \right)/\tilde \gamma \left( 0 \right)} \right) \\ & \ \ \quad \cdot \left[ {{f_{{\rm{down}}}} \!+\! {{\tilde f}_{\rm{I}}}\left( n \right) \!\!-\!\! \Delta {f_{\rm{r}}}\left( n \right) \!\!+\!\! \Delta f\left( 0 \right)\tilde \gamma \left( n \right)/\tilde \gamma \left( 0 \right)} \right] \Big\}\Big\}\end{aligned}$

(13)

因此，利用式(14)可以得到第1个脉冲的中心频率误差的估计值。

$\begin{aligned}\Delta\! \tilde f\left( 0 \right)=& \mathop {{\rm{argmax}}}\limits_{\Delta {f_k}} \sum\limits_{n = 0}^{N - 1} {{s_{{\rm{MF}}}}\left( {{\tau _{{\rm{peak}}}},n} \right)} \\ & \cdot\exp \Big\{ { - {\rm j}{{π}} \Big\{ {\left( {2{\tau _{{\rm{peak}}}} + \Delta {f_k}/\tilde \gamma \left( 0 \right)} \right)} } \\&\cdot\Delta {f_k}+ 2\left( {{\tau _{{\rm{peak}}}} + \Delta {f_k}/\tilde \gamma \left( 0 \right)} \right)\\ &\cdot\left[ {{f_{{\rm{down}}}} + {{\tilde f}_{\rm{I}}}\left( n \right) - \Delta {f_{\rm{r}}}\left( n \right)} \right.\\& {{ + \Delta {f_k}\tilde \gamma \left( n \right)/\tilde \gamma \left( 0 \right)} \Big]} \Big\}\Big\}\end{aligned}\!\!$

(14)

式中，搜索频率 $\Delta {f_k} \in \left[ {\Delta {f_k}_{\min },\Delta {f_k}_{\max }} \right]$ , $\Delta {f_k}_{\min }$ 和 $\Delta {f_k}_{\max }$ 为搜索最小值和最大值。

进而得到直达波到达时间为：

${\tau _{\rm{0}}}{\rm{ = }}{\tau _{{\rm{peak}}}}{\rm{ + }}\Delta {\tau _0}\left( 0 \right) = {\tau _{{\rm{peak}}}}{\rm{ + }}\Delta \tilde f\left( 0 \right)/\tilde \gamma \left( 0 \right)$

(15)

中心频率估计值更新为：

$\begin{aligned}\tilde f_{\rm{I}}'\left( n \right){\rm{ = }}& {{\tilde f}_{\rm{I}}}\left( n \right) - \Delta \tilde f\left( n \right)\\ = & {{\tilde f}_{\rm{I}}}\left( n \right) \!-\! \left( {\Delta \tilde f\left( 0 \right)/\tilde \gamma \left( 0 \right) \!+\! \Delta {\tau _{\rm{r}}}\left( n \right)} \right)\tilde \gamma \left( n \right)\end{aligned}$

(16)

进一步可以得到各脉冲的载频估计值为：

${\tilde f_{\rm{c}}}\left( n \right) = {f_{{\rm{down}}}}{\rm{ + }}\tilde f_{\rm{I}}'\left( n \right)$

(17)

至此，外辐射源发射信号的脉冲重复频率、脉冲宽度、调频率、载频以及直达波达到时间均已获得，为脉间相参积累奠定基础。

3. 长时间积累

通过直达波的处理，可获得外辐射源信号的参数，并将获得的参数估计值代入式(4)可得到重构的脉压参考函数。利用该参考函数，对特定距离门内的回波进行脉压，得到脉压后距离-脉冲维结果。在此基础上，沿脉冲维进行捷变波形的长时间积累技术。目标回波时延为 ${\tau _{\rm{d}}}\left( n \right){\rm{ = 2}}r(n)/{\rm c}$ ，式中，c为光速， $r(n)$ 为目标瞬时斜距，且

$r(n){\rm{ = }}{r_0} + {v_0}\sum\limits_{n' = 0}^n {{T_{\rm{r}}}\left( {n'} \right)} + \frac{{{a_0}}}{2}{\left( {\sum\limits_{n' = 0}^n {{T_{\rm{r}}}\left( {n'} \right)} } \right)^2}$

(18)

式中，r₀为初始斜距，v₀为径向速度，a₀为径向加速度。

结合式(1)，目标回波可以写为：

$\begin{aligned}{s_{{\rm{echo}}}}(\tau ,n) = &{A_2}{\rm{rect}}\left( {\frac{{\tau - {\tau _0}}}{{{T_{\rm{p}}}\left( n \right)}}} \right)\\ &\cdot\exp \left\{ { - {\rm j}2{{π}} {f_{\rm{c}}}\left( n \right){\tau _{\rm{d}}}\left( n \right)} \right\}\\& \cdot\! \exp \!\left\{ {{\rm j}\!{{π}}\! \left[ {2{f_{\rm{I}}}\left( n \right)\tau \!+\! \gamma \left( n \right){{\left( {\tau \!-\! {\tau _{\rm{d}}}\left( n \right)} \right)}^2}} \right]} \right\}\end{aligned}$

(19)

式中，A₂为回波幅度。

用 $\tilde f_{\rm{I}}'\left( n \right)$ 替换式(4)所示参考函数的 ${\tilde f_{\rm{I}}}\left( n \right)$ ，利用更新后的参考函数对式(19)进行脉压，得到

$\begin{aligned}&{s_{{\rm{PC\_echo}}}}(\tau ,n)\\& \ \ = {s_{{\rm{echo}}}}(\tau ,n) \otimes s_{{\rm{ref}}}^*(\tau ,n)\\& \ \ \approx {A_2}{\rm{rect}}\left( {\frac{{\tau - {\tau _{\rm{d}}}\left( n \right)}}{{2{{\tilde T}_{\rm{p}}}\left( n \right)}}} \right)\exp \left\{ { - {\rm j}2{{π}} {\tau _{\rm{d}}}\left( n \right){f_{\rm{c}}}\left( n \right)} \right\}\\&\ \ \quad \cdot {\rm{sinc}}\left\{ {{{π}} \left( {{{\tilde T}_{\rm{p}}}\left( n \right) \!-\! \left| {\tau \!-\! {\tau _{\rm{d}}}\left( n \right)} \right|} \right)\tilde \gamma \left( n \right)\left( {\tau \!-\! {\tau _{\rm{d}}}\left( n \right)} \right)} \right\}\end{aligned}$

(20)

基于式(20)，可得到运动参数空间中，相参积累的结果为：

$\begin{aligned}G\left( {i,j,k} \right) = & \sum\limits_{n = 0}^{N - 1} {{s_{{\rm{PC\_echo}}}}\left( {\tau \left( {{r_i},{v_j},{a_k};n} \right),n} \right)}\\ & \cdot \exp \left\{ {{\rm j}2{{π}} \tau \left( {{r_i},{v_j},{a_k};n} \right){{\tilde f}_{\rm{c}}}\left( n \right)} \right\}\end{aligned}$

(21)

式中，r_i, v_j, a_k分别为搜索的距离、速度和加速度值，且

$\begin{aligned}&\!\!\!\!\tau \left( {{r_i},{v_j},{a_k};n} \right)\\ &\ = \frac{{\rm{2}}}{\rm c}\left[ {{r_i} + {v_i}\sum\limits_{n' = 0}^n {{T_{\rm{r}}}\left( {n'} \right)} + \frac{{{a_i}}}{2}{{\left( {\sum\limits_{n' = 0}^n {{T_{\rm{r}}}\left( {n'} \right)} } \right)}^2}} \right]\end{aligned}$

(22)

式(21)实际上为文献[9–12]中GRFT的一种特殊形式。值得注意的是，式(20)中的脉压输出波形在脉冲间可以是捷变的，因此，其盲速旁瓣^[10]等积累性能与非捷变信号GRFT结果将会有较大区别，篇幅原因，本文不进行讨论。

当 $\tau \left( {{r_i},{v_j},{a_k};n} \right){\rm{ = }}{\tau _{\rm{d}}}\left( n \right)$ 时，式(18)得到峰值：

$G{\left( {i,j,k} \right)_{\max }} = {A_2}N$

(23)

即实现了积累增益随脉冲数增加而线性增加的相参积累结果。将式(21)中得到的积累幅值与恒虚警(Constant False Alarm Rate, CFAR)门限进行比较，即可得到最终检测结果。

对于功率为 ${\sigma ^2}$ 的高斯白噪声背景，N个脉冲相参积累后的噪声输出功率为 $N{\sigma ^2}$ ^[9]，因此，相参积累的输出信噪比为 $NA_2^2/{\sigma ^2}$ ，即为积累前信噪比的N倍。另外，本文方法对输入信噪比的要求与积累脉冲数及检测门限相关。在检测门限一定的条件下，目标回波信噪比越小，所需的积累脉冲数越多，进而所需的波束驻留时间也越长。

4. 数值实验

本文总体算法处理流程如图3所示。

图 3 算法流程

Figure 3. Flowchart of the algorithm

下载: 全尺寸图片幻灯片

本小节将利用仿真和实测数据实验对所提方法进行验证。

仿真实验中，外辐射源雷达工作在P波段，发射波形脉宽、时宽及带宽脉间捷变，场景中存在一个微弱目标，脉压后信噪比为0，运动参数为(120 km, 350 m/s, 20 m/s²)，积累脉冲数为100。部分原始回波实部信号和模值如图4所示。图5为直达波提取结果，该结果验证了本文脉冲提取方法的有效性。图6为目标回波脉压结果。显然，脉压结果中，目标淹没在噪声中。利用动目标检测(Moving Target Detection, MTD)方法，得到结果如图7所示，利用本文长时间积累方法，得到结果如图8所示。显然，通过时间和相位同步，利用本文的波形捷变GRFT方法能够将目标能量有效积累，进而验证了本文算法的有效性。

图 4 原始直达波信号

Figure 4. Original signal of the direct waveform

下载: 全尺寸图片幻灯片

图 5 脉冲提取结果

Figure 5. Pulse extraction result

下载: 全尺寸图片幻灯片

图 6 目标回波脉冲压缩结果

Figure 6. Echo pulse compression result

下载: 全尺寸图片幻灯片

图 7 MTD结果

Figure 7. MTD result

下载: 全尺寸图片幻灯片

图 8 时间相位同步后相参积累结果

Figure 8. Coherent integration result after time and phase synchronization

下载: 全尺寸图片幻灯片

实测数据实验中，外辐射源雷达工作在P波段，发射波形脉宽、时宽及带宽脉间捷变。图9为回波脉压结果，从脉压结果可以看出，场景中存在多个微弱目标，而且在观测时间内，目标存在明显跨距离单元现象。针对图9中弱目标1进行捷变波形GRFT长时间积累，得到目标检测及跟踪结果如图10所示。处理中，每一帧的积累脉冲数为20，平均积累时间约为0.42 s，帧间滑动步长为10个脉冲。图11给出了基于本文算法的相参积累与常规MTD、单脉冲处理的信噪比结果对比。由于目标存在跨距离单元现象，在相同积累脉冲数情况下，本文相参积累性能显著高于常规MTD。

图 9 目标回波脉冲压缩

Figure 9. Echo pulse compression result

下载: 全尺寸图片幻灯片

图 10 弱目标轨迹探测结果

Figure 10. Weak target trace detection result

下载: 全尺寸图片幻灯片

图 11 积累SNR比较

Figure 11. SNR comparisons

下载: 全尺寸图片幻灯片

实验结果表明，利用本文的基于直达波的外辐射源波形估计方法以及捷变波形GRFT相参积累方法，能够将长时间观测的目标回波能量有效积累。

5. 总结

本文针对非合作外辐射源双基地雷达中时间和相位同步问题以及捷变波形的长时间相参积累问题进行了研究。给出了一种直达波脉冲提取方法，并利用直达波参数估计实现了时间和相位同步，在此基础上，进一步给出了波形捷变GRFT长时间相参积累方法，通过仿真和实测数据实验，验证了该方法对微弱目标探测的有效性。

直达波信噪比大小会影响波形参数估计误差，该信噪比越小，波形参数估计误差越大，进而长时间积累性能也越差。估计误差对目标参数测量及相参积累性能的分析将在后续工作中进行定量分析。另外，后续将对双基地目标运动参数的解算、捷变波形GRFT的相参积累性能及其盲速旁瓣特性等方面内容开展研究。

图 1 雷达工作过程示意图

Figure 1. Schematic diagram of radar operating

下载: 全尺寸图片幻灯片

图 2 高速伪随机编码信号产生原理图

Figure 2. Schematic diagram of high-speed pseudo random coded signal

下载: 全尺寸图片幻灯片

图 3 混合采样时序图

Figure 3. Timing diagram of hybrid sampling

下载: 全尺寸图片幻灯片

图 4 脉冲压缩并行运算框架

Figure 4. Structure of pulse compression computed parallelly

下载: 全尺寸图片幻灯片

图 5 多发多收雷达系统原理图

Figure 5. Schematic diagram of MIMO radar system

下载: 全尺寸图片幻灯片

图 6 多发多收雷达系统样机

Figure 6. MIMO radar system

下载: 全尺寸图片幻灯片

图 7 雷达发射信号

Figure 7. Radar transmitting signal

下载: 全尺寸图片幻灯片

图 8 脉冲压缩结果

Figure 8. Results of pulse compression

下载: 全尺寸图片幻灯片

图 9 收发开关级联结构图

Figure 9. Structure of radio frequency switches

下载: 全尺寸图片幻灯片

图 10 各个通道一致性测试结果

Figure 10. Results of consistency of all channels

下载: 全尺寸图片幻灯片

图 11 不同衰减值的闭环测试结果

Figure 11. Results of different attenuation

下载: 全尺寸图片幻灯片

图 12 不同体制雷达的穿墙测试结果

Figure 12. Results of different radar systems though the wall

下载: 全尺寸图片幻灯片

图 13 2发4收对空气中两个目标测试场景

Figure 13. Experimental scene using 4-input 2-output radar to detect two targets in the air

下载: 全尺寸图片幻灯片

图 14 两个目标在不同间距时的测试结果图

Figure 14. Results of different distance of two targets

下载: 全尺寸图片幻灯片

图 15 2发4收穿墙测试场景

Figure 15. Experimental scene using 4-input 2-output radar to detect human through the wall

下载: 全尺寸图片幻灯片

图 16 墙体后目标的测试结果图

Figure 16. Results of human behind the wall

下载: 全尺寸图片幻灯片

表 1 多发多收雷达系统主要参数

Table 1. Main parameter of MIMO radar

参数	数值
通道数量	m发n收(最大4发16收)
码元形式	m序列
码元长度	1023
单个码元宽度(ns)	0.9375
发射信号幅值(Vpp)	5
信号中心频率(GHz)	1.06
–10 dB信号带宽(GHz)	0.2～2.1
实时采样率(MSPS)	200
混合采样率(GSPS)	12.8
平均次数(times)	32
接收机灵敏度(dBm)	–97
脉冲重复频率(KHz)	333.3
采样时窗(ns)	300
扫描率(track/s)	164
成像速率(frame/s)	$\frac{{164}}{{m \cdot n}}$
总功耗(W)	20

下载: 导出CSV

表 2 射频开关详细参数

Table 2. Main parameter of MIMO radar

参数	数值
带宽	0～3 GHz
插损	1.5 dB@0.5～2 GHz
1 dB压缩点	25 dBm@0.5～2 GHz
隔离度	37 dB@0.5～2 GHz
切换速度	25 ns/Turn on; 45 ns/Turn off

下载: 导出CSV

参考文献(18)

[1]	Immoreev I I and Fedotov P G S D V. Ultra wideband radar systems: Advantages and disadvantages[C]. IEEE Conference on Ultra Wideband Systems and Technologies, Baltimore, USA, 2002: 201–205. http://ieeexplore.ieee.org/xpls/icp.jsp?arnumber=1006348
[2]	Parrini F, Fratini M, Pieraccini M,et al.. ULTRA: Wideband ground penetrating radar[C]. 2006 European Radar Conference, UK, 2006: 182–185. http://ieeexplore.ieee.org/document/4058288/
[3]	吴世有, 谭恺, 徐艳云, 等. 超宽带穿墙雷达天线阵列配置分析及运动人体跟踪成像算法[J]. 电子与信息学报, 2012, 34(11): 2601–2607. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX201211009.htm Wu Shi-you, Tan Kai, Xu Yang-yun,et al.. Study on UWB through-wall radar antenna array configuration and moving person tracking and imaging algorithm[J].Journal of Electronics &Information Technology, 2012, 34(11): 2601–2607. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX201211009.htm
[4]	Haddad W S. The Rubble Rescue Radar (RRR): A low power hand-held microwave device for the detection of trapped human personnel[R]. Lawrence Livermore National Laboratory, CA (United States), 1997. http://digital.library.unt.edu/ark:/67531/metadc693400/m2/1/high_res_d/571102.pdf
[5]	Sachs J, Peyerl P, and Zetik R. Stimulation of UWB-sensors: Pulse or maximum sequence[C]. International Workshop on UWB Systems, Oulu, Finland, 2003: 5. http://ieeexplore.ieee.org/document/5117727/authors
[6]	Sachs J, Peyerl P, and Rossberg M. A new UWB-principle for sensor-array application[C]. Instrumentation and Measurement Technology Conference, Venice, Italy, 1999, 3: 1390–1395. doi: 10.1002/9783527651818.ch3/summary
[7]	夏正欢, 张群英, 叶盛波, 等. 一种便携式伪随机编码超宽带人体感知雷达设计[J]. 雷达学报, 2015, 4(5): 527–537. http://radars.ie.ac.cn/CN/abstract/abstract276.shtml Xia Zheng-huan, Zhang Qun-ying, Ye Sheng-bo,et al.. Design of a handheld pseudo random coded UWB radar for human sensing[J].Journal of Radars, 2015, 4(5): 527–537. http://radars.ie.ac.cn/CN/abstract/abstract276.shtml
[8]	张群英, 方广有. 伪随机序列编码脉冲信号在探地雷达中的应用研究[J]. 电子与信息学报, 2011, 33(2): 424–428. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX201102031.htm Zhang Qun-ying and Fang Guang-you. The study of pseudo random sequence’s application to GPR[J].Journal of Electronics&Information Technology, 2011, 33(2): 424–428. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX201102031.htm
[9]	Nicollin F, Barbin Y, Kofman W,et al.. An HF bi-phase shift keying radar: Application to ice sounding in Western Alps and Spitsbergen glaciers[J].IEEE Transactions onGeoscience and Remote Sensing, 1992, 30(5): 1025–1033. doi: 10.1109/36.175337
[10]	Utsi V. Design of a GPR for deep investigations[C]. 2007 4th International Workshop on Advanced Ground Penetrating Radar, Porto Alegre, RS, Brazil, 2007: 222–225.
[11]	Xia Z, Zhang Q, Ye S,et al.. A novel low-frequency coded ground penetrating radar for deep detection[J].IEICE Electronics Express, 2015, 12(11): 20150200. https://www.jstage.jst.go.jp/article/elex/12/11/12_12.20150200/_article
[12]	Xia Z, Fang G, Ye S,et al.. A novel handheld pseudo random coded UWB radar for human sensing applications[J].IEICE Electronics Express, 2014, 11(23): 20140981. https://www.researchgate.net/profile/Fang_Guangyou/publication/278399450_A_novel_handheld_pseudo_random_coded_UWB_radar_for_human_sensing_applications/links/573e80d708ae9ace841136ed.pdf?origin=publication_detail
[13]	Andrieu J, Gallais F, Mallepeyre V,et al.. Land mine detection with an ultra-wideband SAR system[C]. Proceedings SPIE 4742 Detection and Remediation Technologies for Mines and Minelike Targets VII, 2002, 4742: 237–247. http://proceedings.spiedigitallibrary.org/proceeding.aspx?articleid=888254
[14]	Fishler E, Haimovich A, Blum R,et al.. MIMO radar: An idea whose time has come[C]. Proceedings of the 2004 IEEE Radar Conference, Philadelphia PA, USA, 2004: 71–78. http://ieeexplore.ieee.org/iel5/9199/29174/01316398.pdf
[15]	Nag S, Barnes M A, Payment T,et al.. Ultrawideband through-wall radar for detecting the motion of people in real time[C]. Proceedings SPIE 4744: Radar Sensor Technology and Data Visualization, 2002, 4744: 48–57. http://proceedings.spiedigitallibrary.org/mobile/proceeding.aspx?articleid=888676
[16]	Yang Y and Fathy A E. See-through-wall imaging using ultra wideband short-pulse radar system[C]. 2005 IEEE Antennas and Propagation Society International Symposium, Washington DC, USA, 2005, 3: 334–337. http://ieeexplore.ieee.org/document/1552508/
[17]	Yang Y and Fathy A E. Development and implementation of a real-time see-through-wall radar system based on FPGA[J].IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(5): 1270–1280. doi: 10.1109/TGRS.2008.2010251
[18]	Shize G, Shenghe S, and Zhongting Z. A novel equivalent sampling method using in the digital storage oscilloscopes[C]. 1994 IEEE Instrumentation and Measurement Technology Conference, Hamamatsu, Japan, 1994, 2: 530–532. http://ieeexplore.ieee.org/document/351901/

施引文献

期刊类型引用(6)

1.	郁成阳，周婉婷，徐海洲，刘磊，杨君. 非合作式双基地雷达系统设计与验证. 雷达科学与技术. 2024(02): 119-125 . 百度学术
2.	吕敏，孙清洋，刘乐. 基于非合作雷达的外源雷达直达波参数提取方法. 电子技术与软件工程. 2023(04): 64-69 . 百度学术
3.	庄敬敏，宋杰，熊伟，王中训. 非合作辐射源双基地雷达技术综述. 探测与控制学报. 2023(04): 35-45+69 . 百度学术
4.	蒋威，马超，郑园园，曹健. 基于高频外辐射信号的探测效能评估研究. 电子技术与软件工程. 2022(14): 88-93 . 百度学术
5.	裴家正，黄勇，陈宝欣，关键，蔡咪，陈小龙. 联合脉压与Radon傅里叶变换的长时间相参积累方法. 雷达学报. 2021(06): 956-969 . 本站查看
6.	朱拥建，刘远，石林艳，王宏. 非合作侦察定位系统的关键技术. 太赫兹科学与电子信息学报. 2018(03): 452-457 . 百度学术

其他类型引用(4)

资源附件(0)

访问统计

图(16) / 表(2)

计量

文章访问数: 3064
HTML全文浏览量: 621
PDF下载量: 836
被引次数: 10

1. 引言
2. NPBR直达波参数估计
2.1 信号模型
2.2 直达波参数估计
3. 长时间积累
4. 数值实验
5. 总结

一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计

DOI: 10.12000/JR16076

通讯作者:
张晓娟 xjzhang@mail.ie.ac.cn

计量

Design of an Ultra-wideband Pseudo Random Coded MIMO Radar Based on Radio Frequency Switches

1. 引言

2. NPBR直达波参数估计

2.1 信号模型

2.2 直达波参数估计

3. 长时间积累

4. 数值实验

5. 总结

期刊类型引用(6)

其他类型引用(4)

计量

目录

1. 引言

2. NPBR直达波参数估计

2.1 信号模型

2.2 直达波参数估计

3. 长时间积累

4. 数值实验

5. 总结

期刊介绍

联系我们

一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计

DOI: 10.12000/JR16076

通讯作者: 张晓娟 xjzhang@mail.ie.ac.cn

计量

出版历程

Design of an Ultra-wideband Pseudo Random Coded MIMO Radar Based on Radio Frequency Switches

1. 引言

2. NPBR直达波参数估计

2.1 信号模型

2.2 直达波参数估计

3. 长时间积累

4. 数值实验

5. 总结

期刊类型引用(6)

其他类型引用(4)

计量

出版历程

目录

1. 引言

2. NPBR直达波参数估计

2.1 信号模型

2.2 直达波参数估计

3. 长时间积累

4. 数值实验

5. 总结

期刊介绍

联系我们

通讯作者:
张晓娟 xjzhang@mail.ie.ac.cn