CN114710386A - 一种水声ofdm通信非均匀多普勒频偏抑制方法和系统 - Google Patents
一种水声ofdm通信非均匀多普勒频偏抑制方法和系统 Download PDFInfo
- Publication number
- CN114710386A CN114710386A CN202210329421.3A CN202210329421A CN114710386A CN 114710386 A CN114710386 A CN 114710386A CN 202210329421 A CN202210329421 A CN 202210329421A CN 114710386 A CN114710386 A CN 114710386A
- Authority
- CN
- China
- Prior art keywords
- signal
- length
- tail
- frequency offset
- baseband
- 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 32
- 238000004891 communication Methods 0.000 title claims abstract description 26
- 230000001360 synchronised effect Effects 0.000 claims abstract description 22
- 238000012952 Resampling Methods 0.000 claims abstract description 19
- 230000005540 biological transmission Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 abstract description 25
- 230000008859 change Effects 0.000 abstract description 7
- 238000006243 chemical reaction Methods 0.000 description 12
- 230000001427 coherent effect Effects 0.000 description 11
- 230000000694 effects Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 239000000969 carrier Substances 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- 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/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
- H04L27/2691—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation involving interference determination or cancellation
-
- 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/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明属于水声无线通信技术领域,更具体地,涉及一种水声OFDM(OrthogonalFrequency Division Multiplexing,正交频分复用)通信非均匀多普勒频偏抑制方法和系统。
背景技术
OFDM技术通过将高速串行传输的数据流转换为低速并行传输的数据流以及添加循环前缀等操作,可以有效抑制多径时延造成的符号间干扰。此外,OFDM通过子载波间正交性进行解调,因此允许子载波频谱相互重叠,从而可以有效提升系统频谱效率。OFDM技术所具有的这些优势可以克服水声信道的很多恶劣条件,因此被认为是实现高速水声通信的一种重要技术手段。
对于OFDM通信系统,由于其解调依赖于子载波间的正交性,对于频率偏移十分敏感,而水声信道具有强烈的多普勒效应,会导致载波发生频率偏移,从而引起载波间干扰,使得OFDM载波间的正交性被严重破坏,系统信号检测性能大幅下降,甚至无法工作。因此,想要实现数据的精准检测,需要对载波频偏进行抑制。
但是,由于载波频率与信号带宽的比值很小,水声信道本质上呈现出一种宽带的特性,信号经过信道在多普勒效应影响下会产生非均匀频率偏移,即各个子载波频偏大小与子载波自身频率相关,同时,信号在时域会发生展宽或压缩的尺度变化,很难对信号时域尺度变化以及非均匀频率偏移进行统一补偿。因此,一种性能优良的水声OFDM通信非均匀多普勒频偏抑制方法亟待研究。
现有OFDM通信系统均采用块多普勒估计算法来实现多普勒因子估计(例如:CN107426120A、CN107454024A等),其思想是利用原始的同步信号与整个接收信号进行相关运算。采用原始同步信号与整个接收信号进行相关运算,相关运算结果会在首部和尾部产生两个峰值,这两个峰值之间的长度即为接收信号的长度。然而,考虑到水声信道的恶劣条件,信号在经过信道传输之后,会受到严重的多径效应、多普勒效应以及噪声等影响,产生畸变,两端的同步信号会受到相应信道影响,导致与原始的同步信号之间相关性降低,出现相关峰峰值下降、分裂等问题,致使检测性能大幅下降,从而影响多普勒频偏抑制的效果。
发明内容
针对现有技术的缺陷和改进需求,本发明提供了一种水声OFDM通信非均匀多普勒频偏抑制方法和系统,其目的在于利用一种新的信道自适应多普勒估计方法,对信道多普勒因子进行高精准估计,从而确定重采样率,进而通过时域重采样操作对信号时域尺度变化进行消除,同时补偿各子载波的非均匀多普勒频率偏移,使得载波正交性得以恢复,显著提升系统的信号检测性能。
为实现上述目的,按照本发明的第一方面,提供了一种水声OFDM通信非均匀多普勒频偏抑制方法,发送机在发送信号的首部和尾部分别添加LFM同步信号,该方法应用于接收机,包括:
S1.将LFM同步信号长度Tl、多普勒因子最大可能值αmax和原始发送信号长度Tt作为先验信息,对基带接收信号的首部和尾部同步信号分别进行截取;
S2.截取到的首尾两个同步信号进行互相关运算,得到相关峰位置,作为基带接收数据信号长度相对于原始发送信号长度的变化量;
S6.对重采样后的基带接收信号去除掉首部和尾部长度为Tl的信号,得到频偏补偿后的接收数据信号。
优选地,步骤S2中,对截取得到的首尾同步信号进行互相关运算:
其中,C(τ)为互相关运算的结果,τ为信号时延,y(t)为基带接收信号,y*(t)为y(t)的共轭信号。
优选地,步骤S3中,估计出接收数据信号长度:
其中,Δτ为相关峰位置。
为实现上述目的,按照本发明的第二方面,提供了一种水声OFDM通信非均匀多普勒频偏抑制系统,包括:计算机可读存储介质和处理器;
所述计算机可读存储介质用于存储可执行指令;
所述处理器用于读取所述计算机可读存储介质中存储的可执行指令,执行第一方面所述的水声OFDM通信非均匀多普勒频偏抑制方法。
为实现上述目的,按照本发明的第三方面,提供了一种计算机可读存储介质,包括存储的计算机程序;
所述计算机程序被处理器执行时,控制所述计算机可读存储介质所在设备执行第一方面所述的水声OFDM通信非均匀多普勒频偏抑制方法。
总体而言,通过本发明所构思的以上技术方案,能够取得以下有益效果:
本发明利用同步信号长度、发送信号长度以及当前应用场景下最大多普勒因子可能值作为先验信息,截取接收信号的首尾同步信号,截取得到的两个同步信号进行相关运算,根据唯一的相关峰位置,确定多普勒因子估计值。由于截取所得的两段同步信号受到相似的信道特性影响,仍具有良好的相关特性,可以自动匹配不同的水声多径信道,因此本发明相较于传统的块多普勒估计方法,精度更高,自适应能力更强。多普勒因子估计精度的提升,使得OFDM通信系统的非均匀多普勒频偏抑制性能进一步提升,有效补偿了由水声信道多普勒效应所引起的信号时域尺度变化以及载波非均匀频率偏移,使得OFDM通信系统的载波间正交性得以恢复,提升系统的检测性能。
此外,由于接收机与发送机相对运动的不同情况会使得接收信号长度展宽或者收缩,为确保无论在何种情况下都可以截取得到同步信号,因此利用最大多普勒因子可能值这一先验信息对信号进行截取,从而建立统一的算法架构。
最后,本发明只截取首尾两部分信号进行互相关运算,大大降低了算法的计算开销。本发明在高多普勒因子和大子载波数的情况下,系统性能可以得到显著提升。
附图说明
图1是本发明中系统发送数据帧结构示意图;
图2是本发明中截取同步信号操作示意图,包括信号尺度扩展(α<0)和信号尺度压缩(α>0)两种情况;
图3是本发明测试时采用的仿真水声信道路径增益示意图;
图4是本发明与现有方法多普勒因子估计性能的测试结果图;
图5是本发明与块多普勒估计方法以及无重采样处理方法在差分相干检测系统中的信号检测性能的测试结果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供了一种水声OFDM通信非均匀多普勒频偏抑制方法,具体包括以下处理:
步骤一、对于差分相干检测系统,在发射机端,对原始数据进行星座映射、差分编码、串并转换、IFFT、并串转换以及添加LFM(线性调频,Linear Frequency Modulation)同步信号等操作,得到基带时域发送信号,经上变频后得到射频发送信号,由发射机发出。而对于相干检测系统则不需要进行差分编码。发送数据帧如图1所示。
具体地,对于差分相干检测系统,在发送机对原始数据进行QPSK调制,得到原始调制符号bk,然后对原始调制符号进行差分编码,编码规则为:
其中,K为子载波数量,从而到差分编码符号dk。对其进行串并转换后再进行IFFT调制搭载到各个子载波,再由并串转换操作得到发送基带OFDM信号s(t):
其中,T′=T+Tg,T为OFDM信号长度,Tg为保护间隔长度,Δf=1/T为载波间隔。f0为最低载波频率,g(t)=1,ift∈[0,T′],其他区间为0。
进一步地,在生成的OFDM信号两端分别添加长度为Tl的LFM同步信号,构成完整的基带发送信号,经上变频得到射频发送信号。
特别地,在相干检测系统中,不包含差分编码的步骤,其他步骤与上述差分相干检测系统保持一致。
步骤二、信号经过水声信道到达接收机,对接收到的射频信号进行下变频,得到基带接收信号,依据发送时LFM同步信号长度、发送信号长度以及当前应用场景下最大多普勒因子可能的取值等先验信息,对接收信号首部和尾部进行截取,并利用截取所得信号进行互相关运算,如图2所示。
具体地,OFDM射频信号在经过水声信道后到达接收机,对接收到的射频信号进行下变频处理后,得到基带接收信号y(t):
其中,Ap和τp分别表示第p条路径的增益和时延,α表示多普勒因子,n(t)为加性噪声,fk=f0+kΔf为第k个子载波的频率。
依据已知的发送同步信号长度Tl,从接收信号首部截取得到一段同步信号。而在进行尾部同步信号截取时,考虑到由于多普勒因子的取值可能为正或为负,对应接收信号可能发生尺度压缩或者扩展,为保证一定能够截取得到尾部的同步信号,设当前应用场景下最大多普勒因子可能值为αmax,依据这一先验信息,从位置处向后截取,得到尾部同步信号,其中,Tt为原始发送信号长度。然后,对截取得到的首尾同步信号进行互相关运算:
由于截取所得的首部和尾部信号包含相同的LFM同步信号部分,且受到相似的信道特性影响,仍然具有良好的相关特性,因此在运算结果中会产生唯一的一个尖锐的相关峰。
步骤三、根据互相关运算结果得到相关峰位置,从而估计出接收信号的长度,进而估计出多普勒因子的大小。
具体地,截取得到的首尾同步信号互相关运算会产生一个尖锐的相关峰,相关峰位置为Δτ。由图2可知,接收数据信号长度满足:
而根据互相关运算可知:
ΔT=Δτ-Tl
因此接收数据信号长度为:
由于信号受多普勒效应的影响在时域会发生压缩或展宽的尺度变化,接收数据信号与原始发送信号长度之间满足:
因此可以估计出多普勒因子的大小为:
步骤五、对接收数据信号进行串并转换、FFT解调、并串转换、差分检测以及星座解调等操作,恢复出原始的数据符号。
因此,信号的时域尺度变化及非均匀多普勒频率偏移都得到有效补偿。
进一步地,对重采样后得到的信号,去除掉首部和尾部长度为Tl同步信号,接着进行串并转换、FFT解调、并串转换等操作得到解调符号对于差分相干检测系统,利用当前符号与前一个符号之间具有相干性,进行差分检测,得到检测符号:
然后对得到的符号进行符号判决,判决方式为:计算差分检测符号与QPSK星座图上各星座点之间的均方误差,取均方误差最小的星座点作为判决结果。之后再进行QPSK解调,得到原始的数据符号。而对于相干检测系统,则需要将差分检测的公式替换为相干检测系统中的对应公式,其他步骤与上述差分相干检测系统保持一致。
结合以下实验结果对本发明的使用效果进一步说明:
测试时所采用的仿真水声信道,其路径增益如图3所示。在多普勒因子高达0.001,子载波数为1024的条件下,分别采用经典的块多普勒估计算法和本发明提出的方法对信道多普勒因子进行估计,两种方法的估计性能如图4所示。图5展示了在差分相干检测系统中,基于这两种方法进行重采样的差分水声OFDM通信系统的信号检测性能以及不进行重采样处理的信号检测性能的测试结果。
从仿真实验结果可以看出,在高多普勒因子以及大子载波数的条件下,本发明实现的信道自适应多普勒估计方法相比传统的块多普勒估计方法估计精度更高。而基于信道自适应多普勒估计方法进行重采样的差分水声OFDM通信系统相比采用块多普勒估计方法以及不进行重采样操作的系统,系统性能也有显著提升。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种水声OFDM通信非均匀多普勒频偏抑制方法,发送机在发送信号的首部和尾部分别添加LFM同步信号,该方法应用于接收机,其特征在于,包括:
S1.将LFM同步信号长度Tl、多普勒因子最大可能值αmax和原始发送信号长度Tt作为先验信息,对基带接收信号的首部和尾部同步信号分别进行截取;
S2.截取到的首尾两个同步信号进行互相关运算,得到相关峰位置,作为基带接收数据信号长度相对于原始发送信号长度的变化量;
S6.对重采样后的基带接收信号去除掉首部和尾部长度为Tl的信号,得到频偏补偿后的接收数据信号。
5.一种水声OFDM通信非均匀多普勒频偏抑制系统,其特征在于,包括:计算机可读存储介质和处理器;
所述计算机可读存储介质用于存储可执行指令;
所述处理器用于读取所述计算机可读存储介质中存储的可执行指令,执行权利要求1至4任一项所述的水声OFDM通信非均匀多普勒频偏抑制方法。
6.一种计算机可读存储介质,其特征在于,包括存储的计算机程序;
所述计算机程序被处理器执行时,控制所述计算机可读存储介质所在设备执行权利要求1至4任一项所述的水声OFDM通信非均匀多普勒频偏抑制方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210329421.3A CN114710386B (zh) | 2022-03-30 | 2022-03-30 | 一种水声ofdm通信非均匀多普勒频偏抑制方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210329421.3A CN114710386B (zh) | 2022-03-30 | 2022-03-30 | 一种水声ofdm通信非均匀多普勒频偏抑制方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114710386A true CN114710386A (zh) | 2022-07-05 |
CN114710386B CN114710386B (zh) | 2024-03-19 |
Family
ID=82170973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210329421.3A Active CN114710386B (zh) | 2022-03-30 | 2022-03-30 | 一种水声ofdm通信非均匀多普勒频偏抑制方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114710386B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115987741A (zh) * | 2023-03-20 | 2023-04-18 | 四川创智联恒科技有限公司 | 一种基于ofdm系统的多普勒频偏补偿方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106027116A (zh) * | 2016-07-07 | 2016-10-12 | 哈尔滨工程大学 | 一种基于chirp信号的移动水声通信多普勒系数估计方法 |
CN110247867A (zh) * | 2019-05-16 | 2019-09-17 | 中国科学院声学研究所 | 水声多普勒估计方法及装置、水声通信方法及系统 |
US20200366382A1 (en) * | 2017-11-22 | 2020-11-19 | South China University Of Technology | Sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method |
-
2022
- 2022-03-30 CN CN202210329421.3A patent/CN114710386B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106027116A (zh) * | 2016-07-07 | 2016-10-12 | 哈尔滨工程大学 | 一种基于chirp信号的移动水声通信多普勒系数估计方法 |
US20200366382A1 (en) * | 2017-11-22 | 2020-11-19 | South China University Of Technology | Sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method |
CN110247867A (zh) * | 2019-05-16 | 2019-09-17 | 中国科学院声学研究所 | 水声多普勒估计方法及装置、水声通信方法及系统 |
Non-Patent Citations (3)
Title |
---|
LU MA 等: ""Low-Complexity Doppler Compensation Algorithm for Underwater Acoustic OFDM Systems With Nonuniform Doppler Shifts"", 《IEEE COMMUNICATIONS LETTERS》, 28 May 2020 (2020-05-28), pages 2051 - 2054, XP011807707, DOI: 10.1109/LCOMM.2020.2998293 * |
张行 等: ""MIMO-OFDM 移动水声通信系统中 多普勒扩展估计方案"", 《东南大学学报( 自然科学版)》, 31 May 2017 (2017-05-31), pages 215 - 219 * |
杨思琪: ""时变信道下水声单载波频域均衡与同步技术研究"", 《中国优秀硕士学位论文全文数据库》, 15 May 2021 (2021-05-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115987741A (zh) * | 2023-03-20 | 2023-04-18 | 四川创智联恒科技有限公司 | 一种基于ofdm系统的多普勒频偏补偿方法 |
CN115987741B (zh) * | 2023-03-20 | 2023-05-30 | 四川创智联恒科技有限公司 | 一种基于ofdm系统的多普勒频偏补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114710386B (zh) | 2024-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100594084B1 (ko) | 직교 주파수 분할 다중 수신기의 채널 추정 방법 및 채널추정기 | |
CN107547143B (zh) | 一种已知子载波频率的ofdm-mfsk水声通信宽带多普勒估计与补偿方法 | |
CN1988526B (zh) | 一种多输入多输出、正交频分复用无线系统的同步方法 | |
US20090323843A1 (en) | Method and apparatus for estimating symbol timing offset in wireless communication system | |
CN110636024A (zh) | 一种基于索引调制的5g波形系统同步方法 | |
CN110086560B (zh) | 适用于短突发通信时频联合同步的同步帧结构及同步方法 | |
KR20100060456A (ko) | 무선 통신 시스템에서 주파수 추정 방법 및 장치 | |
CN107454031B (zh) | 一种基于分组信噪比置信度的ofdm-mfsk水声通信技术 | |
CN111935050B (zh) | 一种基于相位搜索的单载波频域均衡水声通信系统残余相偏修正方法 | |
CN104836770B (zh) | 一种基于相关平均与加窗的定时估计方法 | |
JP2005079911A (ja) | Ofdm伝送方式における受信装置 | |
CN102377726B (zh) | Ofdm系统的定时同步方法 | |
EP2194666A1 (en) | Radio communication device, and reception quality estimating method | |
CN108768914B (zh) | 联合正交与非正交的高效频分复用传输方法及传输系统 | |
CN114710386A (zh) | 一种水声ofdm通信非均匀多普勒频偏抑制方法和系统 | |
WO2007055469A1 (en) | Method for generating preamble sequence using pn sequence, and method for time synchronization and frequency offset estimation using pn sequence | |
JP2007201523A (ja) | 受信信号等化装置および方法 | |
CN113132276B (zh) | 一种用于远程无人机图传的单载波频域均衡传输系统 | |
CN110071888B (zh) | 一种高速数传中的快速时间同步方法 | |
CN107276654B (zh) | 信号处理方法和系统 | |
CN111884978B (zh) | 一种基于ofdm抗脉冲噪声的符号同步方法 | |
CN107276940B (zh) | 定时同步方法、装置和系统 | |
Nguyen et al. | A multi-stage method for time synchronization in acoustic underwater communications | |
KR20010081771A (ko) | 직교 주파수 분할 다중 신호의 송/수신 장치 및 방법 | |
CN107277913B (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 |