• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 18
    • 下载费用:30 金币  

    重庆时时彩gt娱乐平台: 基于最小描述长度的近场源定位方法、装置及系统.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201310326242.5

    申请日:

    2013.07.30

    公开号:

    CN103424735A

    公开日:

    2013.12.04

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01S 5/14申请日:20130730|||公开
    IPC分类号: G01S5/14 主分类号: G01S5/14
    申请人: 北京邮电大学
    发明人: 张轶凡; 冯志勇; 白杨; 刘尚; 高明菲; 晏潇; 张奇勋
    地址: 100876 北京市海淀区西土城路10号
    优先权:
    专利代理机构: 北京路浩知识产权代理有限公司 11002 代理人: 王莹
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201310326242.5

    授权公告号:

    ||||||

    法律状态公告日:

    2016.06.08|||2013.12.25|||2013.12.04

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本发明提供一种基于最小描述长度的近场源定位方法、装置及系统,涉及移动业务中的信源定位系统技术领域。该方法包含步骤:S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;S2、基于最小描述长度准则,对协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;S4、将步骤S2中的信源数目K与步骤S3中的相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;S5、对所述相位信息进行反解,得到信源方位。本发明既能降低近场源定位的计算复杂度,提高计算效率;又能充分发挥数据预处理的效能,提高定位的精准性与稳定度。

    权利要求书

    权利要求书
    1.  一种基于最小描述长度的近场源定位方法,其特征在于,包含以下步骤:
    S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
    S2、基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
    S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
    S4、将步骤S2中的所述信源数目K与步骤S3中的所述相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
    S5、对所述相位信息进行反解,得到信源方位。

    2.  如权利要求1所述的近场源定位方法,其特征在于,步骤S1中计算协方差矩阵特征向量表达式为:
    Rxx=UH=Σi=1MλieieiH]]>
    式中,U为接收数据协方差矩阵的特征向量矩阵Σ为由特征值组成的对角阵,ei与λi分别是特征向量及与之对应的特征值,M是特征向量数目。

    3.  如权利要求1所述的近场源定位方法,其特征在于,步骤S2中计算信源数目K的表达式为:
    MDL(k)=-log(i=k+1Mli1M-kΣi=k+1MliM-k)(M-k)N+12k(2M-k)log(N)]]>
    式中,k为模型参数,l1>l2>...>lM是样木协方差矩阵,M是特征向量数目。

    4.  如权利要求1所述的近场源定位方法,其特征在于,步骤S4包括步骤:
    S41、对所述两组相关性矩阵与特征向量求内积,得到两组任意两个特征向量间的内积值;
    S42、对所述内积值进行排序,取出所述信源数目K对较大相关性的特征向量对;
    S43、通过所述K对特征向量及其所对应的特征值与信号相位关联关系,得到菲涅尔传播模型中的相位参数。

    5.  如权利要求4所述的近场源定位方法,其特征在于,步骤S43中解算接收信号的相位信息的表达式为:


    式中,N是快拍数,ul(i+1)是对应的第l个特征向量中的第i+1个分量,L是估计结果K,与为菲涅尔模型中的相位参数。

    6.  如权利要求1所述的近场源定位方法,其特征在于,步骤S5中反解信源方位的表达式为:
    θ1~=-arcsin(γ1~λ2πd),l=1,2,...,L]]>
    式中,d为天线间距,为第l个信源的方向角。

    7.  一种基于最小描述长度的近场源定位装置,其特征在于,包含以下部分:
    协方差矩阵???,对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
    信源数目预估计???,基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
    相关性矩阵特征分解???,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
    特征向量匹配???,将信源数目预估计??榧扑愠龅牡男旁词縆与相关性矩阵特征分解??榧扑愠龅南喙匦跃卣筇卣飨蛄拷岷?,进行特征向量匹配,解算出接收信号的相位信息;
    反解信源方位???,对相位信息进行反解,得到信源方位。

    8.  一种基于最小描述长度的近场源定位系统,其特征在于,包含数据接收???、数据采集???,还包含权利要求6所述的近场源定位装置;
    其中,
    数据接收??榘煜哒罅?、滤波器、放大器及A/D转换器,将接收到的信号发送给数据采集???;
    数据采集??榘嗦肥莸缏分战由璞负投嗦反湫诺?,将数据接收??榉⑺偷男藕糯涓≡炊ㄎ蛔爸?。

    9.  如权利要求8所述的近场源定位系统,其特征在于,还包含监测数据库??楹臀恢檬菘饽??;
    所述监测数据库???,对近场源定位装置的测试数据进行存储;
    所述位置数据库???,对近场源定位装置的位置数据进行存储。

    说明书

    说明书基于最小描述长度的近场源定位方法、装置及系统
    技术领域
    本发明涉及移动业务中的信源定位系统技术领域,具体涉及一种基于最小描述长度的近场源定位方法、装置及系统。
    背景技术
    无线电被动定位技术是无线电监测与管理中的重要技术手段,在雷达、声纳、炮兵电声测量等军事领域,石油地震勘测、地下管道泄漏检测定位和故障检测等工业领域,以及移动定位业务、无线电检测等通信领域都有着广泛的应用。
    传统的无线电发射源定位是采用转动接收天线(传感器)的角度来实现的。但这种方法存在测量精度与测量速度的瓶颈,难以应对空间存在多个目标的情况。而通过采用阵列信号处理技术,合理选取阵元数目和阵元间距可以较好地进行多目标方位探测。阵列即按一定方式布置在空间上不同位置的传感器组,利用信号的空域特性增强信号并且有效提取信号的空域信息。因此阵列信号处理也常称为空域信号处理。与传统的单个定向传感器相比,阵列信号处理具有灵活的波束控制、较高的信号增益、极强的抗干扰能力及空间超分辨能力等优点,相关的研究工作不断发展与深入,其应用范围也不断扩大。并且随着微电子、数字信号处理、并行处理等技术的迅猛发展,阵列信号处理技术也得到较快的发展。阵列信号处理最主要的研究内容包括DOA估计和波束形成。较早的DOA估计方法又称为波束形成方法,而该波束形成方法利用了空域维纳滤波的匹配概念,由阵列流形在信号空间中的投影大小判定信号方向,后来随着研究的深入,出现了多重信号分类(MUSIC)基于接收信号矩阵分解的谱估计方法,该方法利用信号空间与噪声空间的正交性,先对接收数据的协方差矩阵进行特征分解, 再对方向向量相关性进行搜索。将相关性的极小值点对应的信源参数作为信源位置的最大似然估计,从而实现了对多目标信源的精确定位。
    然而,以上多目标信源定位的方法,如MUSIC具有较高的计算复杂度。多重信号分类在进行谱峰搜索的过程中需要对信源方位角的候选区间进行全面搜索,对于需要多维空间表征的信号参数,搜索复杂度是指数级增加的。因此在实际应用中,该方法在计算效率上存在严重的缺陷,极大限制了定位系统的时效性。
    发明内容
    (一)解决的技术问题
    针对现有技术的不足,本发明提供一种基于最小描述长度的近场源定位方法、装置及系统,使得近场源定位的计算复杂度降低。
    (二)技术方案
    为实现以上目的,本发明通过以下技术方案予以实现:
    一种基于最小描述长度的近场源定位方法,包含以下步骤:
    S1、对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
    S2、基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目K;
    S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
    S4、将步骤S2中的所述信源数目K与步骤S3中的所述两组相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
    S5、对所述相位信息进行反解,得到信源方位。
    优选的,步骤S1中计算协方差矩阵特征向量表达式为:
    Rxx=UH=Σi=1MλieieiH]]>
    式中,U为接收数据协方差矩阵的特征向量矩阵,∑为由特征值组成的对角阵,ei与λi分别是特征向量及与之对应的特征值,M是特征向量数目。
    优选的,步骤S2中计算信源数目K的表达式为:
    MDL(k)=-log(i=k+1Mli1M-kΣi=k+1MliM-k)(M-k)N+12k(2M-k)log(N)]]>
    式中,k为模型参数,l1>l2>…>lM是样本协方差矩阵,M是特征向量数目。
    优选的,步骤S4包括步骤:
    S41、对所述两组相关性矩阵与特征向量求内积,得到两组任意两个特征向量间的内积值;
    S42、对所述内积值进行排序,取出所述信源数目K对较大相关性的特征向量对;
    S43、通过所述K对特征向量及其所对应的特征值与信号相位关联关系,得到菲涅尔传播模型中的相位参数。
    优选的,步骤S43中解算接收信号的相位信息的表达式为:


    式中,N是快拍数,u1(i+1)是对应的第1个特征向量中的第i+1个分量,L是估计结果K,与为菲涅尔模型中的相位参数。
    优选的,步骤S5中反解信源方位的表达式为:
    θ1~=-arcsin(γ1~λ2πd),l=1,2,...,L]]>
    式中,d为天线间距,为第1个信源的方向角。
    本发明提供了一种基于最小描述长度的近场源定位装置,包含以下部分:
    协方差矩阵???,对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
    信源数目预估计???,基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目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之间存在如下关系:
    ωi=2πdλsinθi---(1)]]>
    式中,d是两个相邻阵元之间的距离,λ为信号波长,阵元间距d满足d≤λ/2,否则相位差ωi可能大于π,进而产生方向模糊,即θi与π+θi都可能是si(n)的波达方向。又由于该阵列在空间上的排布为等距线阵,信号si(n)到达第k个阵元的电波与到达参考阵元的电波之间的相位差为:
    (k-1)ωi=2π(k-1)dλsinθi---(2)]]>
    因此,信号si(n)在第k个阵元上的接收信号为:
    若阵列由M个阵元组成,则信号si(n)到达各阵元的方向向量可表示为:
    a(θi)=[1,e-iωi,..,e-i(M-1)ωi]T=[a1(θi),..,aM(θi)]T---(3)]]>
    如果共有p个信号位于远场(p<M),则在第k个阵元上的观测或接收信号xk(n)为:
    a(θi)=[1,e-iωi,..,e-i(M-1)ωi]T=[a1(θi),..,aM(θi)]T]]>
    式中,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)
    这样,阵元接收数据可以写成如下向量形式:
    x(n)=Σi=1pa(θi)si(n)+ek(n)=A(θ)s(n)+e(n)]]>
    式中,
    A(θ)=[a(θ1),..,a(θp)]]]>
    =1111e-jω1e-jω2..e-jωp........e-j(M-1)ω1e-j(M-1)ω2..e-j(M-1)ωp]]>
    s(n)=[s1(n),s2(n),..,sN(n)]T---(7)]]>
    分别为M×p维方向向量矩阵和N×1维信号向量。具有上述形式范德蒙德矩阵的特点是:若ω1,ω2,..,ωp互不相同,则矩阵的各列相互独立,也即矩阵是列满秩矩阵。
    由阵元接收数据的数学模型,阵列的协方差矩阵
    Rxx=E[x(n)x(n)H]=AE[SSH]AH=ARssAH+RN    (8)
    式中,RN为噪声协方差矩阵,对于空间的理想高斯白噪声,且噪声方差为则有下式
    Rxx=E[x(n)x(n)H]=ARssAH+δn2I---(9)]]>
    进而对Rxx进行特征分解,得到步骤S1中计算协方差矩阵特征向量表达式为:
    Rxx=UH=Σi=1MλieieiH---(10)]]>
    其中,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之间的误差;
    即:Mmdl=argminMi∈M{|Lm(Mi)|+|Lc(D|Mi)|}---(11)]]>
    任何一个观测序列(即对象),可以看成是由两个序列组成,一个是确定性序列,一个是“纯”随机过程的序列(即观测序列与确定性序列的残差序列)。对应于MDL准则,Mi刻画的是确定性序列;模型Mi与观测对象的误差是随机过程序列。对于确定性模型的部分,描述性语言Lc(D|Mi)可以选取如AR模型,多项式模型,插值模型等的描述性模型。此时描述性语言的长度Lm(Mi)为模型中参数的个数。对于不确定性部分,描述性语言Lc(D|Mi)可以选取概率分布模型。信息论专门研究随机性事件。其中的Shannon编码定理表明,用概率分布模型描述随机过程时,其描述长度等于对象概率分布以2为底对数的相反数。已经证明,在大样本数据情况下,最小描述长度准则可有下面的指标度量:
    MDL=-2logf(X|Θ^)+12klogN---(12)]]>
    上式中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)可作如下分解:
    R(k)=Σi=1k(λi-σ2)ViViH+σ2I---(14)]]>
    式中,λ1,λ2,...,λk和V1,V2,...,Vk分别是矩阵R(k)的特征值与特征向量。下面用θ(k)表示模型参数向量,即
    Θ(k)T=(λ1,λ2,...,λk,σ2,V1,V2,...,Vk)---(15)]]>
    参量化后,现在来引出该检测模型的信息论指标。由于假设观测数据是零均值的复高斯随机向量,则它们的联合概率密度有如下形式:
    f(x(t1),x(t2),...,x(tN)|Θ(k))=i=1N1πMdetR(k)e-x(ti)H[R(k)]-1x(ti)---(16)]]>
    对上式取负对数,并且忽略与参数向量θ(k)的无关项,可得如下对数似然函数
    L(Θ(k))=-NlogdetR(k)-tr[R(k)]-1R^---(17)]]>
    式中,是样本协方差矩阵,即
    R^=1NΣi=1Nx(ti)x(ti)H---(18)]]>
    最大似然估计即是求取特定θ(k)使公式(5-10)中的似然函数L(θ(k))取得最大。文献[26]中Anderson给出了该问题的闭式解:
    λl^=li,i=1,2,...,k]]>
    σ^2=1M-kΣi=k+1Mli]]>
    V^1=Ci,i=1,2,...,k---(19)]]>
    这里,l1>l2>…>lM是样本协方差矩阵的特征值,且C1,C2,...,CM为对应的特征向量。将上述最大似然解代入式似然函数,简化可得:
    L(Θ(k))=log(i=k+1Mli1/(M-k)Σi=k+1Mli/(M-k))(M-k)N---(20)]]>
    可以看出,上式中括号内的分式就是M-k个小特征值的几何均值与算术均值的比值。
    这里,可以通过计算θ(k)张成的参数空间的自由度得到θ(k)中自由参数的个数。由于复协方差矩阵的特征值为实数,但其特征向量为复向量,所以θ(k)中实际含有k+1+2Mk个参数。然而这k+1+2Mk个参数的取值是有限制的:特征向量本身相互正交,并且假定特征向量为单位向量(已经过单位化)。这等价于参数空间的自由度因标准化降低了2k, 因正交性降低了则模型自由参数个数为:
    k+1+2Mk-212k(k-1)-2k=k(2M-k)+1---(21)]]>
    进而得到步骤S2中计算信源数目K的表达式为:
    MDL(k)=-log(i=k+1Mli1M-kΣi=k+1MliM-k)(M-k)N+12k(2M-k)log(N)---(22)]]>
    其中,k为模型参数(控制变量),l1>l2>…>lM是样本协方差矩阵,M是特征向量数目。
    S3、与步骤S1并行,对接收信号进行两个相关性矩阵特征分解,得到两组相关性矩阵特征向量;
    定义基于阵元滞后量的二阶统计量矩阵R1。该矩阵为N阶方阵,其m行n列元素定义如下:
    R1(m,n)→defE[xm-n+τ1(k)xm-n-υ1*(k)],1m,nN---(23)]]>
    式中,τ1与υ1共同构成了阵元间距,即阵元滞后量。下面记m-n+τ1为Δ(τ1),记m-n-υ1为Δ(υ1),则当信源与阵元噪声之间两两不相关时



    这里,为第l个信源的信号功率。类似矩阵R1,交换阵元滞后量主键,可定义矩阵R2的m行n列元素如下
    R2(m,n)→defE[xn-m+τz(k)xn-m-υz*(k)],1m,nN]]>
    记n-m+τ2为Δ(τ2),记n-m-υ2为Δ(υ2),则R2(m,n)可表示如下


    现令R1中的τ1+υ1=1,R2中的τ2+υ2=-1,则可将其分别转化为矩阵形式。令Γ=diag{Γs1,Γsz,...,ΓsL},]]>Ω=diag{ejγ1,ez,...,ejL},]]>信号方向向量方向矩阵A=[a1,a2,...,aL],为范德蒙德矩阵。则R1与R2有如下形式
    R1=AΩτ1+v1Λτ1z-υ1zΓAH---(27)]]>
    R2=τz+vzΛτzz-υzzΓAH---(28)]]>
    显见上述统计量是从空间阵元排布的角度考察不同阵元滞后量下阵元接收数据的相关性。下面再从时域的角度,考虑不同时延的相关统计量。由于发射信号为窄带随机过程,所以k时刻的信号s(k)与k+1时刻的信号s(k+1)有相角存在一个时间间隙的差量。由此,可定义如下基于一个时间间隙的二阶统计量:
    R3(m,n)→defE[xm-n+τ3(k+1)xm-n-υ3*(k)],1m,nN---(29)]]>
    假设L个窄带信号的频率分别为ω1,ω2,...,ωL,令Φ=diaθ{ω1,ω2,...,ωL},则经类似R1、R2的化简过程进行处理后,R3(m,n)有如下形式

    将其转化为矩阵形式,可得
    R3=Ωτ3+υ3Λτ32-υ32ΓAH---(31)]]>
    基于上述阵元滞后量的二阶统计特性分析,考虑单个阵元延迟的相关性矩阵。已有约束条件:τ1+υ1=1;τ2+υ2=-1;τ3+υ3=1即可保证延迟量为1,也就是考察单间隔阵元的二阶统计特性。在此基础上,再对延迟量取值做如下限制:
    τ12-υ12=1τ22-υ22=-1τ32-υ32=1---(32)]]>
    结合两组限制条件,可得
    τ1=1,υ1=0τ2=0,υ2=-1τ3=1,υ3=0---(33)]]>
    代入R1、R2、R3,得
    R1=AΩΛΓAHR2=AΩ-1Λ-1AHR3=AΦΩΛΓAH---(34)]]>
    由假设所有信源都具有非零功率,可知对角矩阵Γ对角线元素非零,且为可逆矩阵。此外,由于信源参量及各不相同,因此矩阵A各列向量互不相关,为满秩阵。同时,矩阵R1,R2,R3是秩L的N阶方阵。
    对矩阵R1进行特征值分解,有
    R1=Σl=1Lv1u1u1H=UVUH---(35)]]>
    上式中,U=[u1,u2,...,uL]为R1的L个特征向量构成的矩阵,V=diag{v1,v2,...,vL}为R1的L个特征值构成的对角矩阵。定义矩阵R1的伪逆:
    R1#=Σi=1Lvi-1uiuiH=UV-1UH---(36)]]>
    由式R1、R2、R3的形式,有
    R2R1#A=AΩ-2Λ-2---(37)]]>
    R3R1#A=---(38)]]>
    由上式可见,A中各列向量同为与的特征向量。并且由等式右边的对角矩阵可知矩阵与的公共特征向量ai对应的特征值为和由于这两个特征值是对应于同一个特征向量的,所以特征向量确定后即可将和与进行配对。当然,实际中由于噪声的存在对两个矩阵进行分解后得到的特征向量并不相等。这时考虑采用求取复向量内积的方式提取向量的相关度量,将相关性最大的向量视为“相等”的向量进行配对。
    S4、将步骤S2中的所述信源数目K与步骤S3中的所述相关性矩阵特征向量结合,进行特征向量匹配,解算出接收信号的相位信息;
    优选的,步骤S4包括步骤:
    S41、对所述两组相关性矩阵与特征向量求内积,得到两组任意两个特征向量间的内积值;
    S42、对所述内积值进行排序,取出所述信源数目K对较大相关性的特征向量对;
    所述向量对表达式为:

    S43、通过所述K对特征向量及其所对应的特征值与信号相位关联关系,得到菲涅尔传播模型中的相位参数。
    其中,步骤S43中解算接收信号的相位信息的表达式为:


    其中,N是快拍数,u1(i+1)是对应的第1个特征向量中的第i+1个分量,L是估计结果K,与勾菲涅尔模型中的相位参数。
    S5、对所述相位信息进行反解,得到信源方位。
    步骤S5中反解信源方位的表达式为:
    θ1~=-arcsin(γ1~λ2πd),l=1,2,...,L---(41)]]>
    其中,d为天线间距,为菲涅尔模型中的相位参数,为第1个信源的方向角。
    实施例2:
    如图3所示,本发明提供了一种基于最小描述长度的近场源定位装置,包含以下部分:
    协方差矩阵???,对接收信号进行协方差矩阵特征分解,得到协方差矩阵特征向量;
    信源数目预估计???;基于最小描述长度准则,对所述协方差矩阵特征向量进行信源数目预估计,得到预估计的信源数目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方法相比较,本发明实现了近场源定位效率与定位准确性的双赢。
    需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
    以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

    关 键 词:
    基于 最小 描述 长度 近场 定位 方法 装置 系统
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:基于最小描述长度的近场源定位方法、装置及系统.pdf
    链接地址://www.4mum.com.cn/p-5778381.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03