CN112595902B - 一种阵列天线系统故障诊断方法及装置 - Google Patents

一种阵列天线系统故障诊断方法及装置 Download PDF

Info

Publication number
CN112595902B
CN112595902B CN202011556274.0A CN202011556274A CN112595902B CN 112595902 B CN112595902 B CN 112595902B CN 202011556274 A CN202011556274 A CN 202011556274A CN 112595902 B CN112595902 B CN 112595902B
Authority
CN
China
Prior art keywords
array
fault
antenna
antenna system
array element
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011556274.0A
Other languages
English (en)
Other versions
CN112595902A (zh
Inventor
孙锐
陈力
卫国
张文逸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN202011556274.0A priority Critical patent/CN112595902B/zh
Publication of CN112595902A publication Critical patent/CN112595902A/zh
Application granted granted Critical
Publication of CN112595902B publication Critical patent/CN112595902B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08—Measuring electromagnetic field characteristics
    • G01R29/10—Radiation diagrams of antennas

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本申请提供的一种阵列天线系统故障诊断方法及装置，通过对待诊断阵列天线系统的空域无线信号(均匀线性阵列采样信号或均匀平面阵列采样信号)的采样与处理，根据待诊断阵列天线系统的阵列构型(均匀线性阵列或均匀平面阵列)，确定与阵列构型对应的故障诊断方法；将空域无线信号作为输入，运行故障诊断方法，得到天线阵元的故障响应系数；根据天线阵元的故障响应系数，得到待诊断阵列天线系统的故障天线阵元的索引。本申请通过对待诊断阵列天线系统的空域无线信号的采样与处理，故障诊断能在无参考信号的条件下进行，可以检测并定位阵列天线系统中的故障器件，相比传统的阵列天线系统诊断方法，能够降低故障诊断的复杂度与硬件开销。

Description

一种阵列天线系统故障诊断方法及装置
技术领域
本申请涉及无线通信技术领域，尤其涉及一种阵列天线系统故障诊断方法及装置，可在无参考信号的条件下实施故障诊断。
背景技术
在当前以及未来的无线通信系统中，阵列天线是关键技术之一。近年来，阵列天线系统故障诊断方法得到了广泛研究，目前提出的阵列天线系统故障诊断方法需要依赖于精确的参考信号，参考信号可以通过在微波暗室中特定位置部署参考信号发射源并根据空间坐标计算得到，在微波暗室中可精确控制无线环境，来确保参考信号的精确性。然而，这样的参考信号生成方法不适用于户外场景中的在线诊断。
因此，在不依赖精确参考信号的情况下，如何实现对阵列天线系统中故障天线阵元的检测与定位，是本领域技术人员亟待解决的问题。
发明内容
本申请提供了一种阵列天线系统故障诊断方法及装置，在不依赖精确参考信号的情况下，实现对阵列天线系统中故障天线阵元的检测与定位。
为了实现上述目的，本申请提供了以下技术方案：
一种阵列天线系统故障诊断方法，包括：
获取待诊断阵列天线系统的空域无线信号，得到信号y或者Y，其中，y为均匀线性阵列采样信号，Y为均匀平面阵列采样信号；
根据所述待诊断阵列天线系统的阵列构型，确定与所述阵列构型对应的故障诊断方法，所述陈列构型为均匀线性阵列或均匀平面阵列；
将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数；
根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀线性阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式γ＝[γ1,γ2,...,γG]计算得到离散化角度集合γ，其中，
Figure BDA0002858870250000021
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量；
根据公式Aγ＝[a(γ1)，...，a(γG)]计算得到名义观测矩阵Aγ，其中，
Figure BDA0002858870250000022
为γg对应的导向矢量；
根据公式
Figure BDA0002858870250000023
计算得到校准矩阵B；
根据公式
Figure BDA0002858870250000024
计算得到所述天线阵元的故障响应系数
Figure BDA0002858870250000025
其中，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀平面阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式
Figure BDA0002858870250000026
计算得到天线阵元的故障响应系数
Figure BDA0002858870250000027
其中||·||*表示核范数，λ为正则化参数。
优选的，所述根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引，具体为：
根据公式|fn|＞th确定满足条件的天线阵元索引，所述满足条件的天线阵元索引作为所述待诊断阵列天线系统的故障天线阵元的索引，其中，fn为第n个天线阵元的故障响应系数，th为判决门限。
一种阵列天线系统故障诊断装置，包括：
第一处理单元，用于获取待诊断阵列天线系统的空域无线信号，得到信号y或者Y，其中，y为均匀线性阵列采样信号，Y为均匀平面阵列采样信号；
第二处理单元，用于根据所述待诊断阵列天线系统的阵列构型，确定与所述阵列构型对应的故障诊断方法，所述陈列构型为均匀线性阵列或均匀平面阵列；
第三处理单元，用于将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数；
第四处理单元，用于根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀线性阵列时，所述第三处理单元具体用于：
将所述空域无线信号作为输入，根据公式γ＝[γ1,γ2,...,γG]计算得到离散化角度集合γ，其中，
Figure BDA0002858870250000031
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量；
根据公式Aγ＝[a(γ1)，...，a(γG)]计算得到名义观测矩阵Aγ，其中，
Figure BDA0002858870250000032
为γg对应的导向矢量；
根据公式
Figure BDA0002858870250000033
计算得到校准矩阵B；
根据公式
Figure BDA0002858870250000034
计算得到所述天线阵元的故障响应系数
Figure BDA0002858870250000035
其中，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀平面阵列时，所述第三处理单元具体用于：
将所述空域无线信号作为输入，根据公式
Figure BDA0002858870250000041
计算得到天线阵元的故障响应系数
Figure BDA0002858870250000042
其中||·||*表示核范数，λ为正则化参数。
优选的，所述第四处理单元具体用于：
根据公式|fn|＞th确定满足条件的天线阵元索引，所述满足条件的天线阵元索引作为所述待诊断阵列天线系统的故障天线阵元的索引，其中，fn为第n个天线阵元的故障响应系数，th为判决门限。
一种存储介质，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在的设备执行如上述所述的阵列天线系统故障诊断方法。
一种电子设备，所述电子设备包括至少一个处理器、以及与所述处理器连接的至少一个存储器、总线；其中，所述处理器、所述存储器通过所述总线完成相互间的通信；所述处理器用于调用所述存储器中的程序指令，以执行如上述所述的阵列天线系统故障诊断方法。
本申请提供的一种阵列天线系统故障诊断方法及装置，通过对待诊断阵列天线系统的空域无线信号(均匀线性阵列采样信号或均匀平面阵列采样信号)的采样与处理，根据待诊断阵列天线系统的阵列构型(均匀线性阵列或均匀平面阵列)，确定与阵列构型对应的故障诊断方法；将空域无线信号作为输入，运行故障诊断方法，得到天线阵元的故障响应系数；根据天线阵元的故障响应系数，得到待诊断阵列天线系统的故障天线阵元的索引。本申请通过对待诊断阵列天线系统的空域无线信号的采样与处理，故障诊断能在无参考信号的条件下进行，可以检测并定位阵列天线系统中的故障器件，相比传统的阵列天线系统诊断方法，能够降低故障诊断的复杂度与硬件开销。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的一种阵列天线系统故障诊断方法流程示意图；
图2为本申请实施例提供的阵列天线系统诊断方法的均匀线性阵列的诊断结果展示示意图；
图3为本申请实施例提供的阵列天线系统诊断方法的均匀平面阵列的诊断结果展示示意图；
图4为本申请实施例提供的阵列天线系统诊断方法的诊断成功概率性能示意图；
图5为本申请实施例提供的一种阵列天线系统故障诊断装置结构示意图；
图6为本申请实施例提供的一种电子设备结构示意图。
具体实施方式
申请人发现，随着对系统吞吐量要求的不断增加，阵列天线系统中的器件数量急剧提升，系统结构也日益复杂。然而，阵列天线系统的可靠性问题也随着系统规模的增加而日益突出。例如，集成在单个基板上的多个天线阵元很容易受到外部环境因素(例如沙尘、水渍等)的影响而产生遮挡覆盖，从而导致发射/接收信号功率降低，另外，一些毫米波频段器件较低的热效率和较高的运行频率也直接导致了器件本身的低可靠性。此外，一些新型的阵列天线系统架构(例如混合波束赋形系统)中采用的混合交叉连接结构使得系统复杂性进一步增加，同时也导致器件之间的高度相关性，使得故障器件的影响极易扩散至整个系统，当系统中存在故障器件时，会导致接收信号功率降低、天线方向图失真等一系列问题，最终使得系统性能大幅下降。因此，用于故障器件检测与定位的诊断方法是阵列天线系统维护工作中不可或缺的一部分。
然而，目前提出的阵列天线系统故障诊断方法需要依赖于精确的参考信号，参考信号可以通过在微波暗室中特定位置部署参考信号发射源并根据空间坐标计算得到，在微波暗室中可精确控制无线环境，来确保参考信号的精确性。而对于户外的无线环境而言，户外的无线环境不可能像微波暗室一样可以精确控制，而会受到多径散射等各种衰落因素的影响导致无法根据空间坐标可靠计算参考信号；另外，受限于户外在线诊断场景，参考信号发射源的部署也会极为困难，因此，通过控制无线环境来生成参考信号的方法不适用于户外场景中的在线诊断。
为此，本申请提供一种阵列天线系统故障诊断方法及装置，
本申请的发明目的在于：在不依赖精确参考信号的情况下，实现对阵列天线系统中故障天线阵元的检测与定位。
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
请参见附图2，本申请实施例提供一种阵列天线系统故障诊断方法流程图。如图2所示，本申请实施例提供一种阵列天线系统故障诊断方法，该方法具体包括如下步骤：
S11：获取待诊断阵列天线系统的空域无线信号，得到均匀线性阵列采样信号或均匀平面阵列采样信号。
其中，均匀线性阵列采样信号记作y，均匀平面阵列采样信号记作Y。
S12：根据所述待诊断阵列天线系统的阵列构型，确定与所述阵列构型对应的故障诊断方法，所述陈列构型为均匀线性阵列或均匀平面阵列。
S13：将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数。
需要说明的是，本申请实施例中，分别以均匀线性阵列和均匀平面阵列的诊断为例进行介绍。
本申请实施例中，当所述待诊断阵列天线系统的阵列构型为均匀线性阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式γ＝[γ1,γ2,...,γG]计算得到离散化角度集合γ，其中，
Figure BDA0002858870250000071
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量。
根据公式Aγ＝[a(γ1)，...，a(γG)]计算得到名义观测矩阵Aγ，其中，
Figure BDA0002858870250000072
为γg对应的导向矢量。
根据公式
Figure BDA0002858870250000073
计算得到校准矩阵B。
根据公式
Figure BDA0002858870250000074
计算得到所述天线阵元的故障响应系数
Figure BDA0002858870250000075
其中，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平。
具体的，对于均匀线性阵列：
首先，可以将[-π/2,π/2]的连续角域进行离散化，得到一个包含有G＝cN个离散角度值的离散化角度集合γ＝[γ1,γ2,...,γG]，其中
Figure BDA0002858870250000076
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量，特别地，定义2r＝sin(γg+1-γg)，g＝1，...，G-1为步长，r为步长。
现假设有K个信号入射阵列，由于实际的入射信号角度是在[-π/2,π/2]区间中连续取值的，因此将κ＝[κ1,κ2,...,κG]定义为真实但未知的离散化角度栅格，其中同样包含有G＝cN个离散角度值。则接收到的信号可写为如下形式：
Figure BDA0002858870250000081
其中，xg是第g个格点对应的入射信号的增益，
Figure BDA0002858870250000082
为κg对应的导向矢量；Aκ＝[a(κ1)，...，a(κG)]为真实观测矩阵；x＝[x1,...,xG]为入射信号的增益。
根据一阶泰勒展开，有下式成立：
y′＝Aκx≈(Aγ+B diag(β))x
其中，
Figure BDA0002858870250000083
为校准矩阵；β＝sin(κ)-sin(γ)，由于|sin(κi)-si(nγi)|≤r,≤1i≤G，β中的每个元素取值范围为[-r,r]。
当系统中存在故障天线阵元和观测噪声时，所观测得到的信号可以写为如下形式：
y＝(Aγ+B diag(β))x+f+w
其中，f为故障响应系数，完好天线阵元的故障响应系数为0，而故障天线阵元的故障响应系数不为0；w为观测噪声。
注意到：入射信号的数量是要远小于天线阵元总数N的，因此，x为一个稀疏信号(稀疏度为K)；而由于一般情况下故障阵元数量是要远低于系统中天线阵元总数的，因此，f也是一个稀疏信号。而从公式y＝(Aγ+B diag(β))x+f+w可以看出，信号x和信号f分别在字典(Aγ+B diag(β))和单位阵下呈现出稀疏特性。因此，可以通过求解如下优化问题获得信号x和信号f：
Figure BDA0002858870250000084
其中，s.t.表示“subject to”，即“受限于”，这是数学中常用的缩写，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平；信号
Figure BDA0002858870250000085
中的非零元素标志着故障天线阵元的索引和故障响应系数；
Figure BDA0002858870250000086
为对应变量的估计，同时还可以看出，本申请实施例提供的诊断方法直接采样空域无线信号即可完成诊断，而无需专门的参考信号。
如图2所示，为本申请实施例提供的阵列天线系统诊断方法的均匀线性阵列的诊断结果展示示意图。其中阵列天线的阵元数量为64，故障天线阵元的索引为32和50；入射信号个数为4，其增益服从独立且相同的单位复高斯分布且到达角服从[-π/2,π/2]区间内的均匀分布；信号噪声功率比为30分贝。图(a)为阵列各个阵元的接收信号强度，图(b)为各个阵元的故障响应系数。从图2可以看出，在无参考信号的情况下是无法判断故障阵元的，而本发明所提出的诊断方法可以成功判断出故障阵元的位置并求出故障响应系数。
本申请实施例中，当所述待诊断阵列天线系统的阵列构型为均匀平面阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式
Figure BDA0002858870250000091
计算得到天线阵元的故障响应系数
Figure BDA0002858870250000092
其中||·||*表示核范数，λ为正则化参数。
具体的，对于均匀平面阵列：
对于均匀平面阵列，可以将阵列接收到的信号表示为如下形式：
Figure BDA0002858870250000093
其中，αk表示第k个入射信号的增益，K为入射信号的数量；H和W分别为垂直和水平方向的阵元数量；
Figure BDA0002858870250000094
和
Figure BDA0002858870250000095
分别为垂直方向和水平方向的导向矢量；θk和φk分别为第k个入射信号的水平和垂直到达角；d为相对于波长的阵元间距。
当系统中存在故障天线阵元和观测噪声时，所观测得到的信号可以写为如下形式：
Y＝Y′+F+W
其中，F为故障响应系数，完好天线阵元的故障响应系数为0，而故障天线阵元的故障响应系数不为0；W为观测噪声。
注意到：入射信号的数量是要远小于min(H,W)的，因此Y′为一个低秩矩阵(秩为K)；而由于一般情况下故障阵元数量是要远低于系统中天线阵元总数的，因此，F为一个稀疏矩阵。因此，可以采用稀疏与低秩矩阵分解算法求解如下的优化问题从而得到F：
Figure BDA0002858870250000101
其中，||·||*表示核范数，λ为正则化参数；信号
Figure BDA0002858870250000102
中的非零元素标志着故障阵元的索引和故障响应系数。同时还可以看出，本申请实施例提供的诊断方法直接采样空域无线信号即可完成诊断，而无需专门的参考信号。
如图3所示，为本申请实施例提供的阵列天线系统诊断方法的均匀平面阵列的诊断结果展示示意图。其中阵列天线的阵元数量为4096，故障天线阵元的索引为(16，62)、(14，35)、(29，15)、(51，36)、(47，11)；入射信号个数为4，其增益服从独立且相同的单位复高斯分布且到达角服从[-π/2,π/2]区间内的均匀分布；信号噪声比为30分贝。图(a)为阵列各个阵元的接收信号强度，图(b)为各个阵元的故障响应系数。从图3可以看出，在无参考信号的情况下是无法判断故障阵元的，而本发明所提出的诊断方法可以成功判断出故障阵元的位置并求出故障响应系数。
S14：根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引。
本申请实施例中，上述所述根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引，具体可以为：
根据公式|fn|＞th确定满足条件的天线阵元索引，所述满足条件的天线阵元索引作为所述待诊断阵列天线系统的故障天线阵元的索引，其中，fn为第n个天线阵元的故障响应系数，th为判决门限。
如图4所示，为本申请实施例提供的阵列天线系统诊断方法的诊断成功概率性能示意图。图4以诊断成功概率为指标。成功概率越高说明诊断效果越好。设置不同的故障阵元数为1、2、3(图(a)，均匀线性阵列)以及1、5、10(图(b)，均匀平面阵列)。图3(a)的三条曲线分别用圆形“●”、交叉符号“×”以及空心圆形“○”为曲线的节点作为区分，图3(b)的三条曲线同样分别用圆形“●”、交叉符号“×”以及空心圆形“○”为曲线的节点作为区分。随着信号噪声功率比的增加，成功概率逐渐增加。另外，随着故障阵元数量的增加，相同信号噪声功率比下的成功概率是逐渐降低的。
本申请实施例提供的一种阵列天线系统故障诊断方法，通过对待诊断阵列天线系统的空域无线信号(均匀线性阵列采样信号或均匀平面阵列采样信号)的采样与处理，根据待诊断阵列天线系统的阵列构型(均匀线性阵列或均匀平面阵列)，确定与阵列构型对应的故障诊断方法；将空域无线信号作为输入，运行故障诊断方法，得到天线阵元的故障响应系数；根据天线阵元的故障响应系数，得到待诊断阵列天线系统的故障天线阵元的索引。本申请实施例通过对待诊断阵列天线系统的空域无线信号的采样与处理，故障诊断能在无参考信号的条件下进行，可以检测并定位阵列天线系统中的故障器件，相比传统的阵列天线系统诊断方法，能够降低故障诊断的复杂度与硬件开销。
请参阅图5，基于上述实施例公开的一种阵列天线系统故障诊断方法，本实施例对应公开了一种阵列天线系统故障诊断装置，具体的，该装置包括：
第一处理单元51，用于获取待诊断阵列天线系统的空域无线信号，得到均匀线性阵列采样信号或均匀平面阵列采样信号；
第二处理单元52，用于根据所述待诊断阵列天线系统的阵列构型，确定与所述阵列构型对应的故障诊断方法，所述陈列构型为均匀线性阵列或均匀平面阵列；
第三处理单元53，用于将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数；
第四处理单元54，用于根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀线性阵列时，所述第三处理单元53具体用于：
将所述空域无线信号作为输入，根据公式γ＝[γ1,γ2,...,γG]计算得到离散化角度集合γ，其中，
Figure BDA0002858870250000121
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量；
根据公式Aγ＝[a(γ1)，...，a(γG)]计算得到名义观测矩阵Aγ，其中，
Figure BDA0002858870250000122
为γg对应的导向矢量；
根据公式
Figure BDA0002858870250000123
计算得到校准矩阵B；
根据公式
Figure BDA0002858870250000124
计算得到所述天线阵元的故障响应系数
Figure BDA0002858870250000125
其中，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀平面阵列时，所述第三处理单元53具体用于：
将所述空域无线信号作为输入，根据公式
Figure BDA0002858870250000126
计算得到天线阵元的故障响应系数
Figure BDA0002858870250000127
其中||·||*表示核范数，λ为正则化参数。
优选的，所述第四处理单元54具体用于：
根据公式|fn|＞th确定满足条件的天线阵元索引，所述满足条件的天线阵元索引作为所述待诊断阵列天线系统的故障天线阵元的索引，其中，fn为第n个天线阵元的故障响应系数，th为判决门限。
所述阵列天线系统故障诊断装置包括处理器和存储器，上述第一处理单元、第二处理单元、第三处理单元和第四处理单元等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序单元来实现相应的功能。
处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来达到在不依赖精确参考信号的情况下，实现对阵列天线系统中故障天线阵元的检测与定位，相比传统的阵列天线系统诊断方法，能够降低故障诊断的复杂度与硬件开销。
本申请实施例提供了一种存储介质，其上存储有程序，该程序被处理器执行时实现所述阵列天线系统故障诊断方法。
本申请实施例提供了一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行所述阵列天线系统故障诊断方法。
本申请实施例提供了一种电子设备，如图6所示，该电子设备60包括至少一个处理器601、以及与所述处理器连接的至少一个存储器602、总线603；其中，所述处理器601、所述存储器602通过所述总线603完成相互间的通信；处理器601用于调用所述存储器602中的程序指令，以执行上述的所述阵列天线系统故障诊断方法。
本文中的电子设备可以是服务器、PC、PAD、手机等。
本申请还提供了一种计算机程序产品，当在数据处理设备上执行时，适于执行初始化有如下方法步骤的程序：
获取待诊断阵列天线系统的空域无线信号，得到信号y或者Y，其中，y为均匀线性阵列采样信号，Y为均匀平面阵列采样信号；
根据所述待诊断阵列天线系统的阵列构型，确定与所述阵列构型对应的故障诊断方法，所述陈列构型为均匀线性阵列或均匀平面阵列；
将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数；
根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀线性阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式γ＝[γ1,γ2,...,γG]计算得到离散化角度集合γ，其中，
Figure BDA0002858870250000141
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量；
根据公式Aγ＝[a(γ1)，...，a(γG)]计算得到名义观测矩阵Aγ，其中，
Figure BDA0002858870250000142
为γg对应的导向矢量；
根据公式
Figure BDA0002858870250000143
计算得到校准矩阵B；
根据公式
Figure BDA0002858870250000144
计算得到所述天线阵元的故障响应系数
Figure BDA0002858870250000145
其中，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平。
优选的，当所述待诊断阵列天线系统的阵列构型为均匀平面阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式
Figure BDA0002858870250000146
计算得到天线阵元的故障响应系数
Figure BDA0002858870250000147
其中||·||*表示核范数，λ为正则化参数。
优选的，所述根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引，具体为：
根据公式|fn|＞th确定满足条件的天线阵元索引，所述满足条件的天线阵元索引作为所述待诊断阵列天线系统的故障天线阵元的索引，其中，fn为第n个天线阵元的故障响应系数，th为判决门限。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
在一个典型的配置中，设备包括一个或多个处理器(CPU)、存储器和总线。设备还可以包括输入/输出接口、网络接口等。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)，存储器包括至少一个存储芯片。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (10)

1.一种阵列天线系统故障诊断方法，其特征在于，包括：
获取待诊断阵列天线系统的空域无线信号，得到信号y或者Y，其中，y为均匀线性阵列采样信号，Y为均匀平面阵列采样信号；
根据所述待诊断阵列天线系统的阵列构型，确定与所述阵列构型对应的故障诊断方法，所述阵列构型为均匀线性阵列或均匀平面阵列；
将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，其中，完好天线阵元的故障响应系数为0，故障天线阵元的故障响应系数不为0；
根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引，所述故障天线阵元的索引为用于标识所述故障天线阵元的标识信息。
2.根据权利要求1所述的方法，其特征在于，当所述待诊断阵列天线系统的阵列构型为均匀线性阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式γ＝[γ1,γ2,...,γG]计算得到离散化角度集合γ，其中，
Figure FDA0003604128150000011
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量；
根据公式Aγ＝[a(γ1)，...，a(γG)]计算得到名义观测矩阵Aγ，其中，
Figure FDA0003604128150000012
为γg对应的导向矢量；
根据公式
Figure FDA0003604128150000013
计算得到校准矩阵B；
根据公式
Figure FDA0003604128150000014
s.t.||y-(Aγ+Bdiag(β))x-f||2≤δ计算得到所述天线阵元的故障响应系数
Figure FDA0003604128150000015
其中，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平，
Figure FDA0003604128150000021
和
Figure FDA0003604128150000022
分别为入射信号增益矢量x和校准系数矢量β的估计。
3.根据权利要求1所述的方法，其特征在于，当所述待诊断阵列天线系统的阵列构型为均匀平面阵列时，所述将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，具体为：
将所述空域无线信号作为输入，根据公式
Figure FDA0003604128150000023
s.t.||Y-Y′-F||2≤δ计算得到天线阵元的故障响应系数
Figure FDA0003604128150000024
其中||·||*表示核范数，λ为正则化参数，
Figure FDA0003604128150000025
为对均匀平面阵列无故障接收信号Y′的估计。
4.根据权利要求1所述的方法，其特征在于，所述根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引，具体为：
根据公式|fn|＞th确定满足条件的天线阵元索引，所述满足条件的天线阵元索引作为所述待诊断阵列天线系统的故障天线阵元的索引，其中，fn为第n个天线阵元的故障响应系数，th为判决门限。
5.一种阵列天线系统故障诊断装置，其特征在于，包括：
第一处理单元，用于获取待诊断阵列天线系统的空域无线信号，得到信号y或者Y，其中，y为均匀线性阵列采样信号，Y为均匀平面阵列采样信号；
第二处理单元，用于根据所述待诊断阵列天线系统的阵列构型，确定与所述阵列构型对应的故障诊断方法，所述阵列构型为均匀线性阵列或均匀平面阵列；
第三处理单元，用于将所述空域无线信号作为输入，运行所述故障诊断方法，得到天线阵元的故障响应系数，其中，完好天线阵元的故障响应系数为0，故障天线阵元的故障响应系数不为0；
第四处理单元，用于根据所述天线阵元的故障响应系数，得到所述待诊断阵列天线系统的故障天线阵元的索引，所述故障天线阵元的索引为用于标识所述故障天线阵元的标识信息。
6.根据权利要求5所述的装置，其特征在于，当所述待诊断阵列天线系统的阵列构型为均匀线性阵列时，所述第三处理单元具体用于：
将所述空域无线信号作为输入，根据公式γ＝[γ1,γ2,...,γG]计算得到离散化角度集合γ，其中，
Figure FDA0003604128150000031
为第g个离散角度的取值，d为天线阵元相对于波长的间距，G＝cN为离散化角度的数量，c为过采样系数，N为天线阵元数量；
根据公式Aγ＝[a(γ1)，...，a(γG)]计算得到名义观测矩阵Aγ，其中，
Figure FDA0003604128150000032
为γg对应的导向矢量；
根据公式
Figure FDA0003604128150000033
计算得到校准矩阵B；
根据公式
Figure FDA0003604128150000034
s.t.||y-(Aγ+Bdiag(β))x-f||2≤δ计算得到所述天线阵元的故障响应系数
Figure FDA0003604128150000035
其中，||·||1表示l1范数，||·||2表示l2范数，δ为噪声水平，
Figure FDA0003604128150000036
和
Figure FDA0003604128150000037
分别为入射信号增益矢量x和校准系数矢量β的估计。
7.根据权利要求5所述的装置，其特征在于，当所述待诊断阵列天线系统的阵列构型为均匀平面阵列时，所述第三处理单元具体用于：
将所述空域无线信号作为输入，根据公式
Figure FDA0003604128150000038
s.t.||Y-Y′-F||2≤δ计算得到天线阵元的故障响应系数
Figure FDA0003604128150000039
其中||·||*表示核范数，λ为正则化参数，
Figure FDA00036041281500000310
为对均匀平面阵列无故障接收信号Y′的估计。
8.根据权利要求5所述的装置，其特征在于，所述第四处理单元具体用于：
根据公式|fn|＞th确定满足条件的天线阵元索引，所述满足条件的天线阵元索引作为所述待诊断阵列天线系统的故障天线阵元的索引，其中，fn为第n个天线阵元的故障响应系数，th为判决门限。
9.一种存储介质，其特征在于，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在的设备执行如权利要求1至4中任一项所述的阵列天线系统故障诊断方法。
10.一种电子设备，其特征在于，所述电子设备包括至少一个处理器、以及与所述处理器连接的至少一个存储器、总线；其中，所述处理器、所述存储器通过所述总线完成相互间的通信；所述处理器用于调用所述存储器中的程序指令，以执行如权利要求1至4中任一项所述的阵列天线系统故障诊断方法。
CN202011556274.0A 2020-12-24 2020-12-24 一种阵列天线系统故障诊断方法及装置 Active CN112595902B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011556274.0A CN112595902B (zh) 2020-12-24 2020-12-24 一种阵列天线系统故障诊断方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011556274.0A CN112595902B (zh) 2020-12-24 2020-12-24 一种阵列天线系统故障诊断方法及装置

Publications (2)

Publication Number Publication Date
CN112595902A CN112595902A (zh) 2021-04-02
CN112595902B true CN112595902B (zh) 2022-07-15

Family

ID=75202517

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011556274.0A Active CN112595902B (zh) 2020-12-24 2020-12-24 一种阵列天线系统故障诊断方法及装置

Country Status (1)

Country Link
CN (1) CN112595902B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107315183A (zh) * 2017-06-01 2017-11-03 西南电子技术研究所（中国电子科技集团公司第十研究所） 导航卫星阵列天线接收系统的校准方法
CN111398917A (zh) * 2020-02-27 2020-07-10 四川九洲防控科技有限责任公司 一种天线波束指向测试方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101281997B (zh) * 2007-04-05 2012-05-30 电信科学技术研究院 智能天线故障处理方法及装置
CN102854401A (zh) * 2012-09-11 2013-01-02 电子科技大学 一种阵列天线方向图时域测量方法
CN105277826B (zh) * 2015-10-23 2019-01-25 电子科技大学 阵列天线故障诊断系统
CN105158621B (zh) * 2015-10-23 2018-03-09 电子科技大学 基于分区的阵列天线快速故障诊断方法
CN106053968B (zh) * 2016-06-03 2019-09-10 航天科工空间工程发展有限公司 使用时频分析测量的单道天线远场天线因子估计方法
EP3306327B1 (en) * 2016-10-06 2019-06-12 Rohde & Schwarz GmbH & Co. KG Antenna array, test system and method for testing a device under test
CN109813969B (zh) * 2017-11-21 2021-08-03 深圳市通用测试系统有限公司 一种阵列天线的诊断方法、设备、系统
CN108932381B (zh) * 2018-06-27 2022-10-11 电子科技大学 一种考虑阵列误差的天线阵列故障诊断方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107315183A (zh) * 2017-06-01 2017-11-03 西南电子技术研究所（中国电子科技集团公司第十研究所） 导航卫星阵列天线接收系统的校准方法
CN111398917A (zh) * 2020-02-27 2020-07-10 四川九洲防控科技有限责任公司 一种天线波束指向测试方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"机载雷达天线空域稳定测试系统的设计";齐兴龙;《吉林大学学报》;20120531;第30卷(第3期);第261-264页 *

Also Published As

Publication number Publication date
CN112595902A (zh) 2021-04-02

Similar Documents

Publication Publication Date Title
US9628996B2 (en) Spectrum sensing apparatus and method for cooperative cognitive radio network in non-gaussian noise environment, and fusion center apparatus and cooperative cognitive radio system using the same
Gershman et al. Experimental performance of adaptive beamforming in a sonar environment with a towed array and moving interfering sources
Wan et al. Autonomous vehicle source enumeration exploiting non-cooperative UAV in software defined internet of vehicles
Wang et al. Interference DOA estimation and suppression for GNSS receivers using fully augmentable arrays
CN111814960B (zh) 天气气候系统的识别方法、装置、存储介质及处理器
CN117075154A (zh) 卫星导航阵列抗干扰接收机在给定场景下的载噪比估计方法
CN112595902B (zh) 一种阵列天线系统故障诊断方法及装置
CN110196068B (zh) 一种极区集中滤波组合导航系统残差向量故障检测与隔离方法
CN114254265B (zh) 基于统计流形距离的卫星通信干扰几何分析方法
Singh et al. Discrete interferences optimum beamformer in correlated signal and interfering noise
Tian et al. A low-complexity DOA estimation algorithm for distributed source localization
Gu et al. Robust adaptive beamforming for steering vector uncertainties based on equivalent DOAs method
CN111079591B (zh) 基于改进多尺度主成分分析的不良数据修复方法及系统
CN106443612B (zh) 一种基于信息熵的多雷达混合检测方法
Chen et al. A strong interference suppressor for satellite signals in GNSS receivers
CN108008351B (zh) 一种基于正向衰减逆推发射幅度聚集性的信号定位方法
CN115616477A (zh) 一种非平面波单脉冲测角方法
Yang et al. A modified robust algorithm against large look direction error based on interference-plus-noise covariance matrix reconstruction and steering vector double estimation
CN113965236A (zh) 适用于卫星通信的高鲁棒性自适应波束赋形方法和装置
Richmond Response of sample covariance based MVDR beamformer to imperfect look and inhomogeneities
Sorek et al. Robust Spectrum Sensing Via Probability Measure Transform
Chu et al. A robust algorithm for adaptive interference cancellation
Liang et al. On the Cramér-Rao Bound and the number of resolvable sources in the presence of nonuniform noise for underdetermined DOA estimation
CN115792792B (zh) 极低样本下基于协方差矩阵迭代估计的测向方法与系统
Lu et al. Virtual Reconstruction‐Based Robust Adaptive Beamforming for Distributed Digital Subarray Antennas

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant