亚星游戏官网-www.yaxin868.com

专利名称：一种对传感器天线阵列接收信号的快速处理方法
技术领域：
本发明属于利用传感器天线阵列对无线传播信号的来波方向进行测定的领域，特别是ー种利用天线阵列的旋转不变性(ESPRIT)对单个或多个无线传播信号的来波方向进行快速測定(定向)的方法。采用该方法可在无需为了得到旋转不变关系矩阵进行多次特征值分解的情况下进行快速測定(定向)、并可获得较好的效果(精确度)，以及降低測定成本等优点。
背景技术：
接收无线传播信号的天线阵列信号处理技术已广泛应用。阵列信号处理是将ー组传感器在空间的不同位置按照一定的规则布置，形成传感器阵列，通过多个通道接收，利用一定的空域处理或者空时等多维联合处理的方式，得到所需的信号或信号參数。与传统的单个定向传感器相比，阵列信号处理具有更灵活的波束控制，更高的信号増益和空间超分辨能力，与此相关的研究工作迅速发展，应用范围也不断扩大，现已成功应用于地震、天文、雷达、声纳、医学信号检测、无线通信、地质勘探等领域。信号的来波方向(DOA)是阵列信号处理技术所需要的ー个十分重要的參数，是阵列信号处理领域研究的热点方向。现代超分辨的信号来波方向估计(測定)的基本原理是通过将阵列多个阵元在同一时刻的采样信号排列为一向量，计算多个快拍下的阵列接收信号的自相关矩阵，并建立相应的信号子空间和噪声子空间；利用信号子空间的特性进行来波方向角的估计(測定)。多重信号分类法(MUSIC算法)是超分辨测向技术的代表之一，它利用信号子空间和噪声子空间的正交性原理构造ー个伪谱函数，并在一定的方位内进行谱峰捜索，虽然多重信号分类法能得到很高的角度分辨率，但是其谱峰搜索尤其是在高精度的ニ维搜索时需要的庞大的计算量阻碍了其应用。而以旋转不变技术(ESPRIT)为代表的方法利用局部阵列的平移不变性，无需计算量庞大的谱峰捜索就能得到信号来波方向角的超分辨估计(測定)，其估计(測定)的精度与MUSIC (多重信号分类法)相差无几。因此利用阵列的旋转不变性来估计信号的来波方向角具有重大的研究意义。在ESPRIT处理方法的具体实现方面，典型的有基于不同准则的ESPRIT算法和基于ESPRIT的推广算法。其中，基于最小ニ乘准则的LS-ESPRIT算法和基于总体最小ニ乘准则的TLS-ESPRIT法是两种基于不同准则的ESPRIT算法。LS-ESPRIT算法是用ー个范数平方为最小的扰动去干扰信号子空间，校正信号子空间中的噪声的方法。TLS-ESPRIT算法是同时扰动两个子阵的信号子空间，并使扰动范数的平方保持最�。�同时校正两个子阵的信号子空间中存在的噪声。这两种算法都需要对由两个子阵接收信号组成的组信号向量的协方差矩阵进行特征值分解，亦即要对2倍于子阵的协方差矩阵阶数的矩阵进行特征值分解，因而计算量庞大；而且，在二者构造旋转不变关系矩阵的过程中都需要再次对计算过程中构造的矩阵进行特征值分解；而后者在计算中还需要对一般矩阵进行求逆的计算，这样导致了 ESPRIT算法实施需要的计算量与需要谱峰搜索的MUSIC算法比起来也相差无几；而矩阵束ESPRIT算法则是利用两个子阵接收数据的自相关矩阵和互相关矩阵进行广义特征分解以确定来波方向的方法，该方法虽然可以有效的降低ESPRIT算法的运算量，但是对来波方向測量精度难以满足常规应用的要求。

发明内容
本发明的目的是针对背景技术存在的缺陷，研究设计ー种对传感器天线阵列接收信号的快速处理方法，以克服背景技术存在的缺陷，达到在降低对设备和硬件性能要求的同吋，对目标来波方向进行快速、准确测定，提高精度，简化流程、降低运行成本等目的。本发明的解决方案是首先将均匀直线阵列按照图I所示的划分方法分为两个子阵，将第一个子阵得到的采样信号矩阵与自身进行相关运算得到子阵I的自相关矩阵，再将第一个子阵与第二个子阵得到的信号矩阵进行互相关运算，得到子阵I与子阵2的互相关矩阵；将子阵I的自相关矩阵进行特征值分解得到子阵I的信号子空间矩阵，利用信号子空间矩阵的共轭转置矩阵与子阵I的自相关矩阵相乗、然后将其结果再与子阵I的信号子 空间矩阵相乘，得到ー个对角矩阵，将这个对角矩阵的各对角元均减去子阵I的小特征值的平均值，得到ー个信号特征矩阵，再将所得信号特征矩阵中的各对角元取倒数，建立信号特征矩阵的逆矩阵；将子阵I的信号子空间矩阵的共轭转置矩阵与子阵I与子阵2的互相关矩阵相乘、其结果在再乘以子阵I的信号子空间矩阵，得到一个子阵I与子阵2互相关矩阵的过渡矩阵；该过渡矩阵与上述信号特征矩阵的逆矩阵相乘、即得到旋转不变关系矩阵；最后对旋转不变关系矩阵进行特征值分解、分解为ー组特征值和与该组特征值对应的特征向量，其中所得各特征值即为不同信号来波方向的角度參数，最后经过反正弦运算，即得到不同信号的来波方向角。因而，本发明方法包括步骤A.划分子阵及建立各子阵的信号矩阵将均匀线性阵列(ULA)传感器天线的前M-I个阵元组合成子阵1，再将第2个阵元至第M个阵元组合成子阵2 ;将子阵I和子阵2中各阵元的接收信号分别建立接收信号矩阵、即首先将第一时刻两子阵中各阵元得到的信号分别依次排成一列、作为对应子阵中的第一个列向量，再将下ー时刻两子阵中各阵元得到的信号分别依次排列组成第二列、得到第二个列向量，以此类推按时间先后顺序将各阵元接收到的信号列向量依次排列，分别得到子阵I和子阵2的接收信号矩阵；步骤B.建立自相关矩阵和互相关矩阵根据步骤A所得子阵I和子阵2两个子阵的接收信号矩阵，通过R = Χ.Λ7ν建立自相关矩阵；通过R = XiX211建立互相关矩阵；上式中Art为子阵I接收信号的自相关矩阵，Rvi为子阵I与子阵2接受信号的
互相关矩阵；Xi、X2分别为子阵I、子阵2的接收信号矩阵，(·)H表示矩阵的共轭转置；步骤C.建立信号子空间矩阵对自相关矩阵天通过
R =υΣΙ”!进行特征值分解；其中U为以M-I个正交的特征向量构成的酉矩阵、U=Lu1, U2, ...，u.J，U1, U2,...，Ush为子阵I接收信号的自相关矩阵及Λ的M-I个特征值对应的特征向量；Σ为以M-I个特征值为对角元的对角矩阵、Σ =diag (A1, λ2,…λ ^1), A1, λ2,…λ η为' 的M-I个特征值；(·)Η表示矩阵的共轭转置；以不同信号的个数K为基准，将自相关矩阵h的特征值按照从大到小减值排列，
利用排列后的前K个较大特征值对应的特征向量建立子阵I的信号子空间矩阵；其信号子空间矩阵为Us= [udl, Ud2，…Udk]
将后M-K-I个小特征值取平均值，得到λ E待用；上式中udl，ud2,…Udk依次为前K个大特征值对应的特征向量；步骤D.建立信号特征矩阵及其逆矩阵把子阵I接收信号的自相关矩阵Λνι两边分别乘以步骤C所得的信号子空间矩阵Us的共轭转置矩阵和信号子空间矩阵us，其结果再减去由步骤C提取待用的小特征值平均值λ E构成的KXK阶的对角矩阵；得到ー个KXK阶的信号特征矩阵B=diag ( λ dl_ λ Ε, λ d2_ λ Ε,…λ 亚―λ Ε)其中λdl，λ d2,…λ dK依次为自相关矩阵特征值按照从大到小减值排列的前K个大特征值；建立上述信号特征矩阵的逆矩阵即取该对角矩阵中各对角元的倒数，并以各对角元的倒数建立信号特征矩阵的逆矩阵、逆矩阵为B 1=diag[l/ ( λ dl- λ E)，I/ ( λ d2- λ Ε)，…，I/ ( λ 狀-λ Ε)]步骤Ε·建立过渡矩阵首先将由步骤C中提取待用的小特征值平均值λ Ε乘以(M-I) X (M-I)阶的下三角矩阵Z得到加权矩阵λ Εζ
O O ヽ其中篕=10 .·
OI Oノ利用子阵I与子阵2接收信号矩阵所得的互相关矩阵足112减去上述所得加权矩阵λ ΕΖ ;其结果再两边分别乘以步骤C中所得的信号子空间矩阵Us的共轭转置矩阵和信号子空间矩阵Us，得到过渡矩阵'B，即B^=UZ(Rxix2-XeZ)Us步骤F.建立旋转不变关系矩阵
将步骤E中得到的矩阵A 左乘步骤D中得到的信号特征矩阵的逆矩阵B'即得子阵I的旋转不变关系矩阵Ψ，即ψ=η ·.步骤G.获得不同信号的来波方向角对步骤F中所得的旋转不变关系矩阵Ψ 进行特征值分解，得到所需的K个特征值及与各个特征值对应的特征向量，所得特征值为2…、IK ;将各个特征值经过自然对数和反正弦运算得到对应信号的来波方向角，即:β ニ arcsin(-~—In ζ)
k2-J其中Λ为载波波长，d为阵元间距(一般取为半波长)。本发明采用直接建立旋转不变关系矩阵的方法，只利用矩阵的乘法和加法处理构造阵列旋转不变关系的方法有利于计算机处理、并可通过FPGA实现传感器天线阵列对单个或多个无线传播信号的来波方向进行快速測定(定向)。与现有技术比较，不但处理速度快，在相同的信噪比的情况下来波方向角的測定精度得到显提高，而且本发明方法在满足平移不变性的ニ维阵列中还可运用于测定空间来波方向角；因而本发明具有对设备和硬件性能要求低，对目标来波方向的測定快速、准确，精度高，流程短、运行成本低，应用范围广等特点。


图I为本发明均匀线阵子阵划分和接收信号模型示意图；图2为本发明具体实施方式
流程图。
具体实施例方式本实施方式中无线传播信号在空间高斯白噪声的背景下、各个(来波)信号之间是互不相关的，天线阵列为均匀直线阵列，每个阵元均是全向天线，阵元个数为M=16，阵元间的间隔为d=0. 083米；本实施方式中信号源为来自K=3个不同方向的远场窄带信号、且都为相同频率的载波，其波长为Λ=0. 166米，且3个信号源分别以角度θ1=20°、θ2=30°、θ3=40°入射到传感器均匀直线阵列；相邻两阵元中对同一信号源、其后一个阵元相比于前一个阵元的相位延迟为Φ =2 31 d sin Θ j//Λ = ji sin Θ j ;采样快拍数为N=1024,信噪比SNR=-15dB。步骤A.划分子阵及建立各子阵的信号矩阵前15个阵元组合成子阵1，第2个到第16个阵元组合成子阵2 ;则两个子阵的接收数据矩阵为I1=K^E1 ；X2=A2S+E2=A1 Φ S+E2 ；其中X1为子阵I的接收数据矩阵，X2为子阵2的接收数据矩阵；ん为子阵I的方向矩阵、A2为子阵2的方向矩阵；S为3X1024维的源信号矩阵；Φ为3X3维的旋转不变关系矩阵，E1为子阵I的高斯白噪声矩阵，E2为子阵2的高斯白噪声矩阵；本实施例信号与空间噪声不相关；步骤B.建立自相关矩阵和互相关矩阵根据两个子阵的接收信号矩阵X1和X2，建立自相关矩阵#和和互相关矩阵Λ、.λ，用阵列的导向矢量表示为如下的形式R = A1PAf + σ2Ι= A2PAJ1 +σ Z= ΑΙΡΦΗΑ1Η + σこ Z其中上标(·)Η表示矩阵的共轭转置，P=SSh为信号自相关矩阵，σ 2I和σ 2Z均为15X15阶的方阵、分别为两个子阵接收噪声矩阵的自相关矩阵和互相关矩阵；步骤C.建立信号子空间矩阵对自相关Rvi进行特征值分解，得到以特征值为对角元的对角矩阵Σ =CliagO1, λ 9) ··· λ 15)，以及咯特征向量组成的酉矩阵U=Iiu1, U2，…U15]；将自相关^^的特征值按照降序排列，将前3个较大特征值入”入”入“寸应的特征向量Up u2、U3构成信号子空间，即信号子空间矩阵Us=Iiu1, U2, U3I ；将后12个较小特征值取平均值λ Ε待用；步骤D.建立信号特征矩阵及其逆矩阵把自相关两边分别乘以信号子空间矩阵Us的共轭转置矩阵和信号子空间矩阵队，再减去由自相关矩阵A1中较小特征值的平均值λ “勾成的3X3阶的对角矩阵λΕΙ ；
得到ー个3X3阶的信号特征矩阵B=diag ( λ「λ Ε, λ 2- λ E, λ 3- λ E)然后建立该信号特征矩阵的逆矩阵即取该对角矩阵中各对角元的倒数，并以各对角元的倒数建立信号特征矩阵的逆矩阵、其逆矩阵为B 1=diag[l/ ( λ「λ E)，I/ ( λ 2- λ E)，I/ ( λ 3- λ E)]步骤E.建立过渡矩阵将互相关《_2减去由步骤C中提取待用的较小特征值的均值λ Ε乘以3 X 3阶的下三角矩阵Z得到的加权矩阵λΕΖ，得到的结果再两边分别乘以信号子空间矩阵Us的共轭转置矩阵和信号子空间矩阵Us，得到过渡矩阵B
权利要求
1.一种对传感器天线阵列接收信号的快速处理方法，包括 步骤A.划分子阵及建立各子阵的信号矩阵 将均匀线性阵列传感器天线的前M-I个阵元组合成子阵1，再将第2个阵元至第M个阵元组合成子阵2 ;将子阵I和子阵2中各阵元的接收信号分别建立接收信号矩阵、即首先将第一时刻两子阵中各阵元得到的信号分别依次排成一列、作为对应子阵中的第一个列向量，再将下一时刻两子阵中各阵元得到的信号分别依次排列组成第二列、得到第二个列向量，以此类推按时间先后顺序将各阵元接收到的信号列向量依次排列，分别得到子阵I和子阵2的接收信号矩阵； 步骤B.建立自相关矩阵和互相关矩阵 根据步骤A所得子阵I和子阵2两个子阵的接收信号矩阵，通过
全文摘要
该发明属于对传感器天线阵列接收信号的快速处理方法。包括划分子阵及建立各子阵的信号矩阵，建立自相关矩阵和互相关矩阵，建立信号子空间矩阵，建立信号特征矩阵及其逆矩阵，建立过渡矩阵，建立旋转不变关系矩阵，最后获得不同信号的来波方向角。该发明采用直接建立旋转不变关系矩阵的方法，只利用矩阵的乘法和加法处理构造阵列旋转不变关系的方法实现传感器天线阵列对单个或多个无线传播信号的来波方向进行快速测定；与背景技术相比不但处理速度快，特别是在信噪比较低时对来波方向角的测定精度显提高，而具有对设备和硬件性能要求低，对目标来波方向的测定快速、准确，精度高，流程短、运行成本低，应用范围广等特点。
文档编号G01S5/02GK102662158SQ20121013717
公开日2012年9月12日 申请日期2012年5月4日 优先权日2012年5月4日
发明者张伟, 杨陶柳, 贾海涛, 赵秀粉 申请人:电子科技大学