CN107911133B - 一种移动水声通信的多普勒因子估计及补偿方法 - Google Patents

一种移动水声通信的多普勒因子估计及补偿方法 Download PDF

Info

Publication number
CN107911133B
CN107911133B CN201711148334.3A CN201711148334A CN107911133B CN 107911133 B CN107911133 B CN 107911133B CN 201711148334 A CN201711148334 A CN 201711148334A CN 107911133 B CN107911133 B CN 107911133B
Authority
CN
China
Prior art keywords
doppler
signal
compensation
estimation
underwater sound
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
Application number
CN201711148334.3A
Other languages
English (en)
Other versions
CN107911133A (zh
Inventor
王德清
徐景鑫
胡晓毅
解永军
林志达
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xiamen Tangpu Information Technology Co ltd
Original Assignee
Xiamen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xiamen University filed Critical Xiamen University
Priority to CN201711148334.3A priority Critical patent/CN107911133B/zh
Publication of CN107911133A publication Critical patent/CN107911133A/zh
Application granted granted Critical
Publication of CN107911133B publication Critical patent/CN107911133B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
    • H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
    • H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
    • H04B1/1036—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal with automatic suppression of narrow band noise or interference, e.g. by using tuneable notch filters
    • 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/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2649—Demodulators
    • H04L27/265—Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
    • 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/2662—Symbol synchronisation
    • 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/2666—Acquisition of further OFDM parameters, e.g. bandwidth, subcarrier spacing, or guard interval length
    • 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)
  • Physics & Mathematics (AREA)
  • Discrete Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Radio Relay Systems (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

一种移动水声通信的多普勒因子估计及补偿方法，涉及水声通信。1)大尺度多普勒估计及补偿；2)残余多普勒估计及补偿；3)多普勒相位旋转补偿。为了准确且高效地估计移动水声通信环境下的多普勒因子，从而克服与陆上无线电信道相比更为明显的多普勒效应，消除其对水声OFDM通信系统带来的不良影响，需要一种兼顾准确性和计算复杂度的移动水声通信的多普勒因子估计及补偿方法。在移动水声通信中，尤其应用于移动水声OFDM系统的多普勒因子估计和补偿方法，估计在频域进行，针对水声信道的特点，在传统多普勒估计及补偿的基础之上，更加适用于快速变化的移动水声信道。其估计精度高，同时又适当地降低计算复杂度，实用性优良。

Description

一种移动水声通信的多普勒因子估计及补偿方法
技术领域
本发明涉及水声通信领域，尤其是涉及一种移动水声通信的多普勒因子估计及补偿方法。
背景技术
水下无人航行器(Unmanned Underwater Vehicle,UUV)作为海洋探测和水声通信的重要 载体，已成为各国海洋工程领域的研究热点。UUV具有智能型、隐蔽性、机动性、经济性等 特点，并且重量轻、尺寸小，使用方便，具有广阔的经济和工业价值[1]。利用UUV进行高速 运动过程中的移动水声通信，已成为涵盖海洋技术与信息技术的世界各国急需的高新技术之 一。
不同于陆地上的无线电通信，声波在水下的传播速度远小于电磁波的传播速度，故与陆 空无线通信比较而言，水声信道中存在更为明显的多普勒效应。由于水声信道的可用带宽窄， 水声通信本质上是一个宽带通信系统，因此，水声信道中的多普勒效应将会引起接收信号在 频域上的频率偏移和频谱扩展，在时域上则表现为接收时域波形的压缩或扩展。对于水下无 人航行器，其运动速度一般为1.5～15m/s,此时的多普勒扩展因子可达10-3量级，在无线电通 信系统中唯一能够与之相比的是低轨卫星系统过顶时产生的多普勒效应。为了减小多普勒效 应对水声信号接收和处理的不良影响，有必要对接收到的水声信号进行有效的多普勒扩展因 子估计和补偿，从而能够恢复出原始信号波形。
正交频分复用技术(Orthogonal Frequency Division Multiplexing,OFDM)具有良好的抗频 率选择性衰落和高的带宽利用率，可以用于多径效应严重的水下高速声通信。然而，OFDM 水声通信系统的接收性能对符号同步偏差和频率偏差敏感。相较于陆上无线电信道，在水声 通信中，由于通常工作频率较低，且水声声速(约1500m/s)远低于无线电传输速率，当水 下运动目标之间即使存在很小的相对运动，都会引起较大的多普勒扩展，将在OFDM水声通 信系统接收端引起严重的子载波间干扰(Inter-Carrier Interference,ICI)，导致解调器性能恶 化，甚至不能正确解调。因此，OFDM移动水声通信技术的关键就是如何有效地降低多普勒 扩展引起的子载波间干扰。
当前，水声通信中主要采取以下三类方法来降低多普勒扩展引起的子载波间干扰[2]：
1.均衡的方法。该类方法主要使用时域或频域均衡，通过采用一定的均衡算法来对抗多 普勒扩展，降低子载波间干扰。这类方法对于符号间干扰比较有效，但是对于子载波间干扰， 仅适合小多普勒频移下的高速水声通信或非相干水声通信，对于高速OFDM水声通信，计算 复杂度较高。
2.基于发送端处理的方法。这类方法通过设计发送端数据结构、改变子载波映射方式等 消除子载波间干扰，采用这些方法，通信系统的频带利用率比一般的OFDM要低，除导频及 训练序列外，需要占用额外的带宽。
3.基于估计与补偿的方法。这类方法通过降低多普勒扩展对信号的影响来消除子载波间 干扰，也是目前应用最为广泛的方法。它的思想是：首先估计多普勒因子，然后根据此多普 勒因子对接收信号进行重采样，最后再进行解调。相较于前面两类方法，这类方法效果较好。
公开的移动水声通信多普勒估计方法中多为基于时域相关运算的经典多普勒因子估计算 法，主要包括块多普勒估计法、自循环移位序列估计法、模糊度函数估计法等。其基本思想 为利用首尾的同步信号或是循环前缀进行整体的多普勒因子估计，其中块多普勒估计法对多 普勒因子的估计只能在一次数据包发送完成后进行，不具备实时操作性，而自循环移位序列 估计法和模糊度函数估计法可以直接通过接收端的相关运算进行估计，但对于水声高噪声、 强多径的环境下，训练序列的自相关性往往受到限制，基于时域相关运算的多普勒估计方法 在水声信号处理领域对多普勒因子的估计结果往往不甚理想。
在频域，多普勒因子的估计则可以通过多普勒偏移的检测来完成。经典的基于FFT(Fast Fourier Transformation)的频偏检测算法，由于受到加窗信号FFT变换后的频谱泄漏效应的限 制，导致估计算法的分辨率较低，往往难以满足实际应用的需求。而若通过增加有效数据长 度以提高分辨率，则增大了计算量，降低了有效数据通信速率。为此研究者提出了一系列在 不增加有效信号长度前提下提高频率分辨率的方法。目前已有多种测频算法成功应用于水声 信号测频工程中，诸如复相关、ZOOM-FFT、插值FFT等多普勒估计方法[3][4]，但这些频域 估计方法在计算复杂度和估计精度间无法兼顾。文献[5]提出了一种“基于初始角频率可调FFT (FFTω)的多普勒估计算法”，其特点在于估计精度高，但同时具有较高计算复杂度。
实际应用中用探测信号的估计结果来补偿后续的OFDM符号，有可能会存在较大的偏差。 此外，由于重采样误差、海浪在通信时间内的起伏、收发两端相对运动速度变化等原因，在 大尺度多普勒估计后，仍存在较大的残余多普勒扩展因子。基于空子载波的OFDM频偏估计 算法[6]，是一种通过对空子载波上的能量进行计算，以估计载波频率偏移的大小的方法,但其 以一定频率步长进行预设频带内的估计，其估计精度存在一定偏差。
目前公开的移动水声通信多普勒因子估计方法多可归纳为基于首尾的同步信号或是自循 环移位序列的时域估计方法和基于单频信号的频域估计方法，对接收信号进行整体的多普勒 因子估计。如公开的专利文献中，公开号为CN105282082A的一种基于拷贝相关与空子载波 结合的多普勒估计方法提出在每帧信号首尾插入线性调频信号，进行整体的多普勒因子估计， 该方法仅根据数据帧的整体多普勒扩展作出估计，难以实现对多普勒因子的准确估计；公开 号为CN104901718A的基于直接序列扩频信号载波频率测量的多普勒估计方法，基于直接序 列扩频信号对解扩后的信号测量频率值，对比估计频率值和原载波频率值估计多普勒因子； 公开号为CN103618686A的水声OFDM多普勒因子精确估计方法，设计一种加入带循环前缀 的前同步码和单频信号的OFDM帧，对接收信号进行多普勒估计；公开号为CN102916922A 的文献《水声OFDM自适应搜索多普勒补偿方法》采用CW(ContinuesWave)单频信号作 为训练序列进行多普勒频偏因子粗测，利用FFT对多普勒频偏进行补偿。以上方案均是通过 同步信号或是自循环序列对接收信号进行整体的多普勒因子估计，在具有强时变性的水声信 道环境下，仅以同步信号或自循环序列所在时刻的多普勒因子估计值来补偿后续OFDM数据 符号的多普勒扩展，其估计准确性难以保证。以上方法不能精确地估计出每个OFDM符号的 多普勒因子，从而不能达到对信号进行准确解调的效果。
参考文献
[1]普湛清,王巍,张扬帆等.UUV平台OFDM水声通信时变多普勒跟踪与补偿算法[J]. 仪器仪表学报,2017,38(7):1634-1634.
[2]王彪，支志福，朱志宇等.水声高速OFDM通信的多普勒扩展处理方法[P].中国:201310416303.7,2013.09.12.
[3]Susaki,H.,"Method of high-resolution frequency measurement forpulse-Doppler sonar," Underwater Technology,2002.Proceedings of the2002International Symposium on,vol.,no., pp.39,44,2002.
[4]V.K.Jain,W.L.Collins,and D.C.Davis,“High-accuracy analogmeasurements via interpolated fft,”Instrumentation and Measurement,IEEETransactions on,vol.28,no.2,pp.113 –122,June 1979.
[5]马文翰.移动水声OFDM通信系统多普勒效应估计补偿算法研究[D].厦门大学,2013.
[6]吉磊，李玉柏.OFDM系统的频偏估计算法[J].计算机应用研究，2011，28(4):1480-1485.
发明内容
本发明的目的是提供一种移动水声通信的多普勒因子估计及补偿方法。
本发明包括以下步骤：
1)大尺度多普勒估计及补偿；
在步骤1)中，所述大尺度多普勒估计及补偿，为对接收机接收到的信号进行第一次多 普勒扩展处理，具体包括如下步骤：
(1)接收机接收信号首先进入带通滤波器BPF，滤除带外噪声；
(2)对带通滤波器BPF的输出信号进行同步检测，实现起始信号的定位；
(3)根据发送信号中插入的探测信号，基于“初始角频率可调的FFT变换(FFTω)”的多普勒扩展因子估计算法，结合抛物线拟合算法降低原有多普勒扩展因子估计算法的计算 复杂度，实现对多普勒扩展因子a的粗估计，得到估计值所述探测信号包括三组频率， 三组频率分别为f1,f2,f3的单频信号序列；
(4)以多相滤波器重采样的方法，对带通滤波器BPF输出信号以进行重采样，得到 重采样后信号实现对接收信号的第一次多普勒补偿。
2)残余多普勒估计及补偿；
在步骤2)中，所述残余多普勒估计及补偿的具体方法可为：经过大尺度多普勒估计及补 偿后，系统仍然受残余的多普勒扩展因子的影响，利用空子载波进行残余多普勒估计，相对 于传统的残余频偏估计，其估计精度更高，此为对接收信号进行第二次多普勒扩展处理，具 体步骤如下：
(1)对步骤1)第(4)部分产生的重采样后信号r′(t)，以采样率fs进行采样，采样后得 到接收信号r′(n)；
(2)由已知空子载波个数K、空子载波位置信息k＝[k1,k2,…,kK]以及m个不同残余多 普勒扩展因子组成的集合a′m,m＝1,2,…,M，对接收信号r′(n)进行M次重采样，并进行N点 FFT运算，得到各空子载波位置的幅度：
Rm(i),i＝[0,N-1]；
(3)提取出个空子载波位置处对应的FFT运算结果Ym(j)，进而得该OFDM符号中空子 载波处的能量之和J(m)；
(4)若对空子载波处的数据进行了正确的残余多普勒补偿，则理论上空子载波处不会出 现因子载波干扰引起的从相邻子载波泄露的能量，由此可知步骤(3)得到的J(m)中最小值 所对应补偿的残余频偏最接近真实的残余多普勒，得从而得出估计的残余 多普勒
(5)通过已求的残余多普勒对重采样后的信号r′(n)进行残余多普勒补偿，得到补偿 后信息序列
3)多普勒相位旋转补偿。
在步骤3)中，所述多普勒相位旋转补偿的具体方法可为：当大尺度多普勒补偿后，若各 个符号发送数据不同，或存在符号定时偏差或载波频率偏差，进行残余多普勒补偿后，就会 导致不同符号之间存在相位偏差，即相位旋转，且每个符号的的每个子载波都存在不同的相 位旋转值；当OFDM系统为时域差分系统或块状放置导频的相干系统时，由残余多普勒带来 的相位旋转将会严重影响系统性能；根据所述残余多普勒估计及补偿得到的相位补偿所需偏 差点数对接收信号进行多普勒相位旋转补偿，此为对接收信号进行第三次多普勒扩展处理；
由于多普勒扩展因子是实时发生变化的，对于每一个符号的残余多普勒估计结果也不尽 相同，若M个符号的估计结果分别为假设各个符号的估计结果足够精确，且 在一个符号内残余多普勒扩展因子几乎不变，则此时第m个符号上第k个子载波所需补偿的 相位应为：
综上所述，为了准确且高效地估计移动水声通信环境下的多普勒因子，从而克服与陆上 无线电信道相比更为明显的多普勒效应，消除其对水声OFDM通信系统带来的不良影响，需 要一种兼顾准确性和计算复杂度的移动水声通信的多普勒因子估计及补偿方法。
本发明的有益效果体现在：
(1)构建大尺度多普勒估计及补偿、残余多普勒估计及补偿和多普勒相位旋转补偿的 “三步走”多普勒估计及补偿方案，与现有的“先粗后细”的“两步走”多普勒估计补偿方案相比， 由于增加针对残余多普勒导致的相位旋转补偿，使得多普勒因子估计精度更高，估计更准确。
(2)与现有的多普勒估计及补偿方法相比，首次提出了一种基于FFTω变换的抛物线拟 合多普勒估计算法。与原有FFTω相比，该算法通过抛物线拟合的方式便可得到待测频率的精 确位置，而仅需要进行较少次的FFTω变换，大大减少运算量。与传统的基于FFT变换的抛 物线拟合多普勒估计方法相比，由于FFTω可通过调整变换初始角频率，使得其谱线逐步逼近 待估计的单频信号的谱线位置，最小化频谱泄漏，因此具有更高的估计精度。
(3)通过空子载波上的能量比较，与背景技术中提及的残余频偏估计及补偿相比，本发 明提出的残余多普勒估计及补偿方法能够更加有效地减小残余多普勒带来的ICI，估计精度 更高，补偿效果更好。
(4)本发明提出的多普勒相位旋转补偿方法，根据所述的残余多普勒估计及补偿得到的 相位补偿所需偏差点数，对接收信号每个符号的每个子载波上存在的不同相位旋转值进行多 普勒相位旋转补偿，能够有效地解决时域差分或块状放置导频的OFDM相干系统中残余多普 勒导致的每个符号上的每个子载波存在不同的相位旋转值的问题。
(5)在移动水声通信中，尤其应用于移动水声OFDM系统的多普勒因子估计和补偿方 法，估计在频域进行，针对水声信道的特点，在传统多普勒估计及补偿的基础之上，更加适 用于快速变化的移动水声信道。其估计精度高，同时又适当地降低计算复杂度，实用性优良。
附图说明
图1为应用的OFDM移动水声通信系统框图。
图2为探测信号发送帧格式。
图3为抛物线拟合算法示意图；在图3中，曲线a为原始频谱，b为抛物线拟合。
图4为重采样插值与抽取方法示意图。
图5为多相滤波器结构示意图。
图6为计算机仿真中大尺度多普勒补偿后的数据星座图。
图7为计算机仿真中残余频偏补偿后的数据星座图。
图8为计算机仿真中残余多普勒补偿后的数据星座图。
图9为计算机仿真中多普勒相位旋转补偿后的数据星座图。
图10为海洋试验的行船示意图。
图11为海试试验中大尺度多普勒补偿后的数据星座图。
图12为海试试验中残余频偏补偿后的数据星座图。
图13为海试试验中残余多普勒补偿后的数据星座图。
图14为海试试验中多普勒相位旋转补偿后的数据星座图。
具体实施方式
下面将结合本发明的附图，对本发明的具体实施方式进行清楚、完整地描述：
在本实施例中，采用基带OFDM时域差分系统作为所构建的移动水声通信系统，系统框 图如图1所示：
发送端信源比特流经过串并转换、4DPSK(Differential Phase Shift Keying)符号映射以 及插入空载波(用于接收端的残余多普勒估计)后，经过OFDM调制和串并转换处理，之后 将发射信号插入探测信号和同步chirp信号，用于接收端的大尺度多普勒估计和信号同步。探 测信号为三组频率分别为f1,f2,f3的单频信号序列s1(t),s2(t),s3(t)，该频率均为FFT频率分辨 率的整数倍。如图2所示，在发送帧格式中插入串行和并行两种发送方式的单频探测信号， 以保证接收端单频信号的可靠信噪比。
接收端获取接收信号，根据接收信号与本地拷贝chirp信号作相关运算，由chirp信号强 自相关特性得到相关峰位置，从而确定接收数据信号的起始位置。如图1所示，对多普勒因 子的估计及补偿以“三步走”方案实现：(1)大尺度多普勒估计及补偿、(2)残余多普勒估 计及补偿、(3)多普勒相位旋转补偿。以下对本发明针对多普勒估计及补偿的三步骤方案进 行详细的描述：
1)大尺度多普勒估计及补偿
步骤S1-1：接收机接收信号进入带通滤波器BPF，滤除带外噪声；
步骤S1-2：对带通滤波器BPF的输出信号以本地拷贝chirp信号进行，实现同步检测， 确定起始数据信号的位置；
步骤S1-3：根据发送信号中插入的探测信号(三组频率分别为f1,f2,f3的单频信号序列 s1(t),s2(t),s3(t))，基于FFTω的多普勒扩展因子估计算法，结合抛物线拟合算法降低原有多 普勒扩展因子估计算法的计算复杂度，实现对多普勒扩展因子a的粗估计，得到估计值图3为抛物线拟合算法示意图，图中fd为FFT的采样间隔，(f,Y)为待测频点所在坐标，通 过FFTω变换，求得(f1,Y1)，(f2,Y2)，(f3,Y3)三个采样点坐标，利用抛物线(二次函数)来拟合sinc 函数的顶点(f,Y)的坐标。设该抛物线的表达式为Y＝af2+bf+c，则通过拉格朗日插值多项 式求得表达式各系数a、b、c，进而求得抛物线极值点对应的频率：
Y＝p(f-f2)(f-f3)+q(f-f1)(f-f3)+r(f-f1)(f-f2)
p＝Y1/(f1-f2)(f1-f3)
q＝Y2/(f2-f1)(f2-f3)
r＝Y3/(f3-f1)(f3-f2)
a＝p+q+r
b＝-p(f2+f3)-q(f1+f3)-r(f1+f2)
c＝pf2f3+qf1f3+rf1f2
以此方法分别对每个单频信号序列s1(t),s2(t),s3(t)进行频率估计，估计结果为得到多普勒扩展因子的估计值为：
考虑到个别频点可能由于信道衰落存在较大的误差，对通过求中值得到多普勒扩展因 子的估计值
步骤S1-4：以多相滤波器重采样的方式，对带通滤波器BPF输出信号以进行重采样， 得到重采样后信号实现对接收信号的第一次多普勒补偿。
图4为重采样插值与抽取方法示意图，即先对接收信号进行L倍的插值过采样，经过低 通滤波器之后再进行M倍的抽取降采样，以实现等效L/M因子的采样率转换，从而实现多 普勒因子的补偿。
由于插值和抽取过程均需要经过低通滤波器，而在实际应用中，插值因子L和抽取因子 M通常都很大，因而造成了采样率转换过程中的低通滤波器也具有较高的阶次，因此变采样 的运算量很大，为了降低运算量,同时保证变采样后信号较高的保真度,本发明采用如图5的多 相插值滤波器进行变采样。L/M倍抽样率转换器的输入和输出关系如式所示：
其中h(n)是一个FIR滤波器的单位抽样响应，长度为N，将其分为L组子滤波器，每一 组子滤波器长度为K，<nM>L表示对L取模运算，表示向下取整。
2)残余多普勒估计及补偿
实际应用中用探测信号的估计结果来补偿后续的OFDM符号，由于水声信道的时变性， 有可能会存在较大的偏差，加之在对OFDM系统进行大尺度多普勒效应补偿之后，由于大 尺度多普勒估计误差与重采样误差等原因，系统仍然受残余的多普勒扩展因子的影响。利用 空子载波进行残余多普勒估计，相对于残余频偏估计，其估计精度更高，本发明采用如下方 法加以实现：
步骤S2-1：对步骤S1-4产生的重采样后信号r′(t)，以采样率fs进行采样，采样后得到接 收信号r′(n)，之后将其与步骤S1-2中提及的本地拷贝chirp信号进行再次相关，通过相关峰 位置及相应数据长度得到OFDM数据帧的起始位置；
步骤S2-2：设已知空子载波个数为K，空子载波位置信息k＝[k1,k2,…,kK]且将m个不同 残余多普勒扩展因子组成的集合表示为a′m,m＝1,2,…,M，对接收信号r′(n)进行M次与步骤 S1-4相同的多相滤波重采样，并进行N点FFT运算，得到各空子载波位置的幅度：
其中fc为载波频率。
步骤S2-3：提取出个空子载波位置处对应的FFT运算结果：
Ym(j)＝Rm(kj),j＝1,2,…,K
进而切得该OFDM符号中空子载波处的能量之和：
J(m)＝Ym(j)×Ym(j)H
步骤S2-4：若对空子载波处的数据进行了正确的残余多普勒补偿，则理论上空子载波处 不会出现因子载波干扰引起的从相邻子载波泄露的能量，由此可知步骤S2-3得到的J(m)中 最小值所对应补偿的残余频偏最接近真实的残余多普勒因子
步骤S2-5：通过已求的残余多普勒对步骤S1-4重采样后的信号r′(n)进行残余多普勒 补偿，即以相同的多相滤波器重采样的方式，以再次重采样，得到补偿后信息序列
至此，完成对接收信号残余多普勒效应的消除，实现对接收信号的第二次多普勒补偿。
3)多普勒相位旋转补偿
当大尺度多普勒补偿后，若各个符号发送数据不同，或存在符号定时偏差或载波频率偏 差，进行残余多普勒补偿后，就会导致不同符号之间存在相位偏差，即相位旋转，且每个符 号的的每个子载波都存在不同的相位旋转值。当OFDM系统为时域差分系统或块状放置导频 的相干系统时，有残余多普勒带来的相位旋转将会严重影响系统性能。为此，本发明对接收 信号根据所述残余多普勒估计及补偿得到的相位补偿所需偏差点数对接收信号进行多普勒相 位旋转补偿，此为对接收信号进行第三次多普勒扩展处理。
设发送信号为
其中，m表示第m个OFDM符号，fi＝fc+i/T表示第i个子载波频率(其中 i＝0,1,…,N-1，fc为载波频率)，第m个符号的第i个子载波上的源数据表示为dmi，T表示 符号长度，Tcp表示循环前缀长度。
设水声信道各径具有相同多普勒扩展因子，经过该信道后的接收信号以采样间隔为Ts对 接收信号进行采样后，得到接收序列：
其中，为噪声序列。
设多普勒扩展因子估计结果为且当最大多径时延小于循环前缀长度时ISI(符号间干扰) 约为0，对接收信号的第m个符号进行重采样，则第m个符号可表示为：
其中为多径对第i个子载波带来的相位偏差。
对重采样后的数据进行FFT解调，可得第m个符号上的第k个子载波的解调结果为：
其中，表示噪声于第k个子载波上的干扰。且由于估计误差c表示由频率同步 偏差引起的能量损失，I表示子载波间干扰(ICI)。
由以上公式，对接收信号进行相位补偿：
其中，表示相位补偿后第k个子载波上的残余相位旋转，此时，当忽略噪 声带来的影响，即则通过信道估计或差分解调后，ymk≈dmk。
由于多普勒扩展因子是实时发生变化的，对于每一个符号的残余多普勒估计结果也不尽 相同，若M个符号的估计结果分别为假设各个符号的估计结果足够精确，且 在一个符号内残余多普勒扩展因子几乎不变，则此时第m个符号上第k个子载波所需补偿的 相位应为：
至此，完成对接收信号的多普勒相位旋转补偿，实现对接收信号的第三次多普勒补偿。
为了说明本发明的有效性，通过以下计算机仿真实验和海洋试验进一步说明：
计算机仿真中，人为设定大尺度多普勒估计存在估计误差(即残余多普勒因子)为a′＝1×10-4，图6～9分别为大尺度多普勒补偿后、残余频偏补偿后、残余多普勒补偿后和多普勒相位旋转补偿后的QPSK(正交相移键控，Quadrature Phase Shift Keying)数据星座图的仿 真结果。由此可知，当估计足够精确时，经过3个步骤的补偿后，可消除多普勒对系统的影 响，以实现符号的正确解调。且对比图6～8三图，可知第二步采用残余频偏估计与补偿能有 效减小残余多普勒带来的ICI，但相对于采用残余多普勒估计与补偿，仍然存在一定ICI。
在海洋试验中，行船示意图如图10所示，接收机静止，发射机以5节的速度做直线运动。 由于不同位置下行船方向与收发两端连线的夹角不同，多普勒因子值(a＝vcosθ/c，vcosθ为 相对运动速度，c为声速)也不同，从而为实验提供实时变化的多普勒效应。图11～14分别 为大尺度多普勒补偿后、残余频偏补偿后、残余多普勒补偿后和多普勒相位旋转补偿后的 QPSK数据星座图的海试试验结果。由此可知，在实际应用中，本发明阐述的方法仍可以有 效减小多普勒效应带来的ICI，实现正确的符号解调。

Claims (2)

1.一种移动水声通信的多普勒因子估计及补偿方法，其特征在于包括以下步骤：
1)大尺度多普勒估计及补偿：对接收机接收到的信号进行第一次多普勒扩展处理，具体包括如下步骤：
(1)接收机接收信号首先进入带通滤波器BPF，滤除带外噪声；
(2)对带通滤波器BPF的输出信号进行同步检测，实现起始信号的定位；
(3)根据发送信号中插入的探测信号，基于“初始角频率可调的FFT变换(FFTω)”的多普勒扩展因子估计算法，结合抛物线拟合算法降低原有多普勒扩展因子估计算法的计算复杂度，实现对多普勒扩展因子a的粗估计，得到估计值所述探测信号包括三组频率分别为f1,f2,f3的单频信号序列；
(4)以多相滤波器重采样的方法，对带通滤波器BPF输出信号以进行重采样，得到重采样后信号实现对接收信号的第一次多普勒补偿；
2)残余多普勒估计及补偿：对接收信号进行第二次多普勒扩展处理，具体步骤如下：
(1)对步骤1)第(4)部分产生的重采样后信号r′(t)，以采样率fs进行采样，采样后得到再次重采样后的信号r′(n)；
(2)由已知空子载波个数K、空子载波位置信息k＝[k1,k2,…,kK]以及m个不同残余多普勒扩展因子组成的集合a′m,m＝1,2,…,M，对接收信号r′(n)进行M次重采样，并进行N点FFT运算，得到各空子载波位置的幅度：
Rm(i),i＝[0,N-1]
(3)提取出各空子载波位置处对应的FFT运算结果Ym(j)：
Ym(j)＝Rm(kj),j＝1,2,…,K
进而得到OFDM符号中空子载波处的能量之和J(m)：
J(m)＝Ym(j)×Ym(j)H
(4)步骤2)第(3)部分所得的J(m)中最小值所对应补偿的残余频偏最接近真实的残余多普勒，得从而得出估计的残余多普勒
(5)通过已求的残余多普勒对再次重采样后的信号r′(n)进行残余多普勒补偿，得到补偿后信息序列
3)多普勒相位旋转补偿。
2.如权利要求1所述的一种移动水声通信的多普勒因子估计及补偿方法，其特征在于在步骤3)中，所述多普勒相位旋转补偿的具体方法为：根据所述残余多普勒估计及补偿得到的相位补偿所需偏差点数对接收信号进行多普勒相位旋转补偿，此为对接收信号进行第三次多普勒扩展处理。
CN201711148334.3A 2017-11-17 2017-11-17 一种移动水声通信的多普勒因子估计及补偿方法 Active CN107911133B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711148334.3A CN107911133B (zh) 2017-11-17 2017-11-17 一种移动水声通信的多普勒因子估计及补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711148334.3A CN107911133B (zh) 2017-11-17 2017-11-17 一种移动水声通信的多普勒因子估计及补偿方法

Publications (2)

Publication Number Publication Date
CN107911133A CN107911133A (zh) 2018-04-13
CN107911133B true CN107911133B (zh) 2019-08-23

Family

ID=61846208

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711148334.3A Active CN107911133B (zh) 2017-11-17 2017-11-17 一种移动水声通信的多普勒因子估计及补偿方法

Country Status (1)

Country Link
CN (1) CN107911133B (zh)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108924073B (zh) * 2018-06-28 2019-08-20 中国科学院声学研究所 一种基于伪随机序列同步的快速自适应多普勒估计方法
CN109104387B (zh) * 2018-08-23 2021-01-01 上海交通大学 一种水声OFDM通信中基于Dual-Zadoff-Chu序列的信号同步方法
CN109547372B (zh) * 2018-10-15 2021-06-25 中国人民解放军战略支援部队信息工程大学 正交频分复用水声通信中时变宽带多普勒因子估计方法及装置
CN109143231B (zh) * 2018-10-29 2023-03-28 河海大学 基于循环对消的数字电视无源双基地雷达目标检测方法
CN109743118B (zh) * 2018-12-24 2021-05-11 三亚哈尔滨工程大学南海创新发展基地 一种时变双扩信道条件下的高频谱效率的ofdm水声通信方法
CN109768812B (zh) * 2019-01-30 2020-12-22 杭州电子科技大学 一种基于混沌调频的水声通信多普勒估计与同步方法
CN109714112B (zh) * 2019-02-28 2021-02-26 厦门大学 一种利用移动平台集群的水声通信方法及系统
CN110113279B (zh) * 2019-05-05 2021-09-28 哈尔滨工程大学 一种移动跳频水声通信多普勒因子估计方法
CN110545247B (zh) * 2019-07-31 2021-09-21 成都天奥集团有限公司 低轨卫星多载波通信系统下行链路载波频偏估计与补偿方法
CN110798422B (zh) * 2019-11-25 2022-01-07 成都天奥集团有限公司 低轨卫星多载波通信系统下行链路采样频偏估计与补偿方法
CN111131098B (zh) * 2019-12-31 2022-08-19 哈尔滨哈船海洋信息技术有限公司 基于水声ofdm信号自相关输出峰插值的多普勒因子估计方法
CN112003810B (zh) * 2020-10-28 2021-01-29 鹏城实验室 一种基于ofdm的水声通信同步方法、智能终端及存储介质
CN112152951B (zh) * 2020-11-25 2021-02-23 鹏城实验室 水声通信探测方法、装置、设备及存储介质
CN112187697B (zh) * 2020-11-25 2021-02-23 鹏城实验室 水声通信探测信号生成方法、装置、设备及存储介质
CN112737996B (zh) * 2020-12-25 2022-03-25 中国科学院声学研究所 水声通信时变定时偏移迭代估计和补偿方法及装置
CN113179234B (zh) * 2021-04-26 2022-01-18 哈尔滨工程大学 一种基于分块多普勒补偿的高动态长信号累积方法
CN113259291B (zh) * 2021-05-26 2022-01-18 哈尔滨工程大学 利用水声连续信号动态多普勒跟踪实现的相位补偿方法
CN113939009A (zh) * 2021-10-13 2022-01-14 国网信息通信产业集团有限公司 同步估计方法、装置及接收机和存储介质
CN114124245B (zh) * 2021-11-05 2024-05-03 北京长城电子装备有限责任公司 一种时变信道下的水声高速通信方法
CN114389717B (zh) * 2021-12-29 2022-11-15 北京理工大学 基于正交频分复用ofdm的异步多用户水声通信方法和装置
CN116418636B (zh) * 2023-06-08 2023-08-18 芯迈微半导体(上海)有限公司 一种无线通信的增强的多普勒扩展估计方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103354538A (zh) * 2013-07-20 2013-10-16 华南理工大学 一种对水声通信中的接收信号进行多普勒补偿的方法
CN103491046A (zh) * 2013-09-12 2014-01-01 江苏科技大学 水声高速ofdm通信的多普勒扩展处理方法
CN105227512A (zh) * 2015-10-19 2016-01-06 宁波大学 一种ofdm水声通信系统中的脉冲噪声估计方法
CN105282082A (zh) * 2015-09-28 2016-01-27 哈尔滨工程大学 一种基于拷贝相关与空子载波结合的多普勒估计方法
CN107231176A (zh) * 2017-07-24 2017-10-03 哈尔滨工程大学 一种基于子载波能量的ofdm‑mfsk水声通信宽带多普勒估计与补偿方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103354538A (zh) * 2013-07-20 2013-10-16 华南理工大学 一种对水声通信中的接收信号进行多普勒补偿的方法
CN103491046A (zh) * 2013-09-12 2014-01-01 江苏科技大学 水声高速ofdm通信的多普勒扩展处理方法
CN105282082A (zh) * 2015-09-28 2016-01-27 哈尔滨工程大学 一种基于拷贝相关与空子载波结合的多普勒估计方法
CN105227512A (zh) * 2015-10-19 2016-01-06 宁波大学 一种ofdm水声通信系统中的脉冲噪声估计方法
CN107231176A (zh) * 2017-07-24 2017-10-03 哈尔滨工程大学 一种基于子载波能量的ofdm‑mfsk水声通信宽带多普勒估计与补偿方法

Also Published As

Publication number Publication date
CN107911133A (zh) 2018-04-13

Similar Documents

Publication Publication Date Title
CN107911133B (zh) 一种移动水声通信的多普勒因子估计及补偿方法
CN107231176B (zh) 一种ofdm-mfsk水声通信宽带多普勒估计与补偿方法
CN108199994B (zh) 一种无线综测仪ofdm信号频偏估计方法及系统
CN110798422B (zh) 低轨卫星多载波通信系统下行链路采样频偏估计与补偿方法
US7639750B2 (en) Phase tracking method and device thereof
CN107257324B (zh) 一种ofdm系统中的时频联合同步方法及装置
CN111683034B (zh) 一种基于ofdm的大多普勒无线通信时频同步方法
CN101414990A (zh) 一种单载波频域均衡系统中的载波频偏和时延的捕获方法
CN110177066B (zh) 一种5g nr系统中的大频偏估计方法及装置
CN102137049B (zh) 一种lte系统中的单时隙频偏估计方法和装置
CN105704081B (zh) 多个偏移载波调幅信号的优选和解调方法
CN105187352B (zh) 一种基于ofdm前导的整数频偏估计方法
CN103023831B (zh) 一种适用于突发波形的载波频偏估计方法
Wang et al. Doppler estimation and timing synchronization of underwater acoustic communication based on hyperbolic frequency modulation signal
CN105871765A (zh) 一种基于fft辅助s-pll的无线通信载波跟踪方法
CN111711492B (zh) 用于深潜器的水声通信符号定时估计及补偿方法及装置
CN102571652A (zh) 一种gnss基带信号的评估方法
CN113872909B (zh) 一种低轨卫星物联网短时突发载波捕获方法及系统
CN113253313A (zh) 一种改进pmf-fft系统多普勒频率估计精度的方法
CN116106900A (zh) 基于滤波器组多载波的一体化信号设计及处理方法
CN110278170A (zh) 基于最大似然的短波通信频偏估计方法
Zhang et al. Radar sensing via otfs signaling: A delay doppler signal processing perspective
JP5624527B2 (ja) シングルキャリア受信装置
CN101404633A (zh) 基于块传输的单载波系统的载波跟踪方法
CN113438730B (zh) 一种基于gfdm信号的无线定位方法

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
TR01 Transfer of patent right

Effective date of registration: 20230922

Address after: 361000, Xiamen Torch High tech Zone, Xiamen, Fujian Province, China, Software Park 2, No. 10-3, Wangwanghai Road, 603-1

Patentee after: Xiamen Tangpu Information Technology Co.,Ltd.

Address before: Xiamen City, Fujian Province, 361005 South Siming Road No. 422

Patentee before: XIAMEN University

TR01 Transfer of patent right