• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
    • / 17
    • 下载费用:30 金币  

    重庆时时彩官方网站博彩公司: 一种远场和近场混合信号源的定位方法.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201410174341.0

    申请日:

    2014.04.28

    公开号:

    CN103954931A

    公开日:

    2014.07.30

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01S 5/20申请日:20140428|||公开
    IPC分类号: G01S5/20; G06F19/00(2011.01)I 主分类号: G01S5/20
    申请人: 西安交通大学
    发明人: 左炜亮; 辛景民; 杨紫堃; 郑南宁
    地址: 710049 陕西省西安市咸宁西路28号
    优先权:
    专利代理机构: 西安通大专利代理有限责任公司 61200 代理人: 陆万寿
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410174341.0

    授权公告号:

    ||||||

    法律状态公告日:

    2016.06.08|||2014.08.27|||2014.07.30

    法律状态类型:

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

    摘要

    本发明公开了一种远场和近场混合信号源的定位方法,利用斜投影算子,有效地从混合信号中分离出远场信号,得到仅包含近场信号信息的传感器协方差矩阵,进而计算出近场信号的一维方位信息,另外,对于斜投影算子的计算,提出了一种更有效的估计方法,并引入交替迭代过程来消除近场信源定位中出现的饱和现象,本发明降低了计算复杂度,方法简单有效。

    权利要求书

    权利要求书
    1.  一种远场和近场混合信号源的定位方法,其特征在于,包括以下步骤:利用消除噪声后的混合信号估计远场信号的波达方向角,利用斜投影算子从消除噪声后的混合信号中分离出远场信号,得到仅包含近场信号信息的协方差矩阵,根据仅包含近场信号信息的协方差矩阵估计出近场信号的波达方向角及距离,所述混合信号为入射到对称均匀线阵上的K个非相干窄带信号{sk(n)},所述对称均匀线阵包含2M+1个全向传感器阵元,M的取值范围为M≥K,阵元间距为d,设前K1个入射信号为远场信号,远场信号的方位信息为后K2个入射信号为近场信号,近场信号的方位信息为K=K1+K2,θk表示第k个入射信号的波达方向角,所述波达方向角为第k个入射信号相对于y轴的逆时针夹角,rk是第k个入射信号相对于坐标原点的距离。

    2.  如权利要求1所述一种远场和近场混合信号源的定位方法,其特征在于,所述远场信号的波达方向角的估计方法包括以下步骤:
    1)消除阵列协方差矩阵R中的噪声部分,得到无噪协方差矩阵将分为两部分:
    2)构造代价函数
    gf(θ)=afH(θ)ΠQ^af(θ)---(1)]]>
    式(1)中(·)H表示共轭转置;
    ΠQ^=Q^(Q^HQ^)-1Q^H=Q^(I2M+1-K-P^H(P^P^H+IK)-1P^)Q^H]]>
    Q^=[P^T,-I2M+1-K]T]]>
    P^=(G^HG^)-1G^HH^]]>
    (·)T表示转置,Im表示m×m的单位矩阵;
    3)根据式(1)构造多项式其中p(z)=Δ[zM,...,z,1,z-1,...,z-M]T,z=Δej2πdsinθ/λ,]]>λ表示入射信号波长,j表示单位虚数,j2=-1,通过求多项式的K1个零相位点来估计远场信号的波达方向角。

    3.  如权利要求2所述一种远场和近场混合信号源的定位方法,其特征在于,所述步骤1)具 体包括以下步骤:
    a、根据对称均匀线阵接收的数据求得阵列协方差矩阵R的估计值:
    R=1NΣn=1Nx(n)xH(n)---(2)]]>
    其中,N表示采样数,x(n)表示阵列接收数据;
    b、将阵列协方差矩阵R划分为如下形式

    c、估计噪声方差σ2
    σ2=tr{R22Π}tr{Π}---(4)]]>
    其中tr{·}表示求矩阵的迹,表示广义逆;
    d、无噪协方差矩阵

    4.  如权利要求2所述一种远场和近场混合信号源的定位方法,其特征在于,所述近场信号的波达方向角及距离的估计方法包括以下步骤:

    2.  1)估计斜投影算子

    2.  2)消除无噪协方差矩阵中的远场信号,得到仅包含近场信号信息的协方差矩阵Rn,Rn=(I2M+1-EAf|An)R‾(I2M+1-EAf|An)H;]]>

    2.  3)利用仅包含近场信号信息的协方差矩阵Rn的反对角元素信息构造一个新的协方差矩阵并将新构造的协方差矩阵分块:

    2.  4)构造代价函数
    fn=a‾nH(θ)ΠQ‾^a‾n(θ)---(12)]]>
    式(12)中,ΠQ‾^=Q‾^(Q‾^HQ‾^)-1Q‾^H=Q‾^(IM+1-K2-P‾^H(P‾^P‾^H+IK2)-1P‾^)Q‾^H,P‾^=(G‾^HG‾^)-1G‾^HH‾^,]]>Q‾^=[P‾^T,-IM+1-K2]T;]]>

    2.  5)根据式(12)构造多项式pN(z)=ΔzMpH(z)ΠQ‾^p(z),]]>其中p(z)=Δ[1,z-1,...,z-M]T,]]>z=Δej4πdsinθ/λ,]]>通过求多项式pN(z)=ΔzMpH(z)ΠQ‾^p(z)]]>的K2个零相位点来估计近场信号的波达方向角;

    2.  6)构造代价函数
    f‾N(r)=anH(r,θ^)ΠQ^an(r,θ^)=μnH(θ^,r)DH(θ^)ΠQ^D(θ^)μn(θ^,r)---(13)]]>
    其中μn(θ^,r)=[eM2φ,e(M-1)2φ,...,1]T,φ=Δπd2cos2θ^/(λr),D(θ^)=Δ[Φ(θ^),Φ*(θ^)J]T,]]>是块对角矩阵,Φ(θ^)=diag{e-jω^M,e-jω^(M-1),...,1},]]>(·)*表示取共轭,ω^=Δ-2πdsinθ^/λ,]]>J是一个(M+1)×M的选择矩阵,定义为J=[eM,eM-1,…,e1],em表示单位矩阵IM+1的第m列,m=1,2,…,M;

    2.  7)根据式(13)构造多项式pr(z)=ΔzM2pH(z)DH(θ^)ΠQ^D(θ^)p(z),]]>其中p(z)=Δ[1,z,...,zM2]T,]]>z=Δeπd2cos2θ^/(λr),]]>通过求多项式pr(z)=ΔzM2pH(z)DH(θ^)ΠQ^D(θ^)p(z)]]>的K2个零相位点来估计近场信号的距离;

    2.  8)更新斜投影算子返回步骤2.2),继续顺序执行,如此往复循环,直至满足精度要求时结束。

    5.  如权利要求4所述一种远场和近场混合信号源的定位方法,其特征在于,所述步骤2.1)中斜投影算子的估计方法包括以下步骤:
    a、构造一个(2M+1)×(2M+1)矩阵:
    Ro=ΔR‾ΠAf---(14)]]>
    其中Af是远场信号方向矩阵,由估计出的远场信号的波达方向角计算得出;
    b、对RoΔ进行QR分解得斜投影算子由下述方式获得

    其中是矩阵的后2M+1-K2列,Δ是(2M+1)×(2M+1)的置换矩阵,不改变矩阵Ro的列相关性,表示广义逆。

    6.  如权利要求4所述一种远场和近场混合信号源的定位方法,其特征在于,所述步骤2.3) 具体包括以下步骤:
    仅包含近场信号信息的协方差矩阵Rn的反对角元素rn(i)定义为
    rn(i)=Δ(Rn)i,2M+2-i=Σk=K1+1Krsk2e-j2(M+1-i)ωk---(16)]]>
    i=1,2,…,2M+1,是第k个近场信号的功率;


    7.  如权利要求4所述一种远场和近场混合信号源的定位方法,其特征在于,更新斜投影算子的方法为:
    EAf|An=Af(AfHΠAnAf)-1AfHΠAn---(25)]]>
    Af是远场信号方向矩阵,由估计出的远场信号的波达方向角计算得出。

    说明书

    说明书一种远场和近场混合信号源的定位方法
    技术领域
    本发明属于阵列信号处理技术领域,具体涉及一种远场和近场混合信号源的定位方法(LOFNS)。
    背景技术
    信源定位在雷达、声纳、无线通信等领域具有重要的应用,在一些实际应用,如基于麦克风阵列的话音定位问题中,输入信号为远场和近场混合信号,需要分别估计出远场和近场信号的一维波达方向。尽管已经有很多用于解决远场和近场混合信号源定位问题的方法被提出,如基于高阶统计量的估计方法、基于矩阵差分的估计方法、基于二阶统计量(SOS)的估计方法等。但都存在一定的缺陷和不足,基于高阶统计量的估计方法需要较高的计算复杂度;基于矩阵差分的估计方法通过消除混合信号中的远场信号和噪声,进而估计出近场信号波达方向,此种方法受采样次数的影响较大,仅在大采样数据情况下较为有效;基于SOS的估计方法虽然计算复杂度较低,但是该方法依据主观准则来区分远场和近场信号类型,有一定主观性。
    发明内容
    本发明的目的在于提供一种计算有效的远场和近场混合信号源的定位方法。
    为达到上述目的,本发明采用了以下技术方案。
    下文中,对于任意变量a,表示该变量a的估计值。
    利用消除噪声后的混合信号估计远场信号的波达方向角,利用斜投影算子从消除噪声后的混合信号中分离出远场信号,得到仅包含近场信号信息的协方差矩阵,根据仅包含近场信号信息的协方差矩阵估计出近场信号的波达方向角及距离,所述混合信号为入射到对称均匀线阵上的K个非相干窄带信号{sk(n)},所述对称均匀线阵包含2M+1个全向传感器阵元,M的取值范围为M≥K,阵元间距为d,设前K1个入射信号为远场信号,远场信号的方位信息为后K2个入射信号为近场信号,近场信号的方位信息为K=K1+K2,θk表示第k个入射信号的波达方向角,所述波达方向角为第k个入射信号相对于y轴的逆时针夹角,rk是第k个入射信号相对于坐标原点的距离。
    所述远场信号的波达方向角的估计方法包括以下步骤:
    1)消除阵列协方差矩阵R中的噪声部分,得到无噪协方差矩阵将分为两部分:

    2)构造代价函数
    gf(θ)=afH(θ)ΠQ^af(θ)---(1)]]>
    式(1)中(·)H表示共轭转置;
    ΠQ^=Q^(Q^HQ^)-1Q^H=Q^(I2M+1-K-P^H(P^P^H+IK)-1P^)Q^H]]>
    Q^=[P^T,-I2M+1-K]T]]>
    P^=(G^HG^)-1G^HH^]]>
    (·)T表示转置,Im表示m×m的单位矩阵;
    3)根据式(1)构造多项式其中p(z)=Δ[zM,...,z,1,z-1,...,z-M]T,z=Δej2πdsinθ/λ,]]>λ表示入射信号波长,j表示单位虚数,j2=-1,通过求多项式的K1个零相位点来估计远场信号的波达方向角。
    所述步骤1)具体包括以下步骤:
    a、根据对称均匀线阵接收的数据求得阵列协方差矩阵R的估计值:
    R=1NΣn=1Nx(n)xH(n)---(2)]]>
    其中,N表示采样数,x(n)表示阵列接收数据;
    b、将阵列协方差矩阵R划分为如下形式

    c、估计噪声方差σ2
    σ2=tr{R22Π}tr{Π}---(4)]]>
    其中tr{·}表示求矩阵的迹,表示广义逆;
    d、无噪协方差矩阵
    所述近场信号的波达方向角及距离的估计方法包括以下步骤:
    2.1)估计斜投影算子
    2.2)消除无噪协方差矩阵中的远场信号,得到仅包含近场信号信息的协方差矩阵Rn,Rn=(I2M+1-EAf|An)R‾(I2M+1-EAf|An)H;]]>
    2.3)利用仅包含近场信号信息的协方差矩阵Rn的反对角元素信息构造一个新的协方差矩阵并将新构造的协方差矩阵分块:
    2.4)构造代价函数
    fn=a‾nH(θ)ΠQ‾^a‾n(θ)---(12)]]>
    式(12)中,ΠQ‾^=Q‾^(Q‾^HQ‾^)-1Q‾^H=Q‾^(IM+1-K2-P‾^H(P‾^P‾^H+IK2)-1P‾^)Q‾^H,P‾^=(G‾^HG‾^)-1G‾^HH‾^,]]>Q‾^=[P‾^T,-IM+1-K2]T;]]>
    2.5)根据式(12)构造多项式pN(z)=ΔzMpH(z)ΠQ‾^p(z),]]>其中p(z)=Δ[1,z-1,...,z-M]T,]]>z=Δej4πdsinθ/λ,]]>通过求多项式pN(z)=ΔzMpH(z)ΠQ‾^p(z)]]>的K2个零相位点来估计近场信号的波达方向角;
    2.6)构造代价函数
    f‾N(r)=anH(r,θ^)ΠQ^an(r,θ^)=μnH(θ^,r)DH(θ^)ΠQ^D(θ^)μn(θ^,r)---(13)]]>
    其中μn(θ^,r)=[eM2φ,e(M-1)2φ,...,1]T,φ=Δπd2cos2θ^/(λr),D(θ^)=Δ[Φ(θ^),Φ*(θ^)J]T,]]>是块对角矩阵,Φ(θ^)=diag{e-jω^M,e-jω^(M-1),...,1},]]>(·)*表示取共轭,ω^=Δ-2πdsinθ^/λ,]]>J是一个(M+1)×M的选择矩阵,定义为J=[eM,eM-1,…,e1],em表示单位矩阵IM+1的第m列,m=1,2,…,M;
    2.7)根据式(13)构造多项式pr(z)=ΔzM2pH(z)DH(θ^)ΠQ^D(θ^)p(z),]]>其中p(z)=Δ[1,z,...,zM2]T,z=Δeπd2cos2θ^/(λr),]]>通过求多项式pr(z)=ΔzM2pH(z)DH(θ^)ΠQ^D(θ^)p(z)]]>的 K2个零相位点来估计近场信号的距离;
    2.8)更新斜投影算子返回步骤2.2),继续顺序执行,如此往复循环,直至满足精度要求时结束。
    所述步骤2.1)中斜投影算子的估计方法包括以下步骤:
    a、构造一个(2M+1)×(2M+1)矩阵:
    Ro=ΔR‾ΠAf---(14)]]>
    其中Af是远场信号方向矩阵,由估计出的远场信号的波达方向角计算得出;
    b、对RoΔ进行QR分解得斜投影算子由下述方式获得

    其中是矩阵的后2M+1-K2列,Δ是(2M+1)×(2M+1)的置换矩阵,不改变矩阵Ro的列相关性,表示广义逆。
    所述步骤2.3)具体包括以下步骤:
    仅包含近场信号信息的协方差矩阵Rn的反对角元素rn(i)定义为
    rn(i)=Δ(Rn)i,2M+2-i=Σk=K1+1Krsk2e-j2(M+1-i)ωk---(16)]]>
    i=1,2,…,2M+1,是第k个近场信号的功率;

    更新斜投影算子的方法为:
    EAf|An=Af(AfHΠAnAf)-1AfHΠAn---(25)]]>
    Af是远场信号方向矩阵,由估计出的远场信号的波达方向角计算得出。
    本发明的有益效果是:
    本发明对远场和近场信号方位信息分别进行估计,利用斜投影算子,有效地从混合信号中分离出远场信号,得到仅包含近场信号信息的传感器协方差矩阵,进而计算出近场信号的波达方向角及距离,相比于已有的远场和近场混合信号源定位算法,本发明计算复杂度低,方法简单有效,若以matlab每秒浮点运算次数来描述算法复杂度,本发明的复杂度为每秒运算10(2M+1)2+30(2M+1)3次浮点数。
    另外,本发明提出了一种可以更有效地计算斜投影算子的方法,并引入交替迭代方法来消除近场信源定位中出现的饱和现象。
    附图说明
    图1为计算复杂度随采样数及阵元个数变化曲线:(a)为随采样数,阵元个数设为7,(b)为随阵元数,采样数设为100;虚线:LOFNS(无迭代);实线:EAMA;实线带有“×”:MSLA;虚线加“o”:TSA;点划线:EMPL。
    图2为波达方向角的估计均方根误差随信噪比变化曲线:(a)为远场信号,(b)为近场信号;点线加“Δ”:LOFNS;虚线:LOFNS(无迭代);实线:EAMA;实线带有“×”:MSLA;虚线加“o”:TSA;点划线:EMPL;点线:CRB。
    图3为波达方向角的估计均方根误差随采样数变化曲线:(a)为远场信号,(b)为近场信号;点线加“Δ”:LOFNS;虚线:LOFNS(无迭代);实线:EAMA;实线带有“×”:MSLA;虚线加“o”:TSA;点划线:EMPL;点线:CRB。
    图4为近场信号距离的估计均方根误差随信噪比及采样数的变化曲线:(a)为随信噪比,(b)为随采样数;点线加“Δ”:LOFNS;虚线:LOFNS(无迭代);实线:EAMA;实线带有“×”:MSLA;虚线加“o”:TSA;点划线:EMPL;点线:CRB。
    图5为阵列结构图。
    具体实施方式
    下面结合附图和实施例对本发明做详细描述。
    一种远场和近场混合信号源的定位方法,具体实现步骤概括如下:
    1)计算阵列协方差矩阵R的估计值
    2)消除阵列协方差矩阵中的噪声部分,得到无噪协方差矩阵
    3)由无噪协方差矩阵构造代价函数,通过求解优化问题估计出远场信号的波达方 向角{θk}k=1K1;]]>
    4)计算斜投影算子的初始值;
    5)利用斜投影算子消除无噪协方差矩阵中的远场信号部分,得到仅包含近场信息的新的协方差矩阵
    6)由新的协方差矩阵构造代价函数,通过求解优化问题估计出近场信号的波达方向角{θk}k=K1+1K]]>及距离{rk}k=K1+1K;]]>
    7)将上述求得的远场信号和近场信号的波达方向角以及近场信号的距离作为初始值,重新估计斜投影算子,返回步骤5),反复迭代,直到满足条件。
    下面进行具体描述。
    K个非相干窄带信号{sk(n)}入射到对称均匀线阵ULA上,该对称均匀线阵包含2M+1个全向传感器阵元,M的取值范围为M≥K,阵元间距为d。设前K1个入射信号为远场信号,其方位信息为后K2个入射信号为近场信号,其方位信息为并且K=K1+K2。θk表示第k个入射信号(远场信号或近场信号)相对于y轴的逆时针夹角(波达方向角),rk是第k个入射信号(仅近场信号)相对于坐标原点的距离,参见图5。
    令对称均匀线阵ULA的中心为参考阵元,阵列输出信号为
    x(n)=Afsf(n)+Ansn(n)+w(n)
                       (1)
    =As(n)+w(n)
    其中,Α是阵列响应矩阵,A=Δ[Af,An],Af=Δ[af(θ1),af(θ2),...,af(θK1)],]]>An=Δ[an(rK1+1,θK1+1),an(rK1+2,θK1+2),...,an(rK,θK)],]]>远场和近场的方向向量分别定义为af(θk)=Δ[e-jMωk,...,e-jωk,1,ejωk,...,ejMωk]T,an(rk,θk)=Δ[ejτ-Mk,...,ejτ-k,1,ejτk,...,ejτMk]T,]]>(·)T表示转置。远场信号的相位延迟ωk定义为λ是入射信号波长,j表示单位虚数,j2=-1。当近场距离满足rf∈(0.62(D3/λ)1/2,2D2/λ),其中D是阵列孔径,根据菲涅尔近似,近场信号的相位延迟τmk定义为其中 m=-M,…,-1,0,1,…,M,φk=Δπd2cos2θk/(λrk).]]>
    接收数据的阵列协方差矩阵为
    R=ΔE{x(n)xH(n)}=AfRsfAfH+AnRsnAnH+σ2I2M+1=R‾+σ2I2M+1---(2)]]>
    其中是无噪协方差矩阵,定义为Rsf和Rsn分别是远场和近场信号协方差矩阵,定义为Rsf=ΔE{sf(n)sfH(n)},Rsn=ΔE{sn(n)snH(n)}.]]>另外E{w(n)wH(n)}=σ2I2M+1,blkdiag(·)表示构建分块对角矩阵,(·)H表示矩阵共轭转置,Im是一个m×m单位阵。
    下面描述远场信号波达方向角的估计方法,包括以下步骤:
    1)消除阵列协方差矩阵R中的噪声部分,得到无噪协方差矩阵
    2)无噪协方差矩阵可分为两部分,
    3)构造代价函数
    gf(θ)=afH(θ)ΠQ^af(θ)---(3)]]>
    4)根据上式构造多项式其中p(z)=Δ[zM,...,z,1,z-1,...,z-M]T,z=Δej2πdsinθ/λ,]]>通过求该多项式的K1个零相位点来估计远场信号的波达方向角θk;
    所述步骤1)消除阵列协方差矩阵R中的噪声部分,得到无噪协方差矩阵的过程包括:
    a、根据阵列(对称均匀线阵)接收数据求得阵列协方差矩阵R的估计值:
    R=1NΣn=1Nx(n)xH(n)---(4)]]>
    其中,N表示采样数,x(n)表示阵列接收数据;
    b、阵列协方差矩阵R可以划分为如下形式

    c、估计噪声方差σ2
    σ2=tr{R22Π}tr{Π}---(6)]]>
    其中tr{·}表示求矩阵的迹,表示广义逆;
    d、无噪协方差矩阵
    所述步骤3)中
    ΠQ^=Q^(Q^HQ^)-1Q^H=Q^(I2M+1-K-P^H)(P^P^H+IK)-1P^)Q^H---(7)]]>
    Q^=[P^T,-I2M+1-K]T---(8)]]>
    P^=(G^HG^)-1G^HH^---(9)]]>
    是一个线性变化矩阵,阵列响应矩阵A可以被分为两个子阵,如下

    基于A是满秩矩阵这一基本假设,矩阵A的前K行线性独立,因此,存在一个K×(2M+1-K)的线性变化矩阵P,使得
    A2=PHA1,i.e.QHA=O(2M+1-K)×K  (11)
    其中Op×q是一个p×q的零矩阵,由(11)式可以很容易得到
    QHan(r,θ)=0(2M+1-K)×1  (12)
    QHaf(θ)=0(2M+1-K)×1  (13)
    下面描述近场信号波达方向角及距离估计。
    1)估计斜投影算子
    2)消除无噪协方差矩阵中的远场信号,得到仅包含近场信号信息的协方差矩阵Rn,Rn=(I2M+1-EAf|An)R‾(I2M+1-EAf|An)H;]]>
    3)利用仅包含近场信号信息的协方差矩阵Rn的反对角元素信息,构造新的协方差矩阵
    4)新的协方差矩阵可分为两部分:
    5)构造代价函数
    fn=a‾nH(θ)ΠQ‾^a‾n(θ)---(14)]]>
    6)根据上式构造多项式pN(z)=ΔzMpH(z)ΠQ‾^p(z),]]>其中p(z)=Δ[1,z-1,...,z-M]T,]]>通过求该多项式的K2个零相位点来估计近场信号的波达方向角θk;
    7)构造代价函数
    f‾N(r)=anH(r,θ^)ΠQ^an(r,θ^)=μnH(θ^,r)DH(θ^)ΠQ^D(θ^)μn(θ^,r)---(15)]]>
    其中为步骤6)估计的近场信号波达方向角,φ=Δπd2cos2θ^/(λr),]]>定义为D(θ^)=Δ[Φ(θ^),Φ*(θ^)J]T,]]>是块对角矩阵,Φ(θ^)=diag{e-jω^M,e-jω^(M-1),...,1},ω^=Δ-2πdsinθ^/λ,]]>J是一个(M+1)×M的选择矩阵,定义为J=[eM,eM-1,…,e1],em表示单位矩阵IM+1的第m列,m=1,2,…,M;
    8)根据上式构造多项式其中p(z)=Δ[1,z,...,zM2]T,z=Δeπd2cos2θ^/(λr),]]>通过求该多项式的K2个零相位点来估计近场信号的距离;
    9)更新斜投影算子返回步骤2),继续顺序执行,如此往复循环,直至满足精度要求时结束。
    所述步骤1)中斜投影算子的估计方法包括以下步骤:
    a、构造一个(2M+1)×(2M+1)矩阵如下
    Ro=ΔR‾ΠAf=AnRsnAnHΠAf---(16)]]>
    其中Af是远场信号方向矩阵,由以上估计出的远场信号波达方向角计算得出;
    b、对RoΔ进行QR分解得斜投影算子可以由下述方式获得

    其中是矩阵的后2M+1-K2列,Δ是(2M+1)×(2M+1)的置换矩阵,不改变矩阵Ro的列相关性。
    所述步骤3)具体叙述如下:
    仅包含近场信号信息的协方差矩阵Rn的反对角元素rn(i)定义为
    rn(i)=Δ(Rn)i,2M+2-i=Σk=K1+1Krsk2e-j2(M+1-i)ωk---(18)]]>
    i可取i=1,2,…,2M+1,是第k个近场信号的功率,定义为(·)*表示取共轭。由(18)式可以看到反对角元素rn(i)只包含近场信号波达方向角度信息。类似于远场Toepliz协方差矩阵,近场信号协方差矩阵如下

    将(18)式带入(19)式,经过简单的代数运算,重写协方差矩阵如下式
    R‾n=A‾n(θ)RsnA‾nH(θ)---(20)]]>
    其中A‾n(θ)=Δ[a‾n(θ1),a‾n(θ2),...,a‾n(θK1)],a‾n(θk)=Δ[1,ej2ωk,ej4ωk,...,ej2Mωk]T.]]>
    所述步骤5)中
    ΠQ^=Q^(Q^HQ^)-1Q^H=Q^(I2M+1-K-P^H(P^P^H+IK)-1P^)Q^H---(21)]]>
    Q^=[P^T,-I2M+1-K]T---(22)]]>
    P‾^=(G‾^HG‾^)-1G‾^HH‾^---(23)]]>
    是一个线性变化矩阵,可以被分为两个子阵,如下

    其中是一个满秩矩阵,的行向量可以由的行向量线性表示,因此,存在一个K2×(M+1-K2)的线性算子使得
    A‾n2(θ)=P‾HA‾n1(θ),i.e.Q‾HA‾n(θ)=O(M+1-K2)×K2---(25)]]>
    其中Op×q是一个p×q的零矩阵。
    在步骤8)完成后更新斜投影算子
    EAf|An=Af(AfHΠAnAf)-1AfHΠAn---(26)]]>
    更新斜投影算子的目的在于得到更为精确的的估计值。步骤1)中的估计值有一定误差,因为该估计过程需要利用信号协方差矩阵Rs的块对角结构和仅包含近场信号信息的协方差矩阵Rn的反对角元素。然而,在采样数有限的情况下,Rs并非块对角结构,Rn的反对角元素也并非仅包含近场信号信息,可以由下式看到
    R^o=AnR^snAnHΠAf+AfR^sfnAnHΠAf---(27)]]>
    其中由于的存在,的列空间与An的列空间不完全相等。因此步骤1)中斜投影算子的估计并不准确,并且不受叠加噪声的影响,这意味着尽管信噪比趋于无穷大,的影响仍然存在,估计误差并不随信噪比SNR的增加而减小,称为“饱和现象”。而利用第一次计算得到的远场信号和近场信号波达方向角以及近场信号距离,重新计算斜投影算子,可以进一步优化结果,并且通过迭代,消除“饱和现象“。
    下面通过以下不同情形对上述方法的效果进行说明:
    空间有两个波达方向未知的入射信号,其方向信息分别为(∞,-5°),(1.7λ,25°),对称均匀线阵含有7个阵元,阵元间隔为d=λ/4。仿真中加入基于矩阵差分的算法如EMPL、TSA,基于SOS的算法如EAMA、MSLA与LOFNS对比,同时给出了CRB界。每一个仿真结果都是经由1000次独立重复实验得到的。
    由图1说明LOFNS计算复杂度最低。
    由图2及图4(a)说明对于远场信号波达方向角的估计,LOFNS接近EAMA且优于EAMA,特别是在高信噪比情况下。并且通过迭代,消除了饱和现象。
    由图3及图4(b)说明LOFNS在远场信号和近场信号的估计性能都优于其他方法。

    关 键 词:
    一种 近场 混合 信号源 定位 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:一种远场和近场混合信号源的定位方法.pdf
    链接地址://www.4mum.com.cn/p-6143077.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