专利名称：：强弱信号的波达方向与信源数估计方法
技术领域：
：本发明属于信号处理
技术领域：
，涉及空间谱估计，可用于在有多个强弱信号时，同时估计出强弱信号的波达方向与信号源数目，以进一步提取信号的其他信息或抑制干扰。
背景技术：
：波达方向估计是阵列信号处理的一个重要分支，其应用涉及雷达、通信、导航、测控和电子侦察等众多领域。基于空间谱的波达方向估计方法，如最大似然(ML:Maximumlikelihood)、多重信号分类(MUSIC:MultipleSignalClassification)、旋转不变信号参数估计技术彼PRIT-EstimationofSignalparametersViaRotationalInvarianceTechnique)以及加权子空间拟合(WSF:WeightSubspaceFitting)等方法，由于具有对多信源的测角超分辨能力引起了人们的广泛关注。但是，同时存在功率相差很大的强弱信号时，由于强信号的影响，利用上述方法难以准确估计出弱信号的波达方向。针对强信号背景下弱信号波达方向估计问题，目前已有的方法主要有1.J丄i等发表于1997年IEEETransactiononAerospaceandElectronicSystems上的《AngleandWaveformEstimationViaRELAX》一文中，提出通过迭代方式来估计信号的到达角和波形，此方法中需要反复迭代，因此运算量相当大。2.陈辉，苏海军发表于2006年电子学报上的《强干扰/信号背景下的D0A估计新方法》中提出的JJM方法，利用阻塞矩阵抑制掉强信号，然后再估计出弱信号，但此方法要求强信号的波达方向精确已知，这在实际中很难满足；而且在强弱信号角度很近时，抑制掉强信号的同时也会影响到弱信号。3.张静在2009年硕士学位论文《阵列天线子阵划分设计和测向技术研究》中，提出用特征矢量矩阵对接收数据做变换，然后对变换后各通道数据用最大模值归一化，最后再用MUSIC方法实现强弱信号到达角估计。但该方法需要付出多做一次特征分解和搜索各通道最大值的运算量，且测向性能受数据段最大值波动的影响。以上三种方法都是在假定信源数已知的条件下提出的，但在有强弱信号同时存在时，信源数也很难估计，若估计不准会直接导致波达方向的漏估与过估，影响其他信息的正确提取。
发明内容本发明的目的在于克服上述已有技术的不足，提出了一种强弱信号的波达方向与信源数估计方法，以在信源数未知的情况下，有效地估计出强弱信号的波达方向和信源数，有利于对其他信息的正确提取。为实现上述目的，本发明的估计方法包括如下步骤(1)根据阵列接收数据X(t)，估计其协方差矩阵^^为<formula>formulaseeoriginaldocumentpage3</formula>其中，N为采样快拍数，H表示共轭转置运算；(2)对协方差矩阵Ax做特征分解，得到特征值与相应的特征矢量[U,A]=e/g(Rx)其中，eig()表示特征分解，A为特征值由大到小排列组成的对角矩阵，U=[Ul，u2，，uM]为特征矢量矩阵，Ui(i=1，2，…，M)为相应的特征矢量，M为阵元个数；(3)初始化，i=1;(4)用第i个特征矢量Ui计算第i个特征波束的空间谱Pi(9):P,(6>)=H、"、,aH(0)R;'a(0)其中a(9)为搜索导向矢量，R^为^的逆矩阵；(5)将第i个特征波束的空间谱Pi(e)化为单位为dB的形式PidB(e)，求PidB(9)的最大值Pi—dB^(e)及对应的角度^inlax;求主瓣波束宽度外Pi—dB(9)的平均值^—dB;若P,j^(0)-》,—必，则第i个信号的波达方向为e^，转到步骤(4)，估计第i+1个特征波束的空间谱pi+1(e)，否则，信源数等于i-i，估计结束，其中il为所设定的判断有无信号的门限值。本发明与现有技术相比，具有以下优点1)本发明由于通过比较每个特征波束空间谱的最大值与主瓣波束宽度外空间谱的平均值，能同时估计波达方向与信源数，不需要信源数已知。2)本发明由于在计算特征波束的空间谱时，分子表达式计算简单，分母表达式在所有特征波束空间谱中都相同，能预先计算好再调用，故运算量小。3)本发明相对于现有技术不需要迭代运算，无需强信号波达方向的先验知识。对本发明的目的、特征、优点可通过如下附图和实例详细描述。图l是本发明的流程图；图2是本发明仿真实验所用的天线阵列排布图；图3是本发明的实验场景图；图4是用本发明仿真实验获得的第1个特征波束的空间il图5是用本发明仿真实验获得的第2个特征波束的空间il图6是用本发明仿真实验获得的第3个特征波束的空间il图7是用本发明仿真实验获得的第4个特征波束的空间il具体实施例方式参照图l，本发明的实现步骤如下步骤l，根据天线阵列的接收数据x(t)，估计其协方差矩阵A,为<formula>formulaseeoriginaldocumentpage5</formula><formula>formulaseeoriginaldocumentpage5</formula>其中，t二[1，2，…，N]，N为采样快拍数，H表示共轭转置运算。步骤2:对协方差矩阵Ax进行特征分解，得到特征值与相应的特征矢量[U,A]=e,g(Ax)其中，eig。)表示特征分解，A为特征值由大到小排列组成的对角矩阵，U二[Ul，u2，…，Um]为特征矢量矩阵，Ui(i=1，2，…，M)为相应的特征矢量，M为阵元个数。步骤3:从第一个特征矢量开始，依次估计强弱信号的波达方向，即初始化，i=1。步骤4:用第i个特征矢量Ui计算第i个特征波束的空间谱Pi(9):一、aH(0)u,u"a(。P,(6>)=、",'aH(0)R;'a(0)其中a(9)为搜索导向矢量，6^为Rj勺逆矩阵；分母项aH(6l)&:'a(60与Ui无关，只需在估计第1个特征波束的空间谱P工(e)时计算，计算后存起来，估计第i(l<i<M)个特征波束时直接调用即可。步骤5:估计波达方向与信源数5.l)将第i个特征波束的空间谱Pi(e)化为dB形式Pi亚(e);5.2)求P^(e)的最大值记为Pi^M(e)，最大值Pi^x(e)对应的角度记为9.inmx5.3)求主瓣波束宽度外PidB(e)的平均值记为^j;5.4)设定判断有无信号的门限值IL，本实例的门限值IL设为13.4dB;5.5)判断最大值PidBmax(e)与平均值gdB之差和门限值IL的大小关系，若必HL，则第i个信号的波达方向估计为ei隨，转到步骤4估计Pw(e)，若P,_^max(。-A—必〈IL，则信源数估计为i-l，估计结束。本发明的效果可以通过以下实测数据处理结果进一步说明。1.实验环境天线阵列的排布图如图2所示，24个阵元在yz平面上分两行均匀排列，每行是由12个阵元组成的均匀线阵，阵元的行距和列矩都为80mm，波长设为154mm，由于实验中强弱信号和接受阵列基本处于同一水平面，可不考虑俯仰角，在这里只估计方位角。采用两个螺旋天线在同一个水平面内发射一强一弱两个信号，强信号真实方位角为22°，弱信号真实方位角为-2°，由于周围停放一辆汽车，正好能将强信号的多径反射信号被阵列接收，此信号在传播中有一定的功率损耗，相对原来的强信号来说是弱信号，即共有三个功率不等的强弱信号，实验的场景如图3所示。固定弱信号的发射功率，逐渐增大强信号的发射功率，得到几组阵列接收数据，用本发明进行仿真这些数据，同时估计强弱信号的方位角和信源数。2实验结果(1)对其中一组观察数据进行处理，图4为本发明仿真实验获得的第l个特征波束的空间谱，该空间谱满足&—必，(60-A—必，则其最大值P工9)对应的角度21.8°即为强信号的波达方向；图5为本发明仿真实验获得的第2个特征波束的空间谱，同理，其最大值对应的角度52.4。即为强信号的多径信号的波达方向；图6为本发明仿真实验获得的第3个特征波束的空间谱，同理，其最大值对应的角度-2.3°即为弱信号的波达方向；图7为本发明仿真实验获得的第4个特征波束的空间谱，由图7可看出，最大值P4dBmax(9)与周围角度对应的空间谱没有明显的差异，即&jmax(。-Aj<IL，所以视此波束没有信号，估计出信源数为3，估计结果与实验条件相符合。(2)强弱信号的功率差由17dB变换到29dB，估计如下四组数据强弱信号的方位角和信源数，信源数的估计结果都为3，方位角的估计结果如表1所示。表l不同强弱信号功率差时各个信号的方位角估计结果<table>tableseeoriginaldocumentpage6</column></row><table>由表1可以看出，改变强弱信号之间的功率差，估计结果稳定。综上，在有多个强弱信号时，本发明方法能准确地同时估计出强弱信号的波达方向与信源数。权利要求一种强弱信号的波达方向与信源数估计方法，包括如下步骤(1)根据阵列接收数据x(t)，估计其协方差矩阵为<mrow><msub><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>x</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>其中，N为采样快拍数，H表示共轭转置运算；(2)对协方差矩阵做特征分解，得到特征值与相应的特征矢量<mrow><mo>[</mo><mi>U</mi><mo>,</mo><mi>Λ</mi><mo>]</mo><mo>=</mo><mi>eig</mi><mrow><mo>(</mo><msub><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi></msub><mo>)</mo></mrow></mrow>其中，eig(·)表示特征分解，Λ为特征值由大到小排列组成的对角矩阵，U＝[u1，u2，…，uM]为特征矢量矩阵，ui(i＝1，2，…，M)为相应的特征矢量，M为阵元个数；(3)初始化，i＝1；(4)用第i个特征矢量ui计算第i个特征波束的空间谱Pi(θ)<mrow><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><msub><mi>u</mi><mi>i</mi></msub><msubsup><mi>u</mi><mi>i</mi><mi>H</mi></msubsup><mi>a</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>a</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow></mrow></mfrac></mrow>其中a(θ)为搜索导向矢量，为的逆矩阵；(5)将第i个特征波束的空间谱Pi(θ)化为单位为dB的形式Pi_dB(θ)，求Pi_dB(θ)的最大值Pi_dBmax(θ)及对应的角度θimax；求主瓣波束宽度外Pi_dB(θ)的平均值Pi_dB；若<mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>_</mo><mi>dB</mi><mi>max</mi></mrow></msub><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>P</mi><mo>‾</mo></mover><mrow><mi>i</mi><mo>_</mo><mi>dB</mi></mrow></msub><mo>≥</mo><mi>IL</mi><mo>,</mo></mrow>则第i个信号的波达方向为θimax，转到步骤(4)，估计第i+1个特征波束的空间谱Pi+1(θ)，否则，信源数等于i-1，估计结束，其中IL为所设定的判断有无信号的门限值。FSA00000050717600011.tif,FSA00000050717600013.tif,FSA00000050717600016.tif,FSA00000050717600017.tif2.根据权利要求1中所述的一种强弱信号的波达方向与信源数估计方法，其特征在于判断有无信号的门限值IL，设定为13.4dB。全文摘要本发明公开了一种强弱信号的波达方向与信源数估计方法，主要解决现有方法在信源数未知时不能正确估计强弱信号的波达方向的问题。其实现过程为根据阵列接收数据估计其协方差矩阵；对该协方差矩阵进行特征分解得到由大到小排列的特征值及对应的特征矢量；从第1个特征矢量开始，依次计算各个特征波束的空间谱；通过比较该空间谱的最大值与主瓣波束宽度外的平均值之差和所设定的门限值估计各个信号的波达方向及信源数。本方法简单实用，在有多个强弱信号并存时，能精确估计强弱信号的波达方向和信源数，可用于雷达、通信、导航、测控和电子侦察众多领域中信息的提取或干扰的抑制。文档编号H04L1/06GK101795150SQ20101012078公开日2010年8月4日申请日期2010年3月9日优先权日2010年3月9日发明者廖桂生,徐亮,徐青,曾操,李军,杨志伟,陶海红申请人:西安电子科技大学