[0003] 为解决近场源参数估计问题，在高斯噪声环境下，一般采用基于二阶统计特性的线性预测方法及采用高阶累积量的特征分析方法。上述方法在一定程度上提高了参数的估计精度和复杂度。文献《冲击噪声背景下的近场源二维参数估计方法》[J]，王波，王树勋，电路与系统学报，2005，10(5)：5-9.利用均匀线阵为接收阵列研究了冲击噪声下的近场源二维参数估计问题(方位角和距离)，但该方法在均匀圆阵接收阵列下，无法直接扩展。并且在冲击噪声环境下时，由于SαS噪声当且仅当0
发明内容
[0004] 本发明的目的是针对背景技术的不足之处，改进设计一种冲击噪声环境下的近场源角度和距离计算方法，从而达到目标定位过程中速度快、精度高、运算量小、成本低的目的。
[0005] 本发明的技术方案是一种冲击噪声环境下的近场源角度和距离估计方法，该方法包括：
[0006] 步骤1：在定位环境中规划一个圆形区域，将各天线等角度设置于该圆形区域的边缘上；
[0007] 步骤2:将信号源置于定位环境中,各天线接收到信号源发射的信号并存储；
[0008] 步骤3：将存储的各天线数据采用空间符号函数处理；
[0009] 步骤4：使用步骤3处理后得到的数据采用公式：
[0010]
[0011] 求得第一类协变异矩阵R1(k)，
[0012] 其中，()*表示向量共轭，N表示采样次数，Sk(t)表示第k个阵元的第t次采样，k＝1,2,…,M/2，M表示接收天线总个数，式中 可以为任何非零值；
[0013] 步骤5：使用步骤3处理后得到的数据采用公式：
[0014]
[0015] 求得第二类协变异矩阵R2(k)，其中，k＝1,2,…,3M/4；
[0016] 步骤6：对协变异矩阵R1(k)取相角ωk，采用该相角计算出信号源的波达方向；
[0017] 步骤7：对协变异矩阵R2(k)取相角uk，使用该相角计算出信号源离基站的距离。
[0018] 所述步骤1中接收天线个数一般为8至10个。
[0019] 所述步骤3中的空间符号函数是：
[0020]
[0021] 其中：x(t)表示圆阵所有天线第t次采样的接收信号向量，||·||代表二范数。
[0022] 所述步骤6中对协变异矩阵R1(k)取相角ωk，其中k＝1,2,3,…,M/2，L表示均匀圆阵半径，λ表示波长，θ表示信号入射的方位角变量，γk＝2π(k-1)/M表示第k个阵元的方位角，φ表示信号入射的俯仰角变量，mk表示为整数，[0023] 假设mk＝0，则ωk可以表达为
[0024] 矩阵ω≈γb，其中ω＝[ω1 ω2 … ωM/2]T，则由最小二乘可得到[0025] b＝[b1 b2]T＝(γTγ)-1γTω
[0026] 其中
[0027] 利用下式就可得到角度的估计
[0028]
[0029]
[0030] 所述步骤7中对协变异矩阵R2(k)取相角uk，可得其中k＝1,2,3,…,3M/4；
[0031] 将其表达为矩阵形式有， 其中u＝[u1 u2 … u3M/4]T,,
[0032] v＝sin2φ[cos(2θ-2γ1) cos(2θ-2γ2) … cos(2θ-2γ3M/4)]T[0033]
[0034] 由最小二乘可以得到距离r的估计
[0035]
[0036] 其中()H表示共轭转置，()-1代表求逆，v为包含方位角和俯仰角的向量。
[0037] 本发明是一种冲击噪声环境下的近场源角度和距离计算方法，利用圆形接收天线阵列接收到数据后，先计算各数据的空间符号函数，进而计算得到两类协变异矩阵，通过第一类协变异矩阵的相角计算出信号源的波达方向，通过第二类协变异矩阵的相角计算出信号源离基站的距离，从而冲击噪声环境下目标定位过程中具有速度快、精度高、运算量小、成本低的效果。
附图说明
[0038] 图1为本发明的流程图；
[0039] 图2为本发明实施中均匀圆阵和近场信号源之间的关系图；
[0040] 图3为单个近场信号源入射采用二阶统计量和本发明方法的估计俯仰角误差性能比较图；
[0041] 图4为单个近场信号源入射采用二阶统计量和本发明方法的估计方位角误差性能比较图；
[0042] 图5为单个近场信号源入射采用二阶统计量和本发明方法的估计距离误差性能比较图。
具体实施方式
[0043] 下面结合附图对本发明详细说明。
[0044] 如图2所示半径为L的均匀圆阵为例，阵列具有M＝8个阵元。考虑一个窄带信号s(t)，到达均匀圆阵，与X轴正向夹角为θ，与Z轴正向夹角为 。以均匀圆阵的中心为相位参考点。在t时刻阵列的接收信号向量为
[0045] x(t)＝As(t)+n(t) (1)
[0046] 其中，A＝[a1(r,θ,φ) a2(r,θ,φ) … aM(r,θ,φ)]T，[]T表示向量转置，ω＝2π/λ，λ是信号的波长，n(t)是与信号s(t)独立的冲击噪声。dl(r，θ，φ)为信号源与第l个阵元之间的距离，
[0047]
[0048] 其中，
[0049] ρl(θ，φ)＝sinφcos(θ-(l-1)θ0)，l＝1,…,M,
[0050] 根据泰勒级数展开，可将dl(r，θ，φ)表示为
[0051]
[0052] 则导向矢量可表示为
[0053]
[0054] 用空间符号函数处理接收数据
[0055] 用处理后的接收数据S(t)确定协变异矩阵R1(k)为：
[0056]
[0057] 对协变异矩阵R1(k)取相角为
[0058] 其矩阵形式为：
[0059] 将ωk表达为矩阵形式ω≈γb，则由最小二乘可得到
[0060] b＝[b1 b2]T＝(γTγ)-1γTω
[0061] 其中
[0062] 利用下式就可得到角度的估计
[0063]
[0064]
[0065] 用处理后的接收数据S(t)确定协变异矩阵R2(k)为：
[0066]
[0067] 对协变异矩阵R2(k)取相角uk，可得
[0068] 将其表达为矩阵形式有， 其中u＝[u1 u2 … u3M/4]T,
[0069] v＝sin2φ[cos(2θ-2γ1) cos(2θ-2γ2) … cos(2θ-2γ3M/4)]T,[0070]
[0071] 由最小二乘可以得到距离r的估计
[0072]
[0073] 参见图3、图4、图5。在阵元数为8，方位角为70°，俯仰角为10°，冲击噪声的参数分别为α＝1.3，γ＝1，δ＝0。信噪比(dB)从0到25的情况下，与传统的二阶统计量估计相比，本发明方法(SCM方法)的测向和距离估计方法具有更高的估计精度。从图中可以看出本发明方法受信噪比的影响较小，估计性能较为稳定。
