CN111711584A - 基于成比例归一化最小均方误差的稀疏水声信道估计方法 - Google Patents
基于成比例归一化最小均方误差的稀疏水声信道估计方法 Download PDFInfo
- Publication number
- CN111711584A CN111711584A CN202010444429.5A CN202010444429A CN111711584A CN 111711584 A CN111711584 A CN 111711584A CN 202010444429 A CN202010444429 A CN 202010444429A CN 111711584 A CN111711584 A CN 111711584A
- Authority
- CN
- China
- Prior art keywords
- underwater acoustic
- sparse
- acoustic channel
- channel
- impulse response
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000005316 response function Methods 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000004891 communication Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
Images
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/0202—Channel estimation
-
- 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
-
- 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/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
-
- 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/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0254—Channel estimation channel estimation algorithms using neural network algorithms
-
- 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)
- Power Engineering (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种基于成比例归一化最小均方误差的稀疏水声信道估计方法,用以估计时域簇稀疏水声信道冲激响应函数。首先,针对时变水声信道进行簇稀疏建模,以期获得水声信道结构化稀疏表达框架,在此基础上。采用逐块训练模式,设计一种基于成比例归一化最小均方误差算法迭代寻优得到时域稀疏水声信道冲激响应函数,将所估计得到的时域水声信道信息。该项发明适用于时变水声信道估计、水声通信等。有益效果体现在:本发明基于水声信道的簇稀疏范数约束,有效地减少了冗余信道抽头计算,使得本发明产生的水声信道估计结果具有更高的精度。
Description
技术领域
本发明属于水声通信和水声信号处理等领域,涉及一种基于成比例归一化最小均方误差的稀疏水声信道估计方法,该估计方法将提高对时变水声信道的估计性能。
背景技术
水声信道估计和水声通信等问题都可归结为冲激响应函数的估计优化问题,基于训练序列和接收到的信号用以对时变水声信道稀疏表达估计。目前,对水声信道的估 计方法包括有限冲激响应框架和时域的逐块估计框架。有限冲激响应框架的算法细节 具体参见《New sparse adaptive algorithms based on the natural gradient and the L0-norm》,该文2013年发表于《IEEE Journal of Oceanic Engineering》第38期,起始 页码为323。时域的逐块估计框架详见《Estimation of rapidly time-varying sparsechannels》该文2007年发表于《IEEE Journal of Oceanic Engineering》第32期,起始 页码为927。
由于水声信道的多径扩展和时变特性,导致对水声信道的冲激响应函数估计极为困难,因此,对于有限冲激响应框架下的算法效果欠佳。考虑到水声信道的多径扩展 和时变特性可由水声信道的时域进行刻画。本发明将建立在该模型基础上,对时变水 声信道进行估计。然而待估计的抽头系数繁多,计算量大,所幸的是,稀疏估计策略 可以简化估计的抽头个数,从而有望提高估计效率。考虑到实际中的水声信道冲激响 应函数不是严格意义上的稀疏信号,因此,现行压缩感知算法难以直接应用。本发明 提出在簇稀疏表达的基础上,对稀疏水声信道估计的目标函数施加范数约束,以期提 高水声信道的估计精度。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,有效估计出时变多径的水声信道冲激响应函数,本发明提出一种基于成比例归一化最小均方误差的稀疏水声信道估计方法。是一种基 于时域的迭代寻优方法,该发明利用基于成比例归一化最小均方误差算法的迭代框架, 采用簇稀疏约束的目标函数,从而实现迭代的误差最小化,最终精确地恢复出稀疏水 声信道冲激响应函数。
技术方案
一种基于成比例归一化最小均方误差的稀疏水声信道估计方法,其特征在于步骤如下:
步骤1:设参数L,K和d分别为水声信道冲激响应长度、信道分块个数和信道分 块的长度,它们之间的关系表示为L=Kd;设置迭代步长为μ,簇稀疏约束参数γ,α;
步骤2:设置稀疏光滑度调整参数β=10~12,基于成比例归一化最小均方误差算法调整因子0<δ<<1,初始化水声信道冲激响应函数为h0=0;
步骤3:给定输入训练信号x和输出信号y,进行迭代计算得到信道估计hi,i为 迭代次数,初始化i=0,i=1,2,3,…a,a等于数据长度和信道长度之差:
||hi-1[K]||1表示对第i-1时刻的信道冲激响应函数的第K块向量进行L1范数运算; L1范数运算具体运算为对向量各元素的绝对值求和;
sign(hi-1)表示对第i-1时刻的信道冲激响应函数的符号运算,即当信道冲激响应函 数为负数,符号运算为-1,相反则为正数;
4、生成矩阵:Gi=diag(qi,1,qi,2,...,qi,L),其中:diag()表示将括号中的元素变为对 角矩阵;
所述参数设置μ=0.4。
所述参数设置α=-0.5。
所述参数设置γ=2×10-3。
所述δ=10-5。
有益效果
本发明提出的一种基于成比例归一化最小均方误差的稀疏水声信道估计方法,用以估计时域簇稀疏水声信道冲激响应函数。首先,针对时变水声信道进行簇稀疏建模, 以期获得水声信道结构化稀疏表达框架,在此基础上。采用逐块训练模式,设计一种 基于成比例归一化最小均方误差算法迭代寻优得到时域稀疏水声信道冲激响应函数, 将所估计得到的时域水声信道信息。该项发明适用于时变水声信道估计、水声通信等。 有益效果体现在:本发明基于水声信道的簇稀疏范数约束,有效地减少了冗余信道抽 头计算,使得本发明产生的水声信道估计结果具有更高的精度。
附图说明
图1是本发明方法与传统的成比例归一化最小均方误差算法(Improvedproportionate normalized least-mean-squares,IPNLMS)及其L0范数约束版本算法,最大 相关熵准则(Maximum correntropy criterion,MCC),成比例归一化最大相关熵准则(Proportionate normalized maximum correntropy criterion,PNMCC),混合范数约束的最 大相关熵准则(hybrid-norm constrained MCC,HNCMCC)等方法的对比结果图。
图2是传统的IPNLMS,IPNLMS-L0,MCC,PNMCC,HNCMCC方法与本发明 方法(IPNLMS-BL0)在对稀疏信道从h到-h的估计结果对比图。
图3是基于三种信道估计方法在四种不同信噪比条件下的学习曲线结果对比图。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
现结合附图对本发明作进一步的描述。参照图1,设置稀疏信道冲激响应函数长度为L=100,簇稀疏κ=1,块长度d=5,按照标准正态分布方式产生3000点的随机 高斯信号,算法参数设置β=10,μ=0.4,α=-0.5,δ=10-5,γ=2×10-3,
实施例步骤为:
步骤1:设参数L,K和d分别为水声信道冲激响应长度、信道分块个数和信道分 块的长度,它们之间的关系表示为L=Kd;设置迭代步长为μ,簇稀疏约束参数γ,α;
步骤2:设置稀疏光滑度调整参数β=10~12,基于成比例归一化最小均方误差算法调整因子0<δ<<1,初始化水声信道冲激响应函数为h0=0;
步骤3:给定输入训练信号x和输出信号y,进行迭代计算得到信道估计hi,i为 迭代次数,初始化i=0,i=1,2,3,…a,a等于数据长度和信道长度之差:
||hi-1[K]||1表示对第i-1时刻的信道冲激响应函数的第K块向量进行L1范数运算; L1范数运算具体运算为对向量各元素的绝对值求和;
sign(hi-1)表示对第i-1时刻的信道冲激响应函数的符号运算,即当信道冲激响应函 数为负数,符号运算为-1,相反则为正数;
4、生成矩阵:Gi=diag(qi,1,qi,2,...,qi,L),其中:diag()表示将括号中的元素变为对 角矩阵;
所述参数设置μ=0.4。
所述参数设置α=-0.5。
所述参数设置γ=2×10-3。
所述δ=10-5。
得到的结果如图1所示,从图1可以看出,与传统的IPNLMS,IPNLMS-L0,MCC,PNMCC,HNCMCC方法不同,本发明采用了簇稀疏范数约束结构,提高了对信道稀 疏结构的宽容性,在簇稀疏范数的约束条件下,本发明获得了更加精确的估计结果。
对时变稀疏水声信道的估计结果如图2所示,在算法进行跌到到1500点的时刻,使得信道发生突然变化,用以评估各算法在时变情况下对信道的追踪能力,可以看出, 信道的突变带来算法学习曲线的跳跃,然而,由于本发明能更好地结合水声信道簇稀 疏结构,因此获得更加精确的估计结果。
为进一步考察本接收信号不同信噪比对本发明所产生的不同影响,设置的接收信号信噪比从10dB不断变化到25dB。对比IPNLMS,IPNLMS-L0,以及本发明算法的 学习曲线结果如图3所示,可以看出,随着接收信号信噪比的提高,本发明方法取得 更为明显的信道估计精度。
Claims (5)
1.一种基于成比例归一化最小均方误差的稀疏水声信道估计方法,其特征在于步骤如下:
步骤1:设参数L,K和d分别为水声信道冲激响应长度、信道分块个数和信道分块的长度,它们之间的关系表示为L=Kd;设置迭代步长为μ,簇稀疏约束参数γ,α;
步骤2:设置稀疏光滑度调整参数β=10~12,基于成比例归一化最小均方误差算法调整因子0<δ<<1,初始化水声信道冲激响应函数为h0=0;
步骤3:给定输入训练信号x和输出信号y,进行迭代计算得到信道估计hi,i为迭代次数,初始化i=0,i=1,2,3,…a,a等于数据长度和信道长度之差:
||hi-1[K]||1表示对第i-1时刻的信道冲激响应函数的第K块向量进行L1范数运算;L1范数运算具体运算为对向量各元素的绝对值求和;
sign(hi-1)表示对第i-1时刻的信道冲激响应函数的符号运算,即当信道冲激响应函数为负数,符号运算为-1,相反则为正数;
4、生成矩阵:Gi=diag(qi,1,qi,2,...,qi,L),其中:diag()表示将括号中的元素变为对角矩阵;
2.根据权利要求1所述基于成比例归一化最小均方误差的稀疏水声信道估计方法,其特征在于:所述参数设置μ=0.4。
3.根据权利要求1所述基于成比例归一化最小均方误差的稀疏水声信道估计方法,其特征在于:所述参数设置α=-0.5。
4.根据权利要求1所述基于成比例归一化最小均方误差的稀疏水声信道估计方法,其特征在于:所述参数设置γ=2×10-3。
5.根据权利要求1所述基于成比例归一化最小均方误差的稀疏水声信道估计方法,其特征在于:所述δ=10-5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010444429.5A CN111711584B (zh) | 2020-05-22 | 2020-05-22 | 基于成比例归一化最小均方误差的稀疏水声信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010444429.5A CN111711584B (zh) | 2020-05-22 | 2020-05-22 | 基于成比例归一化最小均方误差的稀疏水声信道估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111711584A true CN111711584A (zh) | 2020-09-25 |
CN111711584B CN111711584B (zh) | 2022-06-07 |
Family
ID=72538276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010444429.5A Expired - Fee Related CN111711584B (zh) | 2020-05-22 | 2020-05-22 | 基于成比例归一化最小均方误差的稀疏水声信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111711584B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114465851A (zh) * | 2021-12-25 | 2022-05-10 | 西北工业大学 | 优化核宽最大箕舌线准则的簇稀疏水声信道估计方法 |
CN114584239A (zh) * | 2022-03-07 | 2022-06-03 | 大连理工大学 | 一种基于学习去噪的otfs水声通信稀疏信道估计方法 |
CN116248206A (zh) * | 2023-02-08 | 2023-06-09 | 中国科学院声学研究所 | 一种水声通信信道质量量化分析方法及系统 |
CN117081894A (zh) * | 2023-07-17 | 2023-11-17 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法 |
CN118764009A (zh) * | 2024-09-06 | 2024-10-11 | 西北工业大学 | 一种用于水声信道估计的自适应滤波方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013138248A1 (en) * | 2012-03-13 | 2013-09-19 | Qualcomm Incorporated | Fast and robust estimation of jointly sparse channels |
US20140098704A1 (en) * | 2012-10-09 | 2014-04-10 | Xiao-an Wang | Channel Estimation By Time-Domain Parameter Extraction |
CN107395536A (zh) * | 2017-07-20 | 2017-11-24 | 西北工业大学 | 对多途环境下的水声信道冲激响应函数进行估计的方法 |
CN107592115A (zh) * | 2017-09-12 | 2018-01-16 | 西北工业大学 | 一种基于非均匀范数约束的稀疏信号恢复方法 |
WO2018179557A1 (en) * | 2017-03-31 | 2018-10-04 | Mitsubishi Electric Corporation | Method for decoding symbol and receiver for receiving and decoding symbol |
CN108833312A (zh) * | 2018-06-08 | 2018-11-16 | 西北工业大学 | 一种基于时延多普勒域的时变稀疏水声信道估计方法 |
CN110071881A (zh) * | 2019-04-26 | 2019-07-30 | 北京理工大学 | 一种自适应开销的活跃用户检测和信道估计方法 |
-
2020
- 2020-05-22 CN CN202010444429.5A patent/CN111711584B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013138248A1 (en) * | 2012-03-13 | 2013-09-19 | Qualcomm Incorporated | Fast and robust estimation of jointly sparse channels |
US20140098704A1 (en) * | 2012-10-09 | 2014-04-10 | Xiao-an Wang | Channel Estimation By Time-Domain Parameter Extraction |
WO2018179557A1 (en) * | 2017-03-31 | 2018-10-04 | Mitsubishi Electric Corporation | Method for decoding symbol and receiver for receiving and decoding symbol |
CN107395536A (zh) * | 2017-07-20 | 2017-11-24 | 西北工业大学 | 对多途环境下的水声信道冲激响应函数进行估计的方法 |
CN107592115A (zh) * | 2017-09-12 | 2018-01-16 | 西北工业大学 | 一种基于非均匀范数约束的稀疏信号恢复方法 |
CN108833312A (zh) * | 2018-06-08 | 2018-11-16 | 西北工业大学 | 一种基于时延多普勒域的时变稀疏水声信道估计方法 |
CN110071881A (zh) * | 2019-04-26 | 2019-07-30 | 北京理工大学 | 一种自适应开销的活跃用户检测和信道估计方法 |
Non-Patent Citations (2)
Title |
---|
CUI ZHI-TAO等: "A New Algorithm to Estimate Mixing-Matrix of Underdetermined Blind", 《2012 EIGHTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL INTELLIGENCE AND SECURITY》, 11 January 2013 (2013-01-11) * |
伍等: "部分范数约束的稀疏恢复算法及其在单载波水声数据遥测中的应用", 《通信学报》, 25 June 2018 (2018-06-25) * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114465851A (zh) * | 2021-12-25 | 2022-05-10 | 西北工业大学 | 优化核宽最大箕舌线准则的簇稀疏水声信道估计方法 |
CN114584239A (zh) * | 2022-03-07 | 2022-06-03 | 大连理工大学 | 一种基于学习去噪的otfs水声通信稀疏信道估计方法 |
CN114584239B (zh) * | 2022-03-07 | 2022-11-18 | 大连理工大学 | 一种基于学习去噪的otfs水声通信稀疏信道估计方法 |
CN116248206A (zh) * | 2023-02-08 | 2023-06-09 | 中国科学院声学研究所 | 一种水声通信信道质量量化分析方法及系统 |
CN116248206B (zh) * | 2023-02-08 | 2023-09-22 | 中国科学院声学研究所 | 一种水声通信信道质量量化分析方法及系统 |
CN117081894A (zh) * | 2023-07-17 | 2023-11-17 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法 |
CN117081894B (zh) * | 2023-07-17 | 2024-05-24 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法及系统 |
CN118764009A (zh) * | 2024-09-06 | 2024-10-11 | 西北工业大学 | 一种用于水声信道估计的自适应滤波方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111711584B (zh) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111711584B (zh) | 基于成比例归一化最小均方误差的稀疏水声信道估计方法 | |
CN109116834B (zh) | 一种基于深度学习的间歇过程故障检测方法 | |
CN111523236A (zh) | 基于Koopman算子的压电陶瓷迟滞模型线性化辨识方法 | |
Ghauri et al. | System identification using LMS, NLMS and RLS | |
CN106875002A (zh) | 基于梯度下降法与广义逆的复值神经网络训练方法 | |
CN108566257A (zh) | 一种基于反向传播神经网络的信号恢复方法 | |
CN111812980B (zh) | 基于未知输入观测器的离散切换系统的鲁棒故障估计方法 | |
CN112561119A (zh) | 一种使用arima-rnn组合模型的云服务器资源性能预测方法 | |
CN103888391B (zh) | 基于双Sigmoid混沌神经网络的信号盲检测方法 | |
CN107395536B (zh) | 对多途环境下的水声信道冲激响应函数进行估计的方法 | |
CN107276561A (zh) | 基于量化核最小均方误差的Hammerstein系统辨识方法 | |
CN107592115B (zh) | 一种基于非均匀范数约束的稀疏信号恢复方法 | |
CN110121184A (zh) | 一种非负约束的分布式在线自适应扩散组合系数优化方法 | |
CN109217844B (zh) | 基于预训练随机傅里叶特征核lms的超参数优化方法 | |
CN114462241B (zh) | 多速率系统的全对称多胞形集员状态估计方法 | |
CN111555994B (zh) | 一种基于最大箕舌线准则算法的簇稀疏信道估计方法 | |
CN113033083B (zh) | 一种基于密度峰值聚类径向基神经网络波达方向估计方法 | |
CN112054973A (zh) | 一种最小均方误差稀疏水声信道估计方法 | |
CN114415110A (zh) | 一种非负稀疏贝叶斯学习的直接定位方法 | |
CN103475608A (zh) | 模拟退火与果蝇混合优化小波广义离散多模盲均衡方法 | |
CN111416595B (zh) | 一种基于多核融合的大数据滤波方法 | |
CN109474258B (zh) | 基于核极化策略的随机傅立叶特征核lms的核参数优化方法 | |
CN109347770A (zh) | 一种基于压缩感知理论的信道估计方法及装置 | |
CN111505576B (zh) | 一种针对tdoa定位的传感器选择方法 | |
CN110598334B (zh) | 基于协同衍化相关熵极限学习机的性能退化趋势预测方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220607 |
|
CF01 | Termination of patent right due to non-payment of annual fee |