作者简介:付永庆(1956-),男,教授.研究方向:通信与信息处理,阵列信号处理.E-mail:fuyongqing@hrbeu.edu.cn
为了实现电磁平面波在不同目标上的选择性聚焦,提出了一种基于波动方程时间反转不变性的时间反转算子分解算法。该算法利用天线阵列发射和接收的目标探测信号获得传输矩阵和时间反转算子,然后对该算子进行特征值分解,其中非零主特征值的个数对应目标个数,而每个非零主特征值所对应的主特征向量则包含了相应目标的方位信息。根据每个目标相对应的特征向量计算目标域中各观测点的目标函数可获得目标域中的目标函数值分布情况,并实现在不同目标处的选择性聚焦。最后,通过仿真实验验证了电磁平面波情况下基于时间反转算子分解的选择性聚焦方法的准确性和有效性。
In order to achieve the selective focusing on different targets in electromagnetic plane wave, a decomposition method of the time reversal operator based on the time reversal invariance of wave equation is proposed. In this method, the transfer matrix and time reversal operator are obtained by emitting and receiving the signal for detecting targets with the antenna array. Then, the time reversal operator is used to perform eigenvalue decomposition. In eigenvalue decomposition, the number of nonzero main eigenvalues is the same as the number of targets, and the main eigenvectors corresponding with these nonzero main eigenvalues contain the direction information of the targets. According to these eigenvectors of targets to calculate the objective function of every detection point in the given target area can be obtained the objective function value distribution of the target area, and achieve the selective focusing on different targets. Finally, the simulations are performed and the results show that the selective focusing method with decomposition of the time reversal operator in electromagnetic plane wave is effective.
时间反转理论由光学中的相位共轭产生, Fink等[1]于1989年验证了时间反转的空间聚焦特性并首次给出了时间反转镜(Time reversal mirror, TRM)的概念。随后, 国内外专家在水声领域做了大量的相关研究及实验并取得了一定的研究成果[2, 3]。到目前为止, 时间反转理论已经在超声碎石[4]、无损探伤[5]、水下目标探测[6]、无线通信[7, 8]和成像[9, 10]等领域得到广泛探讨与研究。时间反转技术最主要的优点是不需要介质性质和阵列分布等先验知识就可以实现自适应聚焦[11, 12], 因此越来越受到人们的重视。对于信号反射强度不同的多个目标, 为了实现在目标上的选择性聚焦, 可通过有限次数的时间反转迭代逐渐削弱反射系数较小的目标处的相对聚焦信号能量, 最终实现在强反射目标上的选择性聚焦[13], 但通过迭代时间反转使信号聚焦于弱反射目标处是不可能实现的。
为了实现在信号反射强度不同的多目标处的选择性聚焦, 本文提出了一种基于电磁平面波波动方程的时间反转算子分解算法, 仿真实验显示, 根据时间反转算子分解的特征值和特征向量可以有效地实现在反射系数不同的目标处的选择性聚焦。
时间反转理论以波动方程的时间反转不变性为基础, 把天线阵列接收到的目标信号进行时间反转后重新发射到空间中, 则时间反转后信号经过反向传输并实现在目标处的重聚焦。
假设电磁场中的标量位函数为φ (r, t), 则在无源场区域满足的标量波动方程为:
式中:μ 和ε 分别为介质的介电常数和磁导率; r为电磁波传播的位移矢量。
考虑平面电磁波, 则式(1)所示标量波动方程的通解为:
式中:f1和f2为任意函数; v=1/
将式(2)右边两项分别用φ 1
可以看出, 式(3)和式(4)分别表示传播方向相反的两类波形, 且它们都是波动方程的解。
对式(3)和式(4)做时间反转处理, 可得:
可见, φ 1
考虑一个由N个天线阵元组成的TRM阵列, 有N× N个交叉阵元脉冲响应, 如图1所示。令klm
式中:τ 为电磁波的传播时延。
从阵元m到阵元l之间的交叉脉冲响应可以表示为:
式中:τ m、τ l分别为目标到第m个和第l个阵元的传播时延; * 表示卷积。
令em
式(9)写成频域形式为:
用矩阵形式表示:
式中:E
根据互易定理, 从阵元m到阵元l之间的交叉脉冲响应与从阵元l到阵元m之间的交叉脉冲响应相同, 因此传输矩阵K
设E0
时域内的时间反转操作等同于频域内的相位共轭, 则天线阵列第二次发射的信号向量为第一次接收信号向量的相位共轭:
经过两次时间反转处理后, 天线阵列的发射信号向量可以表示为:
以此类推, 经过2n次和2n+1次时间反转处理后, 天线阵列的发射信号向量可以分别表示为:
式中:K*
因为K
式中:T表示对矩阵的转置运算。
式(17)说明K*
假设目标域中共有d个目标, 各目标对信号的反射率分别为C1, C2, …, Cd。第i(1≤ i≤ d)个目标与接收阵列第n(1≤ n≤ N)个阵元之间的传输函数为hin
如果点目标理想可分辨, 即在其中一个目标上的时间反转聚焦不会给其他目标提供能量, 则Hi
从式(20)可以看出传输矩阵由3个矩阵组成:描述从阵列到目标的前向传输矩阵, 目标反射系数组成的反射矩阵和描述从目标到阵列的反向传输矩阵。
令天线阵列的发射信号为
则S
由Hi
写成矩阵形式有:
进而可得到:
可见, 时间反转算子对应某目标的特征向量是该目标到天线阵列之间的传输向量的共轭, 那么该特征向量与目标区域中的观测点对应的传输函数的内积在其对应的目标处将趋于最大, 可实现在该目标处的信号能量聚焦。假设目标域中共有M个目标观测点, 于是可得目标域中第j(1≤ j≤ M)个观测点处时间反转算子分解的选择性聚焦目标函数公式为:
式中:Hjn
根据以上分析, 可将基于时间反转算子分解的选择性聚焦算法分为以下几个步骤:
(1)利用天线阵列的第一个阵元发射中心频率为ω 0的目标探测信号, 经目标反射后被天线阵列接收。对接收信号进行傅里叶变换可得第一个阵元与所有阵元之间的交叉脉冲响应。
(2)用同样的信号对天线阵列的其他N-1个阵元进行激励, 重复步骤(1)的操作可获得传输矩阵K
(3)计算时间反转算子K*
(4)确定目标搜索域及各观测点。
(5)根据式(26)并利用目标所对应的主特征向量计算目标域中各观测点的目标函数值, 可取得在该目标处的选择性聚焦。
为了验证基于时间反转算子分解的选择性聚焦特性, 建立如图2所示的均匀圆形天线阵列模型, 其中N个阵元均匀分布在以r为半径的圆周上, 且阵列的几何中心与坐标原点重合。仿真实验中取N=32, r为20 m, 天线阵列发射中心频率为30 MHz的目标探测波。考虑以坐标原点为中心, 半径为20 km的圆周区域为待侦测目标区域, 且在整个圆周域上平均分布着360个观测点, 侦测步长为1° , 目标位于待侦测区域内, 且方位角分别为22° , 45° , 70° 。
为实现时间反转算子分解的选择性聚焦, 首先通过时间反转阵列的发射和接收信号获得传输矩阵, 并计算得到时间反转算子, 然后对时间反转算子进行特征值分解, 其中主特征值的个数为目标个数, 而主特征值对应的主特征向量则包含了对应目标的方位信息。利用主特征向量并根据式(26)计算各观测点的目标函数值, 其中函数值最大点所对应的观测点的方位即为该主特征值所对应目标的方位, 因此利用每个主特征值进行目标域中各观测点的目标函数值计算, 可分别实现在不同目标处的选择性聚焦。
首先取C1∶ C2∶ C3=1∶ 0.5∶ 0.3, 通过仿真计算获得阵元间传输矩阵K
为分析基于时间反转算子分解的选择性聚焦算法在不同信噪比下的聚焦性能, 与图3实验相同, 取C1∶ C2∶ C3=1∶ 0.5∶ 0.3, 3个主特征向量V1、V2、V3分别对应目标1、2、3, 并在仿真时加入信噪比不同的高斯白噪声。利用不同主特征向量实现其对应目标处的选择性聚焦, 并执行100次蒙特卡洛实验, 图4为不同信噪比情况下基于时间反转算子分解的选择性聚焦算法的均方根误差。从图中可以看出:随着信噪比的增加, 该算法选择性聚焦的均方根误差逐渐变小, 而且在信噪比为0 dB的情况下, 不同目标的选择性聚焦误差都在1° 以内, 显示了该算法良好的聚焦性能。
利用与图3同样的方法可以获得C1∶ C2∶ C3=0.5∶ 1∶ 0.3和C1∶ C2∶ C3=0.3∶ 0.5∶ 1条件下的仿真结果, 如图5和图6所示。从图中可以看出:改变目标的反射系数, 时间反转算子分解算法依然可以很好地实现在各不同目标处的选择性聚焦, 只是特征向量V1、V2、V3所对应的目标位置发生了改变。综合图3、图5和图6可以看出:非零特征值的大小与对应目标的信号反射率有关, 目标反射率越大, 对应的特征值也越大, 因此可以根据特征值的大小选择相应的特征向量实现在不同目标处的聚焦。
提出了一种基于平面电磁波波动方程时间反转不变性的时间反转算子分解选择性聚焦方法, 该方法通过对时间反转算子的分解可以获得与目标数目一致的主特征值和与目标方位一一对应的主特征向量。根据每个目标所对应的特征向量计算目标域中所有观测点处的目标函数值可以实现在各目标处的选择性聚焦。仿真实验结果显示该方法具有有效性, 利用弱反射目标所对应的特征向量可实现在弱目标处的聚焦, 利用所有的主特征向量可以实现在不同目标上的选择性聚焦, 另外, 主特征值的大小和各目标所对应的特征向量次序都与目标的反射系数有关。该方法可以用于对多个目标中的弱目标进行检测或者是对多目标进行选择性检测的情况, 具有很好的应用前景。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|