CN113987733A - 一种基于线性预测的信源方向鲁棒定位算法 - Google Patents
一种基于线性预测的信源方向鲁棒定位算法 Download PDFInfo
- Publication number
- CN113987733A CN113987733A CN202011365016.4A CN202011365016A CN113987733A CN 113987733 A CN113987733 A CN 113987733A CN 202011365016 A CN202011365016 A CN 202011365016A CN 113987733 A CN113987733 A CN 113987733A
- Authority
- CN
- China
- Prior art keywords
- linear prediction
- coefficient vector
- matrix
- prediction coefficient
- norm
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000004891 communication Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 8
- 238000009826 distribution Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 239000000203 mixture Substances 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于线性预测的信源方向鲁棒定位算法,先构建数据矩阵,然后通过最小化线性预测误差的加权lp‑范数(1≤p<2),估计线性预测系数矢量a,进而估计信源方向;与传统的加权l2‑范数相比,具有更强的抵抗脉冲尖峰特性;与Bayes‑optimal方法相比,有较高的分辨率;本发明使用迭代重加权最小二乘法求解线性预测系数矢量a,具有不需要初始估计值、计算量小的优点,易于硬件实现。
Description
技术领域
本发明涉及信号处理技术领域,特别涉及一种基于线性预测的信源方向鲁棒定位算法。
背景技术
发展水声通信技术是国家海洋战略及现代海军建设的需要。由于电磁波在水中衰减严重,高频无线电通信无法适用于水下通信的场景。水下航行器与母船/基地的通信联络可采用3-300Hz的超低频频段。实践证明:水声通信的合适载波频率为4kHz左右,可见其BW0.707带宽极低,这与人们熟悉的高频无线电通信有本质区别。此外,海洋环境复杂,噪声干扰大。因此,水声通信技术有相当难度。虽然常规的无线通信技术无法直接应用于水声通信,但在原理上是相通的。一般情况下,我们处理的是基带信号。因此,对于水声通信的信号处理,关键在于传统模型的修正。本专利提出一种基于线性预测的信源方向鲁棒定位算法,用于水下目标物体的定位。
在脉冲噪声环境中,对K个信源定位,信号模型示意图如附图7:
现有N个传感器构成的均匀线性阵列(ULA),用于接收从K个目标发来的波动信号s1(t),s2(t),…,sK(t)。在t时刻,ULA接收到的信号可以建模成:
y(t)=As(t)+n(t),t=1,2,…,T,
其中,
y(t)=[y1(t) … yN(t)]T,
s(t)=[s1(t) … sK(t)]T,
A=[a(θ1) … a(θK)],
这里,噪声向量n(t)建模成脉冲噪声(impulsive noise)。
现有技术中,信源方向鲁棒定位算法有:1.基于字典的鲁棒定位算法;2.利用观测信号的鲁棒定位算法,包括基于鲁棒协方差矩阵(robust covariation matrix)的RCM-MUSIC算法;3.基于分数低阶矩(fractional lower order moment)的FLOM-MUSIC算法,基于符号协方差矩阵(sign covariance matrix)和Kendall’stau协方差矩阵(Kendall’stau covariance matrix)的SCM-MUSIC和TCM-MUSIC算法;4.以鲁棒方式观测信号的二阶统计量并定位,包括基于MM-估计器的算法和lp-MUSIC算法等;5.通过最小化观测信号与信号模型之间的拟合误差定位。
但是,上述信源方向鲁棒定位算法中:1存在不适用于衰减信号,分辨率受字典密度制约而较低,计算量大的问题;2-5也是计算量较大。这些问题导致相应的算法使用受限。
发明内容
针对现实中存在的缺陷,本发明提供一种基于线性预测的信源方向鲁棒定位算法。
为了解决上述技术问题,本发明的技术方案为:
一种基于线性预测的信源方向鲁棒定位算法,包含以下步骤:
首先,构建数据矩阵:
Y=X+N,
其中,Y=[y(1) … y(T)],X=A[s(1) … s(T)]和N=[n(1) … n(T)]分别为干净信号和噪声信号部分。
基于干净信号X的线性预测性质,构建观测信号Y的线性预测误差矩阵:
E=[e(1) … e(T)],
其中,a=[a1 … aK]T为线性预测系数矢量,e(t)=P(t)a-b(t),而
P(t)=Toeplitz([yK(t) … yN-1(t)]T,[yK(t) … y1(t)]),
b(t)=-[yK+1(t) … yN(t)]T,
通过最小化线性预测误差的加权lp-范数(1≤p<2),估计线性预测系数矢量a,进而估计信源方向角:
这里,||e(t)||w,p为e(t)的加权lp-范数,定义为:
其中,W为加权矩阵,定义为:
U(t)=diag([|[P(t)a-b(t)]1|p-2 … |[P(t)a-b(t)]N-K|p-2]T)。
第一步:将U(t)和W设置为单位矩阵,估计线性预测系数矢量a:
第三步:更新线性预测系数矢量的估计:
于是,信源方向角可以通过求解线性预测方程:
求得其估计值:
采用上述技术方案,本法具有以下有益效果:
其一,和l2-范数相比,lp-范数具有较好的抵抗脉冲尖峰的特性。
1)矩阵乘法PH(t)U(t)WP(t)和PH(t)U(t)Wb(t),
2)矩阵求逆运算PH(t)U(t)WP(t),
3)构造矩阵W和U(t)。
附图说明
图1为本发明GMM噪声环境下各个信源定位方法性能比较图。其中,横轴为信噪比,纵轴为均方根误差(度)。
图2为本发明GMM噪声环境下各个信源定位方法性能比较图,双目标,SNR=10dB。其中,横轴为目标间距(度),纵轴为均方根误差(度)。
图3为本发明GGD噪声环境下各个信源定位方法性能比较图。其中,横轴为信噪比,纵轴为均方根误差(度)。
图4为本发明GGD噪声环境下各个信源定位方法性能比较图,双目标,SNR=10dB。其中,横轴为目标间距(度),纵轴为均方根误差(度)。
图5为本发明GMM噪声环境下各个信源定位方法性能表格图。
图6为本发明GGD噪声环境下各个信源定位方法性能表格图。
图7为本发明背景技术中的信号模型示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明。在此需要说明的是,对于这些实施方式的说明用于帮助理解本发明,但并不构成对本发明的限定。此外,下面所描述的本发明各个实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互组合。
实施例一:
一种基于线性预测的信源方向鲁棒定位算法,包含以下步骤:
首先,构建数据矩阵:
Y=X+N,
其中,Y=[y(1) … y(T)],X=A[s(1) … s(T)]和N=[n(1) … n(T)]分别为干净信号和噪声信号部分。
基于干净信号X的线性预测性质,构建观测信号Y的线性预测误差矩阵:
E=[e(1) … e(T)],
其中,a=[a1 … aK]T为线性预测系数矢量,e(t)=P(t)a-b(t),而
P(t)=Toeplitz([yK(t) … yN-1(t)]T,[yK(t) … y1(t)]),
b(t)=-[yK+1(t) … yN(t)]T,
通过最小化线性预测误差的加权lp-范数(1≤p<2),估计线性预测系数矢量a,进而估计信源方向角:
其中,||e(t)||w,p为e(t)的加权lp-范数,定义为:
这里,W为加权矩阵,定义为:
U(t)=diag([|[P(t)a-b(t)]1|p-2 … |[P(t)a-b(t)]N-K|p-2]T)。
第一步:将U(t)和W设置为单位矩阵,估计线性预测系数矢量a:
第三步:更新线性预测系数矢量的估计:
于是,信源方向角可以通过求解线性预测方程:
求得其估计值:
实施例二:
仿真试验结果:
其他参数设置:K=2,N=12,T=20。
目标信号:
θ1在-30°~-20°之间均匀分布,θ2在0°~10°之间均匀分布。
s1(t)和s2(t)服从互为独立的复值高斯分布。
评价指标:
均方根误差(root mean square error,RMSE):
仿真试验结果:
设置:p=1.5,高斯混合模型(GMM)噪声。
运行环境:On a PC with an Intel(R)Core(TM)i7-8750H CPU@2.20 GHz,with16.0GB of installed memory。
参照附图1、附图2和附图5,图中RWLP线条为本发明所述的基于线性预测的信源方向鲁棒定位算法的均方根误差曲线。在高斯混合模型仿真实验下,较传统的算法而言,本发明计算量较低,且定位时不需要初始值,可以为lp-MUSIC算法提供初始值。同时,其分辨率较现有Bayes-optimal方法而言较高。
仿真试验结果:
其他参数设置:K=2,N=12,T=20。
目标信号:
θ1在-30°~一20°之间均匀分布,θ2在0°~10°之间均匀分布。
s1(t)和s2(t)服从互为独立的复值高斯分布。
评价指标:
均方根误差(root mean square error,RMSE):
仿真试验结果:p=1.5,广义高斯分布(GGD)噪声。
运行环境:On a PC with an Intel(R)Core(TM)i7-8750H CPU@2.20GHz,with16.0GB of installed memory。
参照附图3,附图4和附图6本发明所述的基于线性预测的信源方向鲁棒定位算法与传统的鲁棒定位算法之比较:图中RWLP线条为本发明所述的基于线性预测的信源方向鲁棒定位算法的均方根误差曲线。在广义高斯分布仿真实验下,较传统的算法而言,本发明计算量较低,且定位时不需要初始值,可以为lp-MUSIC算法提供初始值。同时,其分辨率较现有Bayes-optimal方法而言较高。
以上结合附图对本发明的实施方式作了详细说明,但本发明不限于所描述的实施方式。对于本领域的技术人员而言,在不脱离本发明原理和精神的情况下,对这些实施方式进行多种变化、修改、替换和变型,仍落入本发明的保护范围内。
Claims (1)
1.一种基于线性预测的信源方向鲁棒定位算法,其特征在于:包含以下步骤:
首先,构建数据矩阵:
Y=X+N,
其中,Y=[y(1) … y(T)],X=A[s(1) … s(T)]和N=[n(1) … n(T)]分别为干净信号和噪声信号部分;
基于干净信号X的线性预测性质,构建观测信号Y的线性预测误差矩阵:
E=[e(1) … e(T)],
其中,a=[a1 … aK]T为线性预测系数矢量,e(t)=P(t)a-b(t),而
P(t)=Toeplitz([yK(t) … yN-1(t)]T,[yK(t) … y1(t)]),
b(t)=-[yK+1(t) … yN(t)]T;
本专利提出,通过最小化线性预测误差的加权lp-范数(1≤p<2),估计线性预测系数矢量a,进而估计信源方向角:
这里,‖e(t)‖W,p为e(t)的加权lp-范数,定义为:
其中,W为加权矩阵,定义为:
U(t)=diag([|[P(t)a-b(t)]1|p-2 … |[P(t)a-b(t)]N-K|p-2]T);
第三步:更新线性预测系数矢量的估计:
于是,信源方向角可以通过求解线性预测方程:
求得其估计值:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011365016.4A CN113987733B (zh) | 2020-11-27 | 2020-11-27 | 一种基于线性预测的信源方向鲁棒定位算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011365016.4A CN113987733B (zh) | 2020-11-27 | 2020-11-27 | 一种基于线性预测的信源方向鲁棒定位算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113987733A true CN113987733A (zh) | 2022-01-28 |
CN113987733B CN113987733B (zh) | 2024-05-14 |
Family
ID=79731634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011365016.4A Active CN113987733B (zh) | 2020-11-27 | 2020-11-27 | 一种基于线性预测的信源方向鲁棒定位算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113987733B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005077376A (ja) * | 2003-09-03 | 2005-03-24 | Furuno Electric Co Ltd | 信号到来方向推定装置、アダプティブアレー、受信信号処理装置および測位装置 |
WO2006067869A1 (ja) * | 2004-12-24 | 2006-06-29 | Fujitsu Limited | 到来方向推定装置及びプログラム |
JP2017216601A (ja) * | 2016-05-31 | 2017-12-07 | 日本電信電話株式会社 | 線形予測係数出力装置、線形予測係数出力方法及び線形予測係数出力プログラム |
CN107544052A (zh) * | 2017-08-07 | 2018-01-05 | 大连大学 | 一种基于矩阵补全的二阶统计量重构doa估计方法 |
US20190050372A1 (en) * | 2017-08-14 | 2019-02-14 | City University Of Hong Kong | Systems and methods for robust low-rank matrix approximation |
CN110412499A (zh) * | 2019-07-16 | 2019-11-05 | 北京工业大学 | 基于压缩感知理论下的rss算法的宽带doa估计算法 |
CN111948599A (zh) * | 2020-08-14 | 2020-11-17 | 电子科技大学 | 一种角度相关互耦影响下相干信号的高分辨率定位方法 |
-
2020
- 2020-11-27 CN CN202011365016.4A patent/CN113987733B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005077376A (ja) * | 2003-09-03 | 2005-03-24 | Furuno Electric Co Ltd | 信号到来方向推定装置、アダプティブアレー、受信信号処理装置および測位装置 |
WO2006067869A1 (ja) * | 2004-12-24 | 2006-06-29 | Fujitsu Limited | 到来方向推定装置及びプログラム |
JP2017216601A (ja) * | 2016-05-31 | 2017-12-07 | 日本電信電話株式会社 | 線形予測係数出力装置、線形予測係数出力方法及び線形予測係数出力プログラム |
CN107544052A (zh) * | 2017-08-07 | 2018-01-05 | 大连大学 | 一种基于矩阵补全的二阶统计量重构doa估计方法 |
US20190050372A1 (en) * | 2017-08-14 | 2019-02-14 | City University Of Hong Kong | Systems and methods for robust low-rank matrix approximation |
CN110412499A (zh) * | 2019-07-16 | 2019-11-05 | 北京工业大学 | 基于压缩感知理论下的rss算法的宽带doa估计算法 |
CN111948599A (zh) * | 2020-08-14 | 2020-11-17 | 电子科技大学 | 一种角度相关互耦影响下相干信号的高分辨率定位方法 |
Non-Patent Citations (2)
Title |
---|
WEN-JUN ZENG: "An ℓp-norm minimization approach to time delay estimation in impulsive noise", 《DIGITAL SIGNAL PROCESSING》, 1 July 2013 (2013-07-01), pages 1247 * |
田野;练秋生;: "基于重加权l_1范数惩罚的远近场混合源定位算法", 电子学报, no. 10, 15 October 2016 (2016-10-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN113987733B (zh) | 2024-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107450047B (zh) | 嵌套阵下基于未知互耦信息的压缩感知doa估计方法 | |
CN110113279B (zh) | 一种移动跳频水声通信多普勒因子估计方法 | |
CN110515038B (zh) | 一种基于无人机-阵列的自适应无源定位装置及实现方法 | |
WO2023165631A1 (zh) | 一种用于智能超表面无线通信的信道参数估计方法 | |
CN107395536B (zh) | 对多途环境下的水声信道冲激响应函数进行估计的方法 | |
CN114745237B (zh) | 一种智能超表面辅助多用户无线通信系统的信道估计方法 | |
CN108562866A (zh) | 基于矩阵填充的双基地mimo雷达角度估算方法 | |
CN105939299A (zh) | 一种基于改进后的sage算法的信道参数估计方法 | |
CN112346014B (zh) | 一种基于信号到达时间差的多基地声纳的定位方法 | |
CN113242191A (zh) | 一种改进的时序多重稀疏贝叶斯学习水声信道估计方法 | |
CN113064147A (zh) | 一种低信噪比下新型匹配场被动定位方法 | |
Bilgehan et al. | Fast detection and DOA estimation of the unknown wideband signal sources | |
CN104808190B (zh) | 改善mimo雷达最差参数估计性能的稳健波形设计方法 | |
CN103323810A (zh) | 一种l阵方位角和俯仰角配对的信号处理方法 | |
CN111665469B (zh) | 一种基于空间时频分布的水下多径信号参数估计方法 | |
CN105656577A (zh) | 面向信道冲激响应的分簇方法和装置 | |
CN113987733B (zh) | 一种基于线性预测的信源方向鲁棒定位算法 | |
CN109188345B (zh) | 基于去预延迟空时结构的相干信号源doa估计方法 | |
CN112034492A (zh) | 一种空时极三维联合导航阵列抗干扰处理方法 | |
CN117240666A (zh) | 一种基于簇自注意力机制的神经网络水声信道估计方法 | |
CN106685507A (zh) | 色噪声环境下基于约束Kalman波束形成方法 | |
CN108051773A (zh) | 基于盖式圆盘准则估计信源数目的epuma方法 | |
CN113992285A (zh) | 一种基于平移不变性的信源数鲁棒检测方法 | |
CN105652264A (zh) | 基于高阶累积量的多路径传播声信号分离方法 | |
Feng et al. | Performance evaluation of acoustic model-based blind channel estimation in ocean waveguides |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |