机载多通道SAR运动目标方位向速度和法向速度联合估计算法

姜文; 牛杰; 吴一戎; 梁兴东

doi:10.11999/JEIT190672

机载多通道SAR运动目标方位向速度和法向速度联合估计算法

doi: 10.11999/JEIT190672

姜文^{1, 2, 3, ,},
牛杰^{1, 2, 3},
吴一戎^{1, 2},
梁兴东^{1, 2}

1.
中国科学院空天信息创新研究院北京 100094
2.
微波成像技术国家重点实验室北京 100190
3.
中国科学院大学北京 100049

详细信息

作者简介:
姜文：男，1990年生，博士生，研究方向为机载多通道SAR地面动目标参数估计

牛杰：男，1993年生，博士生，研究方向为新体制SAR海面纹理增强

吴一戎：男，1963年生，中国科学院院士，研究方向为微波成像理论、微波成像技术、雷达信号处理等

梁兴东：男，1973年生，研究员，研究方向为高分辨率SAR系统、干涉SAR系统、成像处理及应用、实时数字信号处理等

通讯作者:
姜文　jiangwen13@mails.ucas.ac.cn

中图分类号: TN957.52
计量
- 文章访问数: 2784
- HTML全文浏览量: 1081
- PDF下载量: 112
- 被引次数: 2
出版历程
- 收稿日期: 2019-09-02
- 修回日期: 2020-02-14
- 网络出版日期: 2020-03-04
- 刊出日期: 2020-06-22

Joint Estimation Algorithm for Azimuth Velocity and Normal Velocity of Moving Targets in Airborne Multi-channel SAR

Wen JIANG^{1, 2, 3
, ,},
Jie NIU^{1, 2, 3},
Yirong WU^{1, 2},
Xingdong LIANG^{1, 2}

1.
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
2.
National Key Laboratory of Microwave Imaging Technology, Beijing 100190, China
3.
University of Chinese Academy of Sciences, Beijing 100049, China

摘要

摘要:
对运动目标进行SAR成像时，参数估计是必不可少的。现有算法主要针对运动目标的径向速度和方位向速度进行估计，而对3维运动目标的法向速度无法估计。该文利用L型基线的机载多通道SAR系统，提出一种方位向速度和法向速度的联合估计算法。该算法在距离-多普勒域提取运动目标信号，并利用多幅SAR图像之间的相位差进行方位向速度和法向速度的联合估计。该算法不依赖图像配准，不需要解多普勒模糊，因此具有较高的估计精度和鲁棒性，有较强的实际意义和应用价值。
- 机载多通道SAR /
- 运动目标检测 /
- 方位向速度 /
- 法向速度 /
- 速度估计
Abstract:
Parameter estimation is essential for SAR imaging of moving targets. The existing algorithms mainly estimate the radial velocity and azimuth velocity of the moving target, but the normal velocity of the three-dimensional moving target can not be estimated. In this paper, a joint estimation algorithm of azimuth velocity and normal velocity is proposed by using an airborne multi-channel SAR system with L-shaped baseline. The algorithm extracts the moving target signal in Range-Doppler domain, and estimates the azimuth and normal velocity jointly using the phase differences between multiple SAR images. The algorithm does not rely on image registration, does not need to solve Doppler ambiguity. Therefore, the algorithm has high estimation accuracy and robustness, and has strong practical significance and application value.
- Airborne multi-channel SAR /
- Moving target detection /
- Azimuth velocity /
- Normal velocity /
- Velocity estimation

HTML全文

1. 引言

因轨道角动量(Orbital Angular Momentum, OAM)模态间正交性和完整性优势，OAM复用能够在有限信道资源上实现信息高质量和大容量传输，携带OAM的电磁涡旋通信技术能潜在地增加无线通信链路的系统容量和频谱效率^[1-7]，被视为第6代无线通信网络的关键技术之一。构建准确的涡旋信道模型是涡旋通信系统设计、性能评估、网络部署和优化的首要任务。

为了实现电磁涡旋通信系统，亟需在准确刻画OAM信道特性基础上，建立能够复现涡旋信号波前结构和幅值特性的OAM信道模型。建模实质是表征拓扑荷和传播距离等因素对涡旋信号幅度与相位的影响。现有的OAM信道模型主要包括自由空间传播模型^[8]、利用螺旋相位因子的涡旋信道模型^[9]、基于电磁场理论的涡旋信道模型^[10]、基于拉盖尔-高斯模的涡旋信道模型^[11]，以及基于镜面反射的多径信道模型^[12]。文献[8]建立了基于失调均匀同心圆形阵列的自由空间涡旋信道模型，并对多模式复用OAM多输入多输出(Multiple Input Multiple Output, MIMO)系统的频谱效率和能量效率进行了研究，结果表明频谱效率和能量效率优于传统MIMO系统。文献[13]建立了基于拉盖尔-高斯理论的自由空间涡旋信道模型，研究了信道传递函数、时域函数及路径损耗，结果表明与自由空间传播模型相比，提出模型能够更准确地刻画涡旋信道特性。文献[14]建立了基于螺旋相位因子的多天线电磁涡旋信道模型，利用该模型分析了最佳相位条件下的涡旋信道容量，结果表明相比传统MIMO系统，多天线电磁涡旋系统具有更高信道容量。上述研究表明，构建OAM信道是研究涡旋波传播特性以及进行OAM通信系统性能评估的前提^[15]。

现有研究存在以下不足：一是在自由空间传播环境中构建的涡旋信道，仅考虑了路径损耗对涡旋信号的影响，而对涡旋信号相位结构变化的研究还不足；二是OAM通信系统更多是部署在多径环境中，自由空间传播环境下建立的涡旋信道不能准确地刻画出多径环境对涡旋信号的影响。针对多径传播，Yan等人^[16]分析、仿真和测量了两径OAM信道(1条直射径和1条镜面反射径)引起的信道内和信道间串扰，结果表明OAM信道模态数越大，信道内串扰和信道间串扰越严重。基于光学射线理论，文献[17,18]通过计算传播距离，建立了OAM-MIMO通信系统的两径信道模型和峡谷环境下的6径信道模型，结果表明多径效应引起的相位畸变和信道串扰随着模态数增大而增大，模态数较小的涡旋信号可以获得更好的传输增益。Liang等人^[19]构建了混合正交频分复用通信系统的稀疏多径信道模型，该模型包含1条直射径、1条1次反射径、1条2次反射径和1条3次反射径，能够描述简单通信场景下的涡旋信道特性。然而，OAM通信技术被期望应用到“室内、城区、隧道、高速铁路”等热点与特殊场景中，复杂多样的建筑物引起的多径干扰会改变波前相位，破坏涡旋波的形态，致使OAM信道呈现出复杂性、特殊性和时变性等特点。现有的自由空间传播信道模型、两径信道模型、稀疏多径信道模型尚不能准确刻画真实场景下OAM信道多径特性，因此有必要针对实际应用场景中多径传播机制构建准确的OAM信道模型。

信道建模方法主要有确定性建模法和统计性建模法^[20]。现有的OAM信道建模本质上是确定性建模方法，核心思想是在已知传播环境详细电磁参数前提下，利用电波传播理论或光学射线理论来预测传播距离，因此严格依赖于传播场景及建筑物的分布，不适用于复杂的室内环境。在室内环境，毫米波能实现高速通信，但毫米波OAM信道特性尚未明确，尤其是波前相位和幅度随多径传播条件、模态数和传播距离等影响因子的变化规律，毫米波OAM信道统计性模型尚无。

针对以上问题，本文开展走廊场景毫米波OAM信道特性分析和统计建模研究，为室内毫米波OAM通信链路与系统仿真提供涡旋信道基础数据和统计模型。首先，在走廊场景中建立基于均匀圆形天线阵列(Uniform Circle Array, UCA)的OAM辐射传输系统；其次，基于光学射线理论推导视距(Line-of-Sight, LoS)和非视距(Non-LoS, NLoS)传播条件下的信号传播距离，结合UCA辐射理论构建OAM信道理论模型；然后，根据多径传播条件下毫米波OAM信道的波前相位和幅度统计值，结合拟合误差确定统计模型及参数；最后，研究不同传播条件、模态数和传播距离下波前相位与幅度衰落特性。

2. 走廊传播环境

矩形长直狭窄走廊环境由地面、天花板和两面墙壁组成，分别是花岗岩、石膏板和白色瓷釉面砖材料，反射系数用Γ₁, Γ₂和Γ₃表示。地面与天花板间距为3 m，两面墙壁间的距离为2.5 m，长度为3 m。图1(a)和图1(b)分别构建了LoS和NLoS传播条件的OAM辐射传输系统。考虑到天线阵列能够在所有频率范围内产生不同模态的涡旋波，系统模型的发射端是由N个阵元组成的UCA天线阵列，阵元用 $n$ 表示， $0 \le n \le N - 1$ ，接收端使用单个天线。图1(a)中收发端天线之间的第1菲涅尔区内不存在障碍物，存在1条直射径。图1(b)中存在1个几何尺寸略大于发射阵列尺寸的障碍物，导致仅有非视距信号。本文假设发射端UCA的中心 $O$ 与接收端天线的水平距离为 $D$ ，UCA阵列半径为R_t，阵列中心 $O$ 距地面、天花板、墙1和墙2的距离分别为h₁, h₂, h₃和h₄。

图 1 OAM辐射传输系统

下载: 全尺寸图片幻灯片

3. 涡旋信道模型

本节基于光学射线理论建立了走廊场景LoS和NLoS传播条件的多径模型，推导传播路径的距离，基于UCA辐射理论建立OAM信道模型。

3.1 传播模型

为了更加直观地体现多径传播，图2(a)和图2(b)分别给出了狭窄走廊环境下LoS和NLoS传播条件多径传播模型的侧视图和俯视图。涡旋信号在走廊环境中传播时会发生直射、反射和散射等多径传播机制，使得接收端接收到的是来自不同路径涡旋信号的叠加信号。为了便于分析，本文只考虑直射传播和反射传播机制，且因3次及以上反射路径的信号强度过低，只考虑1次与2次反射路径。

图 2 走廊环境多径传播模型

下载: 全尺寸图片幻灯片

图2(a)为LoS传播条件下多径传播模型的侧视图，传播路径包含1条直射路径，1条地面反射路径GR、1条天花板反射路径CR、2条墙壁反射路径WR、1条地面与天花板反射路径GC、2条地墙反射路径GW、1条天花板与地面反射路径CG、2条天花板与墙壁反射路径CW、2条双墙反射路径TW、2条墙地反射路径WG、2条墙壁与天花板反射路径WC等16条反射路径。图2(b)为NLoS传播条件下多径传播模型的俯视图，收发天线间不存在直射路径，只包含了16条反射路径。

以发射端的UCA中心 $O$ 为坐标原点建立直角坐标系。设观测点坐标为 $(x,y,z)$ ，则直射径上发射阵元 $n$ 到观测点 $(x,y,z)$ 的传播距离 $d_n^{\text{D}}$ 表示为

$d_n^{D} = \sqrt {{{(x - {R_{\rm{t}}}\sin {\varphi _n})}^2} + {{(y - {R_{\rm{t}}}\cos {\varphi _n})}^2} + {z^2}}$

(1)

其中， ${\varphi _n} = 2\pi n/N$ 为发射阵元 $n$ 的方位角。

由镜像理论可知，反射路径上的涡旋信号可以看作由镜像发射端发射的反模态涡旋信号^[21]，则发射阵元 $n$ 到观测点 $(x,y,z)$ 的GR路径的传播距离 $d_n^{{\text{GR}}}$ 可以表示为

$\begin{split} d_n^{{\text{GR}}} =& [{x^2} + {y^2} + {z^2} + R_{\text{t}}^{\text{2}} - 2{R_{\text{t}}}(x\sin {\varphi _n} + y\cos {\varphi _n}) \\ &+ 4{h_1}(y + {h_1} - {R_{\text{t}}}\cos {\varphi _n}){]^{\frac{1}{2}}}\\[-10pt] \end{split}$

(2)

以GC和GW为例，求2次反射路径的传播距离。若墙壁为墙1，则发射阵元 $n$ 到观测点 $(x,y,z)$ 的GC和GW2次反射径传播距离 $d_n^{{\text{GC}}}$ 和 $d_n^{{\text{GW}}}$ 分别表示为

$\begin{split} d_n^{{\text{GC}}} =& [{x^2} + {y^2} + {z^2} + R_{\text{t}}^{\text{2}} - 2{R_{\text{t}}}(x\sin {\varphi _n} + y\cos {\varphi _n})\\ &+ 4({h_1} + {h_2})( - y + {h_1} + {h_2} + {R_{\text{t}}}\cos {\varphi _n}){]^{\frac{1}{2}}} \end{split}$

(3)

$\begin{split} d_n^{{\text{GW}}} =& [{x^2} + {y^2} + {z^2} + R_{\text{t}}^{\text{2}} - 2{R_{\text{t}}}(x\sin {\varphi _n} + y\cos {\varphi _n})\\ &+ 4{h_1}(y + {h_1} - {R_{\text{t}}}\cos {\varphi _n}) \\ &- 4{h_3}(x - {h_3} - {R_{\text{t}}}\sin {\varphi _n}){]^{\frac{1}{2}}}\\[-10pt] \end{split}$

(4)

以此类推，可以推导出CR, WR, CG, CW, TW, WG和WC反射路径的传播距离。

3.2 信道模型

收发天线静止不动时传播信道传递函数可以表示为

$h(t) = \sum\limits_{k = 0}^{K - 1} {{a_k}\delta (t - {\tau _k}){{\rm{e}}^{{\text{j}}{\theta _k}}}}$

(5)

其中，K表示多径数，a_k, τ_k, θ_k分别表示幅度、传播时延和相位序列。对于静态传播信道，幅度和相位可以通过电场强度E_k(d_n)获得，d_n为不同路径下第n个阵元到观测点的距离，所以式(5)可以表示为

$h(t) = \sum\limits_{k = 0}^{K - 1} {{E_k}({d_n})\delta (t - {\tau _k})}$

(6)

本文中，LoS传播条件下K为17，NLoS传播条件下K为16，时延τ_k可由各条反射径与直射径的传播距离差除以光速获得。

当发射端的UCA辐射出模态数为l(模态数满足 $\left| l \right| \le N/2$ ^[22])的涡旋波时，需要对每个阵元进行等幅差相馈电，即第 $n$ 个阵元的馈电相位为 $2\pi nl/N$ 。根据UCA辐射理论^[23]，在球坐标系中，直射径上观测点 $(r,\theta ,\varphi )$ 处的电场强度可以表示为

${E_{\text{D}}}(d_n^{\text{D}}) = {E_{\text{D}}}(r,\theta ,\varphi ) = \sum\limits_{n = 0}^{N - 1} {{{\text{e}}^{{\text{j}}2\pi nl/N}}\frac{{{{\text{e}}^{ - {\text{j}}kd_n^{\text{D}}}}}}{{d_n^{\text{D}}}}}$

(7)

其中，r为远场点到发射端UCA中心 $O$ 的距离， $\theta$ 为俯仰角， $\varphi$ 为方位角， $k = 2\pi /\lambda$ 为波数， $\lambda$ 为波长。

由于天线尺寸R_t远小于传输距离 $d_n^{\text{D}}$ ，式(7)的相位项中 $d_n^{\text{D}} \simeq r - {R_{\text{t}}}\sin \theta \cos (\varphi - 2\pi n/N)$ ，幅度项中 $d_n^{\text{D}} \simeq r$ ，进而式(7)可以简化为

$\begin{split} {E_{\text{D}}}(r,\theta ,\varphi ) = &\frac{{{{\rm{e}}^{ - {\text{j}}kr}}}}{r}\sum\limits_{n = 0}^{N - 1} {\exp ( - {\text{j}}\sin \theta \cos (\varphi - 2\pi n/N))}\\ &\cdot \exp ( - {\text{j}}{R_{\text{t}}}2\pi nl/N)\\[-10pt] \end{split}$

(8)

当天线阵元数N足够大时，式(8)简化为

${E_{\text{D}}}(r,\theta ,\varphi ) = \frac{{{{\rm{e}}^{ - {\text{j}}kr}}}}{r}N{{\text{j}}^{ - l}}{{\rm{e}}^{{\text{j}}l\varphi }}{J_l}(k{R_{\text{t}}}\sin \theta )$

(9)

其中， ${J_l}( \cdot )$ 表示第1类l阶的贝塞尔函数。利用公式(9)便可求得各条路径上的电场强度。

将球坐标系转换为直角坐标系，发射端UCA中心O为坐标原点，则观测点 $(x,y,z)$ 处直射径、GR 1次反射径和GC/GW 2次反射径的电场强度分别表示为式(10)—式(12)。

由此，可以依次推导出CR, WR, CG, CW, TW, WG和WC路径上的电场强度，最终获得LoS和NLoS传播条件的OAM信道传递函数。

4. 结果与分析

本节基于OAM信道模型研究不同传播条件下毫米波频段涡旋波波前相位和幅度特性。以相位和幅度的统计结果为拟合目标，使用不同的概率统计模型进行拟合，采用均方根误差(Root Mean Squared Error, RMSE)作为评估指标，验证统计模型的准确性，并且对模型进行统计性分析。具体地，使用均匀分布描述不同传播条件下的OAM信道波前相位；使用莱斯因子确定符合的概率分布模型，利用Nakagami-m模型的衰落因子m确定不同传播条件下的幅度分布。

$\begin{split} {E_{\text{D}}}(x,y,z) =& \frac{{N{{\text{j}}^{- l}}}}{{\sqrt {{x^2} + {y^2} + {z^2}} }}{{\text{e}}^{{{ - {\rm{j}}}}\bigr[l\arctan (\frac{y}{x}){\text{ + }}k\sqrt {{x^2} + {y^2} + {z^2}} \bigr]}}\\ & \cdot {J_l}\left(\frac{{k{R_{\text{t}}}\sqrt {{x^2} + {y^2}} }}{{\sqrt {{x^2} + {y^2} + {z^2}} }}\right)\\[-21pt] \end{split}$

(10)

$\begin{split} &{E_{{\rm{GC}}}}(x,y,z) = \frac{{{\varGamma _1}{\varGamma _2}N{{\text{j}}^{ - l}}}}{{\sqrt {{x^2} + {{(y - 2{h_1} - 2{h_2})}^2} + {z^2}} }}\\ & \quad \cdot{{\text{e}}^{{\text{j}}[l\cdot \arctan (\frac{{y - 2{h_1} - 2{h_2}}}{x}) - k\sqrt {{x^2} + {{(y - 2{h_1} - 2{h_2})}^2} + {z^2}} ]}}\\ & \quad \cdot{J_l}\left(\frac{{k{R_{\text{t}}}\sqrt {{x^2} + {{(y - 2{h_1} - 2{h_2})}^2}} }}{{\sqrt {{x^2} + {{(y - 2{h_1} - 2{h_2})}^2} + {z^2}} }}\right) \end{split}$

(11)

$\begin{split} &{E_{{\rm{GW}}}}(x,y,z) = \frac{{{\varGamma _1}{\varGamma _3}N{{\text{j}}^{ - l}}}}{{\sqrt {{{(x - 2{h_3})}^2} + {{(y + 2{h_1})}^2} + {z^2}} }}\\ & \quad \cdot{{\text{e}}^{{\text{j}}\bigr[l\cdot\arctan (\frac{{y + 2{h_1}}}{{x - 2{h_3}}}) - k\sqrt {{{(x - 2{h_3})}^2} + {{(y + 2{h_1})}^2} + {z^2}} \bigr]}}\\ & \quad \cdot{J_l}\left(\frac{{k{R_{\text{t}}}\sqrt {{{(x - 2{h_3})}^2} + {{(y + 2{h_1})}^2}} }}{{\sqrt {{{(x - 2{h_3})}^2} + {{(y + 2{h_1})}^2} + {z^2}} }}\right)\\[-21pt] \end{split}$

(12)

假设工作频率为30 GHz，阵列半径为0.02 m，由于阵元尺寸和阵列半径的限制，将阵元数设置为8，中心 $O$ 与地面、天花板和两面墙壁的距离分别为1.5 m, 1.5 m, 1.2 m和1.3 m，地面、天花板和墙壁的反射系数分别为0.7, 0.91和0.75^[24]。仿真采用的软件为MATLAB，主要参数如表1所示。

表 1 仿真参数设置

参数	数值
工作频率	30 GHz
发射端UCA的阵元数量	8
发射端UCA的半径	0.02 m
发射端UCA中心距地面的距离	1.5 m
发射端UCA中心距天花板的距离	1.5 m
发射端UCA中心与墙1间的距离	1.2 m
发射端UCA中心与墙2间的距离	1.3 m
地面材料的反射系数	0.70
天花板材料的反射系数	0.91
墙壁材料的反射系数	0.75

下载: 导出CSV

| 显示表格

4.1 波前相位

由文献[13]可知，OAM通信系统适用于短距离传输，因此设两个观测面距发射端UCA的距离分别为0.4 m和1 m，大小均为 ${\text{0}}{\text{.4 m}} \times {\text{0}}{\text{.4 m}}$ ，模态数 $l$ 分别为1, 2, 3。

图3是OAM多径信道的波前相位，图3(a)、图3(c)和图3(e)的观测面距发射端UCA的距离为1 m，图3(b)、图3(d)和图3(f)的观测面为0.4 m，图3(a)和图3(b)是直射径传播，图3(c)和图3(d)是3条传播路径(1条直射径和2条1次反射径)，图3(b)和图3(e)是5条传播路径(1条直射径和4条1次反射径)。结果表明，相比1条直射径传播，多条路径传播的OAM信道相位发生了畸变，模态数越高，畸变越严重；对比不同传播路径，相同传播距离传播路径越多，相位畸变越严重，这是因为多径传播引起的信道间和信道内串扰，致使相位产生畸变，且随着模态数增大和传播路径增多而增大。对比不同传播距离，相同传播条件传播距离越远，相同观测平面内的信道相位畸变越严重，辨别OAM信道模态数的难度越大，这是由OAM波的发散性导致的。

图 3 不同模态的OAM信道波前相位

下载: 全尺寸图片幻灯片

基于图3的结果，对不同模态的OAM相位进行统计，由于1, 2, 3模态的概率密度分布相似，本文只给出模态数为1和2的相位分布，如图4所示，其中直方图为概率分布。结果表明信道相位符合均匀分布，概率分布表示为

图 4 OAM信道相位的统计概率及均匀分布拟合结果

下载: 全尺寸图片幻灯片

$f(x,a,b,{w})=\left\{\begin{aligned} & \frac{{w}}{b-a},a < x < b\\ & 0,其他 \end{aligned}\right.$

(13)

其中，w是调整拟合曲线的参数，a和b分别为数轴上的最小值和最大值。在表征波前相位分布中，a可以取0°，b可以取360°。图4中实线为均匀分布的拟合结果。

不同传播条件下的均匀分布最佳拟合参数以及拟合误差如表2所示。结果表明波前相位分布值与均匀分布间的拟合误差RMSE约为0.5%，不同观测距离间的拟参数w差距在0.5以内，这说明均匀分布可以准确地表征OAM多径信道的波前相位统计特性。

表 2 最佳拟合参数w与拟合误差

传播距离(m)	模态数		直射径	3条路径	5条路径
1	l=1	w	12.46	12.24	12.23
	l=1	RMSE	0.0050	0.0046	0.0047
	l=2	w	12.32	12.42	12.59
	l=2	RMSE	0.0047	0.0048	0.0045
0.4	l=1	w	12.43	12.25	12.67
	l=1	RMSE	0.0039	0.0051	0.0045
	l=2	w	12.56	12.38	12.33
	l=2	RMSE	0.0042	0.0042	0.0044

下载: 导出CSV

| 显示表格

4.2 信道幅度

为了获得多径传播条件下OAM信道幅度特性，图5给出了不同传播条件下不同传播距离的OAM信道幅度的变化规律，其中图5(a)为直射径传播，图5(b)为视距传播(1条直射径、4条1次反射径和12条2次反射径)，图5(c)为非视距传播(4条1次反射径和12条2次反射径)。

图 5 OAM信道幅度

下载: 全尺寸图片幻灯片

从图5可以观察到，直射径和LoS传播条件下，OAM信道幅度随模态数增大而逐渐减小，相同传播距离时模态数越大，信道幅度衰减也越大；当传播距离小于0.4 m时，直射径和LoS传播条件的OAM信道幅度的动态变化范围为0～50 dBμV/m, NLoS传播条件的信道幅度仅仅在0～10 dBμV/m间波动，这是因为NLoS传播条件下，反射径对接收信号的贡献远远低于存在直射径的传播情况；当传播距离大于0.4 m时，直射径传播的信道幅度发生了平滑衰减，LoS和NLoS传播条件的信道幅度在0～10 dBμV/m间剧烈波动，这是因为OAM波的发散性使远距离传输时相位畸变严重，导致幅度剧烈波动，多径效应也导致远距离传输的信道幅度衰减增大。

LoS传播条件下，直射径功率与反射径功率之比定义为莱斯因子，表示为K，不同模态OAM多径信道的K值随传播距离的变化规律如图6所示。从图6可以观察到，当传播距离较小时(小于0.4 m)，K值在0～20 dB之间波动，且由图5(a)和图5(b)可知，LoS传播条件下的信道幅度衰减与只有直射径传播的信道幅度衰减情况相似，此时直射径信号是接收信号中的主导分量，信道幅度分布可能接近于莱斯分布。当传播距离较大时(大于0.4 m)，K值明显小于0 dB，且由图5(b)和图5(c)可知，LoS传播条件下的信道幅度衰减与NLoS传播条件下的信道幅度衰减情况相似，此时反射径信号是接收信号中的主导分量，信道幅度分布可能接近于瑞利分布。因此，本文提出利用Nakagami-m分布模型描述不同传播条件的OAM信道幅度分布。

图 6 莱斯因子随传播距离变化图

下载: 全尺寸图片幻灯片

基于图5(b)和图5(c)的涡旋信号幅度，统计不同模态下较大传输距离(2 m)时OAM多径信道幅度的概率值，使用Nakagami-m分布模型描述信道幅度分布，根据衰落因子m判断信道幅度所服从的概率分布：

(1)m=1时，Nakagami-m分布退化为瑞利分布；

(2)m>1时，Nakagami-m分布与莱斯分布接近。

Nakagami-m分布模型表示为

$f(x,{{c}},m,\varOmega ) = {{c}} \cdot \frac{{2{m^m}{x^{2m - 1}}}}{{\varGamma (m){\varOmega ^m}}}{{\rm{e}}^{ - \frac{{m{x^2}}}{\varOmega }}}$

(15)

其中，c为调整拟合曲线的参数， $\varGamma (m)$ 表示Gamma函数， $\varOmega = E({x^2}) = 2{\sigma ^2}$ 表示平均功率， $\sigma$ 表示多径信号平均功率的标准偏差。

图7给出了LoS和NLoS传播条件下的OAM多径信道幅度分布，直方图为信道模型统计值，实线为Nakagami-m分布拟合值。拟合参数及拟合误差如表3所示。

图 7 OAM多径信道幅度的统计概率及拟合结果

下载: 全尺寸图片幻灯片

表 3 Nakagami-m分布拟合参数和拟合误差

传播条件	参数及误差	模态值
传播条件	参数及误差	l=0	l=1	l=2	l=3
LoS	m	0.864	0.823	1.006	0.904
	c	2.649	2.205	2.275	2.448
	$\sigma$	7.426	3.393	2.754	2.246
	RMSE	0.0103	0.0137	0.0094	0.0061
NLoS	m	1.053	1.140	1.080	0.914
	c	0.494	0.644	0.665	0.521
	$\sigma$	2.846	3.015	3.215	3.319
	RMSE	0.0080	0.0108	0.0054	0.0130

下载: 导出CSV

| 显示表格

从图7和表3可以看出，信道幅度概率分布与Nakagami-m概率分布模型有较好的拟合程度，且衰落因子m都接近1，因此LoS和NLoS两种传播条件下的信道幅度概率分布都接近于瑞利分布。LoS传播条件下， $\sigma$ 随着模态数增大而减小，即多径信号平均功率减小，接收信号中直射径占主导地位；统计值与拟合值间的RMSE随模态数增加而减小，高模态的拟合误差低至0.6%，这说明瑞利分布模型能够较好地刻画出信道条件差的OAM 多径信道幅度变化；与LoS传播条件相反，NLoS 传播条件下 $\sigma$ 随模态数增大而增大，说明信道幅度的衰减随模态数增大而减小，多径传播有利于高模态OAM多径信号接收。此外，信道模态数相同时，LoS传播条件下 $\sigma$ 始终大于NLoS传播条件下的 $\sigma$ ，这说明多径传播使信道幅度衰减更加严重。

统计短距离(0.4 m)LoS传播条件下的OAM多径信道幅度概率值，使用Nakagami-m分布模型描述信道幅度值分布，幅度值概率分布和拟合结果分别如图8所示。直方图为幅度的概率分布，实线为Nakagami-m分布的拟合曲线。拟合参数及其误差如表4所示。

图 8 短距离传播LoS条件下的OAM信道幅度统计概率及拟合结果

下载: 全尺寸图片幻灯片

表 4 LoS传播条件Nakagami-m分布拟合参数和拟合误差

拟合参数	模态数l
拟合参数	0	1	2	3
m	8.703	1.894	1.725	1.534
c	1.843	2.224	2.118	4.396
$\sigma$	12.247	15.811	19.102	22.624
RMSE	0.0105	0.0126	0.0139	0.0107

下载: 导出CSV

| 显示表格

图8和表4表明，信道幅度概率分布与Nakagami-m分布模型有较好的拟合程度，且m都大于1，因此LoS传播条件下，较短距离传播时的信道幅度概率分布更接近于莱斯分布。拟合误差RMSE在1%～1.5%范围内， $\sigma$ 随着模态数增大而增大，说明LoS传播条件下短距离传播时，OAM 信道幅度的衰减随模态数增大而减小。

5. 结束语

本文针对走廊环境下的毫米波OAM通信，提出了一种基于Nakagami-m分布的多径信道建模方法，研究了模态数、传播距离和传播条件对相位结构与信道幅度的影响。结果表明相位波前结构分布与模态数、观察距离和多径传播条件密切相关，且服从均匀分布；较大距离多径传输致使相位畸变严重，信道幅度波动较大，衰减较大，信道幅度分布接近于瑞利分布；较小距离LoS传播时，信道幅度分布接近于莱斯分布，较小模态数的涡旋信号有利于增大接收端的接收功率。本文基于Nakagami-m分布拟合的长距离和短距离OAM多径信道参数，可以用于描述不同模态数、不同传播距离和不同传播条件下的 OAM 多径衰落信道特性，为室内场景毫米波OAM通信系统设计与性能评估提供了有用信息。

图 1 机载SAR 3维运动目标成像模型

下载: 全尺寸图片幻灯片

图 2 机载多通道SAR 3维运动目标成像模型

下载: 全尺寸图片幻灯片

图 3 速度联合估计算法处理流程

下载: 全尺寸图片幻灯片

图 4 运动目标1的多普勒域相位差

下载: 全尺寸图片幻灯片

图 5 运动目标2的多普勒域相位差

下载: 全尺寸图片幻灯片

图 6 运动目标3的多普勒域相位差

下载: 全尺寸图片幻灯片

图 7 运动目标聚焦前后的SAR图像

下载: 全尺寸图片幻灯片

图 8 运动目标方位向剖面图

下载: 全尺寸图片幻灯片

表 1 SAR系统仿真基本参数

中心频率	飞行速度	飞行高度	最近斜距	PRF	信噪比	多普勒带宽	顺轨基线	交轨基线
9.6 GHz	150 m/s	7500 m	20 km	2000	10 dB	500 Hz	1 m	1 m

下载: 导出CSV

表 2 目标3维运动参数(m/s)

	径向向速度	方位向速度	法向速度
运动目标1	0	10	10
运动目标2	0	10	0
运动目标3	0	0	10

下载: 导出CSV

表 3 运动目标速度估计结果及误差(m/s)

	理论方位向速度	理论法向速度	估计方位向速度	估计法向速度	方位向速度估计误差	法向速度估计误差
运动目标1	10.00	10.00	9.85	10.00	0.15	0
运动目标2	10.00	0	9.84	1.15×10^–15	0.16	1.15×10^–15
运动目标3	0	10.00	7.93×10^–2	9.90	7.93×10^–2	0.10

下载: 导出CSV

表 4 SAR平台部分参数

距离向采样点数	PRF	天线间距	飞行速度	飞行高度	带宽
16384	781	1 m	146 m/s	7544 m	300 MHz

下载: 导出CSV

参考文献(19)

PERRY R P, DIPIETRO R C, and FANTE R L. SAR imaging of moving targets[J]. IEEE Transactions on Aerospace and Electronic Systems, 1999, 35(1): 188–200. doi: 10.1109/7.745691

FIENUP J R. Detecting moving targets in SAR imagery by focusing[J]. IEEE Transactions on Aerospace and Electronic Systems, 2001, 37(3): 794–809. doi: 10.1109/7.953237

JAO J K. Theory of synthetic aperture radar imaging of a moving target[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(9): 1984–1992. doi: 10.1109/36.951089

MARQUES P A C and DIAS J M B. Velocity estimation of fast moving targets using a single SAR sensor[J]. IEEE Transactions on Aerospace and Electronic Systems, 2005, 41(5): 75–89. doi: 10.1109/TAES.2005.1413748

BARBAROSSA S and FARINA A. Detection and imaging of moving objects with synthetic aperture radar. 2. Joint time-frequency analysis by wigner-ville distribution[J]. IEE Proceedings F - Radar and Signal Processing, 1992, 139(1): 89–97. doi: 10.1049/ip-f-2.1992.0011

WANG Zhirui, XIA Xianggen, XU Jia, et al. Ground moving target imaging based on 2-D velocity search in high resolution SAR[C]. 2017 IEEE Radar Conference, Seattle, USA, 2017: 68–72.

SHI Hongyin, YANG Xiaoyan, ZHOU Qiuxiao, et al. SAR slow moving target imaging based on over-sampling smooth algorithm[J]. Chinese Journal of Electronics, 2017, 26(4): 876–882. doi: 10.1049/cje.2017.06.005

YANG Wei, CHEN Jie, LIU Wei, et al. Moving target azimuth velocity estimation for the MASA mode based on sequential SAR images[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(6): 2780–2790. doi: 10.1109/JSTARS.2016.2641744

SU Jia, TAO Haihong, WANG Ling, et al. Coherently integrated cubic function based doppler parameters estimation for moving-target imaging[C]. 2017 International Applied Computational Electromagnetics Society Symposium, Suzhou, China, 2017: 1–2.

WANG Hanyun and JIANG Yicheng. Real-time parameter estimation for SAR moving target based on WVD slice and FrFT[J]. Electronics Letters, 2018, 54(1): 47–49. doi: 10.1049/el.2017.1740

LI Zhongyu, WU Junjie, LIU Zhutian, et al. An optimal 2-D spectrum matching method for SAR ground moving target imaging[J]. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(10): 5961–5974. doi: 10.1109/TGRS.2018.2829166

魏翔飞, 王小青, 种劲松. 一种基于局域中心频率的SAR图像舰船方位向速度估计方法[J]. 电子与信息学报, 2018, 40(9): 2242–2249. doi: 10.11999/JEIT170991

WEI Xiangfei, WANG Xiaoqing, and CHONG Jinsong. Ship azimuthal speed estimation method based on local region doppler centroid in SAR images[J]. Journal of Electronics &Information Technology, 2018, 40(9): 2242–2249. doi: 10.11999/JEIT170991

王超, 王岩飞, 王琦, 等. 基于回波序列最小二乘拟合的高分辨率SAR运动目标速度估计[J]. 电子与信息学报, 2019, 41(5): 1055–1062. doi: 10.11999/JEIT180695

WANG Chao, WANG Yanfei, WANG Qi, et al. Velocity estimation of moving targets based on least square fitting of high-resolution SAR echo sequences[J]. Journal of Electronics &Information Technology, 2019, 41(5): 1055–1062. doi: 10.11999/JEIT180695

HE Xiongpeng, LIAO Guisheng, XU Jingwei, et al. Robust radial velocity estimation based on joint-pixel normalized sample covariance matrix and shift vector for moving targets[J]. IEEE Geoscience and Remote Sensing Letters, 2019, 16(2): 221–225. doi: 10.1109/LGRS.2018.2871950

WANG Genyuan, XIA Xianggen, and CHEN V C. Three-dimensional ISAR imaging of maneuvering targets using three receivers[J]. IEEE Transactions on Image Processing, 2001, 10(3): 436–447. doi: 10.1109/83.908519

ZHANG Qun and YEO T S. Three-dimensional SAR imaging of a ground moving target using the InISAR technique[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(9): 1818–1828. doi: 10.1109/TGRS.2004.831863

ZHANG Qun, YEO T S, DU Gan, et al. Estimation of three-dimensional motion parameters in interferometric ISAR imaging[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(2): 292–300. doi: 10.1109/TGRS.2003.815669

汤立波, 李道京, 洪文, 等. 基于顺轨-交轨InSAR技术的运动舰船目标三维成像[J]. 系统工程与电子技术, 2008, 30(9): 1669–1673. doi: 10.3321/j.issn:1001-506X.2008.09.017

TANG Libo, LI Daojing, HONG Wen, et al. Three-dimensional imaging of moving ships with 3D motion based on AT-CT InSAR[J]. Systems Engineering and Electronics, 2008, 30(9): 1669–1673. doi: 10.3321/j.issn:1001-506X.2008.09.017

尹建凤, 李道京, 王爱明, 等. 基于星载毫米波顺轨-交轨InISAR的空间运动目标三维成像技术研究[J]. 宇航学报, 2013, 34(2): 237–245. doi: 10.3873/j.issn.1000-1328.2013.02.013

YIN Jianfeng, LI Daojing, WANG Aiming, et al. Three-dimensional imaging technique of space moving target based on spaceborne along-cross track millimeter-wave In-ISAR[J]. Journal of Astronautics, 2013, 34(2): 237–245. doi: 10.3873/j.issn.1000-1328.2013.02.013

施引文献

期刊类型引用(2)

1.	南作用，钟志刚，陈任翔，王亚. 考虑传输损耗的无线通信毫米波信道建模分析. 计算机仿真. 2025(02): 531-535 . 百度学术
2.	吴杰，胡俊，张忠祥，沙威，黄志祥，吴先良. 具有可重构特征的轨道角动量天线技术研究进展. 电子与信息学报. 2024(04): 1173-1185 . 本站查看

其他类型引用(0)

资源附件(0)

访问统计

图(8) / 表(4)

计量

文章访问数: 2784
HTML全文浏览量: 1081
PDF下载量: 112
被引次数: 2

1. 引言
2. 走廊传播环境
3. 涡旋信道模型
3.1 传播模型
3.2 信道模型
4. 结果与分析
4.1 波前相位
4.2 信道幅度
5. 结束语

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

机载多通道SAR运动目标方位向速度和法向速度联合估计算法

doi: 10.11999/JEIT190672

通讯作者:
姜文　jiangwen13@mails.ucas.ac.cn

计量

Joint Estimation Algorithm for Azimuth Velocity and Normal Velocity of Moving Targets in Airborne Multi-channel SAR

1. 引言

2. 走廊传播环境

3. 涡旋信道模型

3.1 传播模型

3.2 信道模型

4. 结果与分析

4.1 波前相位

4.2 信道幅度

5. 结束语

期刊类型引用(2)

其他类型引用(0)

计量

目录

1. 引言

2. 走廊传播环境

3. 涡旋信道模型

3.1 传播模型

3.2 信道模型

4. 结果与分析

4.1 波前相位

4.2 信道幅度

5. 结束语

留言板

机载多通道SAR运动目标方位向速度和法向速度联合估计算法

doi: 10.11999/JEIT190672

通讯作者: 姜文 jiangwen13@mails.ucas.ac.cn

计量

出版历程

Joint Estimation Algorithm for Azimuth Velocity and Normal Velocity of Moving Targets in Airborne Multi-channel SAR

1. 引言

2. 走廊传播环境

3. 涡旋信道模型

3.1 传播模型

3.2 信道模型

4. 结果与分析

4.1 波前相位

4.2 信道幅度

5. 结束语

期刊类型引用(2)

其他类型引用(0)

计量

出版历程

目录

1. 引言

2. 走廊传播环境

3. 涡旋信道模型

3.1 传播模型

3.2 信道模型

4. 结果与分析

4.1 波前相位

4.2 信道幅度

5. 结束语

通讯作者:
姜文　jiangwen13@mails.ucas.ac.cn