一种传感器列系统最大似然角度分辨率估计方法技术方案

本发明专利技术提出一种传感器阵列系统最大似然角度分辨率估计方法，包括如下步骤：建立多传感器阵元的均匀线性阵列，利用天线的接收信号得到观测值的似然函数，通过对似然函数作两维遍历和峰值搜索得到实际估计的角度值，并根据多次快拍下统计到的角度的后验概率分布来判断阵列传感器系统对相近信源的分辨能力，最终通过幅相信息计算得到理论的分辨率结果。本发明专利技术可以通过后验概率出现的峰的个数判断信源是否可分辨，且本发明专利技术的分辨率精确度优于现有技术得到的分辨率的精确度。

【技术实现步骤摘要】
一种传感器列系统最大似然角度分辨率估计方法
本专利技术涉及本专利技术属于信息传输与处理
，尤其是一种传感器列系统最大似然角度分辨率估计方法。
技术介绍
多年以来，利用传感器阵列来接收空间信号并对其进行信号处理已经在很多领域内得到广泛的研究和发展，例如通信、雷达、声纳和地震探测等。接收空间的信号相当于对空间分布的场信号采样而得到空间离散数据。这样做能够增强所需的信号，抑制无用的干扰和噪声，并提取有用的信号特征和信号信息。传感器阵列具有高信号增益、灵活波束控制和超分辨等特点。在均匀线阵的DOA估计中，已经提出了许多空间谱估计方法，如最大似然(ML)方法及其变型，即MODE方法，主要是通过似然函数的多维优化实现渐近最优性能。最早的超分辨DOA估计方法就是著名的多信号分类(MUSIC)算法和通过旋转不变技术(ESPRIT)算法。这些算法是基于接收信号协方差矩阵特征分解的DOA估计方法，通过数值稳定的奇异值分解精确的确定波达方向。由于利用了特征子空间分解，因此DOA估计可以实现较高的分辨率。传统雷达的分辨率为带宽的倒数，针对传感器阵列系统来说分辨率是阵列孔径的倒数，目前，DOA估计的分辨率主要取决于阵列的长度，阵列长度确定后，其分辨率也就被确定了，这个分辨率成为瑞利限。一般的雷达都会受到瑞利限的限制。
技术实现思路
专利技术目的：本专利技术一种传感器列系统最大似然角度分辨率估计方法，得到的分辨率可以突破瑞利限的限制。技术方案：为实现上述目的，本专利技术提出以下技术方案：一种传感器阵列系统最大似然角度分辨率估计方法，所述传感器阵列为多传感器阵元的均匀线性阵列，且各阵元辐射都是全向性的；该方法对两个信源进行探测，包括步骤：(1)利用所述传感器阵列天线的接收信号得到观测值的似然函数；(2)根据步骤(1)得到的似然函数得到实际角度的估计值，具体步骤为：获得N组接收信号，对每一组接收信号，得到此时的似然函数，然后在二维坐标系中的角度搜索范围内，以角度间隔Δθ进行峰值搜索，在搜索过程中，记录每一个DOA角度向量θ对应的似然函数的变化量，将似然函数的变化量最大的DOA角度向量记为θ0，θ0即为相应实际角度的估计值，θ0在二维坐标系中的坐标记为(θ10，θ20)，θ1、θ2为二维坐标系的两条坐标轴；(3)对步骤(2)得到的N组接收信号的实际角度估计值在二维坐标系中的坐标进行统计，根据每个角度间隔Δθ内估计值坐标出现的次数，统计出估计值的后验概率密度函数p(θ|y)，其中，y为接收信号向量；统计后验概率密度函数p(θ|y)的峰，当出现两个峰时，表示两个信源不相关，信源可分辨，转入步骤(4)；而当仅出现一个峰时，表示两个信源相关，信源不可分辨，则结束本方法；(4)计算接收信号的幅相信息：其中，X表示信源矩阵，Y表示接收信号矩阵，I(X；Y)表示X与Y的互信息。表示两个信源的平均信噪比，λ1表示幅相信息的第一个特征值，λ2表示幅相信息的第二个特征值；N0表示噪声方差，α1表示第一个信源的衰减因子，α2表示第二个信源的衰减因子，E[]表示对信源信号求期望；(5)根据接收信号的幅相信息得到最大似然分辨率Δ：其中，M表示阵列天线的阵元数量，γ2为均方根带宽，d表示相邻阵元之间的位置间距，λ表示波长，θ′表示角度的中线轴，即进一步的，所述似然函数的表达式为：其中，RY为接收信号的方差，RY＝E[YYH]＝A(θ)PAH(θ)+N0U；A(θ)表示传感器阵列的导向矩阵，AH(θ)为A(θ)的共轭转置，U为单位矩阵，P为信源的自相关矩阵，有益效果：与现有技术相比，本专利技术具有以下优势：传统雷达的分辨率为带宽的倒数，针对传感器阵列系统来说分辨率是阵列孔径的倒数，而该方法通过后验概率密度观察峰的位置，可以很明确体现出两个信源能否分辨。最后计算得到的理论分辨率验证了其优于传统分辨率的结论。仿真结果表明，当两个信源相关时，只有一个峰值；当两个信源之间的距离逐渐增大达到分辨率以上时，会出现两个峰值。并且分辨率与两信源的中心轴线存在余弦关系，中心线在0度时分辨率最大。附图说明图1是本专利技术的最大似然角度分辨率的估计流程；图2是为两信源设置在-0.3°和0.3°，间隔为0.6°的后验概率分布图；图3为两信源设置在-0.5°和0.5°，间隔为1°的后验概率分布图；图4是本专利技术的分辨率与信源间距关系图。具体实施方式下面结合附图对本专利技术作更进一步的说明。本专利技术所述系统模型为多传感器阵元的均匀线性阵列，且各阵元辐射都是全向性的，且只针对两个信源探测。阵列天线阵元数为M，各阵元辐射都是全向性的，其接收特性只与其位置相关而与尺寸无关，位置间距为d。远场信源与阵列天线位于同一平面内，信源的反射系数为α，反射系数服从瑞利分布，反射的回波信号到达阵列时视为平面波。同时空间源信号作窄带假设，则接收信号的区别主要体现在到达各阵元的波程差引起的相位差异。阵列中第m个阵元接收信号的时间延时为其中m＝0，1，…，M-1，l＝1...L(本专利技术针对两目标的研究，默认信源数L＝2)，θl表示信源方向，v表示信号的传播速度。忽略空域采样的点数(快拍数N＝1)，则第m阵元上的接收信号为其中：ω0表示载波信号的角频率，xl表示从第l个信源上得到的幅度相位信号，αl为第l个信源的衰减因子，是在[0，2π]中均匀分布的相位变量，wm为第m个阵元上得到的复加性高斯白噪声，各阵元间噪声互不相关，与目标源也不相关，且噪声方差为N0。将所有阵元的接收信号整合为矩阵形式表达式：Y＝A(θ)X+W其中：X为信源矩阵，Y为接收信号矩阵，W为噪声矩阵，阵列的导向矩阵A(θ)表示为其中：τ0(θl)为基准点上阵元接收第l个信源的时间延时，τ1(θl)为基准点旁第一个阵元接收第l个信源的延时，τM-1(θl)为第M-1个阵元接收第l个信源的延时。这里给出第l个信源的信噪比(SignaltoNoiseRatio，SNR)ρl2的定义，E{}为对信源信号求期望。假设反射系数是瑞利分布同时相位是均匀分布的，则信源矩阵X可视为高斯分布的变量。考虑到噪声w矩阵是复高斯变量，因此接收信号矩阵Y也是一个均值为0方差为RY的复高斯向量。RY＝E[YYH]＝A(θ)PAH(θ)+N0U其中：AH(θ)为A(θ)的共轭转置形式，U为单位矩阵，P为信源的自相关矩阵，表示为根据接收到的信号可以确定它的对数似然函数为：其中：y为接收向量，θ为DOA向量。针对一组接收信号，得到此时的似然函数，然后在二维坐标系中的角度搜索范围内，以角度间隔Δθ进行峰值搜索，在搜索过程中，记录每一个DOA角度向量θ对应的似然函数的变化量，将似然函数的变化量最大的DOA角度向量记为θ0，θ0即为相应实际角度的估计值，θ0在二维坐标系中的坐标记为(θ10，θ20)，θ1、θ2为二维坐标系的两条坐标轴，这两个角度值在θ1θ2组成的两维坐标轴上是中心对称出现的。接收多组信号，此时我们可以得每一组接收信号的实际角度的估计值，对这些估计值在二维坐标系中做统计处理，根据每个角度间隔内估计值出现的次数统计出、后验概率密度函数p(θ|y)。当空间中存在两个不相关的信号时，后验概率密度会出现两个峰，此时信源可分辨，可以进行下一步。当空间中的两个信源相关时，后验概率体现为只有一个峰，此时的信源不本文档来自技高网...

【技术保护点】
1.一种传感器阵列系统最大似然角度分辨率估计方法，其特征在于，所述传感器阵列为多传感器阵元的均匀线性阵列，且各阵元辐射都是全向性的；该方法对两个信源进行探测，包括步骤：(1)利用所述传感器阵列天线的接收信号得到观测值的似然函数；(2)根据步骤(1)得到的似然函数得到实际角度的估计值，具体步骤为：获得N组接收信号，对每一组接收信号，得到此时的似然函数，然后在二维坐标系

【技术特征摘要】
1.一种传感器阵列系统最大似然角度分辨率估计方法，其特征在于，所述传感器阵列为多传感器阵元的均匀线性阵列，且各阵元辐射都是全向性的；该方法对两个信源进行探测，包括步骤：(1)利用所述传感器阵列天线的接收信号得到观测值的似然函数；(2)根据步骤(1)得到的似然函数得到实际角度的估计值，具体步骤为：获得N组接收信号，对每一组接收信号，得到此时的似然函数，然后在二维坐标系中的角度搜索范围内，以角度间隔Δθ进行峰值搜索，在搜索过程中，记录每一个DOA角度向量θ对应的似然函数的变化量，将似然函数的变化量最大的DOA角度向量记为θ0，θ0即为相应实际角度的估计值，θ0在二维坐标系中的坐标记为(θ10，θ20)，θ1、θ2为二维坐标系的两条坐标轴；(3)对步骤(2)得到的N组接收信号的实际角度估计值在二维坐标系中的坐标进行统计，根据每个角度间隔Δθ内估计值坐标出现的次数，统计出估计值的后验概率密度函数p(θ|y)，其中，y为接收信号向量；统计后验概率密度函数p(θ|y)的峰，当出...

【专利技术属性】
技术研发人员：闫霄，徐大专，朱思钇，周颖，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32