CN110677318A - 一种基于线性调频z变换的水声信道时延估计方法 - Google Patents
一种基于线性调频z变换的水声信道时延估计方法 Download PDFInfo
- Publication number
- CN110677318A CN110677318A CN201910881611.4A CN201910881611A CN110677318A CN 110677318 A CN110677318 A CN 110677318A CN 201910881611 A CN201910881611 A CN 201910881611A CN 110677318 A CN110677318 A CN 110677318A
- Authority
- CN
- China
- Prior art keywords
- cluster
- channel
- time delay
- estimation
- search
- 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.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L43/00—Arrangements for monitoring or testing data switching networks
- H04L43/08—Monitoring or testing based on specific metrics, e.g. QoS, energy consumption or environmental parameters
- H04L43/0852—Delays
-
- 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
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2602—Signal structure
- H04L27/261—Details of reference signals
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Environmental & Geological Engineering (AREA)
- Noise Elimination (AREA)
Abstract
本发明提供的是一种基于线性调频z变换的水声信道时延估计方法。一、先在簇状信道内进行基于CZT的OMP估计,当前簇的最后一条路径满足预设的簇间跳转幅度门限时,跳转到下一簇内继续搜索,直到搜索出簇内的全部满足门限的路径。二、将上步每簇信道的时域估计范围进行一定程度的扩大,形成新的簇状信道,再将新的簇状信道拼接为一个新的全局信道,运用CZT在这个拼接后的全局信道内搜索剩余幅度较小的路径。本发明的二步搜索策略同样也适用于其他水声信道,只要已知关于信道分布范围的先验信息即可。本发明可在降低计算复杂度的同时提高时延估计精度,对均匀导频与非均匀导频OFDM系统均适用,可用于实现高效高精度低复杂度的信道估计。
Description
技术领域
本发明涉及的是一种水声信道时延估计方法。
背景技术
水声信道通常具有较大的时延扩展与严重的多普勒效应,因此水声信道具有时频域双选性,同时也具有时频域稀疏性。近年来,由于正交频分复用(OFDM)技术的高带宽利用率及低复杂度使其广泛地应用于水声通信系统中,实际应用中大部分水声OFDM系统采用相干接收机,这使得信道估计的精度尤其是信道时延估计的精度,变得十分重要。
作为具有较少计算量的贪婪算法,正交匹配追踪(OMP)成为压缩感知(CS)分支中最受欢迎的算法之一。但是OMP的主要计算量集中在字典矩阵的运算上,若想提高估计精度,则要增加字典矩阵的规模,过大的字典矩阵规模导致OMP算法需要极高的计算复杂度。由于水下传感器节点的成本以及布放操作较为昂贵,能量补充极其困难,这就导致了估计精度和计算复杂度的权衡问题,因此需要研究低复杂度的水声信道估计技术。目前国内外许多学者都对此进行了研究,但大部分技术的计算复杂度仍取决于时延网格精度。
不过,有学者发现水声通信信道的路径分布,通常成簇状分布,许多由反射或折射产生的弱径集中在强径附近,形成簇结构。这意味着在水声信道估计中存在可压缩的成分,协调好可压缩成分与不可压缩成分就可以实现效率与精度上的双赢。
公开号为CN106027445A的专利文件中公开了一种水声块结构稀疏特性的信道估计方法。该方法从降低水声信道的稀疏度入手,通过降低OMP所需的导频数量来降低估计所需的计算复杂度,但没有利用水声信道的簇状特性进行进一步改进。
发明内容
本发明的目的在于提供一种既能提高时延估计精度又能降低计算复杂度的基于线性调频z变换的水声信道时延估计方法。
本发明的目的是这样实现的:
步骤2:输入:信道共有Q条多径,分为NC簇,每簇记为Ci;簇间跳转门限比率χ1,全局搜索迭代停止条件χ2δ2,δ2为噪声功率;
步骤3:分簇搜索:
(3-6)更新变量:l=1,判断是否已搜索完所有预设簇,若未搜完,则i=i+1,继续执行步骤(3-7);
(3-7)判断是否已搜索完所有预设簇,若已搜完,则执行步骤4;若未搜完,则执行步骤3;
步骤4:全局搜索:
(4-5)当q>Q或||rq||2<χ2δ2时,执行步骤5,否则继续执行步骤4;
本发明还可以包括:
1.步骤2中多径信道的分簇划分,是通过在OFDM符号前加入线性调频信号,结合k均值聚类算法对信道进行预处理实现。
2.步骤4中的全局搜索范围是指将步骤3中每簇信道Ci对应的搜索范围按一定比例扩大再按序拼接而成的新的搜索范围。
本发明提出了一种既能提高时延估计精度又能降低计算复杂度的水声信道时延估计方法。与公开号为CN106027445A的专利文件中公开的技术方案相比,本发明直接从OMP算法的原理入手,借助水声信道的内部分布特征,引入线性调频z变换(CZT)来降低OMP所需的计算复杂度。
虽然传统快速傅里叶逆变换(IFFT)算法能够有效降低正交匹配追踪(OMP)信道估计的复杂度,但如果要获得高精度的时延估计,全范围的IFFT仍会带来庞大的计算量,导致资源的浪费。而线性调频z变换(CZT)可以将这种资源的浪费转化为有效范围内更高精度的时延估计。本发明紧密结合水声信道簇状分布的特征,提出基于CZT的两步搜索时延估计模型,第一步,先在簇状信道内进行基于CZT的OMP估计,当前簇的最后一条路径满足预设的簇间跳转幅度门限时,跳转到下一簇内继续搜索,直到搜索出簇内的全部满足门限的路径。虽然簇的划分,本身具有去噪并提升估计精度的作用,但为了达到与OMP相似的性能,在接下来的第二步中,将上步每簇信道的时域估计范围进行一定程度的扩大,形成新的簇状信道,再将新的簇状信道拼接为一个新的全局信道,运用CZT在这个拼接后的全局信道内搜索剩余幅度较小的路径,以达到精度补偿的效果。本发明虽说是基于簇状信道提出的,但是该二步搜索策略同样也适用于其他水声信道,只要已知关于信道分布范围的先验信息即可。本发明可在降低计算复杂度的同时提高时延估计精度,对均匀导频与非均匀导频OFDM系统均适用,可用于实现高效高精度低复杂度的信道估计模型。
本发明的有益效果主要体现在:
1.本发明针对均匀导频与非均匀导频OFDM系统提出一种基于线性调频z变换的OMP内积计算方法,利用该方法并结合水声信道的簇状特性可提高时延估计精度,同时降低了内积计算时的复杂度。
2.本发明提出一种两级搜索策略保证了信道时延估计的完整性,不会遗漏某些不在簇内的弱径。
3.本发明可在降低计算复杂度的同时提高时延估计精度,对非均匀与均匀导频的OFDM系统均适用,可用于实现高效高精度低复杂度的信道估计模型中,为水声通信装置或系统实现节省能源与高效运行上的双赢。
附图说明
图1为水声信道分簇特征的示意图;
图2为基于分簇信道的簇间跳转搜索示意图;
图3为基于两级搜索策略的时延搜索示意图;
图4为基于线性调频z变换的信道估计模型。
具体实施方式
下面举例对本发明做更详细的描述。
A.CP-OFDM系统与水声多径信道模型
本发明所涉及的OFDM系统采用循环前缀模式。假定一个OFDM符号块包含K个子载波,每个子载波上的发送符号为s[k]。一个OFDM符号块时长为T,循环前缀的时长为Tcp,载波频率为fc。因此,待发送的OFDM符号时域表达式如下,
假设信道在一个符号持续时间内保持不变,本文的信道模型为具有L条路径的多径信道,每条径的幅度与时延分别为Al和τl(l=0,1,...,L-1),并假设系统经预处理后只剩余残余多普勒偏移fd。接收端经下变频和滤波处理后的信号为
其中ε=fdT表示归一化多普勒因子。
通常,经多普勒补偿之后,残余的多普勒频移也会被限制,因此(4)式中的矩阵Cε可近似为主对角线附近为非零项的带状矩阵。进一步假设,导频间隔足够大,载波间干扰可被忽略。综上,可得到(3)式中导频项的频域输入输出关系式:
对导频位置的频域观测量进行补偿,得到如下表达式:
B.基于线性调频z变换的时延估计模型
假设簇状信道共有Q条多径,分为NC簇,每簇记为Ci,每簇内的路径数目记为为了充分利用分簇信道的可压缩性,本发明提出簇间跳转搜索模型,传统OMP每次计算内积之后,都会在时域全局搜索内积峰值位置,造成计算资源的浪费。为此,本发明仅在簇内进行峰值搜索,且每次只搜索一簇,极大的节省了计算资源。而且簇的划分,相当于为信道估计添加了去噪窗函数,使得最终的信道估计精度优于无去噪效果的标准OMP。基于分簇信道的内积φ计算公式如下:
为了防止分簇搜索过程中,陷入簇内死循环,本发明设置了簇间跳转门限,如图2所示,如果当前簇内所估计的路径幅度小于该门限,则跳转到下一簇内进行搜索,从而实现了簇间跳转。
利用快速傅里叶逆变换(IFFT)对标准OMP计算内积时,考虑的是全局时域范围,而本发明的时域搜索范围由全局缩小到簇内,若依然选用IFFT作为计算内积的手段则浪费了计算资源,所以本文提出使用线性调频z变换(CZT)仅求解簇内时域范围对应的内积值:
nstart是每簇信道的起始候选点在全局的索引位置,n=0,1,...,MCi,MCi为每簇信道Ci内的时延候选点数量。N表示频域观测向量中参与内积计算的元素数量。
C.两步搜索模型
在OMP中,如果两条路径在时延上非常接近,间隔小于内积函数的主瓣宽度时,OMP将无法通过内积峰值识别出真正的路径时延,而是搜索出几个与真实时延接近且幅度较大的路径,和一些时延与真值偏差较大的且幅度较小的补偿路径,这些小幅度路径分布在真实时延位置左右两侧,有时会落在预设簇信道之外。
因此,为了避免遗漏某些簇外的路径,本发明在时延搜索策略上提出了两步搜索的新思路,如图3所示,第一步,先在簇状信道内进行基于CZT的OMP估计,当前簇的最后一条路径满足预设的簇间跳转幅度门限时,跳转到下一簇内继续搜索,直到搜索出簇内的全部满足门限的路径。为了达到与OMP相似的性能,在接下来的第二步中,将上步每簇信道的时域估计范围进行一定程度的扩大,形成新的簇状信道,再将新的簇状信道拼接为一个新的全局信道,运用CZT在这个拼接后的全局信道内搜索剩余幅度较小的路径,以达到精度补偿的效果。
本发明虽说是基于簇状信道提出的,但是该二步搜索策略也适用于其他水声信道,只要已知关于信道多径分布范围的先验信息即可。
D.基于线性调频Z变换的的两步时延搜索模型
构建流程如图4所示,具体步骤如下:
步骤2:输入:信道共有Q条多径,分为NC簇,每簇记为Ci;簇间跳转门限比率χ1,全局搜索迭代停止条件χ2δ2,δ2为噪声功率;多径信道的分簇划分,是通过在OFDM符号前加入线性调频信号,结合k均值聚类算法(k-means clustering algorithm)对信道进行预处理实现的。
步骤3:分簇搜索模块:
(3-1):利用线性调频z变换估计簇Ci内第l条路径,时延估计结果记为
(3-6):更新变量:l=1,判断是否已搜索完所有预设簇,若未搜完,则i=i+1,继续执行步骤(3-7);
(3-7):判断是否已搜索完所有预设簇,若已搜完,则执行步骤4;若未搜完,则执行步骤3;
步骤4:全局搜索模块,全局搜索范围是指将步骤3中每簇信道Ci对应的搜索范围按一定比例扩大再按序拼接而成的新的搜索范围:
(4-5):当q>Q或||rq||2<χ2δ2时,执行步骤5,否则继续执行步骤4;
步骤5:输出信道估计结果
以上步骤流程仅为说明本发明的技术思想,并不用以限制本发明,凡是在本发明提出的技术方案、技术思想、引进方法上所做的任何改动、改进等,均落入本发明保护范围之内。
Claims (3)
1.一种基于线性调频z变换的水声信道时延估计方法,其特征是:
步骤2:输入:信道共有Q条多径,分为NC簇,每簇记为Ci;簇间跳转门限比率χ1,全局搜索迭代停止条件χ2δ2,δ2为噪声功率;
步骤3:分簇搜索:
(3-6)更新变量:l=1,判断是否已搜索完所有预设簇,若未搜完,则i=i+1,继续执行步骤(3-7);
(3-7)判断是否已搜索完所有预设簇,若已搜完,则执行步骤4;若未搜完,则执行步骤3;
步骤4:全局搜索:
(4-5)当q>Q或||rq||2<χ2δ2时,执行步骤5,否则继续执行步骤4;
2.根据权利要求1所述的基于线性调频z变换的水声信道时延估计方法,其特征是:步骤2中多径信道的分簇划分,是通过在OFDM符号前加入线性调频信号,结合k均值聚类算法对信道进行预处理实现。
3.根据权利要求1或2所述的基于线性调频z变换的水声信道时延估计方法,其特征是:步骤4中的全局搜索范围是指将步骤3中每簇信道Ci对应的搜索范围按一定比例扩大再按序拼接而成的新的搜索范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910881611.4A CN110677318A (zh) | 2019-09-18 | 2019-09-18 | 一种基于线性调频z变换的水声信道时延估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910881611.4A CN110677318A (zh) | 2019-09-18 | 2019-09-18 | 一种基于线性调频z变换的水声信道时延估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110677318A true CN110677318A (zh) | 2020-01-10 |
Family
ID=69078165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910881611.4A Pending CN110677318A (zh) | 2019-09-18 | 2019-09-18 | 一种基于线性调频z变换的水声信道时延估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110677318A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117081894A (zh) * | 2023-07-17 | 2023-11-17 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法 |
CN117240666A (zh) * | 2023-08-21 | 2023-12-15 | 中国科学院声学研究所 | 一种基于簇自注意力机制的神经网络水声信道估计方法 |
CN117081894B (zh) * | 2023-07-17 | 2024-05-24 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103117978A (zh) * | 2013-01-24 | 2013-05-22 | 中国科学院声学研究所 | 基于多普勒频移补偿的水声通信多载波匹配方法及装置 |
US20140279727A1 (en) * | 2013-03-15 | 2014-09-18 | William Marsh Rice University | Sparse Factor Analysis for Analysis of User Content Preferences |
CN106027445A (zh) * | 2016-05-13 | 2016-10-12 | 江苏科技大学 | 一种水声块结构稀疏特性的信道估计方法 |
CN109150772A (zh) * | 2018-07-13 | 2019-01-04 | 哈尔滨工程大学 | 一种基于正交匹配追踪的水声信道时延估计方法 |
-
2019
- 2019-09-18 CN CN201910881611.4A patent/CN110677318A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103117978A (zh) * | 2013-01-24 | 2013-05-22 | 中国科学院声学研究所 | 基于多普勒频移补偿的水声通信多载波匹配方法及装置 |
US20140279727A1 (en) * | 2013-03-15 | 2014-09-18 | William Marsh Rice University | Sparse Factor Analysis for Analysis of User Content Preferences |
CN106027445A (zh) * | 2016-05-13 | 2016-10-12 | 江苏科技大学 | 一种水声块结构稀疏特性的信道估计方法 |
CN109150772A (zh) * | 2018-07-13 | 2019-01-04 | 哈尔滨工程大学 | 一种基于正交匹配追踪的水声信道时延估计方法 |
Non-Patent Citations (1)
Title |
---|
GANG QIAO 等: "Chirp Z-Transform Based Sparse Channel Estimation for Underwater Acoustic OFDM in Clustered Channels", 《IEEE》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117081894A (zh) * | 2023-07-17 | 2023-11-17 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法 |
CN117081894B (zh) * | 2023-07-17 | 2024-05-24 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法及系统 |
CN117240666A (zh) * | 2023-08-21 | 2023-12-15 | 中国科学院声学研究所 | 一种基于簇自注意力机制的神经网络水声信道估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110290581B (zh) | 一种5g系统中的快速时频同步方法及终端 | |
CN109150772B (zh) | 一种基于正交匹配追踪的水声信道时延估计方法 | |
JP2009049996A (ja) | タイミング同期方法およびその装置、並びにプリアンブルおよびその生成方法と装置 | |
CN102075486A (zh) | 一种ofdm系统的同步方法 | |
CN112003809B (zh) | 基于运动平台的多载波水声通信多普勒估计与补偿方法 | |
CN101277288A (zh) | 正交频分复用系统的频率同步方法 | |
WO2021027590A1 (zh) | 突发ofdm数据传输的帧同步方法及系统 | |
CN102377726B (zh) | Ofdm系统的定时同步方法 | |
KR100754840B1 (ko) | 직교주파수 다중분할 방식을 이용한 이동통신 시스템의주파수 동기화 방법 | |
CN110677318A (zh) | 一种基于线性调频z变换的水声信道时延估计方法 | |
CN110059401B (zh) | 一种ofdm系统水声信道冲激响应重构方法 | |
CN110798416A (zh) | OFDM系统中基于局部搜索Capon的CFO估计算法 | |
CN109561042B (zh) | 一种ofdm系统接收机的定时频率同步方法 | |
CN111901266A (zh) | 面向低信噪比和高动态场景的载波多普勒频偏捕获方法 | |
CN111726309A (zh) | 用于移动正交频分复用系统的信道估计方法及其估计装置 | |
CN113259281B (zh) | Dmrs和ptrs联合信道估计方法、装置及接收机 | |
CN111726308B (zh) | 基于频响预插值的正交匹配追踪信道估计方法 | |
CN111953626B (zh) | 一种面向正交切普复用调制的低复杂度频选信道估计方法 | |
KR101501334B1 (ko) | Ofdm 시스템을 기반으로 하는 주파수 오차 추정 방법 및 그 장치 | |
KR100653724B1 (ko) | 직교주파수다중분할 방식을 이용한 이동통신시스템의주파수 옵셋의 추정방법 | |
CN108737313B (zh) | 一种残余载波频偏抑制方法 | |
CN112235216A (zh) | 一种基于禁忌搜索的压缩感知导频设计方法 | |
CN114928520B (zh) | 一种用于补偿频偏的广义频分复用时频同步方法 | |
CN107171776B (zh) | 基于fbmc调制在指定时间区间内生成多载波导频序列的方法 | |
CN111131116A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200110 |