CN103424735B - 基于最小描述长度的近场源定位方法、装置及系统 - Google Patents
基于最小描述长度的近场源定位方法、装置及系统 Download PDFInfo
- Publication number
- CN103424735B CN103424735B CN201310326242.5A CN201310326242A CN103424735B CN 103424735 B CN103424735 B CN 103424735B CN 201310326242 A CN201310326242 A CN 201310326242A CN 103424735 B CN103424735 B CN 103424735B
- Authority
- CN
- China
- Prior art keywords
- characteristic vector
- source
- matrix
- signal
- module
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/74—Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于最小描述长度的近场源定位方法、装置及系统,涉及移动业务中的信源定位系统技术领域。该方法包含步骤:S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;S2、基于最小描述长度准则,对协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;S4、将步骤S2中的信源数目K与步骤S3中的相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;S5、对所述相位信息进行反解,得到信源方位。本发明既能降低近场源定位的计算复杂度,提高计算效率;又能充分发挥数据预处理的效能,提高定位的精准性与稳定度。
Description
技术领域
本发明涉及移动业务中的信源定位系统技术领域,具体涉及一种基于最小描述长度的近场源定位方法、装置及系统。
背景技术
无线电被动定位技术是无线电监测与管理中的重要技术手段,在雷达、声纳、炮兵电声测量等军事领域,石油地震勘测、地下管道泄漏检测定位和故障检测等工业领域,以及移动定位业务、无线电检测等通信领域都有着广泛的应用。
传统的无线电发射源定位是采用转动接收天线(传感器)的角度来实现的。但这种方法存在测量精度与测量速度的瓶颈,难以应对空间存在多个目标的情况。而通过采用阵列信号处理技术,合理选取阵元数目和阵元间距可以较好地进行多目标方位探测。阵列即按一定方式布置在空间上不同位置的传感器组,利用信号的空域特性增强信号并且有效提取信号的空域信息。因此阵列信号处理也常称为空域信号处理。与传统的单个定向传感器相比,阵列信号处理具有灵活的波束控制、较高的信号增益、极强的抗干扰能力及空间超分辨能力等优点,相关的研究工作不断发展与深入,其应用范围也不断扩大。并且随着微电子、数字信号处理、并行处理等技术的迅猛发展,阵列信号处理技术也得到较快的发展。阵列信号处理最主要的研究内容包括DOA估计和波束形成。较早的DOA估计方法又称为波束形成方法,而该波束形成方法利用了空域维纳滤波的匹配概念,由阵列流形在信号空间中的投影大小判定信号方向,后来随着研究的深入,出现了多重信号分类(MUSIC)基于接收信号矩阵分解的谱估计方法,该方法利用信号空间与噪声空间的正交性,先对接收数据的协方差矩阵进行特征分解,再对方向向量相关性进行搜索。将相关性的极小值点对应的信源参数作为信源位置的最大似然估计,从而实现了对多目标信源的精确定位。
然而,以上多目标信源定位的方法,如MUSIC具有较高的计算复杂度。多重信号分类在进行谱峰搜索的过程中需要对信源方位角的候选区间进行全面搜索,对于需要多维空间表征的信号参数,搜索复杂度是指数级增加的。因此在实际应用中,该方法在计算效率上存在严重的缺陷,极大限制了定位系统的时效性。
发明内容
(一)解决的技术问题
针对现有技术的不足,本发明提供一种基于最小描述长度的近场源定位方法、装置及系统,使得近场源定位的计算复杂度降低。
(二)技术方案
为实现以上目的,本发明通过以下技术方案予以实现:
一种基于最小描述长度的近场源定位方法,包含以下步骤:
S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
S2、基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
S4、将步骤S2中的所述信源数目K与步骤S3中的所述两组相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
S5、对所述相位信息进行反解,得到信源方位。
优选的,步骤S1中计算协方差矩阵特征向量表达式为:
式中,U为接收数据协方差矩阵的特征向量矩阵,∑为由特征值组成的对角阵,ei与λi分别是特征向量及与之对应的特征值,M是特征向量数目。
优选的,步骤S2中计算信源数目K的表达式为:
式中,k为模型参数,l1>l2>…>lM是样本协方差矩阵,M是特征向量数目。
优选的,步骤S4包括步骤:
S41、对所述两组相关性矩阵与特征向量求内积,得到两组任意两个特征向量间的内积值;
S42、对所述内积值进行排序,取出所述信源数目K对较大相关性的特征向量对;
S43、通过所述K对特征向量及其所对应的特征值与信号相位关联关系,得到菲涅尔传播模型中的相位参数。
优选的,步骤S43中解算接收信号的相位信息的表达式为:
式中,N是快拍数,u1(i+1)是对应的第1个特征向量中的第i+1个分量,L是估计结果K,与为菲涅尔模型中的相位参数。
优选的,步骤S5中反解信源方位的表达式为:
式中,d为天线间距,为第1个信源的方向角。
本发明提供了一种基于最小描述长度的近场源定位装置,包含以下部分:
协方差矩阵模块,对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
信源数目预估计模块,基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
相关性矩阵特征分解模块,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
特征向量匹配模块,将信源数目预估计模块计算出的的信源数目K与相关性矩阵特征分解模块计算出的相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
反解信源方位模块,对相位信息进行反解,得到信源方位。
本发明还提供了一种基于最小描述长度的近场源定位系统,包含数据接收模块、数据采集模块,还包含权利要求6所述的近场源定位装置;
其中,
数据接收模块包含天线阵列、滤波器、放大器及A/D转换器,将接收到的信号发送给数据采集模块;
数据采集模块包含多路数据电路终接设备和多路传输信道,将数据接收模块发送的信号传输给近场源定位装置。
优选的,该系统近一步包括监测数据库模块和位置数据库模块;
所述监测数据库模块,对近场源定位装置的测试数据进行存储;
所述位置数据库模块,对近场源定位装置的位置数据进行存储。
(三)有益效果
本发明通过提供一种基于最小描述长度的近场源定位方法、装置及系统,通过在信号方向向量匹配前进行信源数目预估计,进而指导方向向量匹配按照估计的参量个数进行,能够把近场源定位的计算复杂度降低,提高了计算效率。
本发明又能充分发挥信源数目预估计的效能,提高定位的精准性与稳定度,实现了系统定位性能的整体提升。
较之传统高阶累积量方法如计算四阶相关性矩阵的平行因子方法,与进行参数空间谱峰搜索的方法如MUSIC方法相比较,本发明实现了近场源定位效率与定位准确性的双赢。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的一种基于最小描述长度的近场源定位方法的流程图;
图2为本发明实施例的线性阵列示意图;
图3为本发明实施例的一种基于最小描述长度的近场源定位装置的结构示意图;
图4为本发明实施例的一种基于最小描述长度的近场源定位系统的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1:
如图1所示,一种基于最小描述长度的近场源定位方法,包含以下步骤:
S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
S2、基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
S4、将步骤S2中的所述信源数目K与步骤S3中的所述相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
S5、对所述相位信息进行反解,得到信源方位。
本发明实施例通过提供通过在信号方向向量匹配前进行信源数目预估计,进而指导方向向量匹配按照估计的参量个数进行,能够把近场源定位的计算复杂度降低,提高了计算效率。
下面对本发明实施例进行详细的说明:
一种基于最小描述长度的近场源定位方法,包含以下步骤:
S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
假设多个窄带信号入射到空间某一均匀线性阵列上,并且所有的阵元是理想的各向同性,不存在通道不一致、互耦等因素的影响。这里的这阵元间距为d,这里的阵元经各自的传输通道传送到对应的处理器。
由于窄带信号的包络变化缓慢,等距线阵各阵元接收到的同一信号的包络相同。对于远场信号,即信源与阵元相距足够远时,电波到达各阵元的波前为平面波,且到达各阵元的方向角θi相同,称为波达方向,定义为信号si(n)到达阵元的直射线与阵列法线方向之间的夹角。这里以阵元1为参考阵元,则空间信号到达其他阵元的时间相对于参考阵元存在延迟(或超前)。令信号si(n)电磁波传播延迟在第二个阵元引起的相位差为ωi,则波达方向θi与相位差之间ωi之间存在如下关系:
式中,d是两个相邻阵元之间的距离,λ为信号波长,阵元间距d满足d≤λ/2,否则相位差ωi可能大于π,进而产生方向模糊,即θi与π+θi都可能是si(n)的波达方向。又由于该阵列在空间上的排布为等距线阵,信号si(n)到达第k个阵元的电波与到达参考阵元的电波之间的相位差为:
因此,信号si(n)在第k个阵元上的接收信号为:
若阵列由M个阵元组成,则信号si(n)到达各阵元的方向向量可表示为:
如果共有p个信号位于远场(p<M),则在第k个阵元上的观测或接收信号xk(n)为:
式中,ek(n)表示第k个阵元上的加性观测噪声,将M个阵元上的观测数据组成M×1维观测数据向量:
x(n)=[x1(n),x2(n),..,xM(n)]T(4)
类似可以定义M×1维噪声向量:
e(n)=[e1(n),e2(n),..,eM(n)]T(5)
这样,阵元接收数据可以写成如下向量形式:
式中,
分别为M×p维方向向量矩阵和N×1维信号向量。具有上述形式范德蒙德矩阵的特点是:若ω1,ω2,..,ωp互不相同,则矩阵的各列相互独立,也即矩阵是列满秩矩阵。
由阵元接收数据的数学模型,阵列的协方差矩阵
Rxx=E[x(n)x(n)H]=AE[SSH]AH=ARssAH+RN(8)
式中,RN为噪声协方差矩阵,对于空间的理想高斯白噪声,且噪声方差为则有下式
进而对Rxx进行特征分解,得到步骤S1中计算协方差矩阵特征向量表达式为:
其中,U为接收数据协方差矩阵的特征向量矩阵,∑为由特征值组成的对角阵,即∑=diag(λ1,λ2,..,λM);ei与λi分别是特征向量及与之对应的特征值,M是特征向量数目。
S2、基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
下面首先介绍下最小描述长度准则:
最小描述长度源于柯氏复杂度。1978年,Rissanen在文献中首先提出了最小描述长度准则(MDL)的概念。MDL准则起源于抽象出Kolmogorov负责度衡量的本质思想,以衡量实际时间的复杂度。信息论相关定理使得MDL准则更加实用和完善。
MDL准则的基本思想是建立可以精确描述对象的数学模型,并在低的描述复杂性和模型精确性之间取得最好的折中。在MDL中,数学模型的表现形式成为描述语言。描述长度越大,模型越复杂;描述长度越短,模型越简单。下面对MDL的思想进行简要介绍。
设模型集合为M,MDL准则选择出来的模型为Mmdl,模型Mmdl的选取标准是最小化一下两项(描述长度)之和:
|Lm(Mi)|:描述模型Mi所需要的位数(可理解为将所有参数编码所需要的编码总和,也称为描述长度)。这里,Mi∈M,Lm(Mi)为描述模型Mi的语言,|*|表述位数(即描述长度)。
|Lc(D|Mi)|:给定模型Mi,描述对象D所需的位数。这里|Lc(C|Mi)|为基于模型Mi描述对象D的语言。给定模型Mi描述对象D,等价于描述数据和模型Mi之间的误差;
即:
任何一个观测序列(即对象),可以看成是由两个序列组成,一个是确定性序列,一个是“纯”随机过程的序列(即观测序列与确定性序列的残差序列)。对应于MDL准则,Mi刻画的是确定性序列;模型Mi与观测对象的误差是随机过程序列。对于确定性模型的部分,描述性语言Lc(D|Mi)可以选取如AR模型,多项式模型,插值模型等的描述性模型。此时描述性语言的长度Lm(Mi)为模型中参数的个数。对于不确定性部分,描述性语言Lc(D|Mi)可以选取概率分布模型。信息论专门研究随机性事件。其中的Shannon编码定理表明,用概率分布模型描述随机过程时,其描述长度等于对象概率分布以2为底对数的相反数。已经证明,在大样本数据情况下,最小描述长度准则可有下面的指标度量:
上式中k是模型独立位置参数的个数,N是观测数据序列的长度,是参数向量θ的最大似然估计,第一项是模型参数最大似然估计的对数。
现假设观测向量为X=[x(t1),x(t2),...,x(tN)],且是零均值统计独立的高斯随机向量,则这里的模型族可以被描述为观测数据的协方差矩阵,即
R(k)=E[XXH]=Ψ(k)+σ2I(13)
这里,k∈{0,1,2,...,M-1}取遍所有可能的信号源数目集合,Ψ(k)(信号向量协方差矩阵)是一个秩为k的半正定矩阵,σ2(噪声功率)为一个未知量。
利用线性代数中的谱表示理论,R(k)可作如下分解:
式中,λ1,λ2,...,λk和V1,V2,...,Vk分别是矩阵R(k)的特征值与特征向量。下面用θ(k)表示模型参数向量,即
参量化后,现在来引出该检测模型的信息论指标。由于假设观测数据是零均值的复高斯随机向量,则它们的联合概率密度有如下形式:
对上式取负对数,并且忽略与参数向量θ(k)的无关项,可得如下对数似然函数
式中,是样本协方差矩阵,即
最大似然估计即是求取特定θ(k)使公式(5-10)中的似然函数L(θ(k))取得最大。文献[26]中Anderson给出了该问题的闭式解:
这里,l1>l2>…>lM是样本协方差矩阵的特征值,且C1,C2,...,CM为对应的特征向量。将上述最大似然解代入式似然函数,简化可得:
可以看出,上式中括号内的分式就是M-k个小特征值的几何均值与算术均值的比值。
这里,可以通过计算θ(k)张成的参数空间的自由度得到θ(k)中自由参数的个数。由于复协方差矩阵的特征值为实数,但其特征向量为复向量,所以θ(k)中实际含有k+1+2Mk个参数。然而这k+1+2Mk个参数的取值是有限制的:特征向量本身相互正交,并且假定特征向量为单位向量(已经过单位化)。这等价于参数空间的自由度因标准化降低了2k,因正交性降低了则模型自由参数个数为:
进而得到步骤S2中计算信源数目K的表达式为:
其中,k为模型参数(控制变量),l1>l2>…>lM是样本协方差矩阵,M是特征向量数目。
S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
定义基于阵元滞后量的二阶统计量矩阵R1。该矩阵为N阶方阵,其m行n列元素定义如下:
式中,τ1与υ1共同构成了阵元间距,即阵元滞后量。下面记m-n+τ1为Δ(τ1),记m-n-υ1为Δ(υ1),则当信源与阵元噪声之间两两不相关时
这里,为第l个信源的信号功率。类似矩阵R1,交换阵元滞后量主键,可定义矩阵R2的m行n列元素如下
记n-m+τ2为Δ(τ2),记n-m-υ2为Δ(υ2),则R2(m,n)可表示如下
现令R1中的τ1+υ1=1,R2中的τ2+υ2=-1,则可将其分别转化为矩阵形式。令 信号方向向量方向矩阵A=[a1,a2,...,aL],为范德蒙德矩阵。则R1与R2有如下形式
显见上述统计量是从空间阵元排布的角度考察不同阵元滞后量下阵元接收数据的相关性。下面再从时域的角度,考虑不同时延的相关统计量。由于发射信号为窄带随机过程,所以k时刻的信号s(k)与k+1时刻的信号s(k+1)有相角存在一个时间间隙的差量。由此,可定义如下基于一个时间间隙的二阶统计量:
假设L个窄带信号的频率分别为ω1,ω2,...,ωL,令Φ=diaθ{ω1,ω2,...,ωL},则经类似R1、R2的化简过程进行处理后,R3(m,n)有如下形式
将其转化为矩阵形式,可得
基于上述阵元滞后量的二阶统计特性分析,考虑单个阵元延迟的相关性矩阵。已有约束条件:τ1+υ1=1;τ2+υ2=-1;τ3+υ3=1即可保证延迟量为1,也就是考察单间隔阵元的二阶统计特性。在此基础上,再对延迟量取值做如下限制:
结合两组限制条件,可得
代入R1、R2、R3,得
由假设所有信源都具有非零功率,可知对角矩阵Γ对角线元素非零,且为可逆矩阵。此外,由于信源参量及各不相同,因此矩阵A各列向量互不相关,为满秩阵。同时,矩阵R1,R2,R3是秩L的N阶方阵。
对矩阵R1进行特征值分解,有
上式中,U=[u1,u2,...,uL]为R1的L个特征向量构成的矩阵,V=diag{v1,v2,...,vL}为R1的L个特征值构成的对角矩阵。定义矩阵R1的伪逆:
由式R1、R2、R3的形式,有
由上式可见,A中各列向量同为与的特征向量。并且由等式右边的对角矩阵可知矩阵与的公共特征向量ai对应的特征值为和由于这两个特征值是对应于同一个特征向量的,所以特征向量确定后即可将和与进行配对。当然,实际中由于噪声的存在对两个矩阵进行分解后得到的特征向量并不相等。这时考虑采用求取复向量内积的方式提取向量的相关度量,将相关性最大的向量视为“相等”的向量进行配对。
S4、将步骤S2中的所述信源数目K与步骤S3中的所述相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
优选的,步骤S4包括步骤:
S41、对所述两组相关性矩阵与特征向量求内积,得到两组任意两个特征向量间的内积值;
S42、对所述内积值进行排序,取出所述信源数目K对较大相关性的特征向量对;
所述向量对表达式为:
S43、通过所述K对特征向量及其所对应的特征值与信号相位关联关系,得到菲涅尔传播模型中的相位参数。
其中,步骤S43中解算接收信号的相位信息的表达式为:
其中,N是快拍数,u1(i+1)是对应的第1个特征向量中的第i+1个分量,L是估计结果K,与勾菲涅尔模型中的相位参数。
S5、对所述相位信息进行反解,得到信源方位。
步骤S5中反解信源方位的表达式为:
其中,d为天线间距,为菲涅尔模型中的相位参数,为第1个信源的方向角。
实施例2:
如图3所示,本发明提供了一种基于最小描述长度的近场源定位装置,包含以下部分:
协方差矩阵模块,对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
信源数目预估计模块;基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
相关性矩阵特征分解模块;对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
特征向量匹配模块,将信源数目预估计模块计算出的的信源数目K与相关性矩阵特征分解模块计算出的相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
反解信源方位模块,对相位信息进行反解,得到信源方位。
实施例3:
如图4所示,本发明实施例还提供了一种基于最小描述长度的近场源定位系统,包含数据接收模块、数据传输模块,还包含权利要求6所述的近场源定位装置;
其中,数据接收模块包含天线阵列、滤波器、放大器及A/D转换器,将接收到的信号发送给数据采集模块;
数据采集模块包含多路数据电路终接设备和多路传输信道,将数据接收模块发送的信号传输给近场源定位装置。
其中,天线阵列基于给定的空间排布,接收指定频段的信号,并根据实际定位场景选择相应增益、方向性系数等参数的天线。对于如图2所示的线性阵元,通过理论推导可以得出其接收数据的相关性。
对于L个独立近场窄带信号入射到阵元间距为d,个数为2N的均匀线阵,其天线接收数据具有如下形式:
式中,ni(k)为第i个阵元在第k时刻的加性噪声干扰;sl(k)是第1个信源在第k时刻的辐射信号;共同构成了菲涅尔模型中的相位延迟项,其中,
λ为信号波长,,θl与rl分别是第l个信号的入射角与距离。
另外,所述数据传输模块,天线采集的数据通过数据传输模块进行传输,多路接收信号进入前端多天线系统后沿多路进行传输,直至到达数据处理单元中的输入缓冲区进行处理。数据传输模块由多路数据电路终接设备(简称DCE)和多路传输信道(线路)构成。
所述数据电路终接设备:数据电路终接设备介于数据终端设备即采集多路信号的阵列天线和传输信道之间,其主要作用是实现信号的变换与编码、解码。在发送端对天线阵元的数据进行数据格式、规范的封装及预处理,将其变换成满足传输信道要求的信号,并对数据信号进行编码,以提高可靠性和有效性;在接收端进行相反变换,解调和解码,还原成所发送的信号。此外,数据电路终接设备还有向数据终端设备发送时钟信号的功能,确保与终端设备信号同步。若传输信道为模拟电路,调制解调器就是一个数据电路终接设备;若传输信道为数字信道,数据电路终接设备就是一种专门的数据服务单元,以实现数字信号的码型转换和电平变换、信道特性的均衡、同步时钟等。
且该系统进一步包括监测数据库模块和位置数据库模块;
监测数据库模块,对近场源定位装置的测试数据进行存储;
位置数据库模块,对近场源定位装置的位置数据进行存储。
所述监测数据库模块部署在基站侧,主要负责存储无线电监测数据,终端的频谱占用情况,结合本地的无线环境,计算出是否存在可用资源,是否存在非法用频情况,备份的监测数据为无线电监管部门进行频谱审计、频谱决策提供依据。数据库存储内容如表1所示。
表1监测数据库模块
监测数据库对区域无线电频谱使用情况进行存档,进而使后期数据分析时实现频域现场重构变得有据可循;同时也方便了后续网络业务评估、频域数据挖掘的进行。
所述位置数据库模块用于接收信源定位结果数据,将数据处理的结果即多信源位置参数存入。以时间标示的信源ID为主键,以信号参数、信源位置为属性进行关系存储;构建一个关系型数据库,即可快速查询监测得到的未知多信源的位置等信息。并且利用时间戳标记,将数据载入,进而能够实现信源定位的时域重现。
综上,本发明实施例通过提供一种基于最小描述长度的近场源定位方法、装置及系统,通过在信号方向向量匹配前进行信源数目预估计,进而指导方向向量匹配按照估计的参量个数进行,能够把近场源定位的计算复杂度降低,提高了计算效率。
在多目标信源定位中,基于接收数据矩阵分解的方法都需要对不同信号的方向向量进行匹配加以区分才能够最终确定所有信源的方位,即进行方向向量的匹配。现有方法均是以向量相关性作为评判指标,通过人为设定阈值加以分类。阈值设定的好坏将直接决定匹配结果的准确性,进而直接决定信源参数估计的精度。即缺少能够缩小方向匹配的解空间的客观准则,缺少能够保证匹配准确性、进而保证多信源参数估计的方法。本发明能充分发挥信源数目预估计的效能,提高定位的精准性与稳定度,实现了系统定位性能的整体提升。
较之传统高阶累积量方法如计算四阶相关性矩阵的平行因子方法,与进行参数空间谱峰搜索的方法如MUSIC方法相比较,本发明实现了近场源定位效率与定位准确性的双赢。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (7)
1.一种基于最小描述长度的近场源定位方法,其特征在于,包含以下步骤:
S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
S2、基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
S4、将步骤S2中的所述信源数目K与步骤S3中的所述相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
S5、对所述相位信息进行反解,得到信源方位;
其中,步骤S3包括:定义基于阵元滞后量的二阶统计量矩阵R1,该矩阵为N阶方阵,其m行n列元素定义如下:
式中,N为阵列天线的天线总数,τ1与v1共同构成了阵元间距,即阵元滞后量;
类似矩阵R1,交换阵元滞后量主键,可定义相关性矩阵R2的m行n列元素如下:
可定义如下基于一个时间间隙的二阶统计量:
结合两组限制条件,可得:
代入R1、R2、R3,得:
由式R1、R2、R3的形式,有:
R2R1 #A=AΩ-2Λ-2
R3R1 #A=AΦ
由上式可见,Α中各列向量同为R2R1 #与R3R1 #的特征向量;
其中,γi为信源参量,rsl为第l个信源的信号功率; Φ=diag{ω1,ω2,…ωL},ω1,ω2,…ωL分别为L个窄带信号的频率;
步骤S4包括步骤:
S41、对所述两组相关性矩阵R2R1 #与R3R1 #特征向量求内积,得到两组任意两个特征向量间的内积值;
S42、对所述内积值进行排序,取出K对较大相关性的特征向量对;
S43、通过所述K对特征向量对及其所对应的特征值与信号相位关联关系,得到菲涅尔传播模型中的相位参数;
步骤S43中解算接收信号的相位信息的表达式为:
式中,P是快拍数,ul(i+1)是R3R1 #对应的第l个特征向量中的第i+1个分量,L是估计出的信源数目K,与为菲涅尔模型中的相位参数。
2.如权利要求1所述的近场源定位方法,其特征在于,步骤S1中计算协方差矩阵特征向量表达式为:
式中,U为接收数据协方差矩阵的特征向量矩阵,Σ为由特征值组成的对角阵,ei与λi分别是特征向量及与之对应的特征值,M是特征向量数目。
3.如权利要求1所述的近场源定位方法,其特征在于,步骤S2中计算信源数目K的表达式为:
式中,k为模型参数,l1>l2>…>lM是样本协方差矩阵的特征值,M是特征向量数目,P为快拍数。
4.如权利要求1所述的近场源定位方法,其特征在于,步骤S5中反解信源方位的表达式为:
式中,d为天线间距,为第l个信源的方向角,λ是信号的波长,L为对信源数目K的估计值,为菲涅尔模型中的相位参数。
5.一种基于最小描述长度的近场源定位装置,其特征在于,包含以下部分:
协方差矩阵模块,对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
信源数目预估计模块,基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
相关性矩阵特征分解模块,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
特征向量匹配模块,将信源数目预估计模块计算出的信源数目K与相关性矩阵特征分解模块计算出的相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
反解信源方位模块,对相位信息进行反解,得到信源方位;
相关性矩阵特征分解模块还用于定义基于阵元滞后量的二阶统计量矩阵R1,该矩阵为N阶方阵,其m行n列元素定义如下:
式中,N为阵列天线的天线总数,τ1与v1共同构成了阵元间距,即阵元滞后量;
类似矩阵R1,交换阵元滞后量主键,可定义相关性矩阵R2的m行n列元素如下:
可定义如下基于一个时间间隙的二阶统计量:
结合两组限制条件,可得:
代入R1、R2、R3,得:
由式R1、R2、R3的形式,有:
R2R1 #Α=ΑΩ-2Λ-2
R3R1 #Α=ΑΦ
由上式可见,Α中各列向量同为R2R1 #与R3R1 #的特征向量;
其中,γi为信源参量,rsl为第l个信源的信号功率; Φ=diag{ω1,ω2,…ωL},ω1,ω2,…ωL分别为L个窄带信号的频率;
特征向量匹配模块,还用于对所述两组相关性矩阵R2R1 #与R3R1 #特征向量求内积,得到两组任意两个特征向量间的内积值,对所述内积值进行排序,取出K对较大相关性的特征向量对,通过所述K对特征向量对及其所对应的特征值与信号相位关联关系,得到菲涅尔传播模型中的相位参数;
其中,解算接收信号的相位信息的表达式为:
式中,P是快拍数,ul(i+1)是R3R1 #对应的第l个特征向量中的第i+1个分量,L是估计出的信源数目K,与为菲涅尔模型中的相位参数。
6.一种基于最小描述长度的近场源定位系统,其特征在于,包含数据接收模块、数据采集模块,还包含权利要求5所述的近场源定位装置;
其中,
数据接收模块包含天线阵列、滤波器、放大器及A/D转换器,将接收到的信号发送给数据采集模块;
数据采集模块包含多路数据电路终接设备和多路传输信道,将数据接收模块发送的信号传输给近场源定位装置。
7.如权利要求6所述的近场源定位系统,其特征在于,还包含监测数据库模块和位置数据库模块;
所述监测数据库模块,对近场源定位装置的测试数据进行存储;
所述位置数据库模块,对近场源定位装置的位置数据进行存储。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310326242.5A CN103424735B (zh) | 2013-07-30 | 2013-07-30 | 基于最小描述长度的近场源定位方法、装置及系统 |
PCT/CN2013/083327 WO2015014007A1 (zh) | 2013-07-30 | 2013-09-11 | 基于最小描述长度的近场源定位方法、装置及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310326242.5A CN103424735B (zh) | 2013-07-30 | 2013-07-30 | 基于最小描述长度的近场源定位方法、装置及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103424735A CN103424735A (zh) | 2013-12-04 |
CN103424735B true CN103424735B (zh) | 2016-06-08 |
Family
ID=49649757
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310326242.5A Active CN103424735B (zh) | 2013-07-30 | 2013-07-30 | 基于最小描述长度的近场源定位方法、装置及系统 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN103424735B (zh) |
WO (1) | WO2015014007A1 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103701539A (zh) * | 2013-12-25 | 2014-04-02 | 陕西福友电器设备技术发展有限公司 | 一种cdma用户数量测定方法 |
CN105182322B (zh) * | 2015-08-31 | 2017-12-15 | 西北大学 | 基于反射信号相位差的被动式定位方法 |
CN106680597B (zh) * | 2016-12-13 | 2019-02-15 | 云南电网有限责任公司电力科学研究院 | 确定最强场源方位的方法及系统 |
CN107071898B (zh) * | 2017-04-14 | 2019-07-19 | 中国人民解放军信息工程大学 | 移动通信信号源数据域直接位置估计方法及其装置 |
CN108966137B (zh) * | 2018-08-10 | 2021-06-11 | 欧普照明股份有限公司 | 一种信源图形化近场识别方法 |
CN110187304B (zh) * | 2019-05-21 | 2021-05-04 | 泰凌微电子(上海)股份有限公司 | 一种信号到达角估计方法及装置 |
CN111323744B (zh) * | 2020-03-19 | 2022-12-13 | 哈尔滨工程大学 | 一种基于mdl准则的目标个数和目标角度估计方法 |
CN111965598B (zh) * | 2020-08-17 | 2023-11-17 | 中国人民解放军国防科技大学 | 一种利用非均匀阵列进行辐射源定位方法 |
CN113030850B (zh) * | 2021-03-04 | 2022-11-11 | 西北工业大学 | 基于复数域表征与学习的近场源智能定位方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1297822C (zh) * | 2003-02-21 | 2007-01-31 | 重庆邮电学院 | 一种基于td-scdma无线定位来波方向的估计方法 |
CN102645646B (zh) * | 2012-05-03 | 2013-06-26 | 西安费斯达自动化工程有限公司 | 一种多信源不确定融合定位方法 |
CN102736063A (zh) * | 2012-07-06 | 2012-10-17 | 绍兴文理学院 | 近场声源定位方法 |
CN102841344B (zh) * | 2012-09-13 | 2015-07-15 | 电子科技大学 | 一种少阵元近场宽带信号源参数估计方法 |
CN102928818A (zh) * | 2012-10-18 | 2013-02-13 | 东南大学 | 一种基于近场波束形成的碰摩声发射源的定位方法 |
-
2013
- 2013-07-30 CN CN201310326242.5A patent/CN103424735B/zh active Active
- 2013-09-11 WO PCT/CN2013/083327 patent/WO2015014007A1/zh active Application Filing
Non-Patent Citations (3)
Title |
---|
A Novel Near-Field Source Localization Algorithm;Xu Qian 等;《IEEE》;20101231;188-191 * |
Detection of Signals by Information Theoretic Criteria;MATI WAX 等;《IEEE TRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING》;19850430;第ASSP-33卷(第2期);387-392 * |
Detection of the Number of Coherent Signals by the MDL Principle;MATI WAX 等;《IEEE TRANSACTIONS ON ACOUSTICS. SPEECH. AND PROCESSING》;19890831;第37卷(第8期);1190-1196 * |
Also Published As
Publication number | Publication date |
---|---|
CN103424735A (zh) | 2013-12-04 |
WO2015014007A1 (zh) | 2015-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103424735B (zh) | 基于最小描述长度的近场源定位方法、装置及系统 | |
Jaranowski et al. | Gravitational-wave data analysis. Formalism and sample applications: the Gaussian case | |
Wan et al. | PD source diagnosis and localization in industrial high-voltage insulation system via multimodal joint sparse representation | |
CN101592721B (zh) | 基于特征值重构的相干信号到达角估计方法 | |
CN102317808B (zh) | 用于确定相干源抵达的方位角和仰角的方法 | |
CA3041907A1 (en) | Direction of arrival estimation | |
CN101779140B (zh) | 在任意传感器网络上使用空间平滑来测量相干源的到达角的方法 | |
Ho et al. | Efficient method for estimating directions-of-arrival of partially polarized signals with electromagnetic vector sensors | |
CN108020812A (zh) | 基于特殊三平行线阵结构的二维doa估计方法 | |
CN108089147A (zh) | 改进的短波单位定位方法 | |
Tadaion et al. | A fast multiple-source detection and localization array signal processing algorithm using the spatial filtering and ML approach | |
Bar-Shalom et al. | Transponder-aided single platform geolocation | |
Li et al. | Multi-TDOA estimation and source direct position determination based on parallel factor analysis | |
CN109239648A (zh) | 基于对称循环频率的谱相关子空间测向方法 | |
Wagner et al. | Compressive MIMO beamforming of data collected in a refractive environment | |
Picard et al. | Localization of multiple emitters by spatial sparsity methods in the presence of fading channels | |
Soares et al. | Environmental inversion using high-resolution matched-field processing | |
Xu et al. | Two‐Dimensional Direction‐of‐Arrival Estimation of Noncircular Signals in Coprime Planar Array with High Degree of Freedom | |
Yang et al. | Service antenna array-based parameter estimation for decimeter level indoor localization | |
Wang et al. | Passive TDOA and DOA based HF geolocation without ionosphere information | |
Zhao et al. | A frequency domain direct localization method based on distributed antenna sensing | |
Lee et al. | Machine Learning-based Indoor Positioning Systems Using Multi-Channel Information. | |
Li et al. | Underwater source localization in the presence of strong interference | |
Chu et al. | Exploration on 2D DOA estimation of linear array motion: Uniform linear motion | |
Nagaraju et al. | Performance Analysis of Sparse Array Using Compressive Sensing in A Closed Multi-Path Environment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |