CN113746567A - 一种零值吸引pmser-df均衡系统及其实现方法 - Google Patents

一种零值吸引pmser-df均衡系统及其实现方法 Download PDF

Info

Publication number
CN113746567A
CN113746567A CN202110548453.8A CN202110548453A CN113746567A CN 113746567 A CN113746567 A CN 113746567A CN 202110548453 A CN202110548453 A CN 202110548453A CN 113746567 A CN113746567 A CN 113746567A
Authority
CN
China
Prior art keywords
filter
feedback
coefficient
sparsity
attraction
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
Application number
CN202110548453.8A
Other languages
English (en)
Other versions
CN113746567B (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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202110548453.8A priority Critical patent/CN113746567B/zh
Publication of CN113746567A publication Critical patent/CN113746567A/zh
Application granted granted Critical
Publication of CN113746567B publication Critical patent/CN113746567B/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/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/03305—Joint sequence estimation and interference removal
    • 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
    • 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
    • H04L2025/03592—Adaptation methods
    • H04L2025/03598—Algorithms
    • H04L2025/03611—Iterative algorithms

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Filters That Use Time-Delay Elements (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

本发明提供一种零值吸引PMSER‑DF均衡系统及其实现方法。所述系统包括均衡模块、误码指示器以及判决反馈单元；所述均衡模块包括滤波器、系数更新单元、步长更新单元、稀疏控制模块以及零值吸引模块。所述方法包括以下步骤：设置滤波器系数、稀疏对角矩阵及控制参数的初始值；计算滤波器输出信号；计算误码指示信号；计算滤波器稀疏度；计算零值吸引算子；更新稀疏对角矩阵；更新滤波器系数；重复上述步骤，直至均衡器收敛。本发明的特点是它充分利用信道稀疏性的先验信息，引入稀疏度控制的零值吸引算子，在稀疏信道下利用较少的训练序列便可达到优良的系统性能。

Description

一种零值吸引PMSER-DF均衡系统及其实现方法
技术领域
本发明涉及无线通信领域，尤其涉及一种零值吸引PMSER-DF均衡系统及其实现方法。
背景技术
近年来随着我国经济的高速发展，对海洋空间和海洋资源的开发越来越迫切。水声通信作为目前实现水下中、远程无线信息传输的唯一有效技术手段，在未来海洋工程领域中的战略意义日益凸显。然而，水声信道环境复杂，水平面反射、散射都会造成严重的多径干扰，实现长距离传输只能使用低频段音频，收发机的相对移动以及水流等影响都会造成严重的多普勒效应，相较于陆地几个符号周期的信道时延，水声通信中的信道时延可达到上百个，码间串扰非常严重，加之水声信道噪声丰富，噪声功率大，信噪比较低，为对抗严重的码间串扰，需要发展出一种在低信噪比条件下的快速有效的均衡算法。
尽管水声信道时延严重，多径时延可到上百个符号周期，但是经过多次反射、折射、衰减后，接收端接收到的信号是只有几条较强路径的信号的叠加，从而水声信道在本质上表现出稀疏性，根据迫零均衡器的设计准则，水声信道的稀疏性造成均衡器也表现出一定的稀疏性。
一种由华南理工大学提供的基于最小误码率准则的自适应信道均衡器及其实现方法(中国发明专利号：CN102916916A)，该专利的特点是它是直接基于最小误码率准则的自适应滤波器，利用误码指示器将滤波后的输出信号与导频信号映射成误码指示信号以调整均衡器的方向，此方法相较于基于最小均方误差准则的均衡器，在信道冲激响应比较短的条件下收敛速度快，但在长信道中性能相较于基于最小均方误差准则的均衡器要差。另一种由华南理工大学提供的一种成比例MSER自适应判决反馈均衡系统及其实现方法(中国发明专利号： CN108306837B),该专利的特点是在每次均衡器系数更新时分配成比例与滤波器抽头值大小的步长，使得系数绝对值大的抽头值分配的步长大，系数绝对值小的抽头值分配的步长小，同时是基于最小误码率准则推导的迭代公式，对于短稀疏信道和长稀疏信道收敛速度相较于基于最小均方误差准则的均衡器收敛速度快，但当均衡器收敛到后期时，那些小的抽头值占主要影响时，算法收敛速度骤降，影响均衡器整体性能。
在通信领域，误码率、算法复杂度和收敛速度是衡量一个均衡算法性能的主要评判准则，但最重要的时误码率。前面的专利中都采用了基于最小误码率准则的均衡算法达到了较低的误码率，但没有充分利用信道稀疏性的先验信息导致收敛速度不够快。
发明内容
针对现有自适应均衡算法的缺陷与不足，本发明公开一种零值吸引 PMSER-DF均衡系统及其实现方法，本发明利用信道稀疏性的先验信息，通过在每次迭代更新时加入零值吸引算子，使得迭代初期均衡器稀疏度较小时，选择那些幅度较小的抽头值进行收缩，当迭代到后期时，小抽头值对均衡算法的收敛性占主要影响时，此时滤波器稀疏度大，零值吸引算子骤降，从而在初期通过零值吸引算子加快均衡算法的收敛速度，后期均匀分配定步长加快后期收敛速度，利用较少的训练序列便可达到良好的系统性能，且适用于短稀疏信道和长稀疏信道。
本发明的目的至少通过如下技术方案之一实现。
一种零值吸引PMSER-DF均衡系统，包括均衡模块、误码指示器以及判决反馈单元；所述均衡模块包括滤波器、系数更新单元、步长更新单元、稀疏控制模块以及零值吸引模块；
所述滤波器用于对当前时刻k接收到的信号序列rk以及判决反馈信号序列
Figure RE-GDA0003308007400000021
进行滤波得到输出信号yk；
所述判决反馈单元将滤波器输出值yk进行判决，得到判决反馈信号序列
Figure RE-GDA0003308007400000022
并输入到均衡模块；
所述误码指示器用于将上一次均衡模块的输出信号yk映射成误码指示信号 Ik，作为均衡模块调整的方向；
所述稀疏控制模块用于对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行处理，得到前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k；
所述零值吸引模块用于对前向滤波器的滤波系数cf,k、反馈滤波器的滤波系数cb,k、前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k进行处理，得到前向滤波器零值吸引算子ρf,k和反馈滤波器零值吸引算子ρb,k；
所述步长更新单元用于计算滤波器抽头值相对大小，并根据前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k和滤波器的滤波系数大小分配步长值，得到前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k；
所述系数更新单元用于依据接收到的复信号rk、判决反馈符号序列
Figure RE-GDA0003308007400000023
误码指示信号Ik、前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k、前向滤波器零值吸引算子ρf,k、反馈滤波器零值吸引算子ρb,k对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行更新。
进一步地，所述滤波器包括前向滤波器和反馈滤波器，对当前时刻接收到的信号序列rk以及判决反馈信号序列
Figure RE-GDA0003308007400000031
进行滤波得到输出信号yk，具体如下：
Figure RE-GDA0003308007400000032
其中cf,k为前向滤波器的滤波系数，下标k为时隙下标，表示当前时刻，cb,k为反馈滤波器的滤波系数，T表示矩阵的转置。
进一步地，所述判决反馈单元将滤波器输出值yk进行判决，得到判决反馈信号序列
Figure RE-GDA0003308007400000033
并输入到均衡模块，具体如下：
Figure RE-GDA0003308007400000034
Figure RE-GDA0003308007400000035
Figure RE-GDA0003308007400000036
其中
Figure RE-GDA0003308007400000037
表示取复数的实部；
Figure RE-GDA0003308007400000038
表示取复数的虚部。
进一步地，误码指示器中，对应关系如下：
Figure RE-GDA0003308007400000039
其中，
Figure RE-GDA00033080074000000310
k为当前时刻，下标D表示均衡模块的输出信号相对于导频信号的延迟；β为用于控制映射关系的常数，用tan(βx)来近似sgn(x)，β的取值范围为[1,5]；dk-D表示导频中的期望信号；Ik用来指示符号错误用以调整均衡器的更新方向。
进一步地，稀疏控制模块中，对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行处理，得到前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k，实现如下：
Figure RE-GDA00033080074000000311
Figure RE-GDA0003308007400000041
其中，Nf和Nb分别为前向滤波器和反馈滤波器的抽头个数；||cf,k||1和||cf,k||2分别表示前向滤波器的L1范数和L2范数；||cb,k||1和||cb,k||2分别表示反馈滤波器的L1范数和L2范数。
进一步地，零值吸引模块中，对前向滤波器的滤波系数cf,k、反馈滤波器的滤波系数cb,k、前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k进行处理，得到前向滤波器零值吸引算子ρf,k和反馈滤波器零值吸引算子ρb,k，实现如下：
ρf,k＝(ρf,k(0),ρf,k(1),…,ρf,k(Nf-1))T；
ρb,k＝(ρb,k(0),ρb,k(1),…,ρb,k(Nb-1))T；
Figure RE-GDA0003308007400000042
Figure RE-GDA0003308007400000043
其中，a为稀疏度放大因子，为一常数，取值范围为(5,20)；b为滤波器抽头绝对值放大因子，为一常数，取值范围为(10,30)；l为零值吸引算子放大因子，为一常数，取值范围为(0.1,1)；|cf,k(n)|和|cb,k(n)|分别表示前向滤波器和反馈滤波器第n+1个抽头系数的绝对值；e是自然常数，取值为2.718；ρf,k(n)表示前向滤波器的第n+1个抽头值上零值吸引算子的大小，ρb,k(n)表示反馈滤波器的第 n+1个抽头值上零值吸引算子的大小。
进一步地，步长更新单元中，计算滤波器抽头值相对大小，并根据前向滤波器稀疏度sf,k、反馈滤波器稀疏度sb,k和滤波器的滤波系数大小分配步长值，得到前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k，具体如下：
Gf,k＝diag{gf,k(0),gf,k(1),…,gf,k(Nf-1)}；
Gb,k＝diag{gb,k(0),gb,k(1),…,gb,k(Nb-1)}；
Figure RE-GDA0003308007400000051
Figure RE-GDA0003308007400000052
其中，α是常数，取值范围为(-1,1)，gf,k(n)和gb,k(n)分别表示分配给前向滤波器和反馈滤波器第n+1个抽头的步长。
进一步地，系数更新单元中，依据接收到的复信号rk、判决反馈符号序列
Figure RE-GDA0003308007400000053
误码指示信号Ik、前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k、前向滤波器零值吸引算子ρf,k、反馈滤波器零值吸引算子ρb,k对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行更新，实现如下：
Figure RE-GDA0003308007400000054
Figure RE-GDA0003308007400000055
其中，μf和μb分别为前向滤波器和反馈滤波器的全局步长；上标H和*分别表示矢量的共轭转置和复数的共轭。
进一步地，接收信号rk为列矢量，rk的元素是从当前时刻开始，按时间递减排列。
一种零值吸引PMSER-DF均衡系统的实现方法，包括以下步骤：
步骤S1、设置前向滤波器和反馈滤波器的抽头值Nf和Nb；设置k＝0时，前向滤波器系数和反馈滤波器系数的初始值为
Figure RE-GDA0003308007400000056
Figure RE-GDA0003308007400000057
设置前向滤波器和反馈滤波器的稀疏对角矩阵的初始值Gf,1和Gb,1，设置控制参数μf、μb、α、a、b、l、λ、β和D的值；
步骤S2、利用当前时刻的滤波系数cf,k和cb,k对接收信号矢量rk进行滤波产生滤波输出信号yk；
步骤S3、利用输出信号yk与导频期望信号dk-D得到误码指示信号Ik；
步骤S4、稀疏控制模块利用当前时刻的滤波系数cf,k、cb,k得到滤波器稀疏度sf,k、sb,k；
步骤S5、依据控制参数α、Nf和Nb，当前滤波系数cf,k和cb,k以及稀疏度sf,k和sb,k将稀疏对角矩阵Gf,k-1和Gb,k-1更新为Gf,k和Gb,k；
步骤S6、利用控制参数a、b和l，当前滤波器系数cf,k和cb,k以及稀疏度sf,k和 sb,k得到零值吸引算子ρf,k和ρb,k；
步骤S7、根据前向滤波器和反馈滤波器的全局步长μf和μb，误码指示信号 Ik，稀疏对角矩阵Gf,k和Gb,k以及零值吸引算子ρf,k和ρb,k分别将前向滤波器系数 cf,k和反馈滤波器系数cb,k更新为cf,k+1和cb,k+1；
步骤S8、重复步骤S2～S7,直至均衡器系数收敛，即||cf,k+1-cf,k||≤ε和 ||cb,k+1-cb,k||≤ε。
本发明相对于现有技术具有如下的优点和效果：
1)本发明利用水声信道稀疏性的先验信息，利用基于最小误码率准则下次梯度投影的方法获得的迭代公式，加入零值吸引算子，前期收敛速度快，利用较少的训练序列便可达到较低的误码率。
2)本发明零值吸引算子中加入了滤波器的稀疏度，利用每次迭代中滤波器的稀疏度控制零值吸引算子的大小，使得在初期滤波器稀疏度较小时零值吸引算子较大，选择那些幅度较小的抽头进行收缩，收敛到后期时滤波器稀疏度较大，零值吸引算子骤降，逐渐失去作用，避免后期小抽头值占主要作用时零值吸引算子对算法收敛性的影响。
附图说明
图1为常见的信道均衡器结构框图；
图2为本发明实现零值吸引PMSER-DF均衡系统的示意图；
图3为本发明在使用4-QAM调制下在matlab平台下仿真时所使用的信道的冲击响应示意图；其中图3a为水声信道脉冲响应的实部示意图；图3b为水声信道脉冲响应的虚部示意图；
图4是几种自适应均衡算法与本发明的零值吸引PMSER-DF均衡系统在不同调制方式下收敛性能的比较结果示意图，其中图4a为在使用BPSK调制时使用图3a信道时的比较结果示意图；图4b为在使用4-QAM调制时使用图3a-图4b 信道时的比较结果示意图。
具体实施方式
以下结合附图和实例对本发明的具体实施作进一步说明，但本发明的实施方式不限于此。
实施例：
一种零值吸引PMSER-DF均衡系统，如图2所示，包括均衡模块、误码指示器以及判决反馈单元；所述均衡模块包括滤波器、系数更新单元、步长更新单元、稀疏控制模块以及零值吸引模块；
所述滤波器包括前向滤波器和反馈滤波器，用于对当前时刻接收到的信号序列rk以及判决反馈信号序列
Figure RE-GDA0003308007400000071
进行滤波得到输出信号yk，具体如下：
Figure RE-GDA0003308007400000072
其中cf,k为前向滤波器的滤波系数，cb,k为反馈滤波器的滤波系数，T表示矩阵的转置；k表示时隙，代表当前时刻；接收信号rk为列矢量，rk的元素是从当前时刻开始，按时间递减排列；
所述判决反馈单元将滤波器输出值yk进行判决，得到判决反馈信号序列
Figure RE-GDA0003308007400000073
并输入到均衡模块，具体如下：
Figure RE-GDA0003308007400000074
Figure RE-GDA0003308007400000075
Figure RE-GDA0003308007400000076
其中
Figure RE-GDA00033080074000000710
表示取复数的实部；
Figure RE-GDA0003308007400000077
表示取复数的虚部。
所述误码指示器用于将上一次均衡模块的输出信号yk映射成误码指示信号 Ik，作为均衡模块调整的方向，对应关系如下：
Figure RE-GDA0003308007400000078
其中，
Figure RE-GDA0003308007400000079
k为当前时刻，下标D表示均衡模块的输出信号相对于导频信号的延迟；β为用于控制映射关系的常数，用tan(βx)来近似sgn(x)，β的取值范围为[1,5]；dk-D表示导频中的期望信号；Ik用来指示符号错误用以调整前向滤波器和反馈滤波器的更新方向。
所述稀疏控制模块用于对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行处理，得到前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k，实现如下：
Figure RE-GDA0003308007400000081
Figure RE-GDA0003308007400000082
其中，Nf和Nb分别为前向滤波器和反馈滤波器的抽头个数；||cf,k||1和||cf,k||2分别表示前向滤波器的L1范数和L2范数；||cb,k||1和||cb,k||2分别表示反馈滤波器的L1范数和L2范数；
所述零值吸引模块用于对前向滤波器的滤波系数cf,k、反馈滤波器的滤波系数cb,k、前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k进行处理，得到前向滤波器零值吸引算子ρf,k和反馈滤波器零值吸引算子ρb,k，实现如下：
ρf,k＝(ρf,k(0),ρf,k(1),…,ρf,k(Nf-1))T；
ρb,k＝(ρb,k(0),ρb,k(1),…,ρb,k(Nb-1))T；
Figure RE-GDA0003308007400000083
Figure RE-GDA0003308007400000084
其中，a为稀疏度放大因子，为一常数，取值范围为(5,20)；b为滤波器抽头绝对值放大因子，为一常数，取值范围为(10,30)；l为零值吸引算子放大因子，为一常数，取值范围为(0.1,1)；|cf,k(n)|和|cb,k(n)|分别表示前向滤波器和反馈滤波器第n+1个抽头系数的绝对值；e是自然常数，取值为2.718；ρf,k(n)表示前向滤波器的第n+1个抽头上零值吸引算子的大小，ρb,k(n)表示反馈滤波器的第n+1 个抽头上零值吸引算子的大小。
所述步长更新单元用于计算滤波器抽头系数相对大小，并根据滤波器稀疏度sf,k、sb,k和滤波器的系数大小分配步长值，得到前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k，具体如下：
Gf,k＝diag{gf,k(0),gf,k(1),…,gf,k(Nf-1)}；
Gb,k＝diag{gb,k(0),gb,k(1),…,gb,k(Nb-1)}；
Figure RE-GDA0003308007400000091
Figure RE-GDA0003308007400000092
其中，α是常数，取值范围为(-1,1)，gf,k(n)和gb,k(n)分别表示分配给前向滤波器和反馈滤波器第n+1个抽头的步长；
所述系数更新单元用于依据接收到的复信号rk、判决反馈符号序列
Figure RE-GDA0003308007400000093
误码指示信号Ik、前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k、前向滤波器零值吸引算子ρf,k、反馈滤波器零值吸引算子ρb,k对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行更新，实现如下：
Figure RE-GDA0003308007400000094
Figure RE-GDA0003308007400000095
其中，μf和μb分别为前向滤波器和反馈滤波器的全局步长；上标H和*分别表示矢量的共轭转置和复数的共轭。
一种零值吸引PMSER-DF均衡系统的实现方法，包括以下步骤：
步骤S1、设置前向滤波器和反馈滤波器的抽头值Nf和Nb；设置k＝0时，前向滤波器系数和反馈滤波器系数的初始值为
Figure RE-GDA0003308007400000096
Figure RE-GDA0003308007400000097
本实施例中，设置前向滤波器和反馈滤波器的稀疏对角矩阵的初始值Gf,1和Gb,1为零矩阵；设置控制参数μf、μb、α、a、b、l、λ、β和D的值；
步骤S2、利用当前时刻的滤波系数cf,k和cb,k对接收信号矢量rk进行滤波产生滤波输出信号yk；
步骤S3、利用输出信号yk与导频期望信号dk-D得到误码指示信号Ik；
步骤S4、稀疏控制模块利用当前时刻的滤波系数cf,k、cb,k得到滤波器稀疏度sf,k、sb,k；
步骤S5、依据控制参数α、Nf和Nb，当前滤波系数cf,k和cb,k以及稀疏度sf,k和sb,k将稀疏对角矩阵Gf,k-1和Gb,k-1更新为Gf,k和Gb,k；
步骤S6、利用控制参数a、b和l，当前滤波器系数cf,k和cb,k以及稀疏度sf,k和 sb,k得到零值吸引算子ρf,k和ρb,k；
步骤S7、根据前向滤波器和反馈滤波器的全局步长μf和μb，误码指示信号 Ik，稀疏对角矩阵Gf,k和Gb,k以及零值吸引算子ρf,k和ρb,k分别将前向滤波器系数 cf,k和反馈滤波器系数cb,k更新为cf,k+1和cb,k+1；
步骤S8、重复步骤S2～S7,直至均衡器系数收敛，即||cf,k+1-cf,k||≤ε和 ||cb,k+1-cb,k||≤ε。
本实施例中，如图1所示为常见的信道均衡器结构框图，其中dk为输入信道中的导频信号，hk是长度为L的信道冲激响应，nk是均值为0，方差为
Figure RE-GDA0003308007400000101
的加性高斯噪声。
其中k时刻接收信号rk可表示为发送的导频信号与信道冲激响应的卷积：
Figure RE-GDA0003308007400000102
均衡模块输入信号矢量为：
rk＝[rk,rk-1,…,rk-Nf+1]T；
前向滤波器系数矢量为：
cf,k＝[cf,k(0),cf,k(1),…,cf,k(Nf-1)]T；
反馈滤波器系数矢量为：
cb,k＝[cb,k(0),cb,k(1),…,cb,k(Nb-1)]T；
判决反馈信号矢量为：
Figure RE-GDA0003308007400000103
当发射端采用QPSK调制方式时，在最下误码率准则下，可建立如下最优化问题：
Figure RE-GDA0003308007400000111
其中
Figure RE-GDA0003308007400000112
表示取复数的实部，
Figure RE-GDA0003308007400000113
表示取复数的虚部；
分别利用接收信号的实部和虚部求解滤波器系数的迭代公式，首先利用滤波器输出的实部：
Figure RE-GDA0003308007400000114
其中
Figure RE-GDA0003308007400000115
cf,k,R表示前向滤波器系数的实部向量，cf,k,I表示前向滤波器系数的虚部向量；
Figure RE-GDA0003308007400000116
cb,k,R表示反馈滤波器系数的实部向量， cb,k,I表示反馈滤波器系数的虚部向量。
滤波器输出结果的实部和虚部互相独立，可将最优化问题重写为：
Figure RE-GDA0003308007400000117
Figure RE-GDA0003308007400000118
构造拉格朗日目标函数：
Figure RE-GDA0003308007400000119
所述目标函数分别对
Figure RE-GDA00033080074000001110
求偏导并令偏导数为零，可得如下迭代公式：
Figure RE-GDA00033080074000001111
Figure RE-GDA00033080074000001112
同理利用均衡器输出的虚部可得如下迭代公式：
Figure RE-GDA0003308007400000121
Figure RE-GDA0003308007400000122
其中：
Figure RE-GDA0003308007400000123
Figure RE-GDA0003308007400000124
Figure RE-GDA0003308007400000125
Figure RE-GDA0003308007400000126
将实部与虚部进行合并可得QPSK迭代公式：
Figure RE-GDA0003308007400000127
Figure RE-GDA0003308007400000128
上面是一种成比例MSER自适应判决反馈算法，该方法在基于最小误码率准则下推导迭代公式，且每次迭代更新时分配成比例与抽头系数绝对值大小的步长，初期收敛速度相对较快，后期当小抽头值占主要影响时，收敛速度慢，在每次更新迭代时没有充分利用信道稀疏性的先验信息。本发明提供一种稀疏控制零值吸引成比例最小误码率判决反馈均衡算法，在以上方法的基础上通过加入稀疏度控制的零值吸引算子在初期进一步加快收敛速度，其中k时刻前向滤波器和反馈滤波器稀疏度计算公式为：
Figure RE-GDA0003308007400000129
Figure RE-GDA0003308007400000131
k时刻零值吸引算子向量计算公式为：
ρf,k＝(ρf,k(0),ρf,k(1),…,ρf,k(Nf-1))T；
ρb,k＝(ρb,k(0),ρb,k(1),…,ρb,k(Nb-1))T；
Figure RE-GDA0003308007400000132
Figure RE-GDA0003308007400000133
其中，a为稀疏度放大因子，为一常数，取值范围为(5,20)；b为滤波器抽头绝对值放大因子，为一常数，取值范围为(10,30)；l为零值吸引算子放大因子，为一常数，取值范围为(0.1,1)；|cf,k(n)|和|cb,k(n)|分别表示前向滤波器和反馈滤波器第n+1个抽头系数的绝对值；e是自然常数，取值为2.718；
前向滤波器和反馈滤波器系数更新公式为：
Figure RE-GDA0003308007400000134
Figure RE-GDA0003308007400000135
上述方法为零值吸引PMSER-DF均衡算法(ZA-PMSER-DFE)，具体实施方式如图2所示。
如图3所示，在使用BPSK调制时所使用的信道具有较强的稀疏性。如图3a 和图3b所示，在使用4-QAM调制时，所使用的信道脉冲响应在实部和虚部均具有较强的稀疏性。在附图4a和4b中，通过matlab仿真给出了在BPSK和 4-QAM调制方式及信噪比SNR＝15dB时几种自适应均衡算法与本发明的稀疏控制零值吸引成比例最小误码率判决反馈均衡算法(SZA-PMSER-DFE)误码率的比较结果，其中β取值为1，α取值为0.5，训练序列长度为2000，接着发送数据符号，均衡器步长取值为相同值，从图中可以看出，无论是BPSK还是4-QAM 调制方式,SZA-PMSER-DFE算法在误码率和收敛速度上均显著优于其他三种算法。
本发明在成比例MSER自适应判决反馈算法的基础上加入了稀疏度控制的零值吸引算子，充分利用了信道稀疏性的先验信息，算法迭代初期，滤波器系数稀疏度较小，零值吸引算子作用于那些系数绝对值较小的抽头值进行收缩，算法收敛到后期小抽头值对算法的收敛性占主要影响时，此时滤波器的稀疏度较大，反比例与稀疏度指数的零值吸引算子骤降，从而后期消除零值吸引算子对算法收敛性的影响，因此在误码率与收敛速度上显著优于成比例MSER自适应判决反馈算法。

Claims (10)

1.一种零值吸引PMSER-DF均衡系统，其特征在于，包括均衡模块、误码指示器以及判决反馈单元；所述均衡模块包括滤波器、系数更新单元、步长更新单元、稀疏控制模块以及零值吸引模块；
所述滤波器用于对当前时刻k接收到的信号序列rk以及判决反馈信号序列dk进行滤波得到输出信号yk；
所述判决反馈单元将滤波器输出值yk进行判决，得到判决反馈信号序列dk并输入到均衡模块；
所述误码指示器用于将上一次均衡模块的输出信号yk映射成误码指示信号Ik，作为均衡模块调整的方向；
所述稀疏控制模块用于对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行处理，得到前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k；
所述零值吸引模块用于对前向滤波器的滤波系数cf,k、反馈滤波器的滤波系数cb,k、前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k进行处理，得到前向滤波器零值吸引算子ρf,k和反馈滤波器零值吸引算子ρb,k；
所述步长更新单元用于计算滤波器抽头值相对大小，并根据前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k和滤波器的滤波系数大小分配步长值，得到前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k；
所述系数更新单元用于依据接收到的复信号rk、判决反馈符号序列dk、误码指示信号Ik、前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k、前向滤波器零值吸引算子ρf,k、反馈滤波器零值吸引算子ρb,k对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行更新。
2.根据权利要求1所述的一种零值吸引PMSER-DF均衡系统，其特征在于，所述滤波器包括前向滤波器和反馈滤波器，对当前时刻接收到的信号序列rk以及判决反馈信号序列dk进行滤波得到输出信号yk，具体如下：
Figure FDA0003074381180000011
其中cf,k为前向滤波器的滤波系数，下标k为时隙下标，表示当前时刻，cb,k为反馈滤波器的滤波系数，T表示矩阵的转置。
3.根据权利要求2所述的一种零值吸引PMSER-DF均衡系统，其特征在于，所述判决反馈单元将滤波器输出值yk进行判决，得到判决反馈信号序列dk并输入到均衡模块，具体如下：
Figure FDA0003074381180000021
Figure FDA0003074381180000022
Figure FDA0003074381180000023
其中
Figure FDA0003074381180000024
表示取复数的实部；
Figure FDA0003074381180000025
表示取复数的虚部。
4.根据权利要求3所述的一种零值吸引PMSER-DF均衡系统，其特征在于，误码指示器中，对应关系如下：
Figure FDA0003074381180000026
其中，
Figure FDA0003074381180000027
Figure FDA0003074381180000028
k为当前时刻，下标D表示均衡模块的输出信号相对于导频信号的延迟；β为用于控制映射关系的常数，用tan(βx)来近似sgn(x)，β的取值范围为[1,5]；dk-D表示导频中的期望信号；Ik用来指示符号错误用以调整均衡器的更新方向。
5.根据权利要求4所述的一种零值吸引PMSER-DF均衡系统，其特征在于，稀疏控制模块中，对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行处理，得到前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k，实现如下：
Figure FDA0003074381180000029
Figure FDA00030743811800000210
其中，Nf和Nb分别为前向滤波器和反馈滤波器的抽头个数；||cf,k||1和||cf,k||2分别表示前向滤波器的L1范数和L2范数；||cb,k||1和||cb,k||2分别表示反馈滤波器的L1范数和L2范数。
6.根据权利要求5所述的一种零值吸引PMSER-DF均衡系统，其特征在于，零值吸引模块中，对前向滤波器的滤波系数cf,k、反馈滤波器的滤波系数cb,k、前向滤波器稀疏度sf,k和反馈滤波器稀疏度sb,k进行处理，得到前向滤波器零值吸引算子ρf,k和反馈滤波器零值吸引算子ρb,k，实现如下：
ρf,k＝(ρf,k(0),ρf,k(1),…,ρf,k(Nf-1))T；
ρb,k＝(ρb,k(0),ρb,k(1),…,ρb,k(Nb-1))T；
Figure FDA0003074381180000031
Figure FDA0003074381180000032
其中，a为稀疏度放大因子，为一常数，取值范围为(5,20)；b为滤波器抽头绝对值放大因子，为一常数，取值范围为(10,30)；l为零值吸引算子放大因子，为一常数，取值范围为(0.1,1)；|cf,k(n)|和|cb,k(n)|分别表示前向滤波器和反馈滤波器第n+1个抽头系数的绝对值；e是自然常数，取值为2.718；ρf,k(n)表示前向滤波器的第n+1个抽头值上零值吸引算子的大小，ρb,k(n)表示反馈滤波器的第n+1个抽头值上零值吸引算子的大小。
7.根据权利要求6所述的一种零值吸引PMSER-DF均衡系统，其特征在于，步长更新单元中，计算滤波器抽头值相对大小，并根据前向滤波器稀疏度sf,k、反馈滤波器稀疏度sb,k和滤波器的滤波系数大小分配步长值，得到前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k，具体如下：
Gf,k＝diag{gf,k(0),gf,k(1),…,gf,k(Nf-1)}；
Gb,k＝diag{gb,k(0),gb,k(1),…,gb,k(Nb-1)}；
Figure FDA0003074381180000033
Figure FDA0003074381180000034
其中，α是常数，取值范围为(-1,1)，gf,k(n)和gb,k(n)分别表示分配给前向滤波器和反馈滤波器第n+1个抽头的步长。
8.根据权利要求7所述的一种零值吸引PMSER-DF均衡系统，其特征在于，系数更新单元中，依据接收到的复信号rk、判决反馈符号序列dk、误码指示信号Ik、前向滤波器稀疏对角矩阵Gf,k和反馈滤波器稀疏对角矩阵Gb,k、前向滤波器零值吸引算子ρf,k、反馈滤波器零值吸引算子ρb,k对前向滤波器的滤波系数cf,k和反馈滤波器的滤波系数cb,k进行更新，实现如下：
Figure FDA0003074381180000041
Figure FDA0003074381180000042
其中，μf和μb分别为前向滤波器和反馈滤波器的全局步长；上标H和*分别表示矢量的共轭转置和复数的共轭。
9.根据权利要求1所述的一种零值吸引PMSER-DF均衡系统，其特征在于，接收信号rk为列矢量，rk的元素是从当前时刻开始，按时间递减排列。
10.一种零值吸引PMSER-DF均衡系统的实现方法，其特征在于，包括以下步骤：
步骤S1、设置前向滤波器和反馈滤波器的抽头值Nf和Nb；设置k＝0时，前向滤波器系数和反馈滤波器系数的初始值为
Figure FDA0003074381180000043
Figure FDA0003074381180000044
设置前向滤波器和反馈滤波器的稀疏对角矩阵的初始值Gf,1和Gb,1，设置控制参数μf、μb、α、a、b、l、λ、β和D的值；
步骤S2、利用当前时刻的滤波系数cf,k和cb,k对接收信号矢量rk进行滤波产生滤波输出信号yk；
步骤S3、利用输出信号yk与导频期望信号dk-D得到误码指示信号Ik；
步骤S4、稀疏控制模块利用当前时刻的滤波系数cf,k、cb,k得到滤波器稀疏度sf,k、sb,k；
步骤S5、依据控制参数α、Nf和Nb，当前滤波系数cf,k和cb,k以及稀疏度sf,k和sb,k将稀疏对角矩阵Gf,k-1和Gb,k-1更新为Gf,k和Gb,k；
步骤S6、利用控制参数a、b和l，当前滤波器系数cf,k和cb,k以及稀疏度sf,k和sb,k得到零值吸引算子ρf,k和ρb,k；
步骤S7、根据前向滤波器和反馈滤波器的全局步长μf和μb，误码指示信号Ik，稀疏对角矩阵Gf,k和Gb,k以及零值吸引算子ρf,k和ρb,k分别将前向滤波器系数cf,k和反馈滤波器系数cb,k更新为cf,k+1和cb,k+1；
步骤S8、重复步骤S2～S7,直至均衡器系数收敛，即||cf,k+1-cf,k||≤ε和||cb,k+1-cb,k||≤ε。
CN202110548453.8A 2021-05-19 2021-05-19 一种零值吸引pmser-df均衡系统及其实现方法 Active CN113746567B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110548453.8A CN113746567B (zh) 2021-05-19 2021-05-19 一种零值吸引pmser-df均衡系统及其实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110548453.8A CN113746567B (zh) 2021-05-19 2021-05-19 一种零值吸引pmser-df均衡系统及其实现方法

Publications (2)

Publication Number Publication Date
CN113746567A true CN113746567A (zh) 2021-12-03
CN113746567B CN113746567B (zh) 2022-10-25

Family

ID=78728329

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110548453.8A Active CN113746567B (zh) 2021-05-19 2021-05-19 一种零值吸引pmser-df均衡系统及其实现方法

Country Status (1)

Country Link
CN (1) CN113746567B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007049841A1 (en) * 2005-10-27 2007-05-03 Electronics And Telecommunications Research Institute Apparatus for equalizing decision feedback based on channel properties change and method therefor
CN103825852A (zh) * 2014-01-28 2014-05-28 华南理工大学 一种双模自适应判决反馈均衡模块及其实现方法
CN103957176A (zh) * 2014-04-30 2014-07-30 华南理工大学 一种自适应rls判决反馈均衡系统及其实现方法
CN108306837A (zh) * 2018-02-09 2018-07-20 华南理工大学 一种成比例mser自适应判决反馈均衡系统及其实现方法
CN111030758A (zh) * 2019-12-16 2020-04-17 哈尔滨工业大学（威海） 带有稀疏约束的自适应零吸引因子盲判决反馈均衡算法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007049841A1 (en) * 2005-10-27 2007-05-03 Electronics And Telecommunications Research Institute Apparatus for equalizing decision feedback based on channel properties change and method therefor
CN103825852A (zh) * 2014-01-28 2014-05-28 华南理工大学 一种双模自适应判决反馈均衡模块及其实现方法
CN103957176A (zh) * 2014-04-30 2014-07-30 华南理工大学 一种自适应rls判决反馈均衡系统及其实现方法
CN108306837A (zh) * 2018-02-09 2018-07-20 华南理工大学 一种成比例mser自适应判决反馈均衡系统及其实现方法
CN111030758A (zh) * 2019-12-16 2020-04-17 哈尔滨工业大学（威海） 带有稀疏约束的自适应零吸引因子盲判决反馈均衡算法

Also Published As

Publication number Publication date
CN113746567B (zh) 2022-10-25

Similar Documents

Publication Publication Date Title
CN103957176B (zh) 一种自适应rls判决反馈均衡系统及其实现方法
CN102123115B (zh) 基于粒子群优化的正交小波盲均衡方法
CN108306837B (zh) 一种成比例mser自适应判决反馈均衡系统及其实现方法
CN107124379B (zh) 一种基于改进狼群优化的正交小波常模盲均衡方法
CN112929308B (zh) 基于最大相关熵和分数阶常模盲均衡的信道畸变消除方法
US20050232347A1 (en) Apparatus and method for noise enhancement reduction in an adaptive equalizer
CN105656823A (zh) 基于最小误码率准则的水下通信Turbo接收系统及方法
CN106411799A (zh) 低轨卫星移动通信单载波频域均衡方法
CN114095320B (zh) 基于动量分数阶多模盲均衡算法的信道均衡方法
CN109818888B (zh) 一种脉冲干扰环境下的组稀疏水声信道估计方法
CN113746567B (zh) 一种零值吸引pmser-df均衡系统及其实现方法
CN1845539A (zh) 一种具有重叠结构的时域自适应均衡器
CN103825852A (zh) 一种双模自适应判决反馈均衡模块及其实现方法
CN107070825B (zh) 一种基于模拟退火狼群优化的小波加权多模盲均衡方法
CN116016061B (zh) 高机动平台短波双选信道双迭代Turbo均衡方法
CN101494625A (zh) 一种线性均衡方法及线性均衡器
CN100364257C (zh) 和nr码结合的时域自适应均衡器
CN102137052B (zh) 一种基于梯度向量的变步长最小均方信道均衡方法
CN112039809B (zh) 基于混合软信息的块迭代均衡器及双向块迭代均衡器
CN107018104B (zh) 一种基于混合猴群优化的小波加权多模盲均衡方法
Feng et al. Research on Single Carrier Frequency Domain Equalization in Wideband High Frequency Communication System
CN103841067A (zh) 一种水声信道通信信号的均衡方法
CN102143104B (zh) 具有双重叠结构的时域自适应判决反馈均衡器
CN115225437B (zh) 用于水声协作通信的联合智能均衡与译码方法
CN115208481B (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