CN114584433B - 一种脉冲噪声环境下多途信道中的同步信号检测方法 - Google Patents
一种脉冲噪声环境下多途信道中的同步信号检测方法 Download PDFInfo
- Publication number
- CN114584433B CN114584433B CN202210170624.2A CN202210170624A CN114584433B CN 114584433 B CN114584433 B CN 114584433B CN 202210170624 A CN202210170624 A CN 202210170624A CN 114584433 B CN114584433 B CN 114584433B
- Authority
- CN
- China
- Prior art keywords
- detection
- residual
- signal
- channel
- received signal
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 58
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 15
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 8
- 238000000034 method Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 7
- 241000238557 Decapoda Species 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 241000143060 Americamysis bahia Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000012885 constant function Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03828—Arrangements for spectral shaping; Arrangements for providing signals with specified spectral properties
- H04L25/03834—Arrangements for spectral shaping; Arrangements for providing signals with specified spectral properties using pulse shaping
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B11/00—Transmission systems employing sonic, ultrasonic or infrasonic waves
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
- H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Power Engineering (AREA)
- Noise Elimination (AREA)
Abstract
本发明提供一种脉冲噪声环境下多途信道中的同步信号检测方法,包括:输入接收信号序列;利用鲁棒正交匹配追踪ROMP估计多途信道及重构接收信号,构造非参数化的伪相关检测器PCD进行检测;判断检测统计量是否大于阈值。若检测统计量大于阈值,则同步信号存在,检测结束;若检测统计量小于阈值,则不存在同步信号,此时返回步骤一更新接收信号序列继续检测。本发明克服脉冲噪声和多途信道对信号检测的影响,可以抑制脉冲噪声的影响同时利用多途信息提高检测概率。
Description
技术领域
本发明涉及一种信号检测方法,属于水声信号处理领域。更确切地说是在脉冲噪声环境下,检测受多途信道影响的同步信号的方法。
背景技术
水声通信中,接收端数据处理的前提条件是准确判定有用信号是否到达及何时到达,信号检测成为一个重要环节。传统的信号检测方法往往假设背景噪声服从高斯分布,但水声信号在传播过程中经常会受到脉冲噪声的干扰。典型的脉冲声源包括极区冰层破裂或撞击产生的噪声和浅海海域鼓虾产生的生物噪声等,脉冲的存在使得这类噪声的概率密度函数不再服从高斯分布。此时,传统的基于二阶矩的检测算法性能严重下降。
受海底海面反射的影响,水声信号在传播过程中的时延可达几十甚至数百毫秒,多途严重。传统检测方法并未考虑多途带来的影响,这在一定程度上造成检测性能损失。因此,为了抑制脉冲噪声的影响,同时充分利用信道多途信息,本发明提出脉冲噪声环境下受多途信道影响的同步信号检测方法。
发明内容
本发明的目的是为克服脉冲噪声和多途信道对信号检测的影响而提供的一种信号检测方法,可以抑制脉冲噪声的影响同时利用多途信息提高检测概率。
本发明的目的是这样实现的:步骤如下:
步骤一:输入接收信号序列;
步骤二:利用鲁棒正交匹配追踪ROMP估计多途信道及重构接收信号,
步骤三:构造非参数化的伪相关检测器PCD进行检测;
步骤四:判断检测统计量是否大于阈值。若检测统计量大于阈值,则同步信号存在,检测结束;若检测统计量小于阈值,则不存在同步信号,此时返回步骤一更新接收信号序列继续检测。
本发明还包括这样一些结构特征:
1.步骤二的具体流程如下:
(a)输入接收信号、字典矩阵、稀疏度,初始化迭代次数、残差、原子索引集合及信道估计值;
(b)计算残差尺度估计值;
(c)将残差尺度化,得到伪残差值;
(d)计算伪残差与各原子的内积并选择内积最大的原子,将它放入原子索引集合中;
(e)计算信道的估计值,更新残差;
(f)若迭代次数未达到预设稀疏度,返回(b)继续迭代;否则,返回重构信号。
2.PCD的具体实现如下:
NP准则下,二元假设检验模型为:
H0:y[n]=w[n]
H1:y[n]=Ah[n]+w[n]
其中,y[n]=[y[n-N+1],...y[n]]T为接收信号,N为当前处理块长度;w[n]=[w[n-N+1],...w[n]]T为脉冲噪声,为发送信号矩阵,D=N-NT为最大时延,A的第i(0≤i≤D-1)列由发送信号s的延时构成,即/>h[n]是长度为D的信道冲激响应;H1表示存在同步信号,H0表示不存在同步信号。
与现有技术相比,本发明的有益效果是:(1)二元假设检验模型中综合考虑了脉冲噪声和多途信道的影响;(2)与正交匹配追踪(OMP)相比,ROMP能抑制脉冲噪声的影响,准确地重构接收信号;(3)PCD利用了信道的多途信息,带来检测性能的提升,同时它操作简单,计算量低;(4)本发明适用于极区通信,能抵抗由冰层破裂或撞击产生的脉冲噪声,在复杂的极区环境中保持鲁棒性。
附图说明
图1为本发明流程图;
图2为ROMP和OMP对比图;
图3为特征指数α=1.5的SαS噪声中不同检测器性能对比图;
图4为特征指数α=1.9的SαS噪声中不同检测器性能对比图;
图5为实测鼓虾噪声中不同检测器性能对比图;
图6为ROMP流程图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
本发明的目的是这样实现的:
步骤一:输入接收信号序列;
步骤二:利用鲁棒正交匹配追踪(ROMP)估计多途信道及重构接收信号,具体流程如下:
(a)输入接收信号、字典矩阵、稀疏度,初始化迭代次数、残差、原子索引集合及信道估计值;
(b)计算残差尺度估计值;
(c)将残差尺度化,得到伪残差值;
(d)计算伪残差与各原子的内积并选择内积最大的原子,将它放入原子索引集合中;
(e)计算信道的估计值,更新残差;
(f)若迭代次数未达到预设稀疏度,返回步骤(b)继续迭代;否则,返回重构信号。
步骤三:构造非参数化的伪相关检测器(PCD)进行检测;
步骤四:判断检测统计量是否大于阈值。若检测统计量大于阈值,则同步信号存在,检测结束;若检测统计量小于阈值,则不存在同步信号,此时返回步骤(1)更新接收信号序列继续检测。
本发明的PCD的具体实现如下:
NP准则下,二元假设检验模型为:
H0:y[n]=w[n]
H1:y[n]=Ah[n]+w[n]
其中,y[n]=[y[n-N+1],...y[n]]T为接收信号,N为当前处理块长度。w[n]=[w[n-N+1],...w[n]]T为脉冲噪声,这里假设服从对称α稳定分布(SαS)。为发送信号矩阵,D=N-NT为最大时延,A的第i(0≤i≤D-1)列由发送信号s的延时构成,即h[n]是长度为D的信道冲激响应。H1表示存在同步信号,H0表示不存在同步信号。
为充分考虑多途传播对检测性能的影响,有必要先估计出信道结构。利用水声信道的稀疏特性,采用基于压缩感知的鲁棒信道估计方法,图6的表为ROMP流程。
输入:接收信号y,字典矩阵A,稀疏度L;
初始化:迭代次数i=0,残差ei=y,原子索引集合信道估计值/>
步骤1:计算尺度估计值
步骤2:将残差尺度化,得到伪残差值
步骤3:计算伪残差与各原子的内积并选择最大的原子,即将i放入集合S中,即S=S∪{i},i=i+1;
步骤4:计算信道M估计更新残差/>
步骤5:若i≤L,返回步骤1;否则,返回重构信号
其中,步骤2中的φ(·)是代价函数f(·)的导数,步骤4中的MFIT(·)表示通过迭代最小二乘方法求信道估计值,w(x)=φ(x)/x为权重函数。ROMP和OMP对比如图2所示,可以看到ROMP的代价函数根据残差的大小变化。对于大残差,f(·)是有界的常函数,对于小误差,f(·)类似于平方函数。这样,与OMP相比,它在保持算法性能的同时有效降低了脉冲噪声的影响。
本发明提出非参数化PCD检测,检测统计量为:
可以看到,与SCD的检测统计量T(y)=sTsign(y)相比,PCD的不同主要体现在(1)将符号函数sign(·)用φ(·)代替以保留更多原始接收信号的信息;(2)考虑了信道多途对检测的影响。实际上,当代价函数f(e)=|e|,ROMP迭代次数为1时,PCD退化为SCD。将检测统计量T(y)与预设阈值γ比较,若T(y)>γ,H1成立,存在同步信号;否则,接收信号中只存在噪声。
结合具体参数给出本发明的仿真研究:
仿真条件:稀疏水声信道的多途个数为20,相邻路径时延差服从均值为1ms的指数分布,平均时延为20ms。每条路径的幅度是独立的瑞利分布,平均能量随时延的增加指数衰减,从第0到第20ms的能量差是15dB。同步信号为线性调频(LFM)信号,频带范围2-6kHz,持续时间100ms。利用能量噪声比ENR=10log(ES/2δ2)dB表征信号和噪声之间的相对关系,其中ES为LFM的能量,2δ2为SαS噪声方差。接收端N=200ms,信道延时长度D=100ms。
设SαS噪声特征指数为α=1.5和α=1.9,虚警概率PFA=1e-4,检测概率随能量噪声比的变化分别如图3和图4所示。其中,“L=1”表示迭代次数为1,此时未利用多途信息;“L=3”表示信道估计迭代次数为3;“LLR”是NP准则下的最优检测,需要已知噪声的概率密度函数,可以作为检测性能的基准;“ACC”表示利用OMP估计信道后构造的累计相关检测。从结果可以看到(1)基于OMP的方法在脉冲噪声中不鲁棒,尤其是α=1.5这种强脉冲环境中。在α=1.9这类弱脉冲噪声中,它的检测性能有所提升,但是仍然逊于ROMP下的检测性能。ROMP能抑制脉冲噪声的干扰,具有强鲁棒性。(2)PCD性能稍逊于LLR最优检测,但它不需要利用噪声分布的先验信息,计算量小。(3)利用信道多途信息后的检测性能(L=3)优于未利用多途信息的性能(L=1)。
在实测鼓虾噪声中比较各检测器的性能。鼓虾噪声的持续时间为1300s,选择前900s数据作为H0假设下的噪声,计算阈值;后400s数据作为H1假设下的噪声,计算检测概率。其中,虚警概率设为PFA=2e-3,检测结果如图5。该结果与仿真SαS中结果类似,证实了传统OMP算法性能的退化,基于ROMP的PCD检测器的鲁棒性以及利用信道多途信息后性能的提升。
Claims (2)
1.一种脉冲噪声环境下多途信道中的同步信号检测方法,其特征在于,步骤如下:
步骤一:输入接收信号序列;
步骤二:利用鲁棒正交匹配追踪ROMP估计多途信道及重构接收信号,
步骤三:构造非参数化的伪相关检测器PCD进行检测;
输入:接收信号y,字典矩阵A,稀疏度L;
初始化:迭代次数i=0,残差ei=y,原子索引集合信道估计值/>
步骤3.1:计算尺度估计值
步骤3.2:将残差尺度化,得到伪残差值φ(·)是代价函数f(·)的导数;
步骤3.3:计算伪残差与各原子的内积并选择最大的原子,即将i放入集合S中,即S=S∪{i},i=i+1;
步骤3.4:计算信道M估计更新残差/>MFIT(·)表示通过迭代最小二乘方法求信道估计值,w(x)=φ(x)/x为权重函数;
步骤3.5:若i≤L,返回步骤3.1;否则,返回重构信号
步骤3.6:非参数化PCD检测,检测统计量为:
将检测统计量T(y)与预设阈值γ比较,若T(y)>γ,存在同步信号;否则,接收信号中只存在噪声;
步骤四:判断检测统计量是否大于阈值;若检测统计量大于阈值,则同步信号存在,检测结束;若检测统计量小于阈值,则不存在同步信号,此时返回步骤一更新接收信号序列继续检测。
2.根据权利要求1所述的一种脉冲噪声环境下多途信道中的同步信号检测方法,其特征在于,步骤二的具体流程如下:
(a)输入接收信号、字典矩阵、稀疏度,初始化迭代次数、残差、原子索引集合及信道估计值;
(b)计算残差尺度估计值;
(c)将残差尺度化,得到伪残差值;
(d)计算伪残差与各原子的内积并选择内积最大的原子,将它放入原子索引集合中;
(e)计算信道的估计值,更新残差;
(f)若迭代次数未达到预设稀疏度,返回(b)继续迭代;否则,返回重构信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210170624.2A CN114584433B (zh) | 2022-02-24 | 2022-02-24 | 一种脉冲噪声环境下多途信道中的同步信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210170624.2A CN114584433B (zh) | 2022-02-24 | 2022-02-24 | 一种脉冲噪声环境下多途信道中的同步信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114584433A CN114584433A (zh) | 2022-06-03 |
CN114584433B true CN114584433B (zh) | 2024-05-28 |
Family
ID=81774970
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210170624.2A Active CN114584433B (zh) | 2022-02-24 | 2022-02-24 | 一种脉冲噪声环境下多途信道中的同步信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114584433B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903630A (zh) * | 2014-03-18 | 2014-07-02 | 北京捷通华声语音技术有限公司 | 一种用于消除稀疏噪声方法及装置 |
US8861655B1 (en) * | 2013-06-11 | 2014-10-14 | King Fahd University Of Petroleum And Minerals | Method of performing structure-based bayesian sparse signal reconstruction |
CN106936514A (zh) * | 2017-02-24 | 2017-07-07 | 哈尔滨工业大学深圳研究生院 | 稀疏信道下基于能量集中的水声前导信号检测方法 |
CN107395536A (zh) * | 2017-07-20 | 2017-11-24 | 西北工业大学 | 对多途环境下的水声信道冲激响应函数进行估计的方法 |
CN109039960A (zh) * | 2018-06-26 | 2018-12-18 | 东南大学 | 一种水声稀疏信道估计变步长稀疏度自适应匹配追踪方法 |
CN109818888A (zh) * | 2019-03-25 | 2019-05-28 | 哈尔滨工程大学 | 一种脉冲干扰环境下的组稀疏水声信道估计方法 |
CN110138459A (zh) * | 2019-04-04 | 2019-08-16 | 中国人民解放军战略支援部队信息工程大学 | 基于基追踪去噪的稀疏水声正交频分复用信道估计方法及装置 |
CN111800362A (zh) * | 2020-06-22 | 2020-10-20 | 哈尔滨工业大学(威海) | 一种改进的romp水声信道估计算法 |
CN112653640A (zh) * | 2020-12-18 | 2021-04-13 | 哈尔滨工程大学青岛船舶科技有限公司 | 一种脉冲噪声抑制水声信道估计方法 |
CN112737702A (zh) * | 2020-12-18 | 2021-04-30 | 哈尔滨工程大学 | 一种稀疏干扰背景下的mimo水声信道估计方法 |
CN113612542A (zh) * | 2021-07-31 | 2021-11-05 | 浙江大学 | 应用于水下无线光通信系统的基于变步长广义正交匹配追踪算法的非线性均衡方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8681902B2 (en) * | 2011-09-16 | 2014-03-25 | Alcatel Lucent | Method and apparatus for low complexity robust reconstruction of noisy signals |
-
2022
- 2022-02-24 CN CN202210170624.2A patent/CN114584433B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8861655B1 (en) * | 2013-06-11 | 2014-10-14 | King Fahd University Of Petroleum And Minerals | Method of performing structure-based bayesian sparse signal reconstruction |
CN103903630A (zh) * | 2014-03-18 | 2014-07-02 | 北京捷通华声语音技术有限公司 | 一种用于消除稀疏噪声方法及装置 |
CN106936514A (zh) * | 2017-02-24 | 2017-07-07 | 哈尔滨工业大学深圳研究生院 | 稀疏信道下基于能量集中的水声前导信号检测方法 |
CN107395536A (zh) * | 2017-07-20 | 2017-11-24 | 西北工业大学 | 对多途环境下的水声信道冲激响应函数进行估计的方法 |
CN109039960A (zh) * | 2018-06-26 | 2018-12-18 | 东南大学 | 一种水声稀疏信道估计变步长稀疏度自适应匹配追踪方法 |
CN109818888A (zh) * | 2019-03-25 | 2019-05-28 | 哈尔滨工程大学 | 一种脉冲干扰环境下的组稀疏水声信道估计方法 |
CN110138459A (zh) * | 2019-04-04 | 2019-08-16 | 中国人民解放军战略支援部队信息工程大学 | 基于基追踪去噪的稀疏水声正交频分复用信道估计方法及装置 |
CN111800362A (zh) * | 2020-06-22 | 2020-10-20 | 哈尔滨工业大学(威海) | 一种改进的romp水声信道估计算法 |
CN112653640A (zh) * | 2020-12-18 | 2021-04-13 | 哈尔滨工程大学青岛船舶科技有限公司 | 一种脉冲噪声抑制水声信道估计方法 |
CN112737702A (zh) * | 2020-12-18 | 2021-04-30 | 哈尔滨工程大学 | 一种稀疏干扰背景下的mimo水声信道估计方法 |
CN113612542A (zh) * | 2021-07-31 | 2021-11-05 | 浙江大学 | 应用于水下无线光通信系统的基于变步长广义正交匹配追踪算法的非线性均衡方法 |
Non-Patent Citations (1)
Title |
---|
基于基追踪去噪的水声正交频分复用稀疏信道估计;尹艳玲;乔钢;刘凇佐;周锋;;物理学报(06);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114584433A (zh) | 2022-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Sequential detection and estimation of multipath channel parameters using belief propagation | |
CN107395536B (zh) | 对多途环境下的水声信道冲激响应函数进行估计的方法 | |
CN107561540B (zh) | 基于混响对称谱特性的声纳泄露目标的检测方法及装置 | |
CN108646249B (zh) | 一种适用于部分均匀混响背景的参数化泄露目标检测方法 | |
Cui et al. | Multiple-input multiple-output radar detectors design in non-Gaussian clutter | |
CN111427042B (zh) | 基于高刷新率声脉冲的水下弱目标回波检测前跟踪方法 | |
Huang et al. | GPS signal detection under multiplicative and additive noise | |
CN113064147A (zh) | 一种低信噪比下新型匹配场被动定位方法 | |
CN111090089B (zh) | 一种基于两类辅助数据的空时自适应检测方法 | |
CN114584433B (zh) | 一种脉冲噪声环境下多途信道中的同步信号检测方法 | |
CN113671473B (zh) | 基于环境约束和黎曼距离的联合匹配场定位方法及系统 | |
Li et al. | Underwater acoustic localization of the long baseline based on track-before-detect | |
Broadhead et al. | Performance of some sparseness criterion blind deconvolution methods in the presence of noise | |
CN108594200B (zh) | 一种被动式mimo雷达的全相参目标检测方法 | |
Tian et al. | Wideband signal detection in multipath environment affected by impulsive noise | |
CN114690175B (zh) | 一种基于被动外辐射源雷达的目标直接检测与跟踪方法 | |
CN114994632B (zh) | 一种基于对称功率谱密度的雷达目标检测方法与系统 | |
De Maio et al. | Estimation of the covariance matrix based on multiple a-priori models | |
Rabaste et al. | Approximate multitarget matched filter for MIMO radar detection via orthogonal matching pursuit | |
Zhao et al. | Parameter estimation of LFM-BPSK composite modulated radar signal based on multi-output regression network | |
KR100681393B1 (ko) | 가상 센서를 이용한 다중경로 검출방법 | |
Zhang et al. | An adaptive receiver with an antenna array for channels with correlated non-Gaussian interference and noise using the SAGE algorithm | |
CN116910453A (zh) | 一种极地脉冲噪声干扰下的时延-多普勒联合估计方法 | |
CN109031239B (zh) | 基于信息融合的压缩感知外辐射源雷达对目标检测方法 | |
Bolisetti et al. | Compressive parametric GLRT detector for airborne MIMO radar |
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 |