买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：宁波大学

摘要：本发明公开了一种基于精确模型的二阶相关近场源三维参数欠定估计方法，其基于L型嵌套阵列的近场源三维空间精确传播模型，联合空时信息计算接收数据的自相关和互相关矩阵，构造虚拟接收数据矩阵，通过特征值分解提取噪声子空间，利用信号子空间与噪声子空间的正交性得到空间谱函数，通过一次二维谱峰搜索和一次一维谱峰搜索实现近场源角度和距离参数的欠定估计；优点是该方法能够实现近场源三维参数的欠定估计，具备良好的参数估计性能，无需参数匹配且复杂度低。

主权项：1.一种基于精确模型的二阶相关近场源三维参数欠定估计方法，其特征在于包括以下步骤：步骤1：建立基于L型嵌套阵列的近场源三维空间精确传播模型，所述模型包括K个窄带近场源、L型嵌套阵列，以及基于窄带近场源和L型嵌套阵列建立的数学模型，L型嵌套阵列由位于三维空间的x轴上和y轴上的两个嵌套阵列构成，对应定义为阵列x和阵列y；第k个窄带近场源的三维参数记为αk,βk,rk，其中，k＝1,2,…,K，K≥1，αk表示第k个窄带近场源的入射方向与x轴的夹角，αk∈[0,π]，βk表示第k个窄带近场源的入射方向与y轴的夹角，βk∈[0,π]，rk表示第k个窄带近场源与三维空间的原点的距离，rk∈0,2D2λ]，D表示L型嵌套阵列的孔径，λ为信号波长；阵列x的总阵元个数为N，阵列x为由两个均匀线阵先后嵌套在一起组成的线阵，组成阵列x的第1个均匀线阵的阵元个数和相邻阵元间距对应为Lx和d、第2个均匀线阵的阵元个数和相邻阵元间距对应为N-Lx和Lx+1d，组成阵列x的第1个均匀线阵的第1个阵元和第2个均匀线阵的最后一个阵元对应作为阵列x的第1个阵元和最后一个阵元，组成阵列x的第1个均匀线阵的最后一个阵元与第2个均匀线阵的第1个阵元之间存在线阵间距，阵列x的第n个阵元在三维空间的横坐标为当n∈[1,Lx]时当n∈[Lx+1,N]时阵列x的第n个阵元与第k个窄带近场源的距离为其中，Lx＞1，N-Lx＞1，n＝1,2,…,N；阵列y的总阵元个数为M，阵列y也为由两个均匀线阵先后嵌套在一起组成的线阵，组成阵列y的第1个均匀线阵的阵元个数和相邻阵元间距对应为Ly和d、第2个均匀线阵的阵元个数和相邻阵元间距对应为M-Ly和Ly+1d，组成阵列y的第1个均匀线阵的第1个阵元和第2个均匀线阵的最后一个阵元对应作为阵列y的第1个阵元和最后一个阵元，组成阵列y的第1个均匀线阵的最后一个阵元与第2个均匀线阵的第1个阵元之间存在线阵间距，阵列y的第m个阵元在三维空间的纵坐标为当m∈[1,Ly]时当m∈[Ly+1,M]时阵列y的第m个阵元与第k个窄带近场源的距离为其中，Ly＞1，M-Ly＞1，m＝1,2,…,M；K个窄带近场源入射到L型嵌套阵列中，将基于窄带近场源和L型嵌套阵列的数学模型表示为：其中，t表示时间变量，xt表示阵列x的接收数据矢量，α表示窄带近场源的入射方向与x轴的夹角，r表示窄带近场源与三维空间的原点的距离，α和r为变量，Axα,r表示阵列x的流形矩阵，Axα,r的维度为N×K，Axα,r＝[axα1,r1,axα2,r2,…,axαK,rK]，axαk,rk为导向矢量，axαk,rk的维度为N×1，表示阵列x的第n个阵元对第k个窄带近场源的空间幅度相位因子，g表示已知的路径损耗衰减指数，e为自然常数，j为虚数单位，λ为信号波长，上标“T”表示矩阵或向量的转置，st表示信号矢量，st的维度为K×1，st＝[s1t,s2t,…,sKt]T，skt表示第k个窄带近场源入射的信号，nxt表示xt中的高斯白噪声矢量，nxt的维度为N×1，yt表示阵列y的接收数据矢量，β表示窄带近场源的入射方向与y轴的夹角，β为变量，Ayβ,r表示阵列y的流形矩阵，Ayβ,r的维度为M×K，Ayβ,r＝[ayβ1,r1,ayβ2,r2,…,ayβK,rK]，ayβk,rk为导向矢量，ayβk,rk的维度为M×1，表示阵列y的第m个阵元对第k个窄带近场源的空间幅度相位因子，nyt表示yt中的高斯白噪声矢量，nyt的维度为M×1；步骤2：在步骤1的基础上，计算阵列x的接收数据矢量xt的延时自相关矩阵Rxxτ，然后对Rxxτ进行矢量化计算，得到Rxxτ的矢量化结果rxxτ， 再取τ＝Ts,2Ts,…,NpTs，得到由阵列x的接收数据矢量xt构造出的虚拟接收数据矩阵Wx，Wx＝[rxxTs,rxx2Ts,…,rxxNpTs]＝BxΡs；其中，τ表示延时变量，E{·}表示数学期望，·H为共轭转置运算，xt+τ表示xt延时τ后的结果，Ax即为Axα,r，Rsτ表示信号矢量st的协方差矩阵，Rsτ＝E{stsHt+τ}，st+τ表示st延时τ后的结果，δ·为狄利克雷函数，RNτ表示xt中的高斯白噪声矢量nxt的协方差矩阵，nxt+τ表示nxt延时τ后的结果，在高斯白噪声假设下当τ＝0时表示噪声功率，IN表示N阶单位阵，当τ≠0时RNτ中的所有元素均为0，vec{·}表示矢量化，⊙为Khatri-Rao积运算，上标“*”表示复共轭，ρsτ＝diag{Rsτ}，diag{·}表示取主对角线元素，Bx表示阵列x的虚拟流形矩阵，Bx的第k列向量为bxαk,rk，bxαk,rk为导向矢量，为Kronecker积运算，Ts表示虚拟采样间隔，Np表示伪快拍数，Ρs表示虚拟信号矩阵，Ρs＝[ρsTs,ρs2Ts,…,ρsNpTs]；步骤3：在步骤2的基础上，计算Wx的协方差矩阵然后对进行特征值分解，并取最小的N2-K个特征值对应的特征向量构成噪声子空间导向矢量bxα,r属于信号子空间，利用信号子空间与噪声子空间UN1的正交性，并根据α的范围α∈[0,π]和r的范围r∈0,2D2λ]，计算关于α和r的空间谱函数Pα,r，再通过对Pα,r进行二维谱峰搜索，得到K个峰值，第k个峰值所在位置对应的两个参数即为αk和rk各自的估计值，对应为和其中，表示Ρs的协方差矩阵，axα,r为axαk,rk中的参数替换成变量后的形式；步骤4：在步骤1的基础上，计算阵列x的接收数据矢量xt与阵列y的接收数据矢量yt的延时互相关矩阵Rxyτ，然后对Rxyτ进行矢量化计算，得到Rxyτ的矢量化结果rxyτ，再取τ＝Ts,2Ts,…,NpTs，得到由阵列x的接收数据矢量xt和阵列y的接收数据矢量yt构造出的虚拟接收数据矩阵Wxy，Wxy＝[rxyTs,rxy2Ts,…,rxyNpTs]＝BΡs；其中，yt+τ表示yt延时τ后的结果，Ay即为Ayβ,r，B表示阵列x和阵列y的虚拟流形矩阵，B的第k列向量为bαk,βk,rk，步骤5：在步骤4的基础上，计算Wxy的协方差矩阵然后对进行特征值分解，并取最小的N×M-K个特征值对应的特征向量构成噪声子空间导向矢量bα,β,r属于信号子空间，利用信号子空间和噪声子空间的正交性，并根据β的范围β∈[0,π]，同时将和代入关于β的空间谱函数中，再通过对Pβ进行一维谱峰搜索，得到一个峰值，这个峰值所在位置对应的参数即为βk的估计值其中，ayβ,r为ayβk,rk中的参数替换成变量后的形式，axα,r为axαk,rk中的参数替换成变量后的形式，为将和代入bα,β,r后的形式。

全文数据：

权利要求：

百度查询： 宁波大学 基于精确模型的二阶相关近场源三维参数欠定估计方法