CN112087266B - 一种基于emd-wfft时变宽带多普勒补偿方法 - Google Patents
一种基于emd-wfft时变宽带多普勒补偿方法 Download PDFInfo
- Publication number
- CN112087266B CN112087266B CN202010849090.7A CN202010849090A CN112087266B CN 112087266 B CN112087266 B CN 112087266B CN 202010849090 A CN202010849090 A CN 202010849090A CN 112087266 B CN112087266 B CN 112087266B
- Authority
- CN
- China
- Prior art keywords
- signal
- channel
- wfft
- function
- carrier
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000004891 communication Methods 0.000 claims abstract description 14
- 230000003044 adaptive effect Effects 0.000 claims abstract description 12
- 238000012952 Resampling Methods 0.000 claims abstract description 9
- 238000012549 training Methods 0.000 claims abstract description 8
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 7
- 238000001514 detection method Methods 0.000 claims description 25
- 230000001427 coherent effect Effects 0.000 claims description 23
- 239000000969 carrier Substances 0.000 claims description 18
- 238000000354 decomposition reaction Methods 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005316 response function Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 2
- 230000033001 locomotion Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- 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/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/2692—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with preamble design, i.e. with negotiation of the synchronisation sequence with transmitter or sequence linked to the algorithm used at the receiver
- H04L27/2694—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with preamble design, i.e. with negotiation of the synchronisation sequence with transmitter or sequence linked to the algorithm used at the receiver adaptive design
-
- 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/2695—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明公开了一种基于EMD‑WFFT时变宽带多普勒补偿方法。步骤1:使用训练序列进行同步;步骤2:将步骤1中的训练序列去掉循环前缀;步骤3:将步骤2得到的信号乘小波函数后再进行WFFT处理;步骤4:将步骤3中经过WFFT处理后的数据再通过EMD处理观测频偏分量;步骤5:将步骤4中EMD处理观测频偏分量后进行重采样,完成对宽带多普勒的粗补偿;步骤6:在经过粗补偿处理后的数据输入自适应梯度算法;步骤7:将数据输入自适应梯度算法后,再通过载波进行信道均衡,去除ICI,对多普勒进行补偿。本发明为了解决水声通信***的信道估计问题。
Description
技术领域
本发明属于水声通信领域;具体涉及一种基于EMD-WFFT时变宽带多普勒补偿方法。
背景技术
水声通信中,起伏的海面,变化的海洋环境,船舶和水下无人潜航器的移动,都会引起多普勒频率偏移。由于在水体中声波的传播速度相对不高,因此即便是比较小的收发端的相对运动都将引起较为严重的多普勒频移。多普勒是影响水声信号检测性能的主要因素,多普勒会破坏正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)载波的正交性,产生载波间的干扰(inter carrier interference,ICI),造成接收信号的频谱畸变,严重影响了数据通信。大量关于多普勒补偿的研究已经展开。研究结果集中于水下平台的高速移动,以及水下宽带信号而产生的宽带多普勒补偿问题。OFDM水声通信中产生的时变多普勒会因子载波频率不同而具有差别,这种宽带多普勒具有非一致性,传统的补偿算法不能抑制信号产生的畸变。
由于传播速度和水下平台运动速度,水声信道中的多普勒要比无线信道大得,主要由于声音在水中传播速度仅为(1500m/s)。水下节点运动时发射的宽带信号导致多普勒在频带内具有非一致性。
水声信道时变宽带多普勒使ICI变得更加复杂,依靠传统FFT后频域均衡的方法,根本无法解决载波间的干扰,而且在进行FFT时,干扰已经混叠在频域内,根本无法消除。传统多普勒抑制采用重采样的方法,这种方法对宽带恒定多普勒有效,也就是通信节点相对速度恒定的情况下。然而实际水下信道中,OFDM等多载波的水声信道是时变的、宽带的,发射节点、接收节点的运动都是随机的。因此原重采样的方法只能对宽带多普勒进行粗补,必然会有残余多普勒无法补偿。
现已提出的水声OFDM信号ICI(inter carrier interference)抑制方法大多专注于后快速傅立叶变换(post-FFT)处理。然而在这样做的过程中,一些有用的信息在FFT解调期间丢失。
发明内容
为了解决水声通信***的信道估计问题,本发明提出一种基于EMD-WFFT时变宽带多普勒补偿方法。
本发明通过以下技术方案实现:
一种基于EMD-WFFT时变宽带多普勒补偿方法,所述补偿方法包括以下步骤:
步骤1:使用训练序列进行同步;
步骤2:将步骤1中的训练序列去掉循环前缀;
步骤3:将步骤2得到的信号乘小波函数后再进行WFFT处理;
步骤4:将步骤3中经过WFFT处理后的数据再通过EMD处理观测频偏分量;
步骤5:将步骤4中EMD处理观测频偏分量后进行重采样,完成对宽带多普勒的粗补偿;
步骤6:在经过粗补偿处理后的数据输入自适应梯度算法;
步骤7:将数据输入自适应梯度算法后,再通过载波进行信道均衡,去除ICI,对多普勒进行补偿。
进一步的,所述步骤2具体为,在水声OFDM通信发射端中,经过串行到并行转换、编码和符号映射的调制过程的发送信号X为:
X=[X0,X1,...,XK-1]
其中,K为子载波个数;
经过IFFT后变成频域信号,通过换能器传输的信号表示为:
其中,K为载波数;N为逆傅里叶变换点数;T为信号OFDM符号的持续时间跨度;j为虚部;k为逆傅里叶变换点数变量去,其取值范围为[0,N-1];t为时间;
接收部分,接收到的信号为:
y(t)=h(t,τ)*x(t)+w(t)
其中,w(t)为噪声;h()*为信道响应共轭;τ为时延;
将接收到的信号的下变频为基带信号,去掉循环前缀后,得到的信号建模为:
其中,W(t)为w(t)的傅里叶变换;Δft为在时间为t时的频偏;Hk(t)为第K个子载波的信道函数;XK为第K个子载波经过调制后的信号。
进一步的,所述步骤3具体为,在最小均方误差准则下,利用一组小波函数逼近信道;
传统的OFDM解调是对接收到的信号进行FFT:
其中,v(t)为信号建模函数;
由于水下信道是一个时变信道,根据增益和延迟,建立信道函数为:
h(t,τ)=A(t)δ(t-τ(t))
其中,A(t)为增益函数;δ(t-τ(t))为冲激响应函数;h(t)为信道响应;H(t)为信道函数;
根据最大似然原理,建立最优接收机,最优接收机前端处理后得到:
信道函数分解到该组基函数上,表示为:
其中,Hk.i为i个小波函数的信道函数;I为小波函数个数;
在WFFT中,通过对Morlet母波进行平移和尺度变换得到一组Morlet小波函数,采用该组小波函数进行信道函数分解,该基函数表达式为:
其中,φi(t)为基函数;ti为第i个小波函数的时间;
则有:
其中,
进一步的,所述步骤3WFFT处理具体为,
定义的输入量为:
yk=[…|yk-1,0…yk-1,I-1|yk,0…yk,I-1|yk+1,0…yk+1,I-1|…]T
此处的y为接收到的信号,k为当前载波,重构部分的输出量xk表示为:
当重构部分的I为1时,WFFT即等效为一个传统的FFT解调结构,它对应一个单抽头解调器;
在接收端,对WFFT信道分解后的信号进行信道重构;差分相干检测技术能有效追踪载波的相位变化,能够抑制OFDM中的ICI。
进一步的,在接收端收到信号后,将截取多普勒粗估计信号,该信号带有频偏信息,将通过相干检测处理后的输入经验模态分解算法,分析频谱成分,进而获得频偏信息;水声信道多普勒频移最终导致时域信号的压缩或扩展,通过重采样技术实现宽带多普勒补偿。
进一步的,在进行EMD-WFFT算法处理后的信号应该按比例进行信号重构,并使用的自适应算法确定近似权重以抑制剩余的ICI,对于未知的信道需要先完成信道估计随后再进行信道均衡。
进一步的,所述信道均衡包括以下步骤:
步骤S1:差分相干检测;
步骤S2:随机梯度算法。
进一步的,所述步骤S1具体为,带宽效率被定义为符号率R与带宽B之比,与1/(1+TgB/K)成比例,其中,Tg是OFDM符号之间的保护间隔,不小于多路径扩展,K为子载波的数量,在给定带宽内,随着载波数量的增加,载波间隔Δf=B/K将减小,因此带宽效率和载波之间的相干性都增加,
由于信道的变化,载波的数量不能无限增加;随着给定带宽内载波的数量增加,信号OFDM符号的持续时间跨度T=1/f增加;
若信道噪声为复基带高斯噪声,接收端接收信号在时域内可表示为:
其中,d(n)为对QPSK信号进行频域差分调制的调制信号序列,并且和均值为0的加性复高斯噪声N(t)不相关,h(t)是包括发射滤波器、接收滤波器和水声信道的全部组合信道的冲激响应,θ(t)为具有载波偏移后的载波相位,d(n)的表达式如下:
d(n)=b(n)d(n-1)
b(n)为QPSK调制信号,
若信号重构部分的抽头数为M,则其输入的采样信号为:
v(n)=[v(nT+T)……v(nT-MT)]T
那么此时均衡器的输出为:
其中a′是a的转置,a为重构部分的抽头系数,
可知输入的b(n)估计如下:
因此获得差分相干检测后的信号与实际信号的误差可表示为:
重构抽头系数的最优解可通过最小化均方误差MMSE来获得,由于误差是一个非线性函数,所以令最小均方误差为0无法得到关于a的随机迭代函数。
其中,xk-1为m阶滤波器的输入为xk-1,...,xk-m,滤波器输出yk是所有输入加权后的线性组合;
为了预防信噪比降低时,对载波相位的跟踪丢失,设置比例梯度为:
其中,gk为随机梯度误差;
根据初始设定的步长递归迭代计算组合器的权值:
其中,μ为步长因子;
当信号中不设置导频符号时,用当前的滤波器的输出和决策值来计算误差ek进而完成重构抽头系数的更新迭代。
本发明的有益效果是:
本发明不需要预先获得信道数据,能跟踪一个OFDM符号内部的多普勒时变。EMD-WFFT相对于传统的FFT,能更好地适应时变信道,多普勒补偿效果明显。信号判决重构过程增加了导频数据的数量,仿真和试验结果表明该算法可行,能够提高***的性能,而且仅使用了很少的导频数据。
附图说明
附图1本发明的流程图。
附图2本发明的OFDM导频与检测顺序。
附图3本发明的WFFT分解框图。
附图4本发明的EMD-WFFT多普勒补偿原理。
附图5本发明的WFFT重构框图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
采取EMD-WFFT方法,为接收器提供了一些pre-FFT信息。同时将这种方法与差分相干检测相结合,以降低全信道估计的复杂性和负担。
一种基于EMD-WFFT时变宽带多普勒补偿方法,所述补偿方法包括以下步骤:
步骤1:使用训练序列即线性调频信号进行同步;
步骤2:将步骤1中的训练序列去掉循环前缀;
步骤3:将步骤2得到的信号乘小波函数后再进行WFFT处理;
步骤4:将步骤3中经过WFFT处理后的数据再通过EMD处理观测频偏分量;
步骤5:将步骤4中EMD处理观测频偏分量后进行重采样,完成对宽带多普勒的粗补偿;
步骤6:在经过粗补偿处理后的数据输入自适应梯度算法;
步骤7:将数据输入自适应梯度算法后,再通过载波进行信道均衡,去除ICI,对多普勒进行补偿。
进一步的,步骤2具体为,在水声OFDM通信发射端中,经过串行到并行转换、编码和符号映射的调制过程的发送信号X为:
X=[X0,X1,...,XK-1]
其中,K为子载波个数;
经过IFFT后变成频域信号,通过换能器传输的信号表示为:
其中,K为载波数;N为逆傅里叶变换点数;T为信号OFDM符号的持续时间跨度;j为虚部;k为逆傅里叶变换点数变量去,其取值范围为[0,N-1];t为时间;
接收部分,接收到的信号为:
y(t)=h(t,τ)*x(t)+w(t)
其中,w(t)为噪声;h()*为信道响应共轭;τ为时延;
将接收到的信号的下变频为基带信号,去掉循环前缀后,得到的信号建模为:
其中,W(t)为w(t)的傅里叶变换;Δft为在时间为t时的频偏;Hk(t)为第K个子载波的信道函数;XK为第K个子载波经过调制后的信号。
进一步的,所述步骤3具体为,在最小均方误差准则下,利用一组小波函数逼近信道;即使用某组选定的基函数逼近接收信号,根据最小均方误差准则求解最优参数值,并使用随机梯度算法更新迭代参数;
传统的OFDM解调是对接收到的信号进行FFT:
其中,v(t)为信号建模函数;
由于水下信道是一个时变信道,根据增益和延迟,建立信道函数为:
h(t,τ)=A(t)δ(t-τ(t))
其中,A(t)为增益函数;δ(t-τ(t))为冲激响应函数;h(t)为信道响应;H(t)为信道函数;
根据最大似然原理,建立最优接收机,最优接收机前端处理后得到:
信道函数分解到该组基函数上,表示为:
其中,Hk.i为i个小波函数的信道函数;I为小波函数个数;
在WFFT中,通过对Morlet母波进行平移和尺度变换得到一组Morlet小波函数,采用该组小波函数进行信道函数分解,该基函数表达式为:
其中,φi(t)为基函数;ti为第i个小波函数的时间;
则有:
其中,
其分解框图如图2。
进一步的,所述步骤3WFFT处理具体为,
定义的输入量为:
yk=[…|yk-1,0…yk-1,I-1|yk,0…yk,I-1|yk+1,0…yk+1,I-1|…]T
此处的y为接收到的信号,k为当前载波,重构部分的输出量xk表示为:
当重构部分的I为1时,WFFT即等效为一个传统的FFT解调结构,它对应一个单抽头解调器;给出重构部分原理框图,如图4;
在接收端,需要对WFFT信道分解后的信号进行信道重构;如果信道冲激响应已知,则可自定义重构抽头系数,然而水声信道是随时间变化的未知信道,接收器必须自适应调整重构系数,这种滤波器自动调整接收信号,称为自适应滤波器,因为它可以跟踪检测并调整到适应***的动态。差分相干检测技术能有效追踪载波的相位变化,能够抑制OFDM中的ICI。
进一步的,在接收端收到信号后,将截取多普勒粗估计信号,该信号带有频偏信息,将通过相干检测处理后的输入经验模态分解(Empirical Mode Decomposition,EMD)算法,分析频谱成分,进而获得频偏信息;水声信道多普勒频移最终导致时域信号的压缩或扩展,通过重采样技术实现宽带多普勒补偿。
进一步的,在进行EMD-WFFT算法处理后的信号应该按比例进行信号重构,并使用的自适应算法确定近似权重以抑制剩余的ICI(Inter Carrier Interference),对于未知的信道需要先完成信道估计随后再进行信道均衡。
进一步的,所述信道均衡包括以下步骤
步骤S1:差分相干检测;
步骤S2:随机梯度算法。
进一步的,所述步骤S1具体为,带宽效率被定义为符号率R与带宽B之比,与1/(1+TgB/K)成比例,其中,Tg是OFDM符号之间的保护间隔,不小于多路径扩展,K为子载波的数量,在给定带宽内,随着载波数量的增加,载波间隔Δf=B/K将减小,因此带宽效率和载波之间的相干性都增加,故本文选择使用频域差分编码;
由于信道的变化,载波的数量不能无限增加;随着给定带宽内载波的数量增加,信号OFDM符号的持续时间跨度T=1/f增加,在这种情况下,载波之间的正交性丢失,由此引起的ICI将引起相干检测和差分相干检测的误差;
水声通信中的多普勒频率变化和水声信道会随时间变化,以及收发端的相对运动等因素会影响载波的检测效果,并最终造成接收机性能的恶化。如何应对载波频率偏移大的局面,已成为水声通信中的一个紧迫问题。目前,相干水声通信中大多数应用的都是具有集成二阶锁相环的决策反馈均衡器的接收机结构,但该结构的性能在大载波频率偏移情况下不是很令人满意。差分相干检测技术可以在QPSK调制中具有天生的载波同步,由于其自身的相位跟踪能力,该技术在水声通信中受到越来越多的重视。差分决策反馈相干检测结构与改进后的LMS算法结合,使***对载波频率的抑制作用进一步提高。相干差分检测自适应算法适应性的结合了线性均衡器与差分相干检测器,以减少平均平方检测误差并成功抑制多路径效应;
若信道噪声为复基带高斯噪声,接收端接收信号在时域内可表示为:
其中,d(n)为对QPSK信号进行频域差分调制的调制信号序列,并且和均值为0的加性复高斯噪声N(t)不相关,h(t)是包括发射滤波器、接收滤波器和水声信道的全部组合信道的冲激响应,θ(t)为具有载波偏移后的载波相位,d(n)的表达式如下:
d(n)=b(n)d(n-1)
b(n)为QPSK调制信号,
若信号重构部分的抽头数为M,则其输入的采样信号为:
v(n)=[v(nT+T)……v(nT-MT)]T
那么此时均衡器的输出为:
其中a′是a的转置,a为重构部分的抽头系数,
可知输入的b(n)估计如下:
因此获得差分相干检测后的信号与实际信号的误差可表示为:
重构抽头系数的最优解可通过最小化均方误差MMSE来获得,由于误差是一个非线性函数,所以令最小均方误差为0无法得到关于a的随机迭代函数;若认为在相邻的两次迭代之间重构抽头系数的变化不大,可采用下文讲述的随机梯度法得到。
进一步的,所述步骤S2具体为,在没有信道信息的情况下,为了获得重构抽头权值ak,将辅助变量结合期望值与导频值之间的误差根据MMSE准则,最终得到单个的权值向量,在实现过程中,算法将把ak和ak-1看作相等时,获得随机梯度误差:
其中,xk-1为m阶滤波器的输入为xk-1,...,xk-m,滤波器输出yk是所有输入加权后的线性组合;
为了预防信噪比降低时,对载波相位的跟踪丢失,设置比例梯度为:
其中,gk为随机梯度误差;
根据初始设定的步长递归迭代计算组合器的权值:
其中,μ为步长因子;
采用了阈值方法,即预先设定误差阈值和梯度阈值,一旦某一值超过设定值,则选择重构抽头系数权值不变,这种阈值方法将有效预防在决策期望值错误时造成重构抽头系数的突然变化,同时这种方法将会进一步提高算法的稳定性能,并防止误差扩大。当信号中不设置导频符号时,用当前的滤波器的输出和决策值来计算误差ek进而完成重构抽头系数的更新迭代。
表10随机梯度算法
Table10Random gradient algorithm
Claims (8)
1.一种基于EMD-WFFT时变宽带多普勒补偿方法,其特征在于,所述补偿方法包括以下步骤:
步骤1:将使用的训练序列进行同步;
步骤2:将步骤1中的训练序列去掉循环前缀;
步骤3:将步骤2得到的信号乘以小波函数后再进行WFFT处理;
步骤4:将步骤3中经过WFFT处理后的数据再通过EMD处理观测频偏分量;
步骤5:将步骤4中EMD处理观测频偏分量后的数据进行重采样,并完成对宽带多普勒的粗补偿;
步骤6:在将经过粗补偿处理后的数据输入自适应梯度算法;
步骤7:将自适应梯度算法计算得到的数据,再通过载波进行信道均衡,去除ICI,对多普勒进行补偿;
所述步骤2具体为,在水声OFDM通信发射端中,经过串行到并行转换、编码和符号映射的调制过程的发送信号X为:
X=[X0,X1,...,XK-1]
其中,K为子载波个数;
经过IFFT后变成时域信号,通过换能器传输的信号表示为:
其中,K为子载波个数;N为逆傅里叶变换点数;T为信号OFDM符号的持续时间跨度;j为虚部;k为当前载波,其取值范围为[0,N-1];t为时间;xk为 第K个被调制的发送信号,
接收部分,接收到的信号为:
y(t)=h(t,τ)*x(t)+w(t)
其中,w(t)为噪声;h()*为信道响应共轭;τ为时延;
接收到的信号的下变频为基带信号,去掉循环前缀后,得到的信号建模为:
其中,W(t)为w(t)的傅里叶变换;Δft为在时间为t时的频偏;Hk(t)为第K个子载波的信道函数;XK为第K个子载波经过调制后的信号。
2.根据权利要求1所述一种基于EMD-WFFT时变宽带多普勒补偿方法,其特征在于,所述步骤3具体为,在最小均方误差准则下,利用一组小波函数逼近信道;
OFDM解调是对接收到的信号进行FFT:
其中,v(t)为信号建模函数;yk为接收到的第k个载波的信号;
由于水下信道是一个时变信道,根据增益和延迟,建立信道函数为:
h(t,τ)=A(t)δ(t-τ(t))
其中,A(t)为增益函数;δ(t-τ(t))为冲激响应函数;h(t)为信道响应;H(t)为信道函数;fk为第k个载波频率;
根据最大似然原理,建立最优接收机,最优接收机前端处理后得到:
信道函数分解到φi(t)基函数上,表示为:
其中,Hk.i为i个小波函数的信道函数;I为小波函数个数;
在WFFT中,通过对Morlet母波进行平移和尺度变换得到一组Morlet小波函数,采用该组小波函数进行信道函数分解,φi(t)基函数表达式为:
其中,φi(t)为基函数;ti为第i个小波函数的时间;
则有:
其中,
4.根据权利要求3所述一种基于EMD-WFFT时变宽带多普勒补偿方法,其特征在于,在接收端收到信号后,将截取多普勒粗估计信号,该信号带有频偏信息,将通过相干检测处理后的输入经验模态分解算法,分析频谱成分,进而获得频偏信息;水声信道多普勒频移最终导致时域信号的压缩或扩展,通过重采样技术实现宽带多普勒补偿。
5.根据权利要求4所述一种基于EMD-WFFT时变宽带多普勒补偿方法,其特征在于,在进行EMD-WFFT算法处理后的信号按比例进行信号重构,并使用的自适应算法确定近似权重以抑制剩余的ICI,对于未知的信道先完成信道估计随后再进行信道均衡。
6.根据权利要求4所述一种基于EMD-WFFT时变宽带多普勒补偿方法,其特征在于,所述信道均衡包括以下步骤:
步骤S1:差分相干检测;
步骤S2:随机梯度算法。
7.根据权利要求6所述一种基于EMD-WFFT时变宽带多普勒补偿方法,其特征在于,所述步骤S1具体为,带宽效率被定义为符号率R与带宽B之比,与1/(1+TgB/K)成比例,其中,Tg是OFDM符号之间的保护间隔,不小于多路径扩展,K为子载波的数量,在给定带宽内,随着载波数量的增加,载波间隔Δf=B/K将减小,因此带宽效率和载波之间的相干性都增加,
由于信道的变化,载波的数量不能无限增加;随着给定带宽内载波的数量增加,信号OFDM符号的持续时间跨度T=1/f增加;
若信道噪声为复基带高斯噪声,接收端接收信号在时域内可表示为:
其中,d(n)为对QPSK信号进行频域差分调制的调制信号序列,并且和均值为0的加性复高斯噪声N(t)不相关,h(t)是包括发射滤波器、接收滤波器和水声信道的全部组合信道的冲激响应,θ(t)为具有载波偏移后的载波相位,d(n)的表达式如下:
d(n)=b(n)d(n-1)
b(n)为QPSK调制信号,
若信号重构部分的抽头数为M,则其输入的采样信号为:
v(n)=[v(nT+T)······v(nT+MT)]
那么此时均衡器的输出为:
其中a′是a的转置,a为重构部分的抽头系数,
可知输入的b(n)估计如下:
因此获得差分相干检测后的信号与实际信号的误差可表示为:
重构抽头系数的最优解可通过最小化均方误差MMSE来获得,由于误差是一个非线性函数,所以令最小均方误差为0无法得到关于a的随机迭代函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010849090.7A CN112087266B (zh) | 2020-08-21 | 2020-08-21 | 一种基于emd-wfft时变宽带多普勒补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010849090.7A CN112087266B (zh) | 2020-08-21 | 2020-08-21 | 一种基于emd-wfft时变宽带多普勒补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112087266A CN112087266A (zh) | 2020-12-15 |
CN112087266B true CN112087266B (zh) | 2022-05-10 |
Family
ID=73729531
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010849090.7A Active CN112087266B (zh) | 2020-08-21 | 2020-08-21 | 一种基于emd-wfft时变宽带多普勒补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112087266B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113098815B (zh) * | 2021-03-05 | 2022-04-01 | 华中科技大学 | 一种水声ofdm载波间干扰的频域抑制方法与*** |
CN114124245B (zh) * | 2021-11-05 | 2024-05-03 | 北京长城电子装备有限责任公司 | 一种时变信道下的水声高速通信方法 |
CN115174330B (zh) * | 2022-06-22 | 2023-08-25 | 苏州大学 | 一种多载波接入网失真信号的补偿方法及非线性均衡器 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4743110A (en) * | 1983-03-29 | 1988-05-10 | Thomson Csf | Laser telemetry and Doppler measurement system with pulse compression |
US6961753B1 (en) * | 2001-04-05 | 2005-11-01 | Osburn Iii Douglas C | Enterprise server for communication for a supervisory control and data acquisition (SCADA) System |
CN102916922A (zh) * | 2012-10-15 | 2013-02-06 | 哈尔滨工程大学 | 水声ofdm自适应搜索多普勒补偿方法 |
CN103095639A (zh) * | 2013-01-15 | 2013-05-08 | 哈尔滨工程大学 | Ofdm水声通信并行迭代ici消除方法 |
CN103888168A (zh) * | 2012-12-21 | 2014-06-25 | 中国科学院声学研究所 | 一种用于合成孔径水声通信的多普勒补偿方法及*** |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8019546B1 (en) * | 2001-04-05 | 2011-09-13 | Dj Inventions, Llc | Tsunami monitoring system |
CN103491046B (zh) * | 2013-09-12 | 2016-06-08 | 江苏科技大学 | 水声高速ofdm通信的多普勒扩展处理方法 |
CN107870005A (zh) * | 2016-09-27 | 2018-04-03 | 重庆邮电大学 | 过采样下经验模态分解的归一化随机共振微弱信号检测 |
CN108566354A (zh) * | 2018-04-03 | 2018-09-21 | 哈尔滨工程大学 | 水声ofdm中dpfft时变宽带多普勒补偿方法 |
CN110456393A (zh) * | 2019-08-21 | 2019-11-15 | 四川航天***工程研究所 | 北斗弱信号快速捕获方法 |
CN110784428B (zh) * | 2019-11-01 | 2022-10-21 | 哈尔滨工程大学 | 水声通信网络中基于morl-FFT的自适应多普勒补偿方法 |
-
2020
- 2020-08-21 CN CN202010849090.7A patent/CN112087266B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4743110A (en) * | 1983-03-29 | 1988-05-10 | Thomson Csf | Laser telemetry and Doppler measurement system with pulse compression |
US6961753B1 (en) * | 2001-04-05 | 2005-11-01 | Osburn Iii Douglas C | Enterprise server for communication for a supervisory control and data acquisition (SCADA) System |
CN102916922A (zh) * | 2012-10-15 | 2013-02-06 | 哈尔滨工程大学 | 水声ofdm自适应搜索多普勒补偿方法 |
CN103888168A (zh) * | 2012-12-21 | 2014-06-25 | 中国科学院声学研究所 | 一种用于合成孔径水声通信的多普勒补偿方法及*** |
CN103095639A (zh) * | 2013-01-15 | 2013-05-08 | 哈尔滨工程大学 | Ofdm水声通信并行迭代ici消除方法 |
Non-Patent Citations (2)
Title |
---|
EMD算法应用于OFDM***频偏估计的新方法;苏亮等;《集美大学学报(自然科学版)》;20160725(第04期);275-281 * |
水声宽带信号的多普勒补偿;黄建纯等;《声学技术》;20090415(第02期);99-103 * |
Also Published As
Publication number | Publication date |
---|---|
CN112087266A (zh) | 2020-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7609789B2 (en) | Phase noise compensation for MIMO WLAN systems | |
CN112087266B (zh) | 一种基于emd-wfft时变宽带多普勒补偿方法 | |
CN106789764B (zh) | 联合加权阈值去噪与均衡判决的变换域二次估计方法 | |
US7912118B2 (en) | Hybrid domain block equalizer | |
JP4958565B2 (ja) | 無線通信装置 | |
Zhang et al. | Frequency-domain turbo equalization with soft successive interference cancellation for single carrier MIMO underwater acoustic communications | |
US20040184550A1 (en) | Receiver which demodulates OFDM symbol | |
JP2006262039A (ja) | 伝搬路推定方法及び伝搬路推定装置 | |
JP3715282B2 (ja) | Ofdm受信装置及びofdm信号の補正方法 | |
Han et al. | Iterative per-vector equalization for orthogonal signal-division multiplexing over time-varying underwater acoustic channels | |
CN107483373B (zh) | 一种抗多径迭代加权的lmmse信道估计方法及装置 | |
US10103909B2 (en) | Channel estimation method | |
JP4311132B2 (ja) | Ofdm伝送方式における受信装置 | |
CN101133580A (zh) | 接收装置 | |
US7340000B1 (en) | Decision feedback equalizer in an OFDM WLAN receiver | |
US8279743B2 (en) | Method for interference estimation for orthogonal pilot patterns | |
Yang et al. | Novel channel estimation method based on PN sequence reconstruction for Chinese DTTB system | |
US20100074346A1 (en) | Channel estimation in ofdm receivers | |
CN110784428B (zh) | 水声通信网络中基于morl-FFT的自适应多普勒补偿方法 | |
WO2018000914A1 (zh) | 一种单载波信道估计方法 | |
KR20060071084A (ko) | Ofdm 통신시스템에서 위상잡음 보상을 위한 등화기 및그 방법 | |
US8619923B2 (en) | System and method for optimizing use of channel state information | |
KR100602518B1 (ko) | 직교주파수분할다중 통신 시스템의 채널 추정 방법 및 장치 | |
CN111245589B (zh) | 一种导频叠加信道估计方法 | |
CN110730144B (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 |