CN112953653B - 一种单载波多用户水声通信方法 - Google Patents

一种单载波多用户水声通信方法 Download PDF

Info

Publication number
CN112953653B
CN112953653B CN202110226022.XA CN202110226022A CN112953653B CN 112953653 B CN112953653 B CN 112953653B CN 202110226022 A CN202110226022 A CN 202110226022A CN 112953653 B CN112953653 B CN 112953653B
Authority
CN
China
Prior art keywords
channel
user
interference
signal
iteration
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
CN202110226022.XA
Other languages
English (en)
Other versions
CN112953653A (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.)
Harbin Xinguang Photoelectric Technology Co ltd
Original Assignee
Harbin Engineering 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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202110226022.XA priority Critical patent/CN112953653B/zh
Publication of CN112953653A publication Critical patent/CN112953653A/zh
Application granted granted Critical
Publication of CN112953653B publication Critical patent/CN112953653B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/0224—Channel estimation using sounding signals
    • 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/0256—Channel estimation using minimum mean square error criteria
    • 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/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03178—Arrangements involving sequence estimation techniques
    • H04L25/03248—Arrangements for operating in conjunction with other apparatus
    • H04L25/03254—Operation with other circuitry for removing intersymbol interference
    • H04L25/03267—Operation with other circuitry for removing intersymbol interference with decision feedback equalisers
    • 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/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03878—Line equalisers; line build-out devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

本发明公开了一种单载波多用户水声通信方法，基于单载波相移键控调制体制，将多个用户的通信信息调制在载波相位上，信源比特信息通过卷积编码器进行信道编码，再采用随机交织器交织，确定接收端通信信号的解码流程，通信接收机采用基于软干扰消除的迭代信道估计器，将接收信号中的干扰用户信号重构并消除，利用产生的信号更新期望用户的信道，被动时间反转镜将消除干扰后的多通道信号合并为单通道信号，在压缩多途信道的同时提高期望用户信噪比，再将单通道信号输入自适应turbo均衡器中，软均衡器输出的符号经过信道解码器解码后输出比特数据，接收机再经过多次迭代实现解码，能够进行用户数较多的数据解码且接收机计算复杂度较低。

Description

一种单载波多用户水声通信方法
技术领域
本发明涉及一种单载波多用户水声通信方法，特别是一种基于迭代软干扰消除的单载波多用户水声通信方法，属于水声通信领域。
背景技术
无人潜航器被广泛应用于海洋资源开发、海洋环境观测、情报收集、装备预置、探测预警等领域，并逐渐形成多无人潜航器协同观测模式。但是，随着任务内容和复杂度的上升，单一的无人潜航器已很难满足日益增长的需求。这使得以多潜航器组成水下集群协同执行任务成为必然的发展方向。
与发展较为成熟的单用户水下通信技术相比，用来支持水下集群通信的多用户通信技术仍然面临较多的困难和挑战，例如水下环境中严重受限的带宽，长达几十毫秒甚至上百毫秒的多途时延信道，相比于无线电通信更加严重的多普勒扩展，信道的快速时变性以及强烈的多址干扰。当水下多用户通信系统的用户数量增多时，多址干扰会逐渐成为影响水下多用户通信机性能最重要的因素。许多技术用来实现水下多用户通信系统，其中直接扩频序列码分多址技术被认为是最具发展前景的技术之一，但是该技术的数据速率过低，无法满足未来场景中各种应用的高速数据传输要求。而交织分多址作为与码分多址相似的技术，被证明在性能上优于基于长码和短码的码分多址系统，但是其仍然面临码率较低且计算复杂度较高等问题。正交频分复用技术作为无线电通信中最关键的技术之一，因为其能够灵活地为各个用户分配子载波和功率而受到广泛的研究。但是正交频分复用技术具有较高的峰均功率比，而且对载波频率偏移十分敏感，这都使得其在水声多用户通信中的应用受到巨大的限制。相比之下，单载波多用户通信系统具有较高的通信速率，并且在峰均功率比和对抗载波频率偏移等方面具有优势，能够满足未来应用对水下多用户通信的需求，但是单载波多用户通信系统面临强多址干扰的巨大挑战。
如何在快速时变的水声信道中且具有强多址干扰的条件下，对单载波多用户通信数据进行高性能解码，成为水声通信技术中重要的技术问题。
发明内容
针对上述现有技术，本发明要解决的技术问题是提供一种基于迭代软干扰消除的单载波多用户水声通信方法，能够进行用户数较多的数据解码且接收机计算复杂度较低。
为解决上述技术问题，本发明的一种单载波多用户水声通信方法，包括以下步骤：
步骤1：输入：第i块接收信号ym,i，m∈[1,M]，M为接收阵元数，全部用户的导频向量Xn，其中，n∈[1,N]，N为用户的数量，即发射天线数，根据信道参数将接收数据划分为NB个数据块、迭代次数Iter，i∈[1,NB]；
步骤2：信道和参数初始化，具体为：根据反馈的软信息，计算反馈符号
Figure BDA0002956183750000021
更新迭代次数Iter，利用MMSE信道估计算法进行初始信道估计，得到
Figure BDA0002956183750000022
其中m∈[1,M]，表示发射阵元n与接收阵元m之间的信道，“1”表示第“1”个数据块；
步骤3：块间干扰消除，消除因为缺少保护间隔而产生的前一个数据块对当前数据块的影响；
步骤4：连续干扰消除，利用上一次迭代反馈的先验符号
Figure BDA0002956183750000023
以及估计的信道
Figure BDA0002956183750000024
重构多址干扰
Figure BDA0002956183750000025
其中l∈[1,N],l≠n，n表示期望用户索引值，l表示干扰用户索引值，并对多址干扰的进行幅度补偿和时间延迟，最后从接收信号中减去，得到仅包含期望信号和噪声的多通道数据
Figure BDA0002956183750000026
步骤5：信道更新，利用消除多址干扰之后的信号
Figure BDA0002956183750000027
更新期望用户信道，使用MMSE信道估计算法，具体为：
Figure BDA0002956183750000028
其中，
Figure BDA0002956183750000029
是噪声功率，IL是维度为L的单位矩阵，(·)H表示矩阵的共轭转置，
Figure BDA00029561837500000210
定义为上一次迭代反馈的符号；
步骤6：时间反转合并，利用被动时间反转合并技术将接收的多通道信号转换成单通道数据
Figure BDA00029561837500000211
具体为：
Figure BDA00029561837500000212
其中，Xn,i和Xl,i分别表示用户n和用户l的第i块发射符号矩阵，qnn,m,i和qnl,m,i为Q函数，
Figure BDA00029561837500000213
为干扰消除尺度因子，zn为噪声项；
步骤7：软判决反馈均衡，使用软判决反馈均衡消除残留的码间干扰，软判决反馈均衡器的输出为：
Figure BDA00029561837500000214
其中，ffk和fbk分别为软判决反馈均衡器的前馈和反馈滤波器，k表示第k个符号，并且使用递归最小二乘自适应算法更新，
Figure BDA00029561837500000215
和
Figure BDA00029561837500000216
分别表示k时刻对应
Figure BDA00029561837500000217
和
Figure BDA00029561837500000218
中的滑动窗向量；
步骤8：判断n＝N是否成立，若成立，则当前数据块所有用户检测完成，执行步骤9，进行下一个数据块的处理；否则，执行步骤4，继续检测当前数据的剩余用户；
步骤9：判断i＝NB是否成立，若成立，则当前迭代的数据块全部处理完成，执行步骤10对检测的符号进行信道解码；否则，执行步骤3，继续检测下一个数据块；
步骤10：计算对数似然比、解交织、信道解码，输出解码结果和反馈软信息，反馈外信息作为下一次迭代的先验信息：
Figure BDA0002956183750000031
其中，Le(cn,j)表示用户n比特cn,j的对数似比，
Figure BDA0002956183750000032
表示以条件概率计算的对数似然比，L(cn,j)表示上一次迭代反馈的对数似然比；
步骤11：判断迭代次数是否达到Iter，若达到最大迭代次数，则解码完成，输出
Figure BDA0002956183750000033
表示用户n比特序列的估计结果，否则执行步骤2继续迭代。
本发明还包括：
1.步骤3中块间干扰消除，消除因为缺少保护间隔而产生的前一个数据块对当前数据块的影响具体为：重构块间干扰，满足：
Figure BDA0002956183750000034
其中，
Figure BDA0002956183750000035
且
Figure BDA0002956183750000036
其中k的取值范围为[Nb+1,Nb+L+1]，Nb是每一个数据块的长度，将重构的块间干扰信号从当前块中减去。
2.步骤4中仅包含期望信号和噪声的多通道数据
Figure BDA0002956183750000037
具体为：
Figure BDA0002956183750000038
其中，
Figure BDA0002956183750000039
为上一次迭代产生的先验符号，在第一次迭代中，为上一个数据块估计结果的硬判决符号；干扰尺度因子
Figure BDA0002956183750000041
其中，
Figure BDA0002956183750000042
是重构干扰信号与目标信号的匹配滤波输出结果，τ是匹配滤波器输出最大峰值的时延值，Pl是发射信号的平均功率，N0是噪声功率。
本发明的有益效果：本发明与传统的多用户水声通信方法不同之处在于通过软干扰消除获得的信号更新期望用户的信道，该信道用来支持高性能的时间反转合并。而时间反转镜将消除多址干扰的多通道信号合并为单通道信号，再将其输入自适应turbo均衡器获得期望用户的符号估计，最后通过信道解码器解码得到比特信息。本发明的优点在于(1)能够有效对抗水下强多址干扰；(2)可降低接收机的计算复杂度；(3)可实现较高的通信速率。(4)可用来实现水下无线网络之间的数据传输。
附图说明
图1是基于迭代软干扰消除的单载波多用户水声通信技术流程图；
图2(a)是试验布置图；
图2(b)是发射接收位置分布阵元深度；
图3用户5迭代处理结果，(a1)-(a4)是均衡器输出星座图(b1)-(b4)为对应解码器输出星座图；
图4是多用户数据处理误码率统计结果。
具体实施方式
下面结合附图对本发明具体实施方式做进一步说明。
本发明通过下述技术方案予以实现：基于单载波相移键控调制体制，将多个用户的通信信息调制在载波相位上。信源比特信息通过卷积编码器进行信道编码，再采用随机交织器交织。确定接收端通信信号的解码流程。通信接收机采用基于软干扰消除的迭代信道估计器，将接收信号中的干扰用户信号重构并消除，利用产生的信号更新期望用户的信道。被动时间反转镜将消除干扰后的多通道信号合并为单通道信号，在压缩多途信道的同时提高期望用户信噪比。再将单通道信号输入自适应turbo均衡器中，软均衡器输出的符号经过信道解码器解码后输出比特数据。接收机再经过多次迭代实现高性能的解码。
本发明的目的是这样实现的：
步骤1：输入第i块接收信号ym,i，其中，m∈[1,M]，M为接收阵元数。全部用户的导频矩阵Xn，其中，n∈[1,N]，N为用户的数量，即发射天线数。根据信道参数将接收数据划分为NB个数据块，接收机迭代次数Iter。
步骤2：信道和参数初始化。根据反馈的软信息，计算反馈符号
Figure BDA0002956183750000051
更新迭代次数Iter。利用MMSE信道估计算法进行初始信道估计，得到
Figure BDA0002956183750000052
表示发射阵元n与接收阵元m之间的信道，“1”表示第“1”个数据。
步骤3：块间干扰消除。消除因为缺少保护间隔而产生的前一个数据块对当前数据块的影响。
步骤4：连续干扰消除。利用上一次迭代反馈的先验符号
Figure BDA0002956183750000053
以及估计的信道
Figure BDA0002956183750000054
重构多址干扰
Figure BDA0002956183750000055
其中l∈[1,N],l≠n，n表示期望用户索引值，l表示干扰用户索引值。并对多址干扰的进行幅度补偿和时间延迟，最后从接收信号中减去，得到仅包含期望信号和噪声的多通道数据
Figure BDA0002956183750000056
步骤5：信道更新。利用消除多址干扰之后的信号
Figure BDA0002956183750000057
更新期望用户信道，使用MMSE信道估计算法。
Figure BDA0002956183750000058
其中，
Figure BDA0002956183750000059
是噪声功率，IL是维度为L的单位矩阵，(·)H表示矩阵的共轭转置。
步骤6：时间反转合并。在软干扰消除的基础上，进一步抑制多址干扰。同时，将多通道信号转换成单通道数据
Figure BDA00029561837500000510
Figure BDA00029561837500000511
其中，Xn,i和Xl,i分别表示用户n和用户l的第i块发射符号矩阵，qnn,m,i和qnl,m,i为Q函数，
Figure BDA00029561837500000512
为干扰消除尺度因子，zn为噪声项。
步骤7：软判决反馈均衡。使用软判决反馈均衡消除残留的码间干扰。
Figure BDA00029561837500000513
其中，ffk和fbk分别为软判决反馈均衡器的前馈和反馈滤波器，k表示第k个符号，
Figure BDA00029561837500000514
和
Figure BDA00029561837500000515
分别表示k时刻对应
Figure BDA00029561837500000516
和
Figure BDA00029561837500000517
中的滑动窗向量。
步骤8：判断：n＝N，若当前数据块所有用户检测完成，则进行下一个数据块的处理。否则，继续检测当前数据的剩余用户。
步骤9：判断：i＝NB，若当前迭代的数据块全部处理完成，则对检测的符号进行信道解码，否则，继续检测下一个数据块。
步骤10：计算对数似然比、解交织、信道解码，反馈外信息作为下一次迭代的先验信息。
Figure BDA0002956183750000061
其中,Le(cn,j)表示用户n比特cn,j的对数似比，
Figure BDA0002956183750000062
表示以条件概率计算的对数似然比，L(cn,j)表示上一次迭代反馈的对数似然比。
步骤11：判断：迭代次数是否达到Iter，若达到，则解码完成，输出
Figure BDA0002956183750000063
表示用户n比特序列的估计结果,否则继续迭代。
结合图1，本发明具体实施方式具体为：
1、水下多用户通信技术的具体实现：
(1)搭建多用户水声通信系统，将多用户系统建模为多个单输入多输出系统的叠加。接收信号模型为
Figure BDA0002956183750000064
其中，
Figure BDA0002956183750000065
是第m个水听器的接收的第i个数据块，Nb是每一个数据块的长度，
Figure BDA0002956183750000066
为第n个发射换能器与第m个接收水听器之间的信道，信道长度为L，wm,i为第m个水听器上采集的均值为0方差σ2的加性高斯白噪声，Xn,i定义为
Figure BDA0002956183750000067
其中的xn,k为发射端发送的符号。
(2)信道和参数初始化。根据反馈的软信息，计算先验符号
Figure BDA0002956183750000068
更新迭代次数Iter。利用MMSE信道估计算法进行初始信道估计，得到
Figure BDA0002956183750000069
其中m∈[1,M]。
(3)块间干扰消除。利用前一个数据块的检测符号和估计的信道，重构块间干扰。并且，块间干扰消除在每一个数据块开始处理之前进行。重构的块间干扰符号为
Figure BDA00029561837500000610
其中，
Figure BDA0002956183750000071
且
Figure BDA0002956183750000072
其中k的取值范围为[Nb+1,Nb+L+1]。将重构的块间干扰信号从当前块中减去。
(4)连续干扰消除。
(a)基于信干噪比的排序。计算所有用户的信干噪比SINR，多用户接收机的解码顺序就由各个用户的SINR大小顺序决定，SINR大的用户先解码
Figure BDA0002956183750000073
其中，Pn为期望用户n的平均发射功率，Pl为干扰用户l的平均发射功率，N0是噪声功率。
(b)多址干扰重构。
Figure BDA0002956183750000074
其中，
Figure BDA0002956183750000075
为上一次迭代产生的用户l的先验符号，在第一次迭代中，为上一个数据块估计结果的硬判决符号。
(c)干扰尺度因子。
Figure BDA0002956183750000076
其中，
Figure BDA0002956183750000077
是重构干扰信号与目标信号的匹配滤波输出结果，τ是匹配滤波器输出最大峰值的时延值。
(e)干扰消除。将重构的多址干扰从接收信号中减去，结果表示为
Figure BDA0002956183750000078
其中，
Figure BDA0002956183750000079
(5)信道更新。利用消除多址干扰之后的信号更新期望用户信道，使用MMSE信道估计算法。
Figure BDA00029561837500000710
其中，
Figure BDA00029561837500000711
是噪声功率，IL是维度为L的单位矩阵，(·)H表示矩阵的共轭转置，
Figure BDA00029561837500000712
定义为
Figure BDA0002956183750000081
(6)时间反转合并。利用被动时间反转合并技术将接收的多通道信号合并成单通道信号
Figure BDA0002956183750000082
其中，qnn,m,i和qnl,m,i为Q函数，Xn,i和Xl,i分别表示用户n和用户l的第i块发射符号矩阵，zn为噪声项。
(7)软判决反馈均衡。软判决反馈均衡器的输出为
Figure BDA0002956183750000083
其中，ffk和fbk分别为软判决反馈均衡器的前馈和反馈滤波器，并且使用递归最小二乘自适应算法更新。k表示第k个符号，
Figure BDA0002956183750000084
和
Figure BDA0002956183750000085
分别表示k时刻对应
Figure BDA0002956183750000086
和
Figure BDA0002956183750000087
中的滑动窗向量。
Figure BDA0002956183750000088
的时间滑动窗范围为[k-K2,k+K1]，
Figure BDA0002956183750000089
的滑动窗范围为[k-K4,k+K3]，其中，K1,K2,K3和K4表示均衡器抽头系数长度。
(8)判断：n＝N，若当前数据块所有用户检测完成，则进行下一个数据块的处理。否则，继续检测当前数据的剩余用户。
(9)判断：i＝NB，若当前迭代的数据块全部处理完成，则对检测的符号进行信道解码，否则，继续检测下一个数据块。
(10)计算对数似然比、解交织、信道解码，输出解码结果和反馈软信息。
Figure BDA00029561837500000810
其中,Le(cn,j)表示用户n比特cn,j的对数似比，
Figure BDA00029561837500000811
表示以条件概率计算的对数似然比，L(cn,j)表示上一次迭代反馈的对数似然比。
(11)判断：迭代次数是否达到Iter，若达到，则解码完成，输出
Figure BDA00029561837500000812
表示用户n比特序列的估计结果,否则继续迭代。
2、试验研究
试验条件：
2019年松花湖试验布置如图2(a)所示，其中R0为32元垂直阵的布放位置，P1和P2为两个分开的发射位置，分别距离R0位置568m和1440m。试验中收发阵元的布放如图2(b)所示，实验地点水深约为70m，R0处的32元阵间距为1m，其1号阵元H1距离水面2m。P1和P2位置均使用间距为2m的4元阵发射多用户信号，4元阵分布布置在10-16m的水深范围内。试验中发射的多用户数据均为单载波相移键控调制信号，系统采样率为100kHz，载波频率为10kHz，根升余弦滤波器滚降因子为1。其中，P1点发射3路同频同时的信号，符号周期为1ms，故带宽为2kHz。P2点发射4路通频同时的信号，符号周期为0.5ms，所以带宽为4kHz。另外，在P2位置单独发射了两组单载波信号，系统采样率为96kHz，载波频率分别为12kHz和15kHz，带宽均为4kHz，使用QPSK调制方式。其发射深度均为10m。将采集到的4用户数据、3用户数据和两组单用户数据在基带叠加成9用户数据。使用30个阵元的接收数据进行解码。
图3是用户5数据4次迭代处理结果中均衡器和解码器输出散点图。从图中可以看出，第一次迭代均衡器输出的星座图信噪比较低，导致后续的信道解码器输出结果也存在误码。而经过两次迭代之后接收机达到收敛，信噪比性能具有很大提升。
表1展示了9×30多用户通信数据的误码率结果。从图中可以看到，相对于均衡器的误码率结果，信道解码器利用纠错码能够实现极大的误码率性能提升。用户4、5和6因为具有较低的功率，解码结果性能低于高功率用户。经过3次迭代，均衡器输出误码率性能达到收敛，无法进一步提升。而信道解码器在第二次迭代完成时，就实现了所有用户数据的无误码解码。
表1
Figure BDA0002956183750000091
图4比较了传统接收机和本发明接收机的误码率性能。其中，I表示MMSE信道估计器和基于被动时间反转镜的单通道判决反馈均衡器，II表示硬判决干扰消除MMSE信道估计器和基于被动时间反转镜的单通道判决反馈均衡器，III表示软判决干扰消除MMSE信道估计器和基于被动时间反转镜的单通道软判决反馈均衡器，IV表示软判决干扰消除MMSE信道估计器和基于被动时间反转镜的单通道自适应turbo均衡器。从图中可以看出，经过4次迭代，传统的多用户接收机无法实现高性能的数据解码，而本发明的接收机经过2次迭代就实现了采集数据的无误码解码。

Claims (1)

1.一种单载波多用户水声通信方法，其特征在于，包括以下步骤：
搭建多用户水声通信系统，将多用户系统建模为多个单输入多输出系统的叠加，接收信号模型为：
Figure FDA0003636215440000011
其中，
Figure FDA0003636215440000012
是第m个水听器的接收的第i个数据块，Nb是每一个数据块的长度，
Figure FDA0003636215440000013
为第n个发射换能器与第m个接收水听器之间的信道，信道长度为L，wm,i为第m个水听器上采集的均值为0方差σ2的加性高斯白噪声，Xn,i定义为
Figure FDA0003636215440000014
其中的xn,k为发射端发送的符号；
步骤1：输入：第i块接收信号ym,i，m∈[1,M]，M为接收阵元数，全部用户的导频向量Xn，其中，n∈[1,N]，N为用户的数量，即发射天线数，根据信道参数将接收数据划分为NB个数据块、迭代次数Iter，i∈[1,NB]；
步骤2：信道和参数初始化，具体为：根据反馈的软信息，计算反馈符号
Figure FDA0003636215440000015
更新迭代次数Iter，利用MMSE信道估计算法进行初始信道估计，得到
Figure FDA0003636215440000016
其中m∈[1,M]，表示发射阵元n与接收阵元m之间的信道，“1”表示第“1”个数据块；
步骤3：块间干扰消除，消除因为缺少保护间隔而产生的前一个数据块对当前数据块的影响，具体为：块间干扰消除在每一个数据块开始处理之前进行，重构的块间干扰符号为：
Figure FDA0003636215440000017
其中，
Figure FDA0003636215440000018
且
Figure FDA0003636215440000019
其中k的取值范围为[Nb+1,Nb+L+1]，将重构的块间干扰信号从当前块中减去；
步骤4：连续干扰消除，利用上一次迭代反馈的先验符号
Figure FDA00036362154400000110
以及估计的信道
Figure FDA00036362154400000111
重构多址干扰
Figure FDA0003636215440000021
其中l∈[1,N],l≠n，n表示期望用户索引值，l表示干扰用户索引值，并对多址干扰的进行幅度补偿和时间延迟，最后从接收信号中减去，得到仅包含期望信号和噪声的多通道数据
Figure FDA0003636215440000022
具体为：多址干扰重构：
Figure FDA0003636215440000023
其中，
Figure FDA0003636215440000024
为上一次迭代产生的用户l的先验符号，在第一次迭代中，为上一个数据块估计结果的硬判决符号；
干扰尺度因子：
Figure FDA0003636215440000025
其中，
Figure FDA0003636215440000026
是重构干扰信号与目标信号的匹配滤波输出结果，τ是匹配滤波器输出最大峰值的时延值，Pl是发射信号的平均功率，N0是噪声功率；
将重构的多址干扰从接收信号中减去，结果表示为：
Figure FDA0003636215440000027
其中，
Figure FDA0003636215440000028
步骤5：信道更新，利用消除多址干扰之后的信号
Figure FDA0003636215440000029
更新期望用户信道，使用MMSE信道估计算法，具体为：
Figure FDA00036362154400000210
其中，
Figure FDA00036362154400000211
是噪声功率，IL是维度为L的单位矩阵，(·)H表示矩阵的共轭转置，
Figure FDA00036362154400000212
定义为上一次迭代反馈的符号；
步骤6：时间反转合并，利用被动时间反转合并技术将接收的多通道信号转换成单通道数据
Figure FDA00036362154400000213
具体为：
Figure FDA00036362154400000214
其中，Xn,i和Xl,i分别表示用户n和用户l的第i块发射符号矩阵，qnn,m,i和qnl,m,i为Q函数，
Figure FDA00036362154400000215
为干扰消除尺度因子，zn为噪声项；
步骤7：软判决反馈均衡，使用软判决反馈均衡消除残留的码间干扰，软判决反馈均衡器的输出为：
Figure FDA0003636215440000031
其中，ffk和fbk分别为软判决反馈均衡器的前馈和反馈滤波器，k表示第k个符号，并且使用递归最小二乘自适应算法更新，
Figure FDA0003636215440000032
和
Figure FDA0003636215440000033
分别表示k时刻对应
Figure FDA0003636215440000034
和
Figure FDA0003636215440000035
中的滑动窗向量，
Figure FDA0003636215440000036
的时间滑动窗范围为[k-K2,k+K1]，
Figure FDA0003636215440000037
的滑动窗范围为[k-K4,k+K3]，其中，K1,K2,K3和K4表示均衡器抽头系数长度；
步骤8：判断n＝N是否成立，若成立，则当前数据块所有用户检测完成，执行步骤9，进行下一个数据块的处理；否则，执行步骤4，继续检测当前数据的剩余用户；
步骤9：判断i＝NB是否成立，若成立，则当前迭代的数据块全部处理完成，执行步骤10对检测的符号进行信道解码；否则，执行步骤3，继续检测下一个数据块；
步骤10：计算对数似然比、解交织、信道解码，输出解码结果和反馈软信息，反馈外信息作为下一次迭代的先验信息：
Figure FDA0003636215440000038
其中，Le(cn,j)表示用户n比特cn,j的对数似比，
Figure FDA0003636215440000039
表示以条件概率计算的对数似然比，L(cn,j)表示上一次迭代反馈的对数似然比；
步骤11：判断迭代次数是否达到Iter，若达到最大迭代次数，则解码完成，输出
Figure FDA00036362154400000310
Figure FDA00036362154400000311
表示用户n比特序列的估计结果，否则执行步骤2继续迭代。
CN202110226022.XA 2021-03-01 2021-03-01 一种单载波多用户水声通信方法 Active CN112953653B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110226022.XA CN112953653B (zh) 2021-03-01 2021-03-01 一种单载波多用户水声通信方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110226022.XA CN112953653B (zh) 2021-03-01 2021-03-01 一种单载波多用户水声通信方法

Publications (2)

Publication Number Publication Date
CN112953653A CN112953653A (zh) 2021-06-11
CN112953653B true CN112953653B (zh) 2022-07-15

Family

ID=76246969

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110226022.XA Active CN112953653B (zh) 2021-03-01 2021-03-01 一种单载波多用户水声通信方法

Country Status (1)

Country Link
CN (1) CN112953653B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115208482B (zh) * 2022-06-30 2023-02-03 哈尔滨工程大学 一种极地脉冲干扰下的水声通信方法
CN115208480B (zh) * 2022-06-30 2023-05-12 哈尔滨工程大学 一种基于联合消息传递的冰下水声通信方法
CN115694578B (zh) * 2022-09-19 2023-11-24 华工未来科技(江苏)有限公司 一种多智能反射面同步控制方法、系统、装置及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0822926D0 (en) * 2008-12-17 2009-01-21 Yeo Hong K Multiple output multiple user methods and/or systems of underwater acoustic communication
CN109462443A (zh) * 2018-12-05 2019-03-12 唐山照澜海洋科技有限公司 一种5g多载波水声通信方法
CN111901011A (zh) * 2020-06-19 2020-11-06 中国船舶重工集团公司第七一五研究所 一种基于单矢量水听器的多用户水声通信方法
CN112134822A (zh) * 2020-09-22 2020-12-25 天津工业大学 一种基于ofdm+fdma技术提高海洋用感应耦合锚系链多用户传输性能的方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ATE526732T1 (de) * 1999-11-29 2011-10-15 Commw Scient Ind Res Org Kommunikationssystem
US8116411B2 (en) * 2007-08-06 2012-02-14 Indian Institute Of Science Method to detect data transmitted from multiple antennas and system thereof
CN110768719B (zh) * 2019-09-25 2020-09-11 北京科技大学 一种双路水下可见光通信系统及通信方法
CN111147157B (zh) * 2019-12-30 2021-04-23 浙江大学 一种用于水声信道上正交信号分复用的软干扰消除Turbo均衡方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0822926D0 (en) * 2008-12-17 2009-01-21 Yeo Hong K Multiple output multiple user methods and/or systems of underwater acoustic communication
CN109462443A (zh) * 2018-12-05 2019-03-12 唐山照澜海洋科技有限公司 一种5g多载波水声通信方法
CN111901011A (zh) * 2020-06-19 2020-11-06 中国船舶重工集团公司第七一五研究所 一种基于单矢量水听器的多用户水声通信方法
CN112134822A (zh) * 2020-09-22 2020-12-25 天津工业大学 一种基于ofdm+fdma技术提高海洋用感应耦合锚系链多用户传输性能的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
IDMA水声多用户下行通信链路迭代接收算法;贾宁等;《信号处理》;20190930;第35卷(第9期);全文 *
Underwater Acoustic CommuniTime Delay Shift Coding Schemecation Based on Pattern;殷敬伟;《China Ocean Engineering》;20061231;第20卷(第3期);全文 *
浅海环境下单载波时域均衡水声通信关键技术研究;韩笑;《中国博士学位论文全文数据库 信息科技辑》;20180615(第6期);全文 *
高速水声MIMO通信中的迭代信道估计与均衡技术研究;黄福朋;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20190115(第1期);全文 *

Also Published As

Publication number Publication date
CN112953653A (zh) 2021-06-11

Similar Documents

Publication Publication Date Title
CN112953653B (zh) 一种单载波多用户水声通信方法
CN108833311B (zh) 联合时域聚类去噪与均衡判决的变换域二次估计方法
CN106789764B (zh) 联合加权阈值去噪与均衡判决的变换域二次估计方法
CN101621327B (zh) 单载波发送系统中的无线电通信方法和装置
EP1233565B1 (en) Turbo-reception method and turbo-receiver for a MIMO system
CN108900443B (zh) 一种水声通信中的水声信道干扰消除方法
US8325588B2 (en) Iterative reception method and iterative receiver
US20040132416A1 (en) Equalisation apparatus and methods
CN100395975C (zh) 多载波扩频信号的接收
EP1861918A2 (en) Iterative mimo receiver using group-wise demapping
CN101138212B (zh) 用于复值符号的软解映射的方法
WO2016030758A2 (en) Multiple input multiple output communications over nonlinear channels using orthogonal frequency division multiplexing
Qin et al. Bayesian iterative channel estimation and turbo equalization for multiple-input–multiple-output underwater acoustic communications
CN113381951A (zh) 时变频选衰落信道下mftn联合信道估计与均衡方法
Yin et al. Iterative channel estimation-based soft successive interference cancellation for multiuser underwater acoustic communications
WO2018000914A1 (zh) 一种单载波信道估计方法
CN102790746B (zh) Ofdm系统的信道估计方法
Şahin et al. Spectrally efficient iterative MU-MIMO receiver for SC-FDMA based on EP
CN116707662A (zh) 一种mimo-fmt-tr水声通信系统及方法
Yang et al. Joint Multiple Turbo Equalization for Harsh Time-Varying Underwater Acoustic Channels
JP4463852B2 (ja) 通信路伝達関数を反復的に推定する装置及び方法
CN112039809B (zh) 基于混合软信息的块迭代均衡器及双向块迭代均衡器
Sharma et al. Channel Estimation and Equalization Using FIM for MIMO-OFDM on Doubly Selective Faded Noisy Channels
Jing et al. A Learned Denoising-Based Sparse Adaptive Channel Estimation for OTFS Underwater Acoustic Communications
CN118089741B (zh) 基于延迟多普勒域Turbo均衡的导航数据处理方法

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

Effective date of registration: 20231007

Address after: 1294 Chuangxin Road, Songbei District, Harbin City, Heilongjiang Province, 150001

Patentee after: HARBIN XINGUANG PHOTOELECTRIC TECHNOLOGY Co.,Ltd.

Address before: 150001 Intellectual Property Office, Harbin Engineering University science and technology office, 145 Nantong Avenue, Nangang District, Harbin, Heilongjiang

Patentee before: HARBIN ENGINEERING University