技术领域
[0001] 本
发明属于
信号处理技术领域，涉及空间谱估计，可用于在有多个强弱信号时，同时估计出强弱信号的波达方向与信号源数目，以进一步提取信号的其他信息或抑制干扰。
背景技术
[0002] 波达方向估计是阵列信号处理的一个重要分支，其应用涉及雷达、通信、导航、测控和
电子侦察等众多领域。基于空间谱的波达方向估计方法，如最大似然(ML：Maximum likelihood)、多重信号分类(MUSIC：Multiple Signal Classification)、旋转不变信号参数估计技术(ESPRIT：Estimation of Signal parameters Via RotationalInvariance Technique)以及加权子空间拟合(WSF：Weight Subspace Fitting)等方法，由于具有对多信源的测
角超分辨能
力引起了人们的广泛关注。但是，同时存在功率相差很大的强弱信号时，由于强信号的影响，利用上述方法难以准确估计出弱信号的波达方向。
[0003] 针对强信号背景下弱信号波达方向估计问题，目前已有的方法主要有：
[0004] 1.J.Li等发表于1997年IEEE Transaction on Aerospace and Electronic Systems上的《Angle and Waveform Estimation Via RELAX》一文中，提出通过
迭代方式来估计信号的到达角和
波形，此方法中需要反复迭代，因此运算量相当大。
[0005] 2.陈辉，苏海军发表于2006年电子学报上的《强干扰/信号背景下的DOA估计新方法》中提出的JJM方法，利用阻塞矩阵抑制掉强信号，然后再估计出弱信号，但此方法要求强信号的波达方向精确已知，这在实际中很难满足；而且在强弱信号角度很近时，抑制掉强信号的同时也会影响到弱信号。
[0006] 3.张静在2009年硕士学位论文《阵列天线子阵划分设计和测向技术研究》中，提出用特征矢量矩阵对接收数据做变换，然后对变换后各通道数据用最大模值归一化，最后再用MUSIC方法实现强弱信号到达角估计。但该方法需要付出多做一次特征分解和搜索各通道最大值的运算量，且测向性能受数据段最大值
波动的影响。
[0007] 以上三种方法都是在假定信源数已知的条件下提出的，但在有强弱信号同时存在时，信源数也很难估计，若估计不准会直接导致波达方向的漏估与过估，影响其他信息的正确提取。
发明内容
[0008] 本发明的目的在于克服上述已有技术的不足，提出了一种强弱信号的波达方向与信源数估计方法，以在信源数未知的情况下，有效地估计出强弱信号的波达方向和信源数，有利于对其他信息的正确提取。
[0009] 为实现上述目的，本发明的估计方法包括如下步骤：
[0010] (1)根据阵列接收数据x(t)，估计其协方差矩阵 为：
[0011]
[0012] 其中，N为
采样快拍数，H表示共轭转置运算；
[0013] (2)对协方差矩阵 做特征分解，得到特征值与相应的特征矢量：
[0014]
[0015] 其中，eig(·)表示特征分解，Λ为特征值由大到小排列组成的对角矩阵，[0016] U＝[u1，u2，…，uM]为特征矢量矩阵，ui(i＝1，2，…，M)为相应的特征矢量，M为阵元个数；
[0017] (3)初始化，i＝1；
[0018] (4)用第i个特征矢量ui计算第i个特征波束的空间谱Pi(θ)：
[0019]
[0020] 其中：a(θ)为搜索导向矢量，
[0021] 为 的逆矩阵；
[0022] (5)将第i个特征波束的空间谱Pi(θ)化为单位为dB的形式Pi_dB(θ)，求Pi_dB(θ)的最大值Pi_dBmax(θ)及对应的角度θimax；求主瓣波束宽度外Pi_dB(θ)的平均值Pi_dB；若 则第i个信号的波达方向为θimax，转到步骤(4)，估计第i+1个特征波束的空间谱Pi+1(θ)，否则，信源数等于i-1，估计结束，其中IL为所设定的判断有无信号的
门限值。
[0023] 本发明与
现有技术相比，具有以下优点：
[0024] 1)本发明由于通过比较每个特征波束空间谱的最大值与主瓣波束宽度外空间谱的平均值，能同时估计波达方向与信源数，不需要信源数已知。
[0025] 2)本发明由于在计算特征波束的空间谱时，分子表达式计算简单，分母表达式在所有特征波束空间谱中都相同，能预先计算好再调用，故运算量小。
[0026] 3)本发明相对于现有技术不需要迭代运算，无需强信号波达方向的先验知识。
[0027] 对本发明的目的、特征、优点可通过如下
附图和实例详细描述。
附图说明
[0029] 图2是本发明仿真实验所用的天线阵列排布图；
[0030] 图3是本发明的实验场景图；
[0031] 图4是用本发明仿真实验获得的第1个特征波束的空间谱；
[0032] 图5是用本发明仿真实验获得的第2个特征波束的空间谱；
[0033] 图6是用本发明仿真实验获得的第3个特征波束的空间谱；
[0034] 图7是用本发明仿真实验获得的第4个特征波束的空间谱。
具体实施方式
[0035] 参照图1，本发明的实现步骤如下：
[0036] 步骤1，根据天线阵列的接收数据x(t)，估计其协方差矩阵 为：
[0037]
[0038] 其中，t＝[1，2，…，N]，N为采样快拍数，H表示共轭转置运算。
[0039] 步骤2：对协方差矩阵 进行特征分解，得到特征值与相应的特征矢量：
[0040]
[0041] 其中，eig(·)表示特征分解，Λ为特征值由大到小排列组成的对角矩阵，U＝[u1，u2，…，uM]为特征矢量矩阵，ui(i＝1，2，…，M)为相应的特征矢量，M为阵元个数。
[0042] 步骤3：从第一个特征矢量开始，依次估计强弱信号的波达方向，即初始化，i＝1。
[0043] 步骤4：用第i个特征矢量ui计算第i个特征波束的空间谱Pi(θ)：
[0044]
[0045] 其中：a(θ)为搜索导向矢量，
[0046] 为 的逆矩阵；
[0047] 分母项 与ui无关，只需在估计第1个特征波束的空间谱P1(θ)时计算，计算后存起来，估计第i(1＜i＜M)个特征波束时直接调用即可。
[0048] 步骤5：估计波达方向与信源数
[0049] 5.1)将第i个特征波束的空间谱Pi(θ)化为dB形式Pi_dB(θ)；
[0050] 5.2)求Pi_dB(θ)的最大值记为Pi_dBmax(θ)，最大值Pi_dBmax(θ)对应的角度记为θimax；
[0051] 5.3)求主瓣波束宽度外Pi_dB(θ)的平均值记为Pi_dB；
[0052] 5.4)设定判断有无信号的门限值IL，本实例的门限值IL设为13.4dB；
[0053] 5.5)判断最大值Pi_dBmax(θ)与平均值Pi_dB之差和门限值IL的大小关系，若则第i个信号的波达方向估计为θimax，转到步骤4估计Pi+1(θ)，若则信源数估计为i-1，估计结束。
[0054] 本发明的效果可以通过以下实测
数据处理结果进一步说明。
[0055] 1.实验环境
[0056] 天线阵列的排布图如图2所示，24个阵元在yz平面上分两行均匀排列，每行是由12个阵元组成的均匀线阵，阵元的行距和列矩都为80mm，
波长设为154mm，由于实验中强弱信号和接受阵列基本处于同一
水平面，可不考虑
俯仰角，在这里只估计方位角。
[0057] 采用两个螺旋天线在同一个水平面内发射一强一弱两个信号，强信号真实方位角为22°，弱信号真实方位角为-2°，由于周围停放一辆
汽车，正好能将强信号的多径反射信号被阵列接收，此信号在传播中有一定的功率损耗，相对原来的强信号来说是弱信号，即共有三个功率不等的强弱信号，实验的场景如图3所示。
[0058] 固定弱信号的发射功率，逐渐增大强信号的发射功率，得到几组阵列接收数据，用本发明进行仿真这些数据，同时估计强弱信号的方位角和信源数。
[0059] 2实验结果
[0060] (1)对其中一组观察数据进行处理，图4为本发明仿真实验获得的第1个特征波束的空间谱，该空间谱满足 则其最大值P1_dBmax(θ)对应的角度21.8°即为强信号的波达方向；图5为本发明仿真实验获得的第2个特征波束的空间谱，同理，其最大值对应的角度52.4°即为强信号的多径信号的波达方向；图6为本发明仿真实验获得的第3个特征波束的空间谱，同理，其最大值对应的角度-2.3°即为弱信号的波达方向；图
7为本发明仿真实验获得的第4个特征波束的空间谱，由图7可看出，最大值P4_dBmax(θ)与周围角度对应的空间谱没有明显的差异，即 所以视此波束没有信号，
估计出信源数为3，估计结果与实验条件相符合。
[0061] (2)强弱信号的功率差由17dB变换到29dB，估计如下四组数据强弱信号的方位角和信源数，信源数的估计结果都为3，方位角的估计结果如表1所示。
[0062] 表1 不同强弱信号功率差时各个信号的方位角估计结果
[0063]强弱功率差/dB 强信号/° 弱信号/° 强信号多径/°
17 21.7 -2.1 51.8
20 21.9 -2.6 52.3
23 21.9 -2.5 52
29 21.8 -2.3 52.4
[0064] 由表1可以看出，改变强弱信号之间的功率差，估计结果稳定。
[0065] 综上，在有多个强弱信号时，本发明方法能准确地同时估计出强弱信号的波达方向与信源数。
