一种基于线性预测的信源方向鲁棒定位算法的制作方法

1.本发明涉及信号处理技术领域，特别涉及一种基于线性预测的信源方向鲁棒定位算法。


背景技术：

2.发展水声通信技术是国家海洋战略及现代海军建设的需要。由于电磁波在水中衰减严重，高频无线电通信无法适用于水下通信的场景。水下航行器与母船/基地的通信联络可采用3-300hz的超低频频段。实践证明：水声通信的合适载波频率为4khz左右，可见其bw0.707带宽极低，这与人们熟悉的高频无线电通信有本质区别。此外，海洋环境复杂，噪声干扰大。因此，水声通信技术有相当难度。虽然常规的无线通信技术无法直接应用于水声通信，但在原理上是相通的。一般情况下，我们处理的是基带信号。因此，对于水声通信的信号处理，关键在于传统模型的修正。本专利提出一种基于线性预测的信源方向鲁棒定位算法，用于水下目标物体的定位。
3.在脉冲噪声环境中，对k个信源定位，信号模型示意图如附图7：
4.现有n个传感器构成的均匀线性阵列(ula)，用于接收从k个目标发来的波动信号s1(t)，s2(t)，
…
，sk(t)。在t时刻，ula接收到的信号可以建模成：
5.y(t)＝as(t) n(t)，t＝1，2，
…
，t，
6.其中，
7.y(t)＝[y1(t)
ꢀ…ꢀyn
(t)]
t
，
[0008]
s(t)＝[s1(t)
ꢀ…ꢀ
sk(t)]
t
，
[0009]
a＝[a(θ1)
ꢀ…ꢀ
a(θk)]，
[0010][0011]
λ为信号波长，d为阵元间距，θk为第k个目标的方向角，而n(t)＝[n1(t)
ꢀ…ꢀnn
(t)]
t
为t时刻各阵元上的噪声。
[0012]
这里，噪声向量n(t)建模成脉冲噪声(impulsive noise)。
[0013]
现有技术中，信源方向鲁棒定位算法有：1.基于字典的鲁棒定位算法；2.利用观测信号的鲁棒定位算法，包括基于鲁棒协方差矩阵(robust covariation matrix)的rcm-music算法；3.基于分数低阶矩(fractional lower order moment)的flom-music算法，基于符号协方差矩阵(sign covariance matrix)和kendall’stau协方差矩阵(kendall’s tau covariance matrix)的scm-music和tcm-music算法；4.以鲁棒方式观测信号的二阶统计量并定位，包括基于mm-估计器的算法和l
p-music算法等；5.通过最小化观测信号与信号模型之间的拟合误差定位。
[0014]
但是，上述信源方向鲁棒定位算法中：1存在不适用于衰减信号，分辨率受字典密度制约而较低，计算量大的问题；2-5也是计算量较大。这些问题导致相应的算法使用受限。


技术实现要素：

[0015]
针对现实中存在的缺陷，本发明提供一种基于线性预测的信源方向鲁棒定位算法。
[0016]
为了解决上述技术问题，本发明的技术方案为：
[0017]
一种基于线性预测的信源方向鲁棒定位算法，包含以下步骤：
[0018]
首先，构建数据矩阵：
[0019]
y＝x n，
[0020]
其中，y＝[y(1)
ꢀ…ꢀ
y(t)]，x＝a[s(1)
ꢀ…ꢀ
s(t)]和n＝[n(1)
ꢀ…ꢀ
n(t)]分别为干净信号和噪声信号部分。
[0021]
基于干净信号x的线性预测性质，构建观测信号y的线性预测误差矩阵：
[0022]
e＝[e(1)
ꢀ…ꢀ
e(t)]，
[0023]
其中，a＝[a1ꢀ…ꢀak
]
t
为线性预测系数矢量，e(t)＝p(t)a-b(t)，而
[0024]
p(t)＝toeplitz([yk(t)
ꢀ…ꢀyn-1
(t)]
t
，[yk(t)
ꢀ…ꢀ
y1(t)])，
[0025]
b(t)＝-[y
k 1
(t)
ꢀ…ꢀyn
(t)]
t
，
[0026]
表示第一列为c1，第一行为的toeplitz矩阵。
[0027]
通过最小化线性预测误差的加权l
p-范数(1≤p＜2)，估计线性预测系数矢量a，进而估计信源方向角：
[0028][0029]
这里，||e(t)||
w，p
为e(t)的加权l
p-范数，定义为：
[0030][0031]
其中，w为加权矩阵，定义为：
[0032][0033]
为脉冲噪声的方差，为脉冲噪声的方差，表示第一列为c1，最后一行为的hankel矩阵；矩阵u(t)定义为：
[0034]
u(t)＝diag([|[p(t)a-b(t)]1|
p-2
ꢀ…ꢀ
|[p(t)a-b(t)]
n-k
|
p-2
]
t
)。
[0035]
使用迭代重加权最小二乘法求解
[0036]
第一步：将u(t)和w设置为单位矩阵，估计线性预测系数矢量a：
[0037][0038]
第二步：基于构造矩阵u(t)和w；
[0039]
第三步：更新线性预测系数矢量的估计：
[0040][0041]
第四步：循环重复第二步和第三步，直到线性预测系数矢量在相邻两次迭代中的估计之差的l
2-范数小于10-13
，输出作为线性预测系数矢量的估计。
[0042]
于是，信源方向角可以通过求解线性预测方程：
[0043]
的根
[0044]
求得其估计值：
[0045]
的相位，k＝1，2，
…
，k。
[0046]
采用上述技术方案，本法具有以下有益效果：
[0047]
其一，和l
2-范数相比，l
p-范数具有较好的抵抗脉冲尖峰的特性。
[0048]
其二，求解线性预测系数矢量的估计主要计算量在于：
[0049]
1)矩阵乘法ph(t)u(t)wp(t)和ph(t)u(t)wb(t)，
[0050]
2)矩阵求逆运算ph(t)u(t)wp(t)，
[0051]
3)构造矩阵w和u(t)。
[0052]
因此，总计算量约为较现有技术而言，计算量低。
附图说明
[0053]
图1为本发明gmm噪声环境下各个信源定位方法性能比较图。其中，横轴为信噪比，纵轴为均方根误差(度)。
[0054]
图2为本发明gmm噪声环境下各个信源定位方法性能比较图，双目标，snr＝10db。其中，横轴为目标间距(度)，纵轴为均方根误差(度)。
[0055]
图3为本发明ggd噪声环境下各个信源定位方法性能比较图。其中，横轴为信噪比，纵轴为均方根误差(度)。
[0056]
图4为本发明ggd噪声环境下各个信源定位方法性能比较图，双目标，snr＝10db。其中，横轴为目标间距(度)，纵轴为均方根误差(度)。
[0057]
图5为本发明gmm噪声环境下各个信源定位方法性能表格图。
[0058]
图6为本发明ggd噪声环境下各个信源定位方法性能表格图。
[0059]
图7为本发明背景技术中的信号模型示意图。
具体实施方式
[0060]
下面结合附图对本发明的具体实施方式作进一步说明。在此需要说明的是，对于这些实施方式的说明用于帮助理解本发明，但并不构成对本发明的限定。此外，下面所描述的本发明各个实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互组合。
[0061]
实施例一：
[0062]
一种基于线性预测的信源方向鲁棒定位算法，包含以下步骤：
[0063]
首先，构建数据矩阵：
[0064]
y＝x n，
[0065]
其中，y＝[y(1)
ꢀ…ꢀ
y(t)]，x＝a[s(1)
ꢀ…ꢀ
s(t)]和n＝[n(1)
ꢀ…ꢀ
n(t)]分别为干净信号和噪声信号部分。
[0066]
基于干净信号x的线性预测性质，构建观测信号y的线性预测误差矩阵：
[0067]
e＝[e(1)
ꢀ…ꢀ
e(t)]，
[0068]
其中，a＝[a1ꢀ…ꢀak
]
t
为线性预测系数矢量，e(t)＝p(t)a-b(t)，而
[0069]
p(t)＝toeplitz([yk(t)
ꢀ…ꢀyn-1
(t)]
t
，[yk(t)
ꢀ…ꢀ
y1(t)])，
[0070]
b(t)＝-[y
k 1
(t)
ꢀ…ꢀyn
(t)]
t
，
[0071]
表示第一列为c1，第一行为的toeplitz矩阵。
[0072]
通过最小化线性预测误差的加权l
p-范数(1≤p＜2)，估计线性预测系数矢量a，进而估计信源方向角：
[0073][0074]
其中，||e(t)||
w，p
为e(t)的加权l
p-范数，定义为：
[0075][0076]
这里，w为加权矩阵，定义为：
[0077][0078]
为脉冲噪声的方差，为脉冲噪声的方差，表示第一列为c1，最后一行为的hankel矩阵；矩阵u(t)定义为：
[0079]
u(t)＝diag([|[p(t)a-b(t)]1|
p-2
ꢀ…ꢀ
|[p(t)a-b(t)]
n-k
|
p-2
]
t
)。
[0080]
使用迭代重加权最小二乘法求解
[0081]
第一步：将u(t)和w设置为单位矩阵，估计线性预测系数矢量a：
[0082][0083]
第二步：基于构造矩阵u(t)和w；
[0084]
第三步：更新线性预测系数矢量的估计：
[0085][0086]
第四步：循环重复第二步和第三步，直到线性预测系数矢量在相邻两次迭代中的估计之差的l
2-范数小于10-13
，输出作为线性预测系数矢量的估计。
[0087]
于是，信源方向角可以通过求解线性预测方程：
[0088]
的根
[0089]
求得其估计值：
[0090]
的相位，k＝1，2，
…
，k。
[0091]
实施例二：
[0092]
仿真试验结果：
[0093]
噪声模型：高斯混合模型(gaussian mixture model，gmm)——c1＝0.1，c2＝0.9，
[0094]
其他参数设置：k＝2，n＝12，t＝20。
[0095]
目标信号：
[0096]
θ1在-30
°
～-20
°
之间均匀分布，θ2在0
°
～10
°
之间均匀分布。
[0097]
s1(t)和s2(t)服从互为独立的复值高斯分布。
[0098]
评价指标：
[0099]
均方根误差(root mean square error，rmse)：
[0100][0101]
其中s为试验次数，为第s次试验θk的估计值。
[0102]
仿真试验结果：
[0103]
设置：p＝1.5，高斯混合模型(gmm)噪声。
[0104]
运行环境：on a pc with an intel(r)core(tm)i7-8750h cpu@2.20 ghz，with 16.0gb of installed memory。
[0105]
参照附图1、附图2和附图5，图中rwlp线条为本发明所述的基于线性预测的信源方向鲁棒定位算法的均方根误差曲线。在高斯混合模型仿真实验下，较传统的算法而言，本发明计算量较低，且定位时不需要初始值，可以为l
p-music算法提供初始值。同时，其分辨率较现有bayes-optimal方法而言较高。
[0106]
仿真试验结果：
[0107]
噪声模型：广义高斯分布(generalized gaussian distribution，ggd)——
[0108]
其他参数设置：k＝2，n＝12，t＝20。
[0109]
目标信号：
[0110]
θ1在-30
°
～一20
°
之间均匀分布，θ2在0
°
～10
°
之间均匀分布。
[0111]
s1(t)和s2(t)服从互为独立的复值高斯分布。
[0112]
评价指标：
[0113]
均方根误差(root mean square error，rmse)：
[0114][0115]
其中s为试验次数，为第s次试验θk的估计值。
[0116]
仿真试验结果：p＝1.5，广义高斯分布(ggd)噪声。
[0117]
运行环境：on a pc with an intel(r)core(tm)i7-8750h cpu@2.20ghz，with 16.0gb of installed memory。
[0118]
参照附图3，附图4和附图6本发明所述的基于线性预测的信源方向鲁棒定位算法与传统的鲁棒定位算法之比较：图中rwlp线条为本发明所述的基于线性预测的信源方向鲁棒定位算法的均方根误差曲线。在广义高斯分布仿真实验下，较传统的算法而言，本发明计算量较低，且定位时不需要初始值，可以为l
p-music算法提供初始值。同时，其分辨率较现有bayes-optimal方法而言较高。
[0119]
以上结合附图对本发明的实施方式作了详细说明，但本发明不限于所描述的实施方式。对于本领域的技术人员而言，在不脱离本发明原理和精神的情况下，对这些实施方式进行多种变化、修改、替换和变型，仍落入本发明的保护范围内。