一种基于双通道接收机的无线信号到达角估计方法及系统与流程

1.本发明涉及无线测向领域，具体涉及一种基于双通道接收机的无线信号到达角估计方法及系统。


背景技术：

2.现有的天线阵列信号处理算法中，最具代表性的是music算法及一系列衍生算法。该算法通过计算采样周期内接收数据的平均相关矩阵，基于数据的统计特征，对无线信号的到达角做出估计。但是该算法具有以下缺点：1)其估计结果依赖于数据的统计特性，需要在信道参数稳定的状态下进行高频次的采样，对设备成本要求较高；2)难以处理复杂环境，无法对多径进行辨识，需要先假定多径数量后再进行估计功耗较大。因此，亟需提出一种低成本、低功耗的无线信号的到达角估计方法及系统。


技术实现要素：

3.本技术提供一种基于双通道接收机的无线信号到达角估计方法及系统，以至少解决相关技术中无线信号的到达角估计中成本及功耗高的技术问题。
4.本技术第一方面实施例提出一种基于双通道接收机的无线信号到达角估计方法，所述方法包括：
5.利用双通道接收机接收多天线阵列预设时段内各轮天线信号数据；
6.基于所述预设时段内各轮天线信号数据和稀疏重构的最优化算法对接收的多天线阵列对应的天线信号入射强度谱进行估计，得到所述预设时段内各轮天线信号数据对应的天线信号入射强度谱；
7.对所述预设时段内各轮天线信号数据对应的天线信号入射强度谱的谱峰对应的角度进行筛选，并将筛选后的角度按照降序进行排序形成角度序列；
8.将所述角度序列作为无线信号到达角的估计值。
9.本技术第二方面实施例提出一种基于双通道接收机的无线信号到达角估计系统，其特征在于，所述系统包括：
10.接收模块，用于利用双通道接收机接收多天线阵列预设时段内各轮天线信号数据；
11.估计模块，用于基于所述预设时段内各轮天线信号数据和稀疏重构的最优化算法对接收的多天线阵列对应的天线信号入射强度谱进行估计，得到所述预设时段内各轮天线信号数据对应的天线信号入射强度谱；
12.筛选模块，用于对所述预设时段内各轮天线信号数据对应的天线信号入射强度谱的谱峰对应的角度进行筛选，并将筛选后的角度按照降序进行排序形成角度序列；
13.确定模块，用于将所述角度序列作为无线信号到达角的估计值。
14.本技术第三方面实施例提出一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时，实现如本技术第一方
面的预测方法。
15.本技术第四方面实施例提出一种非临时性计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现如本技术第一方面的预测方法。
16.本技术的实施例提供的技术方案至少带来以下有益效果：
17.本发明提供的一种基于双通道接收机的无线信号到达角估计方法及系统，所述方法包括：利用双通道接收机接收多天线阵列预设时段内各轮天线信号数据；基于所述预设时段内各轮天线信号数据和稀疏重构的最优化算法对接收的多天线阵列对应的天线信号入射强度谱进行估计，得到所述预设时段内各轮天线信号数据对应的天线信号入射强度谱；对所述预设时段内各轮天线信号数据对应的天线信号入射强度谱的谱峰对应的角度进行筛选，并将筛选后的角度按照降序进行排序形成角度序列；将所述角度序列作为无线信号到达角的估计值。本发明提供的技术方案，通过稀疏重构的方式对无线信号的到达角进行估计，可以减少所需的接收机数量，有效降低设备的成本与功耗。
18.本技术附加的方面以及优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本技术的实践了解到。
附图说明
19.本技术上述的和/或附加的方面以及优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
20.图1是根据本技术一个实施例提供的一种基于双通道接收机的无线信号到达角估计方法的流程图；
21.图2是根据本技术一个实施例提供的一种基于双通道接收机的无线信号到达角估计系统的结构图；
22.图3是根据本技术一个实施例提供的一种基于双通道接收机的无线信号到达角估计系统中估计模块的结构图；
23.图4是根据本技术一个实施例提供的一种基于双通道接收机的无线信号到达角估计系统中筛选模块的结构图。
具体实施方式
24.下面详细描述本技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本技术，而不能理解为对本技术的限制。
25.本技术提出的一种基于双通道接收机的无线信号到达角估计方法及系统，所述方法包括：利用双通道接收机接收多天线阵列预设时段内各轮天线信号数据；基于所述预设时段内各轮天线信号数据和稀疏重构的最优化算法对接收的多天线阵列对应的天线信号入射强度谱进行估计，得到所述预设时段内各轮天线信号数据对应的天线信号入射强度谱；对所述预设时段内各轮天线信号数据对应的天线信号入射强度谱的谱峰对应的角度进行筛选，并将筛选后的角度按照降序进行排序形成角度序列；将所述角度序列作为无线信号到达角的估计值。本发明提供的技术方案，通过稀疏重构的方式，对无线信号的到达角进行估计，可以减少所需的接收机数量，有效降低设备的成本与功耗。
26.实施例1
27.图1为本公开实施例提供的一种基于双通道接收机的无线信号到达角估计方法的流程图，如图1所示，所述方法包括：
28.步骤1：利用双通道接收机接收多天线阵列预设时段内各轮天线信号数据；
29.步骤2：基于所述预设时段内各轮天线信号数据和稀疏重构的最优化算法对接收的多天线阵列对应的天线信号入射强度谱进行估计，得到所述预设时段内各轮天线信号数据对应的天线信号入射强度谱；
30.在本公开实施例中，所述基于所述预设时段内各轮天线信号数据和稀疏重构的最优化算法对接收的多天线阵列对应的天线信号入射强度谱进行估计，得到所述预设时段内各轮天线信号数据对应的天线信号入射强度谱，包括：
31.步骤2-1：利用所述预设时段内各轮天线信号数据确定所述天线阵列预设时段对应的相关矩阵及阵列流形；
32.步骤2-2：基于所述天线阵列预设时段对应的相关矩阵及阵列流形分别确定所述相关矩阵对应的一维向量及所述天线阵列预设时段对应的阵列流形与单位阵的kronecker积；
33.步骤2-3：基于所述相关矩阵对应的一维向量及所述天线阵列预设时段对应的阵列流形与单位阵的kronecker积和稀疏重构的最优化算法构建目标函数
34.步骤2-4：采用牛顿梯度下降法对所述目标函数进行迭代优化，得到所述天线阵列预设时段对应的天线信号入射强度谱；
35.其中，为相关矩阵对应的一维向量，p为所述天线阵列预设时段对应的阵列流形sv与单位阵的kronecker积，并将对角线元素对应的行的元素值减半，得到所述天线阵列预设时段对应的矩阵，为预设时段内所述多天线的强度投影对应的复数矩阵b的元素沿列向量方向展开得到复数矩阵b对应的一维向量，g为所述复数矩阵b所有行向量的模的累和，λ为惩罚系数，为所述相关矩阵进行稀疏重构后的重构误差，其中，为了实现更好的收敛性能，经过大量数据验证，最终惩罚系数λ取值0.1，并且在生成p时，对角线上的元素的权重只有50％。
36.需要说明的是，所述相关矩阵对应的一维向量的获取过程包括：
37.步骤f1：获取预设时段内各轮天线信号数据中两两天线间的互相关和单个天线间的自相关；
38.步骤f2：根据预设时段内各轮天线信号数据中两两天线间的互相关和单个天线间的自相关确定预设时段内各轮天线信号数据对应的相关矩阵；
39.示例的，若某时刻(即某轮)获得天线数据a＝[a0,a1,a2,a3]
t
,则该时刻对应的相关矩阵为：
[0040][0041]
其中，ah表示a的共轭转置，表示ai的共轭。即为天线i的自相关项，即为天线i和天线j的互相关项。
[0042]
在本实施例中，由于使用了双通道接收机，因此只能同时获得两根天线的接收数据，对于4天线阵列，一组完整的测量数据应该包含次测量，分别为a
0,1
,a
0,2
,a
0,3
,a
1,3
,a
2,3
，其中a
i,j
＝[ai,aj]
t
，同时有因此根据6次不同天线组合的测量结果，在信道短时间保持稳定的近似假设下，可以获得该时刻对应的相关矩阵。
[0043]
步骤f3：对预设时段内各轮天线信号数据对应的相关矩阵求平均值，得到所述天线阵列预设时段对应的相关矩阵ry；
[0044]
步骤f4：将所述天线阵列预设时段对应的相关矩阵ry的元素沿列向量方向展开，并将对角线对应元素减半，得到所述相关矩阵对应的一维向量
[0045]
需要说明的是，所述天线阵列预设时段对应的阵列流形的获取过程包括：
[0046]
基于预设的空间分辨率和多天线阵列的天线阵元确定所述天线阵列预设时段对应的阵列流形。
[0047]
步骤3：对所述预设时段内各轮天线信号数据对应的天线信号入射强度谱的谱峰对应的角度进行筛选，并将筛选后的角度按照降序进行排序形成角度序列；
[0048]
在本公开实施例中，所述对所述预设时段内各轮天线信号数据对应的天线信号入射强度谱的谱峰对应的角度进行筛选，并将筛选后的角度按照降序进行排序形成角度序列，包括：
[0049]
筛选出在预设角度范围内的所述天线信号入射强度谱的谱峰对应的角度；
[0050]
将所述筛选出的角度按照降序进行排序形成角度序列。
[0051]
步骤4：将所述角度序列作为无线信号到达角的估计值。
[0052]
结合上述基于深度学习的发电设备异常检测方法对本技术的具体方法进行举例说明：
[0053]
以下以产品所用的四天线阵列 双通道接收机为例。该设备在一个测量周期内共进行6轮测量，每轮遍历6种天线组合。将同一轮次测量的数据近似视为同时测量，获得的数据可以表示为：d＝{φu(te)|u＝0..5，e＝0..5}，u∈[1,...,n-1],e∈[1,...,k-1],n为每轮的天线组合数，k为一个测量周期内测量的轮数。对于每一个时刻对应的数据，可以分别计算出两两天线之间的互相关以及单个天线的自相关，得到对应的相关矩阵，再对6个时刻的相关矩阵求平均，得到最后的ry。
[0054]
进一步的，空间分辨率选择以1
°
为间隔等分[-90
°
,90
°
]，即θ＝[-90
°
，-89
°
，-88
°
..90
°
]，天线阵元n＝4，载波频率记为f，则阵列流形表示为sv＝exp(1i*[0,1,2,3]
t
*sin
θ)，将ry的元素沿列向量方向展开，并将对角线对应元素减半，得到
[0055]
计算sv矩阵与大小为4的单位阵的kronecker积，并将对角线元素对应的行(1,6,11,16)的元素值减半，得到矩阵p。
[0056]
用复数矩阵b[181][4]表示181个方向上4根天线的强度投影，则即为矩阵b的元素沿列向量方向展开后组成的一维向量，g即为矩阵b所有行向量的模的累和。
[0057]
至此上述目标函数中的参量与变量都构建完成，使用随机初值对矩阵b初始化后，即可采用牛顿梯度下降法对目标函数进行优化。
[0058]
至优化收敛后，由矩阵b所有行向量的模组成的向量spec[181]，即为所求的信号入射强度谱。根据已知的先验经验，对强度谱进行后处理，如：强度低于一定阈值的角度视为无信号；超出一定范围的角度对应的强度置为0等。
[0059]
最后遍历处理后的强度谱，将谱峰对应的角度按照峰值降序排序，即为计算出的对到达角的估计值。
[0060]
综上所述，本实施例提供的一种基于双通道接收机的无线信号到达角估计方法，通过稀疏重构的方式，对无线信号的到达角进行估计，可以减少所需的接收机数量，有效降低设备的成本与功耗。
[0061]
实施例2
[0062]
图2为本公开实施例提供的一种基于双通道接收机的无线信号到达角估计系统的结构图，如图2所示，所述系统包括：
[0063]
接收模块，用于利用双通道接收机接收多天线阵列预设时段内各轮天线信号数据；
[0064]
估计模块，用于基于所述预设时段内各轮天线信号数据和稀疏重构的最优化算法对接收的多天线阵列对应的天线信号入射强度谱进行估计，得到所述预设时段内各轮天线信号数据对应的天线信号入射强度谱；
[0065]
筛选模块，用于对所述预设时段内各轮天线信号数据对应的天线信号入射强度谱的谱峰对应的角度进行筛选，并将筛选后的角度按照降序进行排序形成角度序列；
[0066]
确定模块，用于将所述角度序列作为无线信号到达角的估计值。
[0067]
在本公开实施例中，如图3所示，所述估计模块，包括：
[0068]
第一确定单元，用于利用所述预设时段内各轮天线信号数据确定所述天线阵列预设时段对应的相关矩阵及阵列流形；
[0069]
第二确定单元，用于基于所述天线阵列预设时段对应的相关矩阵及阵列流形分别确定所述相关矩阵对应的一维向量及所述天线阵列预设时段对应的阵列流形与单位阵的kronecker积；
[0070]
构建单元，用于基于所述相关矩阵对应的一维向量及所述天线阵列预设时段对应的阵列流形与单位阵的kronecker积和稀疏重构的最优化算法构建目标函数
[0071]
优化单元，用于采用牛顿梯度下降法对所述目标函数进行迭代优化，得到所述天线阵列预设时段对应的天线信号入射强度谱；
[0072]
其中，为相关矩阵对应的一维向量，p为所述天线阵列预设时段对应的阵列流形sv与单位阵的kronecker积，并将对角线元素对应的行的元素值减半，得到所述天线阵列预设时段对应的矩阵，为预设时段内所述多天线的强度投影对应的复数矩阵b的元素沿列向量方向展开得到复数矩阵b对应的一维向量，g为所述复数矩阵b所有行向量的模的累和，λ为惩罚系数，为所述相关矩阵进行稀疏重构后的重构误差。
[0073]
需要说明的是，所述相关矩阵对应的一维向量的获取过程包括：
[0074]
获取预设时段内各轮天线信号数据中两两天线间的互相关和单个天线间的自相关；
[0075]
根据预设时段内各轮天线信号数据中两两天线间的互相关和单个天线间的自相关确定预设时段内各轮天线信号数据对应的相关矩阵；
[0076]
对预设时段内各轮天线信号数据对应的相关矩阵求平均值，得到所述天线阵列预设时段对应的相关矩阵ry；
[0077]
将所述天线阵列预设时段对应的相关矩阵ry的元素沿列向量方向展开，并将对角线对应元素减半，得到所述相关矩阵对应的一维向量
[0078]
需要说明的是，所述天线阵列预设时段对应的阵列流形的获取过程包括：
[0079]
基于预设的空间分辨率和多天线阵列的天线阵元确定所述天线阵列预设时段对应的阵列流形。
[0080]
在本公开实施例中，如图4所示，所述筛选模块，包括：
[0081]
筛选单元，用于筛选出在预设角度范围内的所述天线信号入射强度谱的谱峰对应的角度；
[0082]
排序单元，用于将所述筛选出的角度按照降序进行排序形成角度序列。
[0083]
综上所述，本实施例提供的一种基于双通道接收机的无线信号到达角估计系统，通过稀疏重构的方式，对无线信号的到达角进行估计，可以减少所需的接收机数量，有效降低设备的成本与功耗。
[0084]
实施例3
[0085]
为了实现上述实施例，本公开还提出一种计算机设备。
[0086]
本实施例提供的计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时，实现实施例1中的方法。
[0087]
实施例4
[0088]
为了实现上述实施例，本公开还提出一种非临时性计算机可读存储介质。
[0089]
本实施例提供的计算机设备，其上存储有计算机程序，计算机程序被处理器执行时实现实施例1中的方法。
[0090]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本技术的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结
合和组合。
[0091]
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本技术的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本技术的实施例所属技术领域的技术人员所理解。
[0092]
尽管上面已经示出和描述了本技术的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本技术的限制，本领域的普通技术人员在本技术的范围内可以对上述实施例进行变化、修改、替换和变型。