基于微波光子阵列的超宽带信号波达方向角估计方法与流程

本发明属于光电通信技术领域，涉及一种超宽带信号波达方向角估计方法,特别涉及一种基于微波光子阵列的超宽带信号波达方向角估计方法,可用于目标检测与无源定位。

背景技术：

超宽带信号波达方向角估计是利用空间中的天线阵列接收多个不同方向的信号源发出的信号，运用现代信号处理方法快速准确的获得信号源方向的技术，在雷达、声纳、无线通信等领域具有重要应用价值。针对该问题构建的模型中，出现较早且应用较为广泛的是基于子空间的模型，之前的大部分超宽带信号波达方向角估计都是利用该模型生成的。

目前，最新超宽带信号波达方向角估计方法是用微波光子估计超宽带信号波达方向角，这种方法旨在采用光电器件和光电学方法估计超宽带信号波达方向角，与传统超宽带信号波达方向角估计方法相比，具有估计带宽大、损耗小、抗干扰能力强、系统小型等优点，能够克服电子瓶颈，适应复杂的电磁环境的优点。例如申请公布号为cn107528638a,名称为“基于微波光子滤波的宽带微波信号到达角估计方法”的专利申请中，公开了一种基于微波光子滤波的宽带微波信号到达角估计方法，该方法利用偏振复用马赫增德尔调制器，结合差分延时模块，构建双抽头微波光子滤波器，对输入电信号进行陷波滤波并观察结果，但是由于该方法是基于两阵元实现的，因此一次只能估计一个角度，且估计效率和估计精度较低。

技术实现要素：

本发明的目的在于克服上述现有技术存在的不足，提出了一种基于微波光子阵列的超宽带信号波达方向角估计方法，用于解决现有技术中存在的超宽带信号波达方向角估计效率和估计精度较低的技术问题。

为克服上述现有技术的不足，本发明采取的技术方案包括如下步骤：

(1)获取一阶边带信号e0(t)：

将天线0接收到的微波信号v0(t)和激光信号源输出的光载波信号e0输入到马赫增德尔调制器mzm-0中进行调制，得到多边带信号，并对多边带信号进行窄带滤波，得到一阶边带信号e0(t)：

(2)构建微波光子阵列d：

构建由m个并联的天线单元组成的天线阵列，相邻天线单元之间的距离为d，m＝1…m，m≥2，第m个天线单元上级联有马赫增德尔调制器mzm-m、滤波器f-m和光功率计opm-m，得到微波光子阵元dm，m个微波光子阵元组成微波光子阵列d；

(3)获取微波光子阵列d输出的超宽带信号em的光功率pm：

(3a)将第m个天线单元接收到的微波信号vm(t)和mzm-0输出的一阶边带信号e0(t)输入到马赫增德尔调制器mzm-m中进行调制，并通过滤波器f-m对调制得到的第m个调制信号进行窄带滤波，得到第m个超宽带信号em：

其中φk为第m个天线单元接收到的微波信号vm(t)的相位差，φkm表示m个天线单元接收到的微波信号的相位差之和，lm为马赫增德尔调制器mzm-m的插入损耗；

(3b)将第m个超宽带信号em输入到光功率计opm-m中，计算超宽带信号em的光功率pm：

(4)定义光功率比值向量q：

(4a)计算pm的相对光功率比值qm：

(4b)定义m×1维的光功率比值向量：q＝[q1,q2,…,qm…,qm]^t，其中[·]^t表示矩阵[·]的转置；

(5)定义m×nθ维矩阵

(5a)根据微波信号源的空域稀疏特性，采用空间网格划分方法，将可搜索的空间角度域等间隔划分成nθ份,得到超宽带微波信号em的相位的取值空间为：

(5b)定义m×nθ维的矩阵

(6)获取微波光子阵列d输出的超宽带信号em的相位

(6a)定义一个nθ×1维的稀疏表示系数向量：该向量为未知向量；

(6b)通过改进的l1-sracv算法模型，构建基于稀疏重构的约束优化方程式：

其中||·||1和||·||2分别表示l1范数和l2范数，s.t.表示约束关系，ε表示误差的允许值；

(6c)估计稀疏重构的约束优化方程中的最优解，并用的最优解求解超宽带信号em的相位

(7)计算微波光子阵列d输出的超宽带信号em的波达方向角θm：

(7a)计算超宽带信号em的相位与相邻天线单元接收到的微波信号vm(t)的时差δt：

其中，ωm为第m个天线单元接收到的微波信号vm(t)的角频率，ρ是一个整数值，可根据天线阵列的物理参数获取，c表示光速，l为两个相邻的马赫增德尔调制器之间的光纤长度；

(7b)计算超宽带信号em的波达方向角θm：

其中，d表示相邻天线单元之间的距离。

本发明与现有技术相比，具有以下优点：

1)本发明在超宽带信号的获取中采用微波光子阵列，无需事先估计入射信号的个数，同时可用于处理相干信号源和非相干信号源，避免了传统超宽带信号的波达角度估计中角度搜索或角度匹配过程，降低了运算量，与现有技术相比，提高了超宽带信号波达方向角的估计效率。

2)本发明在超宽带信号波达方向角的估计中采用稀疏表示的思想将波达方向角度估计问题转化成稀疏重构问题，利用信号源在空域稀疏的特点进行建模，突破了微波光子的瑞丽限，与现有技术相比，提高了超宽带信号波达方向角的估计精度。

附图说明

图1是实现本发明采用的估计系统的结构示意图；

图2是本发明的实现流程图。

具体实施方式

下面结合附图和具体实施例，对本发明作进一步详细描述：

参照图1、实现本发明采用的估计系统，包括微波光子阵列，一阶边带信号生成器。在空间等间隔放置m个天线，m＝1…m，m≥2，分别定义为天线1，天线2，…，天线m，…，天线m，相邻天线之间的距离为d，m个并联的天线组成的微波光子阵列，第m个天线上级联有马赫增德尔调制器mzm-m、滤波器f-m和光功率计opm-m。一阶边带信号生成器是将激光信号源天线0、和一个滤波器f-0并联到马赫增德尔调制器mzm-0上。一阶边带信号生成器中的滤波器f-0与微波光子阵列d中每一个天线相连接，微波微波光子阵列d中的每一个光功率计分别级联1个光功率比值器，m个光功率比值器并联到稀疏表示模型中。

参照图2、一种基于微波光子阵列的超宽带信号波达方向角估计方法，包括如下步骤：

步骤1)获取一阶边带光电信号e0(t)：

将天线0接收到的微波信号v0(t)和激光信号源输出的光载波信号e0输入到马赫增德尔调制器mzm-0中进行调制成为光电信号，得到多边带光电信号，并对多边带光电信号进行窄带滤波，滤除光载频以外的信号频率分量，得到一阶边带光电信号e0(t)：

其中，ω表示光载波信号e0的载波角频率，p表示光载波信号e0的载波功率，k为调制系数，l0为mzm-0的插入损耗,ω0为微波信号v0(t)的角频率，j1(·)表示一阶贝赛尔函数。

马赫增德尔调制器是将输入的光载波信号和微波信号调制成为光电信号的光电器件，将光信号和电信号调制成为光电信号，即微波光子。

步骤2)构建微波光子阵列d：

在空间放置m个天线单元，m＝1…m，m≥2，分别定义为天线单元1，天线单元2，…，天线单元m，天线单元m，这m个天线单元等间隔放置，相邻天线单元之间的距离为d，构建由m个并联的天线单元组成的天线阵列，第m个天线单元上级联有马赫增德尔调制器mzm-m、滤波器f-m和光功率计opm-m，得到微波光子阵元dm，这m个微波光子阵元组成微波光子阵列d，将一阶边带光电信号e0(t)输入到这m个微波光子阵元组成微波光子阵列d中，可同时得到m个超宽带信号，对m个超宽带信号进行波达方向角的估计，这样可以避免传统超宽带信号的波达角度估计中搜索或角度匹配问题，降低了运算量，提高超宽带信号波达方向角的估计效率，具体实施例中，m＝50。

步骤3)获取微波光子阵列d输出的超宽带信号em的光功率pm：

步骤3a)将第m个天线单元接收到的微波信号vm(t)和mzm-0输出的一阶边带光电信号e0(t)输入到马赫增德尔调制器mzm-m中进行调制，并通过滤波器f-m对调制得到的第m个调制信号进行窄带滤波，得到第m个超宽带信号em：

其中φk为第m个天线单元接收到的微波信号vm(t)的相位差，φkm表示m个天线单元接收到的微波信号的相位差之和，lm为马赫增德尔调制器mzm-m的插入损耗；

步骤3b)将第m个超宽带信号em输入到光功率计opm-m中，计算超宽带信号em的光功率pm：

步骤4)定义光功率比值向量q：

步骤4a)计算pm的相对光功率比值qm：

(4b)定义m×1维的光功率比值向量：q＝[q1，q2,…,qm…,qm]^t，其中[·]^t表示矩阵[·]的转置；

步骤5)定义与阵列流型向量相关的m×nθ维矩阵

步骤5a)根据微波信号源的空域稀疏特性，采用空间网格划分方法，将可搜索的空间角度域[-90°,90°]等间隔划分成nθ份，得到超宽带微波信号em的相位的取值空间为：nθ＞＞m，网格划分间隔的取值nθ根据期望达到的角度估计精度进行设定，网格划分间隔越小，则最终得到的角度估计值精度越高，具体实施例中，nθ＝100；

步骤5b)定义m×nθ维的矩阵

步骤6)获取微波光子阵列d输出的超宽带信号em的相位

步骤6a)定义一个nθ×1维的稀疏表示系数向量：该向量为未知向量；

步骤6b)信号源在空域有稀疏特性，应用于超宽带信号的波达方向角可以突破微波光子阵列分辨率的瑞丽限，进一步提高超宽带信号波达方向角精度。通过改进的l1-sracv算法模型，可得稀疏重构约束化方程为：

其中||·||1和||·||2分别表示l1范数和l2范数，s.t.表示约束关系，ε表示误差的允许值；

改进的l1-sracv算法模型是利用阵列接收数据奇异值分解(svd)得到的信号子空间构造稀疏表示模型，重构稀疏约束化方程。

步骤6c)估计稀疏重构的约束优化方程中的最优解，并用的最优解求解超宽带信号em的相位

采用凸优化方法求解稀疏重构约束化方程，凸优化是一种比较特殊的优化，是指目标函数和约束函数均为凸函数的优化问题，凸优化问题有一套非常完备的解决算法，在此采用现有针对凸优化问题的软件包cvx(grantm,boyds.cvx:matlabsoftwarefordisciplinedconvexprogramming[j].2008[online]available:http://stanfordedu/～boyd/cvx)来求解，通过该方法能够快速地得到空域稀疏向量的最优解。

以超宽带信号em的相位取值范围的值为x轴坐标，以稀疏表示系数向量的最优解的幅度值为y轴坐标，绘制幅度谱图，在该幅度谱图中按照从高到低的顺序寻找幅值较大的前m个谱峰，第m谱峰的峰值点所对应的x轴坐标即为所求第m个超宽带信号em的相位

步骤7)计算微波光子阵列d输出的超宽带信号em的波达方向角θm：

步骤7a)计算超宽带信号em的相位与相邻天线单元接收到的微波信号vm(t)的时差δt：

其中，ωm为第m个天线单元接收到的微波信号vm(t)的角频率，ρ是一个整数值，可根据天线阵列的物理参数获取，c表示光速，l为两个相邻的天线单元之间的光纤长度，具体实施例中，ρ＝1，l＝1米；

步骤7b)计算超宽带信号em的波达方向角θm：

其中，d表示相邻天线单元之间的距离，在具体实施例中d＝1米。