CN109379308B - 一种基于稀疏感知的自适应空时均衡方法 - Google Patents

一种基于稀疏感知的自适应空时均衡方法 Download PDF

Info

Publication number
CN109379308B
CN109379308B CN201811519562.1A CN201811519562A CN109379308B CN 109379308 B CN109379308 B CN 109379308B CN 201811519562 A CN201811519562 A CN 201811519562A CN 109379308 B CN109379308 B CN 109379308B
Authority
CN
China
Prior art keywords
tap
zero
time
cross
pilot
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
CN201811519562.1A
Other languages
English (en)
Other versions
CN109379308A (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201811519562.1A priority Critical patent/CN109379308B/zh
Publication of CN109379308A publication Critical patent/CN109379308A/zh
Application granted granted Critical
Publication of CN109379308B publication Critical patent/CN109379308B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03012Arrangements for removing intersymbol interference operating in the time domain
    • H04L25/03019Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03891Spatial equalizers

Landscapes

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

Abstract

一种基于稀疏感知的自适应空时均衡方法,在发送端将待发射的数据分块,然后在每个数据块的末尾加上独特字;发送端发送信号到达接收端并在接收端经过定时同步预处理后,截取接收到的导频;然后利用贪婪算法或者互相关矢量算法来确定时域滤波器的非零抽头位置,再计算选中位置的接收导频块的自相关矩阵与互相关矢量,从而计算自适应权值;接收端将接收到的数据进行空间合并,再进行去独特字操作,实现信道的均衡。本发明提出的方法能够有效降低多径效应和频率选择性衰落对无线通信过程的影响,通过合理的时域滤波器稀疏抽头设置,有效改善无线通信系统的性能,对抗信道的时延扩展。

Description

一种基于稀疏感知的自适应空时均衡方法
技术领域
本发明属于无线通信领域,涉及一种基于稀疏感知的自适应空时均衡方法。
背景技术
无线信号传播过程中,信号传播环境复杂,由于电磁波的传播路径不同,使得信号能量不断消耗,信号幅度、相位及时间不断变化,引起原始信号失真或者出错,造成多径效应。此外,电磁波的传播除了直射波、反射波以外,不可避免地会受到障碍物、其他无线信号等环境因素的影响,使得信号幅度迅速下降。当信号带宽大于相干带宽时,还会引起频率选择性衰落。因此,无线信道相对有线信道而言,具有很大的随机性,信道估计是实现无线通信的一项关键技术,信道估计的精度将直接影响到整个系统的性能。
离散采样信道大部分的采样值都为零,只有极少数采样点的值非零,通常称这样的信道为稀疏信道。随着研究的不断深入,大量文献表明,无线信道本身往往呈现出稀疏性。当信道具有稀疏性时,时域滤波的抽头也应相应具有稀疏性。一方面,信道稀疏时,冗余的抽头分量带来的可能都是无用的符号间干扰,对提升系统性能改善不大。另一方面,在给定非零抽头个数情况下,时域滤波可以通过合理的稀疏抽头设置对抗信道的大时延扩展。
发明内容
本发明的目的在于克服现有技术的缺点,提供一种基于稀疏感知的自适应空时均衡方法。
为达到上述目的,本发明采用以下技术方案予以实现:
一种基于稀疏感知的自适应空时均衡方法,包括以下步骤:
步骤一:在发送端将待发射的数据分块,然后在每个数据块的末尾加上独特字;
步骤二:发送端发送信号到达接收端并在接收端经过定时同步预处理后,截取接收到的导频;然后在空时均衡的过程中采用贪婪算法或互相关矢量算法自适应地寻找最优时域滤波器非零抽头位置,然后计算自适应权值;
步骤三:接收端将接收到的数据进行分块,再利用步骤二得到的自适应权值对每根天线上的接收数据块进行合并,并进行去独特字操作,从而实现信道的均衡。
本发明进一步的改进在于,步骤二中采用贪婪算法自适应地寻找最优时域滤波器非零抽头位置的具体过程如下:
假设接收机有M根天线,Δ为时域滤波器可能的最大阶数,Q为限定的时域非零抽头系数的个数,Q≤Δ+1;接收机多天线集中放置,多天线多径信道的时延特征相同,空时滤波时,多天线时域滤波非零抽头的位置相同;
采用N表示计算空时权值的导频长度,L表示数据信号的长度,定义以下符号:
s:N×1维理想导频信号;
r:L×1维的数据信号部分;
xk,m:第m根天线上延迟k位的N×1维接收导频信号,m=1,…,M;k=0,…,Δ;
zk,m:第m根天线上延迟k位的L×1维接收数据信号,m=1,…,M;k=0,…,Δ;
Xk:所有天线上延迟k位的N×M维接收导频时间扩展矩阵,Xk=[xk,1,xk,2,…,xk,M];
Zk:所有天线上延迟k位的L×M维接收数据时间扩展矩阵,Zk=[zk,1,zk,2,…,zk,M];
wq:第q个滤波系数矢量,由所有天线的第q个滤波系数构成,M×1维;
①采用Π={0,1,2,…,Δ}来表示待选择的抽头位置集合,采用Ω={}表示已经选择的非零抽头位置,
Figure BDA0001902877800000021
表示已经选择的非零抽头位置对应的导频块;
②从时域滤波器非零抽头系数只有一个的情况开始,这时待选择的抽头共有T=Δ+1个,基于最小均方误差准则准则,将这T=Δ+1个可能的抽头遍历,从中选择一个,即
Figure BDA0001902877800000022
其中,
q从0遍历到Δ,对应的接收导频块为X0,X1,……,XΔ,其中
Figure BDA0001902877800000031
为权值测试矢量,由下式求得:
Figure BDA0001902877800000032
其中,自相关矩阵
Figure BDA0001902877800000033
互相关矩阵
Figure BDA0001902877800000034
③遍历之后选择第一个非零抽头,然后从Π中剔除掉选中的
Figure BDA0001902877800000035
将其并入到Ω中,
Figure BDA0001902877800000036
中并入第
Figure BDA0001902877800000037
个抽头位置对应的导频块
Figure BDA0001902877800000038
即:
Figure BDA0001902877800000039
④在确定了第一个抽头位置之后,再遵循最小均方误差准则,从剩余的T-1位置中选择下一个非零抽头,前两个非零位置确定后,再类似从剩余的T-2个位置中挑选出一个,依次类推,直到挑选出所有Q个非零抽头位置。
本发明进一步的改进在于,步骤二中采用互相关矢量算法自适应地寻找最优时域滤波器非零抽头位置的具体过程如下:
①接收端收到空时接收导频扩展矩阵Y=[X0 X1 … XΔ]之后,计算其互相关矢量p=YHs;其中,s为N×1维理想导频信号;
②将互相关矢量p表示为
Figure BDA00019028778000000310
其中
Figure BDA00019028778000000311
表示对应第q个位置的接收导频信号互相关矢量;
③计算||pq||2,q=0,1,…,Δ,从这Δ个值中挑选出最大的Q个值,即
Figure BDA00019028778000000312
将其对应的索引作为需要选择的Q个非零抽头位置。
本发明进一步的改进在于,计算自适应权值的具体过程为:根据贪婪算法或互相关矢量算法得到的非零抽头位置,确定对应的接收端导频块,再利用最小均方误差准则准则,求出维纳解,得到自适应权值
Figure BDA00019028778000000313
Figure BDA00019028778000000314
其中,R为自相关矩阵,p为互相关矢量。
本发明进一步的改进在于,采用贪婪算法选择的抽头位置集合为Ω,选中的接收导频块集合为
Figure BDA0001902877800000041
则自相关矢量
Figure BDA0001902877800000042
互相关矢量
Figure BDA0001902877800000043
本发明进一步的改进在于,步骤三的具体过程如下:
利用步骤二得到的自适应权值
Figure BDA0001902877800000044
对各数据块进行均衡:
Figure BDA0001902877800000045
其中,
Figure BDA0001902877800000046
为接收之后恢复出来的数据部分,
Figure BDA0001902877800000047
为选中抽头位置对应的接收数据块,去掉恢复之后各数据块末尾的独特字,完成信道的均衡。
与现有技术相比,本发明具有以下有益效果:
本发明中在发送端将待发射的数据分块,然后在每个数据块的末尾加上独特字;发送端发送信号到达接收端并在接收端经过定时同步预处理后,截取接收到的导频;然后利用贪婪算法或者互相关矢量算法来确定时域滤波器的非零抽头位置,再计算选中位置的接收导频块的自相关矩阵与互相关矢量,从而计算自适应权值;接收端将接收到的数据进行空间合并,再进行去独特字操作,实现信道的均衡。本发明提出的方法能够有效降低多径效应和频率选择性衰落对无线通信过程的影响,通过合理的时域滤波器稀疏抽头设置,有效改善无线通信系统的性能,对抗信道的时延扩展。与传统空时均衡方法相比,在相同信道时延扩展情况下,本发明能够有效减少时域滤波的非零抽头个数,从而极大地降低了复杂度。在复杂度给定的条件下,本发明通过自适应优化时域滤波器非零抽头位置,能够对抗更大的信道时延扩展;在时域抽头位置确定后,本发明的空时权值依据最小均方误差(MMSE)准则获得。因此,本发明具备一定抗干扰能力,有效提升无线通信链路的顽存性。
进一步的,本发明提出的贪婪算法和互相关矢量算法,从时域滤波器的所有抽头中借助相关参数选取最为匹配的Q个抽头,避免了传统时空均衡方法直接选择前Q个抽头时剩余抽头带来的干扰。
附图说明
图1为本发明采用的系统模型框图;
图2为本发明采用的信号帧结构图;
图3为本发明采用的空时均衡处理的过程示意图;
图4为本发明同对比方案在不同非零抽头个数且无干扰情况下的BER性能示意图;
图5为本发明在有干扰情况下的BER性能示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
本发明的方法包括以下步骤:
步骤一:发射端数据分块加独特字(Unique word,记为:UW):在发送端将待发射的数据分块,然后在每个数据块的末尾加上独特字;具体过程如下:
信号帧包括集中式导频和数据两部分,导频部分用于计算自适应权值;在发送端将待发射的数据进行分块并在每个数据块的末尾添加独特字,且一个完整数据块的长度与导频的长度相同。
步骤二:自适应权值计算:发送端发送信号到达接收端并在接收端经过定时同步预处理后,截取接收到的导频;然后在空时均衡的过程中自适应地寻找最优时域滤波器非零抽头位置,具体利用贪婪算法或互相关矢量算法来具体确定多天线时域滤波器的非零抽头位置,进一步计算出自适应权值;具体过程如下:
假设接收机有M根天线,Δ为时域滤波器可能的最大阶数,Q为限定的时域非零抽头系数的个数,Q≤Δ+1。考虑接收机多天线集中放置,多天线多径信道的时延特征相同,即多天线信道具有相同的非零系数位置。基于此,空时滤波时,多天线时域滤波非零抽头的位置相同。
采用N表示计算空时权值的导频长度,L表示数据信号的长度,定义以下符号:
s:N×1维理想导频信号;
r:L×1维的数据信号部分;
xk,m:第m根天线上延迟k位的N×1维接收导频信号,m=1,…,M;k=0,…,Δ;
zk,m:第m根天线上延迟k位的L×1维接收数据信号,m=1,…,M;k=0,…,Δ;
Xk:所有天线上延迟k位的N×M维接收导频时间扩展矩阵,Xk=[xk,1,xk,2,…,xk,M];
Zk:所有天线上延迟k位的L×M维接收数据时间扩展矩阵,Zk=[zk,1,zk,2,…,zk,M];
wq:第q个滤波系数矢量,由所有天线的第q个滤波系数构成,M×1维。
利用最小均方误差(MMSE)准则,权值计算可表示为:
Figure BDA0001902877800000061
其中
Figure BDA0001902877800000062
T=Δ+1表示空时接收导频扩展矩阵,
Figure BDA0001902877800000063
由于有稀疏性要求,wq矢量中只有Q个非零,其余全为零矢量。
稀疏性时域滤波器抽头确定的核心在于寻找最优的非零权值的位置。以下方案用于在T个可能的位置中挑选出Q个非零抽头系数。
采用贪婪算法来具体确定多天线时域滤波器的非零抽头位置,具体过程如下:
①用Π={0,1,2,…,Δ}来表示待选择的抽头位置集合,用Ω={}表示已经选择的非零抽头位置,
Figure BDA0001902877800000064
表示已经选择的非零抽头位置对应的导频块;
②从时域滤波器非零抽头系数只有一个的情况开始考虑。这时待选择的抽头共有T=Δ+1个,基于MMSE准则,将这T=Δ+1个可能的抽头遍历,从中选择出一个。即
Figure BDA0001902877800000065
q从0遍历到Δ,对应的接收导频块为X0,X1,……,XΔ,其中
Figure BDA0001902877800000071
为权值测试矢量,它由下式求得:
Figure BDA0001902877800000072
其中,自相关矩阵
Figure BDA0001902877800000073
互相关矩阵
Figure BDA0001902877800000074
③遍历之后选择出第一个非零抽头,然后从Π中剔除掉选中的
Figure BDA0001902877800000075
将其并入到Ω中,
Figure BDA0001902877800000076
中并入第
Figure BDA0001902877800000077
个抽头位置对应的导频块
Figure BDA0001902877800000078
即:
Figure BDA0001902877800000079
④在确定了第一个抽头位置之后,再按照相同的步骤,遵循最小均方误差(MMSE)准则,从剩余的T-1位置中选择下一个非零抽头,前两个非零位置确定后,再类似从剩余的T-2个位置中挑选出一个。依次类推,直到挑选出所有Q个非零抽头位置。
最终,选择的抽头位置集合为Ω,选中的导频块集合为
Figure BDA00019028778000000710
对应的自适应权矢量为
Figure BDA00019028778000000711
采用互相关矢量算法来具体确定多天线时域滤波器的非零抽头位置的具体过程如下:
①接收端收到空时接收导频扩展矩阵Y=[X0 X1 … XΔ]之后,计算其互相关矢量p=YHs;
②将互相关矢量p表示为
Figure BDA00019028778000000712
其中
Figure BDA00019028778000000713
表示对应第q个位置的接收导频信号互相关矢量。
③计算||pq||2,q=0,1,…,Δ,从这Δ个值中挑选出最大的Q个值,即
Figure BDA00019028778000000714
将其对应的索引作为需要选择的Q个非零抽头位置。
根据贪婪算法或互相关矢量算法得到的非零抽头位置,确定对应的接收端导频块,再利用MMSE准则,求出维纳解即可得到自适应权值
Figure BDA0001902877800000081
Figure BDA0001902877800000082
若将选中的接收端导频块的集合用
Figure BDA0001902877800000083
表示,则自相关矢量
Figure BDA0001902877800000084
互相关矢量
Figure BDA0001902877800000085
步骤三:信号恢复:接收端将接收到的数据进行分块,再利用步骤二得到的自适应权值对每根天线上的接收数据块进行合并,并进行去独特字操作,从而实现信道的均衡;具体过程如下:
利用步骤二得到的自适应权值
Figure BDA0001902877800000086
对各数据块进行均衡:
Figure BDA0001902877800000087
其中,
Figure BDA0001902877800000088
为接收之后恢复出来的数据部分,
Figure BDA0001902877800000089
为选中抽头位置对应的接收数据块,去掉恢复之后各数据块末尾的独特字,完成信道的均衡。
下面通过一个实施例对本发明进行详细描述。
本发明适用的系统模型参见图1。在本发明的无线通信系统模型中,发射端采用单天线,接收端是由M根天线组成的天线阵列。无线信道模型是服从瑞利衰落的多径延时信道。
图2为本发明采用的信号帧结构。信号帧包括集中式导频和数据两部分,导频部分用于计算自适应权值;在信号帧结构的基础上,将待发射的数据进行分块并在数据块的末尾添加Unique word,且一个完整数据块的长度应与导频长度相同。后续仿真过程中,导频部分采用BPSK调制,数据部分采用QPSK调制。
图3为本发明空时处理的过程示意图。与传统空时处理的不同之处在于,需要从时域滤波器的所有抽头中选择非零抽头以及如何确定这些非零抽头,进而计算出自适应权值,进行信号合并。其计算方式如下:
假设接收机有M根天线,Δ为时域滤波器可能的最大阶数,Q为限定的时域非零抽头系数的个数,Q≤Δ+1。考虑接收机多天线集中放置,多天线多径信道的时延特征相同,即多天线信道具有相同的非零系数位置。基于此,空时滤波时,多天线时域滤波非零抽头的位置相同。
采用N表示计算空时权值的导频长度,L表示数据信号的长度,定义以下符号:
s:N×1维理想导频信号;
r:L×1维的数据信号部分;
xk,m:第m根天线上延迟k位的N×1维接收导频信号,m=1,…,M;k=0,…,Δ;
zk,m:第m根天线上延迟k位的L×1维接收数据信号,m=1,…,M;k=0,…,Δ;
Xk:所有天线上延迟k位的N×M维接收导频时间扩展矩阵,Xk=[xk,1,xk,2,…,xk,M];
Zk:所有天线上延迟k位的L×M维接收数据时间扩展矩阵,Zk=[zk,1,zk,2,…,zk,M];
wq:第q个滤波系数矢量,由所有天线的第q个滤波系数构成,M×1维。
利用最小均方误差(MMSE)准则,权值计算可表示为:
Figure BDA0001902877800000091
其中
Figure BDA0001902877800000092
T=Δ+1表示空时接收导频扩展矩阵,
Figure BDA0001902877800000093
由于有稀疏性要求,wq矢量中只有Q个非零,其余全为零矢量。
稀疏性时域滤波器抽头确定的核心在于寻找最优的非零权值的位置。以下方案用于在T个可能的位置中挑选出Q个非零抽头系数。
方案一:贪婪算法
①从时域滤波器非零抽头系数只有一个的情况开始考虑。这时待选择的抽头共有T=Δ+1个,用Π={0,1,2,…,Δ}来表示待选择的抽头位置集合,用Ω={}表示已经选择的非零抽头位置,
Figure BDA0001902877800000094
表示已经选择的非零抽头位置对应的导频块;
②基于最小均方误差(MMSE)准则,将这T=Δ+1个可能的抽头遍历,从中选择出一个。即
Figure BDA0001902877800000101
q从0遍历到Δ,对应的接收导频块为X0,X1,……,XΔ,其中
Figure BDA0001902877800000102
为权值测试矢量,它由下式求得:
Figure BDA0001902877800000103
其中,自相关矩阵
Figure BDA0001902877800000104
互相关矩阵
Figure BDA0001902877800000105
③遍历之后选择出第一个非零抽头,然后从Π中剔除掉选中的
Figure BDA0001902877800000106
将其并入到Ω中,
Figure BDA0001902877800000107
中并入第
Figure BDA0001902877800000108
个抽头位置对应的导频块
Figure BDA0001902877800000109
即:
Figure BDA00019028778000001010
④在确定了第一个抽头位置之后,再按照相同的步骤,遵循最小均方误差(MMSE)准则,从剩余的T-1个位置中选择下一个非零抽头,前两个非零位置确定后,再类似从剩余的T-2个位置中挑选出一个。依次类推,直到挑选出所有Q个非零抽头系数。最终,选择的抽头位置集合为Ω,对应的自适应权矢量为
Figure BDA00019028778000001011
方案二:互相关矢量算法
①接收端收到空时接收导频扩展矩阵Y=[X0 X1 … XΔ]之后,计算其互相关矢量p=YHs;
②将p表示为
Figure BDA00019028778000001012
其中
Figure BDA00019028778000001013
表示对应第q个位置的接收导频信号互相关矢量。
③计算||pq||2,q=0,1,…,Δ,从这Δ个值中挑选出最大的Q个值,即
Figure BDA00019028778000001014
将其对应的索引作为需要选择的Q个非零抽头的位置,进一步确定对应的导频块。
④利用最小均方误差(MMSE)准则,求出维纳解即可得到自适应权矢量
Figure BDA00019028778000001015
Figure BDA0001902877800000111
若将选中的导频块的集合用
Figure BDA0001902877800000112
表示,则自相关矢量
Figure BDA0001902877800000113
互相关矢量
Figure BDA0001902877800000114
利用步骤二得到的自适应权值
Figure BDA0001902877800000115
对各数据块进行均衡:
Figure BDA0001902877800000116
其中,
Figure BDA0001902877800000117
为接收之后恢复出来的数据部分,
Figure BDA0001902877800000118
为选中抽头位置对应的接收数据块,去掉恢复之后各数据块末尾的独特字,完成信道的均衡。
图4为本发明提出的贪婪算法、互相关矢量算法同对比方案(即非自适应算法,在时域滤波器的T个抽头中选前Q个作为非零抽头)的BER性能示意图。参数设置为:发送端发送的导频长度为512,每个数据块512个符号,其中UW长度为64;基带等效信道长度为40,信道的非零抽头位置随机确定,非零抽头个数为10,均服从等功率瑞利衰落;接收端的天线数M=4,时域滤波器共有T=40个抽头,其中非零抽头个数分别为Q=10,20,30。此外还给出了Q=T=40的仿真结果作为参照。从仿真结果可以看出:随着Q取值增加,三种方案的BER性能都呈优化趋势。但是不论Q=10,20,30三种取值中的哪一种,本发明提出的贪婪算法和互相关矢量算法的BER性能基本上优于传统的非自适应算法。
图5为本发明在有干扰情况下的BER性能示意图。参数设置与图4基本相同,时域滤波器的非零抽头个数Q=20,30,信干比SIR=0dB,类似地,给出Q=T=40的仿真结果作为参照。仿真表明,本发明所提出的贪婪算法和互相关矢量算法对外部干扰造成的BER性能下降也能在极大程度上优化。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (2)

1.一种基于稀疏感知的自适应空时均衡方法,其特征在于,包括以下步骤:
步骤一:在发送端将待发射的数据分块,然后在每个数据块的末尾加上独特字;
步骤二:发送端发送信号到达接收端并在接收端经过定时同步预处理后,截取接收到的导频;然后在空时均衡的过程中采用贪婪算法或互相关矢量算法自适应地寻找最优时域滤波器非零抽头位置,然后计算自适应权值;
步骤三:接收端将接收到的数据进行分块,再利用步骤二得到的自适应权值对每根天线上的接收数据块进行合并,并进行去独特字操作,从而实现信道的均衡;
步骤二中采用贪婪算法自适应地寻找最优时域滤波器非零抽头位置的具体过程如下:
假设接收机有M根天线,Δ为时域滤波器可能的最大阶数,Q为限定的时域非零抽头系数的个数,Q≤Δ+1;接收机多天线集中放置,多天线多径信道的时延特征相同,空时滤波时,多天线时域滤波非零抽头的位置相同;
采用N表示计算空时权值的导频长度,L表示数据信号的长度,定义以下符号:
s:N×1维理想导频信号;
r:L×1维的数据信号部分;
xk,m:第m根天线上延迟k位的N×1维接收导频信号,m=1,…,M;k=0,…,Δ;
zk,m:第m根天线上延迟k位的L×1维接收数据信号,m=1,…,M;k=0,…,Δ;
Xk:所有天线上延迟k位的N×M维接收导频时间扩展矩阵,Xk=[xk,1,xk,2,…,xk,M];
Zk:所有天线上延迟k位的L×M维接收数据时间扩展矩阵,Zk=[zk,1,zk,2,…,zk,M];
wq:第q个滤波系数矢量,由所有天线的第q个滤波系数构成,M×1维;
①采用Π={0,1,2,…,Δ}来表示待选择的抽头位置集合,采用Ω={}表示已经选择的非零抽头位置,
Figure FDA0002341679400000011
表示已经选择的非零抽头位置对应的导频块;
②从时域滤波器非零抽头系数只有一个的情况开始,这时待选择的抽头共有T=Δ+1个,基于最小均方误差准则准则,将这T=Δ+1个可能的抽头遍历,从中选择一个,即
Figure FDA0002341679400000021
其中,
q从0遍历到Δ,对应的接收导频块为X0,X1,……,XΔ,其中
Figure FDA0002341679400000022
为权值测试矢量,由下式求得:
Figure FDA0002341679400000023
其中,自相关矩阵
Figure FDA0002341679400000024
互相关矩阵
Figure FDA0002341679400000025
③遍历之后选择第一个非零抽头,然后从Π中剔除掉选中的
Figure FDA0002341679400000026
将其并入到Ω中,
Figure FDA0002341679400000027
中并入第
Figure FDA0002341679400000028
个抽头位置对应的导频块
Figure FDA0002341679400000029
即:
Figure FDA00023416794000000210
④在确定了第一个抽头位置之后,再遵循最小均方误差准则,从剩余的T-1位置中选择下一个非零抽头,前两个非零位置确定后,再类似从剩余的T-2个位置中挑选出一个,依次类推,直到挑选出所有Q个非零抽头位置;
步骤二中采用互相关矢量算法自适应地寻找最优时域滤波器非零抽头位置的具体过程如下:
①接收端收到空时接收导频扩展矩阵Y=[X0 X1…XΔ]之后,计算其互相关矢量p=YHs;其中,s为N×1维理想导频信号;
②将互相关矢量p表示为
Figure FDA00023416794000000211
其中
Figure FDA00023416794000000212
表示对应第q个位置的接收导频信号互相关矢量;
③计算||pq||2,q=0,1,…,Δ,从这Δ个值中挑选出最大的Q个值,即
Figure FDA00023416794000000213
将其对应的索引作为需要选择的Q个非零抽头位置;
计算自适应权值的具体过程为:根据贪婪算法或互相关矢量算法得到的非零抽头位置,确定对应的接收端导频块,再利用最小均方误差准则准则,求出维纳解,得到自适应权值
Figure FDA0002341679400000031
Figure FDA0002341679400000032
其中,R为自相关矩阵,p为互相关矢量;
采用贪婪算法选择的抽头位置集合为Ω,选中的接收导频块集合为
Figure FDA0002341679400000033
则自相关矢量
Figure FDA0002341679400000034
互相关矢量
Figure FDA0002341679400000035
2.根据权利要求1所述的一种基于稀疏感知的自适应空时均衡方法,其特征在于,步骤三的具体过程如下:
利用步骤二得到的自适应权值
Figure FDA0002341679400000036
对各数据块进行均衡:
Figure FDA0002341679400000037
其中,
Figure FDA0002341679400000038
为接收之后恢复出来的数据部分,
Figure FDA0002341679400000039
为选中抽头位置对应的接收数据块,去掉恢复之后各数据块末尾的独特字,完成信道的均衡。
CN201811519562.1A 2018-12-12 2018-12-12 一种基于稀疏感知的自适应空时均衡方法 Active CN109379308B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811519562.1A CN109379308B (zh) 2018-12-12 2018-12-12 一种基于稀疏感知的自适应空时均衡方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811519562.1A CN109379308B (zh) 2018-12-12 2018-12-12 一种基于稀疏感知的自适应空时均衡方法

Publications (2)

Publication Number Publication Date
CN109379308A CN109379308A (zh) 2019-02-22
CN109379308B true CN109379308B (zh) 2020-03-31

Family

ID=65373322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811519562.1A Active CN109379308B (zh) 2018-12-12 2018-12-12 一种基于稀疏感知的自适应空时均衡方法

Country Status (1)

Country Link
CN (1) CN109379308B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111917522B (zh) * 2020-06-30 2021-05-28 西安交通大学 一种抗干扰宽带单载波空时编码传输方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103401824A (zh) * 2013-07-22 2013-11-20 西安电子科技大学 基于修正牛顿法的频率选择性mimo系统空时盲均衡器方法
CN108650007A (zh) * 2018-05-03 2018-10-12 西安交通大学 一种基于空频自适应滤波的高可靠信道均衡方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7440988B2 (en) * 2004-04-08 2008-10-21 Raytheon Company System and method for dynamic weight processing
US8275019B2 (en) * 2009-06-16 2012-09-25 New Jersey Institute Of Technology Transmitting and receiving schemes for multiuser single-carrier space time spreading with frequency domain equalization
CN105553898A (zh) * 2015-12-18 2016-05-04 中国电子科技集团公司第三研究所 均衡器及反馈均衡方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103401824A (zh) * 2013-07-22 2013-11-20 西安电子科技大学 基于修正牛顿法的频率选择性mimo系统空时盲均衡器方法
CN108650007A (zh) * 2018-05-03 2018-10-12 西安交通大学 一种基于空频自适应滤波的高可靠信道均衡方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Burst-Interference Suppression Based on space time processing";WANG HUILIANG;《IEEE ACCESS》;20171218;全文 *
"On Adaptive Length of Temporal Filter for Space-Time Equalization With Cochannel Interference";Ge Yinghao;《IEEE SIGNAL PROCESSING LETTERS》;20180731;全文 *

Also Published As

Publication number Publication date
CN109379308A (zh) 2019-02-22

Similar Documents

Publication Publication Date Title
US6151487A (en) Demodulation structure for fast fading cellular channels
US8369438B2 (en) Iterative tree search-based precoding technique for multiuser MIMO communication system
JP2001517399A (ja) 自己同期等化方法及びシステム
CN101150343A (zh) 一种mimo移动通信方法及系统
CN108650007B (zh) 一种基于空频自适应滤波的高可靠信道均衡方法
CN104702390A (zh) 分布式压缩感知信道估计中的导频分配方法
CN111212005A (zh) 一种基于重定时同步和干扰抵消的信号检测方法
CN109379308B (zh) 一种基于稀疏感知的自适应空时均衡方法
CN114665930A (zh) 去蜂窝大规模mimo系统的下行链路盲信道估计方法
CN104468432B (zh) 一种短波信道下单载波频域均衡的信道估计去噪声方法
CN108832978B (zh) 一种包含直传链路的多用户mimo中继系统的联合预编码方法
CN104253638B (zh) 基于Stiefel流形上共轭梯度法的MIMO干扰对齐算法
EP1311094A2 (en) Calculation and transmission of coefficients for beamforming or equalisation
CN110350985B (zh) 一种基于主动时间反转的水声并行传输方法
CN101958862A (zh) 基于叠加结构的自适应判决反馈均衡器
CN107171990A (zh) 基于分数间隔与迭代算法的单载波信号时域均衡方法
CN101951277B (zh) 基于最大似然检测的二维块扩频系统中的残留多用户干扰抑制方法
CN111181879B (zh) 一种基于自适应基展开的抗干扰时变信道均衡方法
CN111917522B (zh) 一种抗干扰宽带单载波空时编码传输方法
CN107154813B (zh) 自适应Rake接收机及接收方法
CN105262524B (zh) 基于反馈补偿的mimo系统的多用户调度方法
JP2008072403A (ja) 周波数選択性通信路に於ける固有モード伝送方法
CN101262265A (zh) 一种时分双工无线通信系统的延时分集发射与接收的方法
CN110535796B (zh) 一种低复杂度空频自适应抗干扰信道均衡方法
CN110912590A (zh) 基于信道反转技术的大规模衰落mimo系统的干扰抑制预编码方法

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