基于压缩感知的外辐射源雷达分数阶杂波估计

杨鹏程; 吕晓德; 柴致海; 张丹; 岳琦; 杨璟茂

doi:10.11999/JEIT170079

基于压缩感知的外辐射源雷达分数阶杂波估计

doi: 10.11999/JEIT170079

计量
- 文章访问数: 1431
- HTML全文浏览量: 116
- PDF下载量: 198
- 被引次数: 2
出版历程
- 收稿日期: 2017-01-22
- 修回日期: 2017-04-14
- 刊出日期: 2017-11-19

Fractional-delay Clutter Estimation for Passive Radar Based on Compressive Sensing

摘要

摘要: 参考信号信杂比(SCR)是评估外辐射源雷达积累增益损失的重要参数。利用数字电视地面广播(DTTB)辐射源中的PN序列进行信杂比估计时，将出现接收信号相对于本地帧头的分数阶时延问题，导致信杂比估计出现严重偏差。针对该问题，该文利用信号在时延维的稀疏性，提出基于压缩感知的信杂比估计算法。仿真结果表明，该算法对不同强弱的信号都能得到精确的时延和强度估计，从而保证了信杂比估计的准确性。实测数据处理结果表明，实测数据的参考信号信杂比较高，积累损失较小，约为0.5 dB；信杂比与基线距离有关，基线越长信杂比越低，积累损失越大。
- 外辐射源雷达 /
- 压缩感知 /
- 参考信号 /
- 信杂比 /
- 分数阶时延 /
- PN序列
Abstract: Signal-Clutter-Ratio (SCR) in reference channel is an important parameter for the evaluation of integration loss of passive radar. When PN sequence in the Digital Terrestrial Television Broadcast (DTTB) illuminator is utilized for SCR estimation, there will be the problem of fractional-delay received signal relative to local PN sequence, which leads to a severe deviation of SCR estimation. For this problem, a novel algorithm based on compressive sensing is proposed by exploiting the sparsity of signal in delay dimension. Simulation demonstrates that accurate estimation of delay and strength for signals of different strength can be obtained by the proposed algorithm, which ensures the accuracy of SCR estimation. The processing of experimental data shows that the whole SCR in the received data is relatively high and the integration loss is trivial, about 0.5 dB. In addition, SCR decreases with the length of baseline, which leads to the increase of integration loss.
- Passive radar /
- Compressive Sensing (CS) /
- Reference signal /
- Signal-Clutter-Ratio (SCR) /
- Fractional-delay /
- PN sequence

HTML全文

1. 引言

电波在大气中传播时，大气折射效应会导致电波轨迹发生弯曲及电波传播速度不能按光速传播。该现象会引入大气折射误差，严重影响各个无线电系统探测定位性能，比如：卫星导航系统GNSS(Global Navigation Satellite System)定位精度^[1–4]、雷达探测系统对目标的跟踪精度以及通信系统时间同步精度^[5]。此外，大气折射误差导致的信号时延还可用于反演大气折射率剖面^[6]及监测电离层电子浓度(Total Electronic Contents, TEC)^[7]。因此，大气折射误差的研究对于无线电物理、地球物理、大地测量学、气象学、大气遥感等学科都具有非常重要的意义。为了提高各个无线电系统的测量精度，需要了解大气折射误差的修正方法。

现有关于折射误差的修正方法主要有：射线描迹法或射线追踪法^[8–12]、简化修正法^[13–17]和新方法^[18,19]。射线描迹法从射线方程出发，严格推导得到各个物理量的折射误差计算公式，精度较高。多篇文献对该方法的效果进行了验证^[4,20,21]，并提出了一些基于射线描迹法的快速算法^[22,23]。相比于射线描迹法，其它两种方法做了近似和简化处理，计算复杂度低。简化修正法主要分为两类：一类是经验模型，比如ITU-R P.834-8中的大气折射距离误差修正模型^[12]；另一类是基于天顶延迟(Zenith Tropospheric Delay, ZTD)的映射函数转换斜延迟法^[14–17]。这种简化方法较多地应用于高仰角条件下的大气折射误差修正，尤其是第2类简化方法。ZTD数据一方面可以通过全球GPS站和IGS提供，另一方面可以利用气象资料计算，为斜延迟的预测建立良好基础^[17]。新方法主要有微波辐射计法^[18]和针对电离层提出的实时修正方法^[19]，该方法的优点是实时性好，缺点是会引入其它测量系统误差。综上，虽然大气折射误差修正已经开展了大量的理论和实验研究，并取得了一些研究成果，但这些方法皆是针对标准大气环境提出的，无法直接应用于非标准大气环境。因此，异常大气(比如：大气湍流、对流层波导)环境的大气折射误差问题亟待解决。

为解决该问题，本文开展任意大气环境折射误差修正方法研究。首先，在第2节对大气折射误差涉及的折射率模型和传统射线描迹法进行了简单介绍，并基于几何光学原理给出大气折射误差修正方法中涉及到的物理量的计算公式；针对该方法的局限性，本文在第3节基于传播射线描迹法推导了射线描迹法微分形式，并基于此提出了适用于任意大气的折射误差修正方法；此外，利用测试数据对该方法的可行性进行了验证；最后在第4节以蒸发波导为例开展了异常大气折射误差的仿真分析。

2. 基于传统射线描迹法大气折射误差修正方法

电波在非自由空间中传播时，会产生大气折射效应。该效应是由于实际大气球面分层和水平非均匀特性造成。大气层分为对流层、平流层、电离层和磁层。其中对流层和平流层为60 km以下的大气层，被称为中性大气层。不同大气层折射率空间分布不同。考虑到大气折射率是影响折射误差的一个重要参数，关于折射率获取方法的研究从未中断。利用观测站长期积累的气象数据建立折射率模型是获取折射率的一种普遍使用方法。本文研究过程中即选取典型的折射率模型进行大气环境模拟。由于异常大气通常指对流层大气中异常折射率分布，因此本文主要探讨对流层大气中折射误差修正。

2.1 折射率模型

现有中性大气层的折射率模型主要有：线性模型、指数模型、ITU-R指数模型^[13]、双指数模型、分段模型和Hopfield模型^[4]。其中使用较广的是分段模型和Hopfield模型。分段模型可以更全面地描述大气层折射率 $N$ 特性，它与高度的关系可以表示为

$N\left( h \right) = \left\{ {\begin{aligned} & {{N_s} + \left( {h - {h_s}} \right)\Delta N,}\\ &\quad\quad\quad\quad\quad\quad\, {{{h_s} \le h \le {h_s} + \;{\rm{1\ km}}} }\\& {{N_{\rm{1}}}{e^{ - c\left( {h - {h_s} - {\rm{1}}} \right)}},}\\ &\quad\quad\quad\quad\quad\quad\, {{h_s} + {\rm{1}} \le h \le {\rm{9 \;km}}} \\& {{N_{\rm{9}}}{e^{ - {c_{\rm{9}}}\left( {h - {\rm{9}}} \right)}},}\;\;\;\;{{\rm{ 9\;km}} \le h \le {\rm{60\;km}}} \end{aligned}} \right.$

(1)

其中， ${N_s}$ 为地表面折射率， ${h_s}$ 为地表面海拔高度， $\Delta N$ 为1 km以内折射率梯度， ${N_{\rm{1}}}$ 和 ${N_{\rm{9}}}$ 分别为海拔1 km和9 km处的修正折射率， ${c_{\rm{1}}}$ 为地面以上1 km至海拔9 km的指数衰减率， ${c_{\rm{9}}}$ 为海拔9 km至60 km的指数衰减率。

再由大气折射率和折射指数 $n\left( h \right)$ 的关系，可计算得到大气折射指数：

$N\left( h \right) = \left[ {n\left( h \right) - {\rm{1}}} \right] \times {\rm{1}}{{\rm{0}}^{\rm{6}}}$

(2)

Hopfield模型将折射率分为湿项和干项部分^[4]，研究过程中利用该模型的优点是可以分别考虑它们对电波传播的影响。该模型和分段模型都能较好地描述实际大气折射率特性，本文在研究过程中选用分段模型。

2.2 基于传统射线描迹法大气折射误差修正方法

大气折射误差是指观测站测得的目标视在物理量和真实物理量的差值，其中包括：距离误差、高度误差和角度误差。图1给出了电波传播轨迹示意图，A和B分别为观测站和目标位置，C为地心。A到B沿着射线轨迹 $s$ 的视在距离 ${R_a}$ (单位为km)为

${R_a} = {\rm{c}}\int_{{{{s}}_{}}}^{} {\frac{{ds}}{{{{\rm{c}} / n}}}} = \int_{{h_{\rm{A}}}}^{{h_{\rm{B}}}} {n\csc \theta {\rm{d}}h}$

(3)

其中， ${\rm{c}}$ 为光速， $\theta$ (单位为 ° )为传播路径高度 $h$ 处的本地仰角， $n$ 为电波传播路径上任意高度 $h$ 处的折射指数， ${h_{\rm{A}}}$ (单位为km)和 ${h_{\rm{B}}}$ (单位为km)分别为观测站和目标的海拔高度。

图 1 射线传播轨迹示意图

下载: 全尺寸图片幻灯片

球坐标系下的射线描迹法为

$n\left( {{r_{\rm{0}}} + h} \right)\cos \theta = {n_{\rm{A}}}\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)\cos {\theta _{\rm{A}}}$

(4)

其中， ${n_{\rm{A}}}$ 为高度 ${h_{\rm{A}}}$ 处的折射指数， ${r_{\rm{0}}} = {\rm{6370 \;km}}$ 为地球半径， ${\theta _{\rm{A}}}$ (单位为 ° )为视在仰角。

利用式(4)，式(3)可重新写为

${R_a} =\!\!\! \int_{{h_{\rm{A}}}}^{{h_{\rm{B}}}} \!\!\!\!\!{\frac{{{n^2}({r_0} + h)}}{{\sqrt {{n^2}{{({r_0} + h)}^2} - n_{\rm{A}}^2{{({r_0} + {h_{\rm{A}}})}^2}{{\cos }^2}{\theta _{\rm{A}}}} }}} {\rm{d}}h$

(5)

由于目标的真实位置往往是未知的，可通过视在距离反求目标的真实高度 ${h_{\rm{B}}}$ 。

A到B的真实距离 $R$ (单位为km)为

$R = \sqrt {{{\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)}^{\rm{2}}} + {{\left( {{r_{\rm{0}}} + {h_{\rm{B}}}} \right)}^{\rm{2}}} - {\rm{2}}\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)\left( {{r_{\rm{0}}} + {h_{\rm{B}}}} \right)\cos {\phi _{\rm{T}}}}$

(6)

${\phi _{\rm{T}}}$ 为观测站和目标的地心角。

根据射线描迹法和图1中射线传播轨迹， ${\phi _{\rm{T}}}$ 为

${\phi _{\rm{T}}} = \int_{{h_{\rm{A}}}}^{{h_{\rm{B}}}} {\frac{{\cot \theta }}{{{r_{\rm{0}}} + h}}} {\rm d}h = {n_{\rm{A}}}\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)\cos {\theta _{\rm{A}}}\int_{{h_{\rm{A}}}}^{{h_{\rm{B}}}} {\frac{{\operatorname{d}\! h}}{{\left( {{r_{\rm{0}}} + h} \right)\sqrt {{n^{\rm{2}}}{{\left( {{r_{\rm{0}}} + h} \right)}^{\rm{2}}} - n_{\rm{A}}^{\rm{2}}{{\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)}^{\rm{2}}}{{\cos }^{\rm{2}}}{\theta _{\rm{A}}}} }}}$

(7)

距离误差 $\Delta R$ 为

$\Delta R{\rm{ = }}{R_a} - R = {R_a} - \sqrt {{{\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)}^{\rm{2}}} + {{\left( {{r_{\rm{0}}} + {h_{\rm{B}}}} \right)}^{\rm{2}}} - {\rm{2}}\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)\left( {{r_{\rm{0}}} + {h_{\rm{B}}}} \right)\cos {\phi _{\rm{T}}}}$

(8)

A到B的高度误差 $\Delta h$ (单位为km)为

$\Delta h{\rm{ = }}{h_a} - {h_{\rm{B}}} = \sqrt {{{\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)}^{\rm{2}}} + R_{{a}}^{\rm{2}} + {\rm{2}}\left( {{r_{\rm{0}}} + {h_{\rm{A}}}} \right)R_{{a}}^{}\sin {\theta _{\rm{A}}}} - {r_{\rm{0}}} - {h_{\rm{B}}}$

(9)

A到B的仰角误差 $\Delta \theta$ (单位为 ° )为

$\Delta \theta {\rm{ = }}{\theta _a} - {\theta _{\rm{T}}} = {\theta _a} - \arctan \left( {\cot {\phi _{\rm{T}}} - \frac{{{r_{\rm{0}}} + {h_{\rm{A}}}}}{{{r_{\rm{0}}} + {h_{\rm{B}}}}}\csc {\phi _{\rm{T}}}} \right)$

(10)

其中， ${\theta _{\rm{T}}}$ 为真实仰角，它可根据图1中射线几何示意图得到。

以上即为基于传统射线描迹法的大气折射误差修正方法，该方法利用几何光学原理，有严格的数学推导，不涉及任何近似，有多篇文献验证了该方法的有效性。但该方法仅适用于标准大气环境下折射误差修正，对于异常大气，比如对流层波导，传播路径上射线的高度并不总是随射线轨迹长度增大而增大^[24]，视在距离及地心仰角不能简单地通过积分求得，使得式(5)和式(7)失效。因此，针对该问题，在下一节推导射线描迹法微分形式下大气折射误差，并提出利用迭代的方法计算大气折射误差。

3. 基于射线描迹法微分形式的大气折射误差修正方法

3.1 基于传统射线描迹法大气折射误差修正方法

基于第2节中传统射线描迹法(式(4))，推导射线描迹微分形式的大气折射误差修正方法。

式(4)两边分别对高度 $h$ 求微分可得

$\begin{align}\frac{{{\rm{d}}n}}{{{\rm{d}}h}}&\left( {{r_{\rm{0}}} + h} \right)\cos \theta + n\cos \theta - n\left( {{r_{\rm{0}}} + h} \right)\sin \theta \frac{{{\rm{d}}\theta }}{{{\rm{d}}h}} \\&= {\rm{0}}\end{align}$

(11)

则

$\frac{{{\rm{d}}\theta }}{{{\rm{d}}h}}=\frac{{{\rm{d}}n}}{{{\rm{d}}h}}\frac{{\cot \theta }}{n} + \frac{{\cot \theta }}{{{r_{\rm{0}}} + h}}$

(12)

式(12)即为射线描迹法的微分形式。

基于式(12)和图1中射线的几何关系：

$\frac{{{\rm{d}}h}}{{{\rm{d}}s}} = \sin \theta \hspace{10pt}$

(13)

$\frac{{{\rm{d}}\phi }}{{{\rm{d}}s}} = \frac{{\cos \theta }}{{{r_{\rm{0}}} + h}}$

(14)

大气折射误差修正可通过以下步骤步进迭代求解。

步骤1 初始化各物理量。初始位置处视在距离 ${R_{a,{\rm{0}}}} = {\rm{0}}$ ；本地仰角 ${\theta _{\rm{0}}} = {\theta _{\rm{A}}}$ ；地心角 ${\phi _{\rm{0}}} = {\rm{0}}$ ；高度 ${h_{\rm{0}}} = {h_{\rm{A}}}$ 。

步骤2 计算第 $k$ ( $k = {\rm{1,2,}} ·\!· · ,K$ )步各物理量的值。固定射线轨迹长度步长 $\Delta s$ ，则由式(3)可得

${R_{a,k}} = {R_{a,k - {\rm{1}}}} + \Delta {R_{a,k}} = {R_{a,k - {\rm{1}}}} + {n_{k - {\rm{1}}}} \Delta s$

(15)

由式(13)可得

${h_k} = {h_{k - {\rm{1}}}} + \Delta {h_k} = {h_{k - {\rm{1}}}} + \Delta s \sin {\theta _{k - {\rm{1}}}}$

(16)

由式(14)可得

${\phi _k} = {\phi _{k - {\rm{1}}}} + \Delta {\phi _k} = {\phi _{k - {\rm{1}}}} + \frac{{\cos {\theta _{k - {\rm{1}}}}}}{{{h_{k - {\rm{1}}}} + {r_{\rm{0}}}}} \Delta s$

(17)

由式(12)—式(14)可得

$\begin{align}{\theta _k}& = {\theta _{k - {\rm{1}}}} + \Delta {\theta _k} \\&= {\theta _{k - {\rm{1}}}} + \left[ {{\rm{1}} + \frac{{{h_{k - {\rm{1}}}} + {r_{\rm{0}}}}}{{{n_{k - {\rm{1}}}}}} {{\left( {\frac{{{\rm d}n}}{{{\rm d}h}}} \right)}_{{h_{k - {\rm{1}}}}}}} \right] \Delta {\phi _k}\end{align}$

(18)

步骤3 设置迭代结束的目标函数：

$K = \mathop {\arg \min }\limits_k {\rm{ }}\varphi \left\{ {{R_{a,k}},{R_a}} \right\}$

(19)

利用2阶范数作为该问题的目标函数：

$\varphi = {\Bigl\| {{R_{a,k}} - {R_a}} \Bigl\|^2}$

(20)

步骤4 调整步骤2中 $k$ 的值，利用步骤2计算 ${R_{a,k}}$ ，直至目标函数式(20)达到最优。

步骤5 计算大气折射误差。当目标函数式(20)达到最优时，迭代步数 $k = K$ ，此时目标真实高度 ${h_{\rm{B}}} = {h_{K}}$ ，真实仰角 ${\theta _{\rm{T}}} = {h_{K}}$ ，地心仰角 ${\phi _{\rm{T}}} = {\phi _{K}}$ 。由此，可计算得到大气折射误差，并对其进行修正。

以上步骤即为本文基于射线描迹法微分形式提出的大气折射误差修正方法。相对第2节的传统射线描迹法，该方法不仅适用于标准大气环境的折射误差修正，也适用于异常大气环境的折射误差修正，在下一节对该方法进行验证。

3.2 试验验证

利用文献[21]中的测试数据，对本文所提大气折射误差修正方法进行验证。该数据来源于地面测站对某高空目标的跟踪，角度范围是3°～26°，通过对测试数据和精确轨道数据的处理，可得到每个仰角条件下测站与目标的真实距离、真实仰角及大气折射引入的误差值。信号传播的大气环境对中性大气层采用折射率分段模型(式(1))进行拟合，模型参数为 ${N_s} = {\rm{333}}{\rm{.23}}$ , $\Delta N = - 60.17$ km, $c = 0.1227$ km, ${c_9} = 0.1432$ km，得到60 km以下的大气折射率剖面如图2(a)所示；电离层采用测站位置处实测的电子浓度反演大气折射率，得到的折射率剖面如图2(b)所示。

图 2 大气折射率剖面

下载: 全尺寸图片幻灯片

基于以上的折射率剖面，分别利用第2.2节(式(3)—式(10))和第3.1节中的大气折射误差修正方法详细步骤计算距离和仰角误差，计算结果如图3所示。利用图3的计算结果对测试数据中折射误差(图4中实线)进行修正，结果如图4所示。图3(a)和图3(b)分别给出了距离和仰角的误差，两种方法的计算结果具有较好的一致性。图4(a)和图4(b)分别给出了距离和仰角误差修正后的结果，距离误差修正后的残差在10 m左右，仰角误差修正后的残差在0° 左右。分析以上结果可以得出：基于射线描迹法微分形式的大气折射误差修正方法与基于传统射线描迹法的大气折射误差修正方法具有较好的一致性。

图 3 大气折射误差分布图

下载: 全尺寸图片幻灯片

图 4 大气折射误差修正结果

下载: 全尺寸图片幻灯片

此外，针对传统射线描迹法无法修正异常大气折射误差的问题，本文利用第3.1节的方法对蒸发波导条件下大气折射误差进行仿真分析。

4. 蒸发波导环境大气折射误差仿真分析

蒸发波导是海上对流层波导中出现概率最高的波导类型，本文以其为例探讨海上波导环境对大气折射误差的影响。

4.1 蒸发波导

Paulus-Jeske提出的单参数蒸发波导模型^[25]，简称为PJ模型，该模型应用最为广泛，其修正折射率 $M$ 对应的解析式为

$M(h) = {M_0} + 0.13\left( {h - {h_d} \ln \left(\frac{{h + {z_0}}}{{{z_0}}}\right)} \right)$

(21)

其中， ${M_0}$ 为地表面折射率， ${h_d}$ 为蒸发波导高度， ${z_0} = 0.00015$ m为粗糙度因子。

利用式(22)可以将大气修正折射率转化为折射率^[26]：

$N(h) = M(h) - 0.157h$

(22)

4.2 仿真分析

基于蒸发波导环境折射率模型(式(21))，且考虑到全球统计资料给出的蒸发波导高度变化范围0～40 m，本文分别选取两个波导高度 ${h_d} = 20$ m和 ${h_d} = 30$ m，探讨蒸发波导环境下的大气折射误差，并将其与标准大气环境下的仿真结果作对比。标准大气的折射率模型参见式(1)，其中各个参数取文献[8]给出的全国平均值： $\Delta N = - 39.4$ km, $c = 0.1258$ km, ${c_9} = 0.1434$ km。仿真中雷达天线高度设置为 ${h_{\rm{A}}} = 15$ m，位于蒸发波导内部。由于波导环境对陷获角度的限制，应用于该环境电波传播特性研究的天线波束宽度都比较小。为确保电波射线仰角在1°以内，本文研究过程选取仰角分别为–0.15°, 0°和0.15° 3种情况。

根据以上条件，利用第3.1节中式(15)–式(18)，绘制固定视在距离100 km时，蒸发波导和标准环境条件下，在直角坐标系中的电波传播轨迹图，结果如图5所示。由该图可以直观看出，大气环境和仰角对电波传播轨迹的影响是非常显著的。相比标准大气环境，蒸发波导环境中的电波传播轨迹更加多样化，并不简单随高度增大电波轨迹背离地球表面传播；当仰角为0°时(虚线)，电波被陷获在波导内部向前传播，且波导高度越高，陷获越强。这种电波传播特性的不同必然会导致两种大气环境下大气折射误差结果的差异。因此，为进一步定量分析大气环境对大气折射误差的影响，利用第3.1节中的方法步骤及第2节中折射误差的定义，图6—图8分别仿真了固定视在距离，两种环境、3个仰角条件下的距离误差、高度误差和仰角误差。这3个图中的横坐标均表示视在距离。此外，图9对视在距离100 km时的结果进行了统计。

图 5 电波传播轨迹

下载: 全尺寸图片幻灯片

图 6 距离误差

下载: 全尺寸图片幻灯片

图 7 高度误差

下载: 全尺寸图片幻灯片

由图6可以得出：两种大气环境均表现出随视在距离增大，距离误差增大的规律；相同条件下，标准大气引入的距离误差大于蒸发波导(图9(a))，这是由于电波受波导影响，射线在波导内变化比较平缓导致的(如图5所示)。此外，蒸发波导环境中，电波被陷获越强，距离误差越大。

相比距离误差，大气环境和仰角变化对高度和仰角误差的影响更为显著和复杂。分析图7的仿真结果：当仰角大于等于0°时，高度误差随视在距离的增大而增大，蒸发波导环境大气折射导致的高度误差远大于标准大气。比如图9(b)中对于蒸发波导环境，高度误差在0°仰角时近800 m，而标准大气为200 m左右。仰角小于0°时，在视在距离较小的情况下，视在高度小于目标真实高度，随视在距离增大，视在高度大于真实高度，相应的高度误差也反向增大，蒸发波导环境这种变化规律比较明显。

通过分析图8可以得出：当仰角大于等于0°时，仰角误差随视在距离的增大而增大，蒸发波导环境大气折射导致的仰角误差远大于标准大气。比如图9(c)中对于蒸发波导环境，仰角误差在0°时为0.45°，而标准大气为0.11°。仰角小于0°时，随视在距离的增大，受目标高度变化的影响，仰角误差的变化也呈现与高度误差相同的规律。

图 8 仰角误差

下载: 全尺寸图片幻灯片

图 9 视在距离100 km时大气折射误差统计图

下载: 全尺寸图片幻灯片

综上，由图6—图9可以得出：完全陷获在蒸发波导内部的电波引入的大气折射误差是最显著的；考虑到影响蒸发波导大气折射误差的因素多，相比标准大气，该环境下的大气折射误差修正问题更为严峻。

5. 结束语

为解决传统射线描迹法无法处理异常大气折射误差的问题，本文推导了射线描迹法的微分形式，并基于此提出了一种可以计算任意大气环境折射误差修正方法。该方法通过数值迭代和优化的过程，给出目标真实位置相关的物理量，进而对目标的视在位置进行修正。利用某地面测控站的试验数据对该方法进行了验证，结果显示：距离误差修正后的残差在10 m左右，仰角误差修正后的残差在0°左右。

此外，本文还以蒸发波导为例探讨了异常大气对大气折射误差的影响。通过仿真分析，得出：完全陷获在蒸发波导内部的电波引入的大气折射误差是最显著的；考虑到影响蒸发波导大气折射误差的因素多，相比标准大气，该环境下的大气折射误差修正问题更为严峻。

综上，本文所提方法能够较好地解决任意大气折射误差的修正问题，是现有大气折射误差修正方法的一个补充，为提高无线电系统在任意大气环境下的测量精度提供一种新的技术支持。

参考文献(24)

高志文, 陶然. 外辐射源雷达的相干积累增益计算及性能分析[J]. 电子学报, 2008, 36(6): 1227-1230. doi: 10.3321/j.issn: 0372-2112.2008.06.038.

GRIFFITHS H and BAKER C. Passive coherent location radar systems. Part 1: Performance prediction[J]. IEE Proceedings-Radar, Sonar and Navigation, 2005, 152(3): 153-159. doi: 10.1049/ip-rsn:20045082.

GAO Zhiwen and TAO Ran. Computation and analysis of coherent accumulation gain for passive radar[J]. Acta Electronica Sinica, 2008, 36(6): 1227-1230. doi: 10.3321/ j.issn:0372-2112.2008.06.038.

赵永科. 无源相参雷达参考信号影响及重构方法研究[D]. [硕士论文], 中国科学院大学, 2014: 37-39.

ZHAO Yongke. Research on effects and reconstructing algorithm of reference signal for passive coherent radar[D]. [Master dissertation], University of Chinese Academy of Sciences, 2014: 37-39.

冯远. 数字电视辐射源雷达参考信号获取及干扰抑制算法研究 [D]. [博士论文], 北京理工大学, 2014: 11-13.

FENG Yuan. Research on algorithms of reference signal acquisition and interference suppression in DTTB based passive radar[D]. [Ph.D. dissertation], Beijing Institute of Technology, 2014: 11-13.

BACZYK M K, KULPA K, SAMCZYNSKI P, et al. The impact of reference channel SNR on targets detection by passive radars using DVB-T signals[C]. 2015 IEEE Radar Conference, Arlington, USA, 2015: 708-712. doi: 10.1109/ RADAR.2015.7131088.

万显荣, 岑博, 易建新, 等. 中国移动多媒体广播外辐射源雷达参考信号获取方法研究[J]. 电子与信息学报, 2012, 34(2): 338-343. doi: 10.3724/SP.J.1146.2011.00572.

WAN Xianrong, CEN Bo, YI Jianxin, et al. Reference signal extraction methods for CMMB-based passive bistatic radar [J]. Journal of Electronics Information Technology, 2012, 34(2): 338-343. doi: 10.3724/SP.J.1146.2011.00572.

WANG F. Direct signal recovery and masking effect removal exploiting sparsity for passive bistatic radar[C]. IET International Radar Conference, Hangzhou, China, 2015: 1-6. doi: 10.1049/cp.2015.0969.

JIANG L, FENG T, ZHANG W, et al. The direct wave purifying based on WIFI signal for passive radar[J]. Wireless Communications, Networking and Applications, 2016, 348: 223-234. doi: 10.1007/978-81-322-2580-5_22.

朱家兵, 洪一. 基于复倒谱技术的无源雷达直达波提纯方法 [J]. 现代雷达, 2007, 29(8): 75-78. doi: 10.3969/j.issn.1004- 7859.2007.08.021.

ZHU Jiabing and HONG Yi. A direct-path-wave purified approach to passive radar based on cepstrum technique[J]. Modern Radar, 2007, 29(8): 75-78. doi: 10.3969/j.issn.1004- 7859.2007.08.021.

中国国家标准化管理委员会. 国标: 20600-2006. 数字电视地面广播传输系统帧结构, 信道编码和调制[S]. 2006.

SAC(Standardization Administration of the Peoples Republic of China). GB: 20600-2006. Framing Structure, Channel Coding and Modulation for Digital Television Terrestrial Broadcasting System[S]. 2006.

DONOHO D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289-1306. doi: 10.1109/ TIT.2006.871582.

王海涛, 王伟, 符渭波. 基于压缩感知的无源雷达目标检测和方向估计[J]. 空军预警学院学报, 2015, 29(4): 246-250. doi: 10.3969/j.issn.2095-5839.2015.04.004.

WANG Haitao, WANG Wei, and FU Weibo. Target detection and DOA estimation for based compressed sensing-based passive radar[J]. Journal of Air Force Early Warning Academy, 2015, 29(4): 246-250. doi: 10.3969/j.issn.2095-5839.2015. 04.004.

LU Xiaode, YANG Pengcheng, LI Daojing, et al. Interference cancellation based on compressive sensing for passive coherent radar (PACOR)[C]. 2015 IEEE Radar Conference, Arlington, USA, 2015. doi: 10.1109/RADAR.2015.7131055.

杨鹏程, 吕晓德, 李纪传, 等. 基于压缩感知的外辐射源雷达目标检测[J]. 雷达科学与技术, 2015, 13(4): 384-389. doi: 10.3969/j.issn.1672-2337.2015.04.009.

YANG Pengcheng, LU Xiaode, LI Jichuan, et al. Target detection for passive radar using compressive sensing[J]. Radar Science and Technology, 2015, 13(4): 384-389. doi: 10.3969/j.issn.1672-2337.2015.04.009.

CHEN S S, DONOHO D L, and SAUNDERS M A. Atomic decomposition by basis pursuit[J]. SIAM Journal on Scientific Computing, 1998, 20(1): 33-61.

TROPP J A and GILBERT A C. Signal recovery from random measurements via orthogonal matching pursuit[J]. IEEE Transactions on Information Theory, 2007, 53(12): 4655-4666.

施引文献

期刊类型引用(1)
1. 罗春兰，林文，祝晓霞. (47, 24, 11)QR码的快速代数译码算法. 三明学院学报. 2022(06): 60-66+78 . 百度学术
其他类型引用(1)

资源附件(0)

访问统计

计量

文章访问数: 1431
HTML全文浏览量: 116
PDF下载量: 198
被引次数: 2

1. 引言
2. 基于传统射线描迹法大气折射误差修正方法
2.1 折射率模型
2.2 基于传统射线描迹法大气折射误差修正方法
3. 基于射线描迹法微分形式的大气折射误差修正方法
3.1 基于传统射线描迹法大气折射误差修正方法
3.2 试验验证
4. 蒸发波导环境大气折射误差仿真分析
4.1 蒸发波导
4.2 仿真分析
5. 结束语

1. 引言
2. 基于传统射线描迹法大气折射误差修正方法
2.1 折射率模型
2.2 基于传统射线描迹法大气折射误差修正方法
3. 基于射线描迹法微分形式的大气折射误差修正方法
3.1 基于传统射线描迹法大气折射误差修正方法
3.2 试验验证
4. 蒸发波导环境大气折射误差仿真分析
4.1 蒸发波导
4.2 仿真分析
5. 结束语

参考文献(24)

施引文献

资源附件(0)

访问统计

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

留言板

基于压缩感知的外辐射源雷达分数阶杂波估计

doi: 10.11999/JEIT170079

计量

Fractional-delay Clutter Estimation for Passive Radar Based on Compressive Sensing

1. 引言

2. 基于传统射线描迹法大气折射误差修正方法

2.1 折射率模型

2.2 基于传统射线描迹法大气折射误差修正方法

3. 基于射线描迹法微分形式的大气折射误差修正方法

3.1 基于传统射线描迹法大气折射误差修正方法

3.2 试验验证

4. 蒸发波导环境大气折射误差仿真分析

4.1 蒸发波导

4.2 仿真分析

5. 结束语

期刊类型引用(1)

其他类型引用(1)

计量

目录

1. 引言

2. 基于传统射线描迹法大气折射误差修正方法

2.1 折射率模型

2.2 基于传统射线描迹法大气折射误差修正方法

3. 基于射线描迹法微分形式的大气折射误差修正方法

3.1 基于传统射线描迹法大气折射误差修正方法

3.2 试验验证

4. 蒸发波导环境大气折射误差仿真分析

4.1 蒸发波导

4.2 仿真分析

5. 结束语

留言板

基于压缩感知的外辐射源雷达分数阶杂波估计

doi: 10.11999/JEIT170079

计量

出版历程

Fractional-delay Clutter Estimation for Passive Radar Based on Compressive Sensing

1. 引言

2. 基于传统射线描迹法大气折射误差修正方法

2.1 折射率模型

2.2 基于传统射线描迹法大气折射误差修正方法

3. 基于射线描迹法微分形式的大气折射误差修正方法

3.1 基于传统射线描迹法大气折射误差修正方法

3.2 试验验证

4. 蒸发波导环境大气折射误差仿真分析

4.1 蒸发波导

4.2 仿真分析

5. 结束语

期刊类型引用(1)

其他类型引用(1)

计量

出版历程

目录

1. 引言

2. 基于传统射线描迹法大气折射误差修正方法

2.1 折射率模型

2.2 基于传统射线描迹法大气折射误差修正方法

3. 基于射线描迹法微分形式的大气折射误差修正方法

3.1 基于传统射线描迹法大气折射误差修正方法

3.2 试验验证

4. 蒸发波导环境大气折射误差仿真分析

4.1 蒸发波导

4.2 仿真分析

5. 结束语