CN114785645B - 基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法 - Google Patents
基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法 Download PDFInfo
- Publication number
- CN114785645B CN114785645B CN202210209866.8A CN202210209866A CN114785645B CN 114785645 B CN114785645 B CN 114785645B CN 202210209866 A CN202210209866 A CN 202210209866A CN 114785645 B CN114785645 B CN 114785645B
- Authority
- CN
- China
- Prior art keywords
- signal
- transmitted
- ofdm
- representing
- antenna
- 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 92
- 230000005540 biological transmission Effects 0.000 claims abstract description 190
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 77
- 230000000737 periodic effect Effects 0.000 claims abstract description 56
- 238000013507 mapping Methods 0.000 claims abstract description 55
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 24
- 238000005562 fading Methods 0.000 claims description 113
- 238000005070 sampling Methods 0.000 claims description 49
- 230000008569 process Effects 0.000 claims description 25
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 230000006835 compression Effects 0.000 claims description 15
- 238000007906 compression Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 14
- 230000003111 delayed effect Effects 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 4
- 108010074864 Factor XI Proteins 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000001934 delay Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 abstract description 5
- 239000013256 coordination polymer Substances 0.000 abstract 1
- 238000004891 communication Methods 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 7
- 238000010295 mobile communication Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008054 signal transmission Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 239000007822 coupling agent Substances 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- 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/0238—Channel estimation using blind estimation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
-
- 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/0242—Channel estimation channel estimation algorithms using matrix methods
-
- 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
-
- 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/2602—Signal structure
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Radio Transmission System (AREA)
Abstract
本发明公开了一种基于CP‑Free MIMO‑OFDM信号循环平稳特性的盲信噪比估计方法,其通过使用OFDM子载波集映射,使OFDM信号具有循环平稳特性;设计一种匹配信号来消除去掉CP而产生的码间干扰,并与OFDM信号相加得到发送信号;通过改变子载波映射集的长度,选择不同的循环频率来排除其他信道的干扰,将MIMO信道转换为多个单输入单输出信道;利用发送信号的周期自相关函数能量分布规律,实现盲信噪比估计;优点是其信噪比估计准确度高,且其使得OFDM信号具有循环平稳特性,降低了噪声的影响并能提高系统的频谱效率和降低系统的传输时延。
Description
技术领域
本发明涉及一种信噪比估计技术,尤其是涉及一种基于CP-Free(Cyclic PrefixFree,无循环前缀)MIMO-OFDM(Multiple Input Multiple Output,多入多出;OrthogonalFrequency Divided Multiplexing,正交频分复用)信号循环平稳特性的盲信噪比估计方法。
背景技术
20世纪80年代中期以后,出现了第一代模拟移动通信系统,其采用模拟信号通信方法,使得人们可以通过移动电话实现远程实时通信。随后,第二代移动通信系统采用了数字信号传输方式,提高了无线通信中信号传输的稳定性和抗噪声能力,人们的基本语音通话需求已经得到满足,期盼着能够通过移动通信系统对图片、视频等媒体文件进行传输和处理,实现看电影、视频通话等功能。于是,采用码分多址技术的第三代移动通信系统应运而生。随着移动互联网和电子信息技术的进一步发展,移动通信用户数量以及人们对于数据流量的需求出现爆发式增长,第四代移动通信技术OFDM(Orthogonal FrequencyDivided Multiplexing,正交频分复用)技术和MIMO(Multiple Input Multiple Output,多入多出)技术,大大提高了移动通信系统的频带利用率和系统信道容量。MIMO技术与OFDM技术相结合的MIMO-OFDM技术被广泛应用于军事、交通、物联网等领域中。但是OFDM技术通过插入循环前缀(CP,Cyclic Prefix)来消除多径效应带来的码间干扰,会造成频谱效率降低,传输时延增加,不符合5G低时延、高速率的特点。因此CP-Free MIMO-OFDM的研究显得尤为必要。
MIMO技术和OFDM技术虽然实现了信号传输速率和系统信道容量的提高,推进了无线通信系统的发展,但是由于无线通信系统中天线数、传输信号和用户数量的增加,因此也增加了信道参量估计、信号检测的难度,带来了导频干扰、信号功率分配等新的问题。在无线通信系统中,信噪比值的高低反映了通信质量的好坏,准确的信噪比值是实现高速、高可靠性传输的重要指标。在通信信号处理的很多研究领域,功能强大的信号处理算法在运行过程中都需要通信信道的信噪比值,若无法获得精确的信噪比信息,那么许多信号处理算法的性能就会明显下降,无法满足无线通信系统的性能要求。由于受大型建筑物和天气等因素的影响,发射信号容易出现绕射、反射等现象,同时由于传输环境的不确定,容易产生多径衰落、频率偏移、信号衰减、时延等,空气中存在的热噪声、宇宙噪声等也会对信号造成干扰。因此,如何在信号接收端对大量的信号进行准确的信噪比估计已成为无线通信系统中的研究重点。
目前,多用户无CP的CP-Free MIMO-OFDM系统的盲信噪比估计还未有相关报道,而单个用户的信噪比估计方法大致可以分为两类:一类是基于数据辅助的,即在各个符号中插入导频(训练序列);另一类是基于非数据辅助的,即不依赖发送端的已知信号,只根据接收信号来进行处理,也就是盲信噪比估计方法。这两类信噪比估计方法各有优缺点,基于数据辅助的信噪比估计方法的精确度较高,但是需要发送端发送导频序列与接收端完成同步或在接收端完成接收判决才能进行估计;而基于非数据辅助的盲信噪比估计方法不需要发送导频信息,这样大大提高了频谱利用率,也不需要与接收端同步和完成判决,但其计算复杂度较高。
盲信噪比估计方法主要是基于信号本身的统计特性,如分离符号累计量估计方法(SSME方法)、二阶四阶矩方法(M2M4方法)、信号方差比方法(SVR方法)、平方信号噪声方差方法(SNV方法)等。但是这些方法在使用中都存在不同的问题,例如SSME方法只能在高斯白噪声信道下二进制相移键控(Binary Phase Shift Keying,BPSK)调制时才能进行信噪比估计,而SNV方法在低信噪比时性能较差等。对于CP-Free MIMO-OFDM系统,目前并未有对其信噪比进行估计的相关报道。在传统的MIMO-OFDM系统中,为消除符号间干扰(ISI,Inter-Symbol Interference)而加入的循环前缀(CP,Cyclic Prefix),使得OFDM信号具有循环平稳特性,但当CP过度使用时,会造成系统频谱效率降低,传输时延增加。
发明内容
本发明所要解决的技术问题是提供一种基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其信噪比估计准确度高,且其使得OFDM信号具有循环平稳特性,降低了噪声的影响并能够提高系统的频谱效率和降低系统的传输时延。
本发明解决上述技术问题所采用的技术方案为:一种基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于包括以下步骤:
步骤1:设定CP-Free MIMO-OFDM系统中的各个信道采用多径衰落模型,即CP-FreeMIMO-OFDM系统中的各个信道为多径衰落信道;设定CP-Free MIMO-OFDM系统的发送端具有MT个发射天线且接收端具有MR个接收天线;其中,MT表示发射天线的个数,MT>1,MR表示接收天线的个数,MR>1;
步骤2:在CP-Free MIMO-OFDM系统的发送端,将接收的频域数据信号变换为MT个时域数据信号,每个时域数据信号为由多个长度一致的OFDM符号构成的OFDM信号,所有时域数据信号各自包含的OFDM符号的个数相同;然后对每个OFDM信号中的各个OFDM符号进行子载波映射,得到对应的子载波映射后的OFDM信号;接着对每个子载波映射后的OFDM信号先后进行数模转换和射频模块处理,得到对应的射频模块处理后的OFDM信号;再在每个射频模块处理后的OFDM信号上叠加一个用于消除码间干扰的匹配信号,得到对应的发送信号;
步骤3:CP-Free MIMO-OFDM系统的发送端将MT个发送信号各通过一个发射天线相互平行地传输给CP-Free MIMO-OFDM系统的接收端;
步骤4:在CP-Free MIMO-OFDM系统的接收端,每个接收天线接收到的接收信号由一个发送信号通过一个发射天线发送后经过MT个多径衰落信道后形成的具有不同时延的MT路信号组成,将通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被第r个接收天线接收到的接收信号记为yr(n),其中,1≤i≤MT,ti(n)表示通过第i个发射天线发送的发送信号,1≤r≤MR,1≤m≤Lh,Lh表示多径衰落信道的多径阶数,Si,r,m表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的信号功率因子,e为自然基数,j表示复数中的虚数单位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的到达相位,τi,r,m表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的时延,ti(n-τi,r,m)表示通过第i个发射天线发送的发送信号ti(n)经过hi,r的第m径的时延后的信号,N表示高斯白噪声的噪声功率因子,vi,r(n)表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r上的平稳高斯白噪声;
步骤5:在CP-Free MIMO-OFDM系统的接收端,根据自相关函数的定义,获取yr(n)的自相关函数,记为 其中,τ表示延时变量,E()表示数学期望,yr(n+τ)表示yr(n)延时τ后的信号,表示yr(n+τ)的共扼,为yr(n)的自相关函数,1≤m1≤Lh,1≤m2≤Lh,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的信号功率因子,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的信号功率因子,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的到达相位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的到达相位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的时延,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的时延,表示通过第i个发射天线发送的发送信号ti(n)经过hi,r的第m1径的时延后的信号,表示的共扼,表示ti(n)延时τ后的信号ti(n+τ)经过hi,r的第m2径的时延后的信号,表示ti(n)的自相关函数延时后的值,δ(τ)为vi,r(n)的自相关函数,δ(τ)是延时变量τ的冲击函数;
步骤6:在CP-Free MIMO-OFDM系统的接收端,对作傅立叶级数展开,得到yr(n)的周期自相关函数,记为 其中,k为循环频率,P表示OFDM信号的循环周期,亦为一个OFDM符号的长度,P=M+2L,M表示一个OFDM符号的数据长度,即表示OFDM信号的子载波数,Λ为引入的中间符号,表示ti(n)的周期自相关函数延时后的值,δ(τ)δ(k)为vi,r(n)的周期自相关函数,δ(k)是变量k的冲击函数;
步骤7:在CP-Free MIMO-OFDM系统的接收端,根据通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径和第m2径的时间间隔相等与否这两种情况,将转化为 然后根据发送信号的周期自相关函数的能量分布,令得到其中,表示ti(n)的周期自相关函数;
步骤8:在CP-Free MIMO-OFDM系统的接收端,假设需要估计通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比,那么通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数为然后根据yb(n)是否有通过指定的第a个发射天线发送的发送信号ta(n),将转化为其中,a∈[1,MT],b∈[1,MR],Si,b,m表示通过第i个发射天线发送的发送信号ti(n)到指定的第b个接收天线接收所经过的多径衰落信道hi,b的第m径的信号功率因子,τi,b,m表示通过第i个发射天线发送的发送信号ti(n)到指定的第b个接收天线接收所经过的多径衰落信道hi,b的第m径的时延,1≤J≤MT,且J≠a,SJ,b,m表示通过第J个发射天线发送的发送信号tJ(n)到指定的第b个接收天线接收所经过的多径衰落信道hJ,b的第m径的信号功率因子,Sa,b,m表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的第m径的信号功率因子,表示通过第J个发射天线发送的发送信号tJ(n)的周期自相关函数,表示通过指定的第a个发射天线发送的发送信号ta(n)的周期自相关函数,τJ,b,m表示通过第J个发射天线发送的发送信号tJ(n)到指定的第b个接收天线接收所经过的多径衰落信道hJ,b的第m径的时延,τa,b,m表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的第m径的时延;
步骤9:在CP-Free MIMO-OFDM系统的接收端,对通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比进行估计,具体过程为:
步骤9_1:在的基础上,令ta(n)的循环频率等于ka且ka≠0,在排除ha,b的干扰信道的影响后,得到ta(n)的循环频率等于ka时yb(n)的周期自相关函数,记为 然后对的等号两边同时取绝对值,得到其中,表示循环频率k等于ka时通过指定的第a个发射天线发送的发送信号ta(n)的周期自相关函数,符号“||”为取绝对值符号,ξa,b表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的信号功率因子;
步骤9_2:分析当k=0且τ=0时,转化为然后对的等号两边同时取绝对值,得到其中,表示k=0且τ=0时通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数,表示k=0且τ=0时通过第i个发射天线发送的发送信号ti(n)的周期自相关函数,ξ1,b表示通过第1个发射天线发送的发送信号t1(n)到指定的第b个接收天线接收所经过的多径衰落信道h1,b的信号功率因子,表示k=0且τ=0时通过第1个发射天线发送的发送信号t1(n)的周期自相关函数,ξ2,b表示通过第2个发射天线发送的发送信号t2(n)到指定的第b个接收天线接收所经过的多径衰落信道h2,b的信号功率因子,表示k=0且τ=0时通过第2个发射天线发送的发送信号t2(n)的周期自相关函数,表示通过第MT个发射天线发送的发送信号到指定的第b个接收天线接收所经过的多径衰落信道的信号功率因子,表示k=0且τ=0时通过第MT个发射天线发送的发送信号的周期自相关函数;
步骤9_3:对通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数进行统计量估计,将统计量估计后的值记为 其中,Nnum表示每个OFDM信号中包含的OFDM符号的个数,wb(n)表示通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线实际接收到的接收信号,wb(n+τ)表示wb(n)延时τ后的信号,wb *(n+τ)为wb(n+τ)的共轭;
步骤9_4:估计通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的信号功率因子ξa,b的值和高斯白噪声的噪声功率因子N的值,将ξa,b的估计值和N的估计值对应记为和 然后计算通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比的估计值,记为ρa,b,其中,表示对进行统计量估计后得到的值,表示对进行统计量估计后得到的值,和通过计算得到,表示ξ1,b的估计值,表示ξ2,b的估计值,表示的估计值,和通过计算得到。
所述的步骤2中,将接收的频域数据信号变换为MT个时域数据信号的具体过程为:
步骤2_1:对接收的频域数据信号进行串并转化,分成MT个子符号流;
步骤2_2:对每个子符号流进行无失真压缩编码,并在无失真压缩编码过程中在该子符号流中加入用于克服该子符号流在多径衰落信道中受到的干扰和噪声影响的冗余信息,得到对应的无失真压缩编码后的频域数据信号;
步骤2_3:对每个无失真压缩编码后的频域数据信号进行空时调制,得到对应的空时调制后的频域数据信号;
步骤2_4:对每个空时调制后的频域数据信号先后进行傅里叶逆变换和正交频分复用调制处理,得到对应的时域数据信号,即共得到MT个时域数据信号。
所述的步骤2_2中,采用信道编码技术对每个子符号流进行无失真压缩编码。
所述的步骤2_3中,利用调制器对每个无失真压缩编码后的频域数据信号进行空时调制。
所述的步骤2中,对每个OFDM信号中的各个OFDM符号进行子载波映射的具体过程为:对于任一个OFDM符号,将该OFDM符号中的第1个采样点至第L个采样点按序构成第一映射集,并将该OFDM符号中的第L+1个采样点至第2L个采样点按序构成第二映射集;然后将该OFDM符号中的第1个采样点作为第一起始点,间隔p个采样点后将第一映射集中的各个采样点按序替换该OFDM符号中的第p+1个采样点至第p+L个采样点;并将该OFDM符号中的第L+1个采样点作为第二起始点,间隔p个采样点后将第二映射集中的各个采样点按序替换该OFDM符号中的第p+L+1个采样点至第p+2L个采样点,经过两次采样点替换操作后实现了该OFDM符号的子载波映射过程;其中,L表示第一映射集和第二映射集各自中包含的采样点的个数,即为第一映射集和第二映射集各自的长度,M表示一个OFDM符号的数据长度,即表示OFDM信号的子载波数,p=M。
所述的步骤2中,在每个射频模块处理后的OFDM信号上叠加一个用于消除码间干扰的匹配信号,得到对应的发送信号的具体过程为:将在第i'个射频模块处理后的OFDM信号上叠加匹配信号后得到的发送信号记为ti'(n),ti'(n)=si'(n)+ci'(n);其中,1≤i'≤MT,n表示离散时间点,ti'(n)为离散时间信号,si'(n)表示第i'个射频模块处理后的OFDM信号,ci'(n)表示在si'(n)上叠加的用于消除码间干扰的匹配信号,C表示信道矩阵,Cinter为C的干扰部分,CH为C的厄米特矩阵。
所述的步骤9_1中,排除ha,b的干扰信道的影响的具体过程为:令通过指定的第a个发射天线发送的发送信号ta(n)中的OFDM符号进行子载波映射时的第一映射集和第二映射集的长度为La,并令通过第J个发射天线发送的发送信号tJ(n)中的OFDM符号进行子载波映射时的第一映射集和第二映射集的长度为LJ,其中,a∈[1,MT],1≤J≤MT,且J≠a,La∈[1,L],LJ∈[1,L];然后令再使得通过第J个发射天线发送的发送信号tJ(n)的周期自相关函数处于零点,即使得且从而使得转化为排除了ha,b的干扰信道的影响。
所述的步骤9执行完毕后,重复执行步骤2至步骤9共Number次,得到ha,b上的盲信噪比的Number个估计值,然后将Number个估计值的平均值作为ha,b上的盲信噪比的最终估计值;其中,Number∈[500,1000]。
与现有技术相比,本发明的优点在于:
去掉CP以后,子载波全部用来传输数据,提高了MIMO-OFDM系统的频谱效率,降低了传输时延,并且设计匹配信号来消除码间干扰。通过利用OFDM信号子载波映射的方式使其具有循环平稳特性,降低了噪声对MIMO-OFDM系统的影响。通过选取适当映射集长度以及发送信号的周期自相关函数不同的零点,可以有效地把MIMO信道转化成多个单输入单输出信道,然后利用二阶统计量可准确估计出多径环境下各个信道的信噪比值,提高了MIMO-OFDM系统的频带利用率,仿真结果表明其信噪比估计准确性更高。
附图说明
图1为本发明方法(发射天线和接收天线的个数均为2的情况下)的总体实现框图;
图2为多径衰落模型的示意图;
图3为发送信号的周期自相关函数的能量分布图(幅度图);
图4为在每个OFDM信号中包含的OFDM符号的个数Nnum=1000,实际信噪比值分别为-10dB、-5dB、0dB、5dB、10dB、15dB的条件下,实际信噪比值、利用现有的基于CP OFDM的信噪比估计方法估计得到的信噪比值、利用本发明方法估计得到的盲信噪比值的比较示意图;
图5为利用现有的基于CP OFDM的信噪比估计方法估计得到的信噪比值和利用本发明方法估计得到的盲信噪比值各自的归一化均方误差随信噪比的变化曲线图;
图6为利用现有的基于CP OFDM的信噪比估计方法估计得到的信噪比值和利用本发明方法估计得到的盲信噪比值各自的归一化偏差随信噪比的变化曲线图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明提出的一种基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其总体实现框图如图1所示,其包括以下步骤:
步骤1:设定CP-Free MIMO-OFDM系统中的各个信道采用多径衰落模型,即CP-FreeMIMO-OFDM系统中的各个信道为多径衰落信道,图2给出了多径衰落模型的示意图;设定CP-Free MIMO-OFDM系统的发送端具有MT个发射天线且接收端具有MR个接收天线;其中,MT表示发射天线的个数,MT>1,如MT=2,MR表示接收天线的个数,MR>1,如MR=2。
步骤2:在CP-Free MIMO-OFDM系统的发送端,将接收的频域数据信号变换为MT个时域数据信号,每个时域数据信号为由多个长度一致的OFDM符号构成的OFDM信号,所有时域数据信号各自包含的OFDM符号的个数相同;然后对每个OFDM信号中的各个OFDM符号进行子载波映射,得到对应的子载波映射后的OFDM信号;接着对每个子载波映射后的OFDM信号先后进行数模转换和射频模块处理,得到对应的射频模块处理后的OFDM信号;由于循环前缀CP的去除使得多径效应带来的码间干扰(ISI)无法被消除,因此再在每个射频模块处理后的OFDM信号上叠加一个用于消除码间干扰的匹配信号,得到对应的发送信号。
在此,数模转换和射频模块处理均采用现有技术,射频模块处理主要起到信号调制(把信号加到频率较高的频带上)、功率放大等作用。
在本实施例中,步骤2中,将接收的频域数据信号变换为MT个时域数据信号的具体过程为:
步骤2_1:对接收的频域数据信号进行串并转化,分成MT个子符号流。
步骤2_2:对每个子符号流进行无失真压缩编码,并在无失真压缩编码过程中在该子符号流中加入用于克服该子符号流在多径衰落信道中受到的干扰和噪声影响的冗余信息,得到对应的无失真压缩编码后的频域数据信号。在此,采用信道编码技术对每个子符号流进行无失真压缩编码。
步骤2_3:对每个无失真压缩编码后的频域数据信号进行空时调制,得到对应的空时调制后的频域数据信号。在此,利用调制器对每个无失真压缩编码后的频域数据信号进行空时调制。
步骤2_4:对每个空时调制后的频域数据信号先后进行傅里叶逆变换和正交频分复用调制处理,得到对应的时域数据信号,即共得到MT个时域数据信号。
上述,串并转化、傅里叶逆变换、正交频分复用调制处理均采用现有技术。
在本实施例中,步骤2中,对每个OFDM信号中的各个OFDM符号进行子载波映射的具体过程为:对于任一个OFDM符号,将该OFDM符号中的第1个采样点至第L个采样点按序构成第一映射集,并将该OFDM符号中的第L+1个采样点至第2L个采样点按序构成第二映射集;然后将该OFDM符号中的第1个采样点作为第一起始点,间隔p个采样点后将第一映射集中的各个采样点按序替换该OFDM符号中的第p+1个采样点至第p+L个采样点;并将该OFDM符号中的第L+1个采样点作为第二起始点,间隔p个采样点后将第二映射集中的各个采样点按序替换该OFDM符号中的第p+L+1个采样点至第p+2L个采样点,经过两次采样点替换操作后实现了该OFDM符号的子载波映射过程;其中,L表示第一映射集和第二映射集各自中包含的采样点的个数,即为第一映射集和第二映射集各自的长度,在本实施例中取M表示一个OFDM符号的数据长度,即表示OFDM信号的子载波数,在本实施例中取M=32,实际应用过程中还可取16、或32、或64、或128等,p=M。
在本实施例中,步骤2中,在每个射频模块处理后的OFDM信号上叠加一个用于消除码间干扰的匹配信号,得到对应的发送信号的具体过程为:将在第i'个射频模块处理后的OFDM信号上叠加匹配信号后得到的发送信号记为ti'(n),ti'(n)=si'(n)+ci'(n);其中,1≤i'≤MT,n表示离散时间点,ti'(n)为离散时间信号,si'(n)表示第i'个射频模块处理后的OFDM信号,ci'(n)表示在si'(n)上叠加的用于消除码间干扰的匹配信号,C表示信道矩阵,Cinter为C的干扰部分,CH为C的厄米特矩阵。
在循环前缀CP去除后,第r个接收信号yr(n)为yr(n)=hi,r*ti(n)+vi,r(n)=hi,r*(si(n)+ci(n))+vi,r(n),然后将上式转换为矩阵相乘的形式,描述为: 为码间干扰部分,为了消除码间干扰,令由此可得到匹配信号消除码间干扰后的接收信号yr(n)为其中,1≤i≤MT,1≤r≤MR,hi,r表示第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道,符号“*”为卷积运算符号,vi,r(n)表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r上的平稳高斯白噪声,Cu表示C中无干扰的信道卷积矩阵,CX表示信道循环矩阵,C表示信道矩阵,Cinter为C的干扰部分,si(n)表示第i个射频模块处理后的OFDM信号,ci(n)表示在si(n)上叠加的用于消除码间干扰的匹配信号。
步骤3:CP-Free MIMO-OFDM系统的发送端将MT个发送信号各通过一个发射天线相互平行(即同一时刻)地传输给CP-Free MIMO-OFDM系统的接收端。
步骤4:在CP-Free MIMO-OFDM系统的接收端,每个接收天线接收到的接收信号由一个发送信号通过一个发射天线发送后经过MT个多径衰落信道后形成的具有不同时延的MT路信号组成,将通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被第r个接收天线接收到的接收信号记为yr(n),其中,1≤i≤MT,ti(n)表示通过第i个发射天线发送的发送信号,1≤r≤MR,1≤m≤Lh,Lh表示多径衰落信道的多径阶数,在本实施例中将多径衰落信道等效成4阶FIR(Finite Impulse Response)滤波器,即Lh=4,Si,r,m表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的信号功率因子,e为自然基数,e=2.71…,j表示复数中的虚数单位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的到达相位,τi,r,m表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的时延,ti(n-τi,r,m)表示通过第i个发射天线发送的发送信号ti(n)经过hi,r的第m径的时延后的信号,N表示高斯白噪声的噪声功率因子,vi,r(n)表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r上的平稳高斯白噪声,对于平稳高斯白噪声,其数学期望和方差与时间无关,自相关函数只与时间间隔有关。
;其中,τ表示延时变量,E()表示数学期望,yr(n+τ)表示yr(n)延时τ后的信号,表示yr(n+τ)的共扼,为yr(n)的自相关函数,1≤m1≤Lh,1≤m2≤Lh,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的信号功率因子,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的信号功率因子,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的到达相位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的到达相位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的时延,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的时延,表示通过第i个发射天线发送的发送信号ti(n)经过hi,r的第m1径的时延后的信号,表示的共扼,表示ti(n)延时τ后的信号ti(n+τ)经过hi,r的第m2径的时延后的信号,表示ti(n)的自相关函数延时后的值,δ(τ)为vi,r(n)的自相关函数,δ(τ)是延时变量τ的冲击函数。
步骤6:在CP-Free MIMO-OFDM系统的接收端,对作傅立叶级数展开,得到yr(n)的周期自相关函数,记为 其中,k为循环频率,P表示OFDM信号的循环周期,亦为一个OFDM符号的长度,P=M+2L,M表示一个OFDM符号的数据长度,即表示OFDM信号的子载波数,在本实施例中取M=32,实际应用过程中还可取16、或32、或64、或128等,在本实施例中取π=3.14…,Λ为引入的中间符号,表示ti(n)的周期自相关函数延时后的值,δ(τ)δ(k)为vi,r(n)的周期自相关函数,δ(k)是变量k的冲击函数。
步骤7:在CP-Free MIMO-OFDM系统的接收端,根据通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径和第m2径的时间间隔相等与否这两种情况,将转化为 然后根据发送信号的周期自相关函数的能量分布(图3给出了发送信号的周期自相关函数的能量分布图),令得到其中,表示ti(n)的周期自相关函数。
步骤8:分析ti(n)的自相关函数 其中,ti(n-τ)表示ti(n)延时τ后的信号,ti(n-τ)H为ti(n-τ)的转置共轭,E(ti(n)ti(n-τ)H)表示ti(n)的自相关函数,表示ti(n)的方差,δ(τ)是变量τ的冲击函数,δ(τ-p)表示δ(τ)(向右)延时p后的值,δ(n-l)表示δ(n)延时l后的值,δ(n)是变量n的冲击函数,δ(τ+p)表示δ(τ)(向左)延时p后的值,在本实施例中取M表示一个OFDM符号的数据长度,即表示OFDM信号的子载波数,在本实施例中取M=32,从ti(n)的自相关函数中可以看出只有当τ=0,±p时,ti(n)的自相关函数才不为零,当τ=p时,此时,ti(n)的周期自相关函数为抽样脉冲函数,当k=P/L的整数倍时存在零点,为ti(n)的方差。在CP-Free MIMO-OFDM系统的接收端,假设需要估计通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比,那么通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数为然后根据yb(n)是否有通过指定的第a个发射天线发送的发送信号ta(n),将转化为其中,a∈[1,MT],b∈[1,MR],Si,b,m表示通过第i个发射天线发送的发送信号ti(n)到指定的第b个接收天线接收所经过的多径衰落信道hi,b的第m径的信号功率因子,τi,b,m表示通过第i个发射天线发送的发送信号ti(n)到指定的第b个接收天线接收所经过的多径衰落信道hi,b的第m径的时延,1≤J≤MT,且J≠a,SJ,b,m表示通过第J个发射天线发送的发送信号tJ(n)到指定的第b个接收天线接收所经过的多径衰落信道hJ,b的第m径的信号功率因子,Sa,b,m表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的第m径的信号功率因子,表示通过第J个发射天线发送的发送信号tJ(n)的周期自相关函数,表示通过指定的第a个发射天线发送的发送信号ta(n)的周期自相关函数,τJ,b,m表示通过第J个发射天线发送的发送信号tJ(n)到指定的第b个接收天线接收所经过的多径衰落信道hJ,b的第m径的时延,τa,b,m表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的第m径的时延。
步骤9:在CP-Free MIMO-OFDM系统的接收端,对通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比进行估计,具体过程为:
步骤9_1:要估计通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比,就必须考虑其干扰信道的影响,因此在的基础上,令ta(n)的循环频率等于ka且ka≠0,在排除ha,b的干扰信道的影响后,得到ta(n)的循环频率等于ka时yb(n)的周期自相关函数,记为 然后对的等号两边同时取绝对值,得到其中,表示循环频率k等于ka时通过指定的第a个发射天线发送的发送信号ta(n)的周期自相关函数,符号“||”为取绝对值符号,ξa,b表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的信号功率因子。
在本实施例中,步骤9_1中,排除ha,b的干扰信道的影响的具体过程为:令通过指定的第a个发射天线发送的发送信号ta(n)中的OFDM符号进行子载波映射时的第一映射集和第二映射集的长度为La,并令通过第J个发射天线发送的发送信号tJ(n)中的OFDM符号进行子载波映射时的第一映射集和第二映射集的长度为LJ,其中,a∈[1,MT],1≤J≤MT,且J≠a,La∈[1,L],LJ∈[1,L];然后令再使得通过第J个发射天线发送的发送信号tJ(n)的周期自相关函数处于零点,即使得且这样就只剩下ha,b本身的影响(若需要估计其它信道,可以选择不同的循环频率使得干扰信号的周期自相关函数处于零点),从而使得转化为排除了ha,b的干扰信道的影响,只得到所需的单输入单输出的ha,b。
步骤9_2:分析当k=0且τ=0时,转化为然后对的等号两边同时取绝对值,得到其中,表示k=0且τ=0时通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数,表示k=0且τ=0时通过第i个发射天线发送的发送信号ti(n)的周期自相关函数,ξ1,b表示通过第1个发射天线发送的发送信号t1(n)到指定的第b个接收天线接收所经过的多径衰落信道h1,b的信号功率因子,表示k=0且τ=0时通过第1个发射天线发送的发送信号t1(n)的周期自相关函数,ξ2,b表示通过第2个发射天线发送的发送信号t2(n)到指定的第b个接收天线接收所经过的多径衰落信道h2,b的信号功率因子,表示k=0且τ=0时通过第2个发射天线发送的发送信号t2(n)的周期自相关函数,表示通过第MT个发射天线发送的发送信号到指定的第b个接收天线接收所经过的多径衰落信道的信号功率因子,表示k=0且τ=0时通过第MT个发射天线发送的发送信号的周期自相关函数。
步骤9_3:对通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数进行统计量估计,将统计量估计后的值记为 其中,Nnum表示每个OFDM信号中包含的OFDM符号的个数,wb(n)表示通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线实际接收到的接收信号,wb(n+τ)表示wb(n)延时τ后的信号,wb *(n+τ)为wb(n+τ)的共轭。
步骤9_4:估计通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的信号功率因子ξa,b的值和高斯白噪声的噪声功率因子N的值,将ξa,b的估计值和N的估计值对应记为和 然后计算通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比的估计值,记为ρa,b,其中,表示对进行统计量估计后得到的值,表示对进行统计量估计后得到的值,和通过计算得到,表示ξ1,b的估计值,表示ξ2,b的估计值,表示的估计值,和通过计算得到。
在本实施例中,步骤9执行完毕后,重复执行步骤2至步骤9共Number次,得到ha,b上的盲信噪比的Number个估计值,然后将Number个估计值的平均值作为ha,b上的盲信噪比的最终估计值;其中,Number∈[500,1000]。
为了进一步说明本发明方法的性能,对本发明方法进行实验。
在实验中选用的对比方法为现有的基于CP OFDM的信噪比估计方法,其引用自Shunli Hong,Youming Li,Yao Gao,Jianding Yu.Blind SNR estimation for MIMO-OFDMsystem[C].CCWMC2011,ISBN:978-1-84919-505-8,2011:6-9.(洪顺利等,MIMO OFDM系统的盲信噪比估计,IET国际无线移动与计算通讯会议)。
图4给出了在每个OFDM信号中包含的OFDM符号的个数Nnum=1000,实际信噪比值分别为-10dB、-5dB、0dB、5dB、10dB、15dB的条件下,实际信噪比值、利用现有的基于CP OFDM的信噪比估计方法估计得到的信噪比值、利用本发明方法估计得到的盲信噪比值的比较示意图。图4中示出了通过指定的第1个发射天线发送的发送信号t1(n)到指定的第1个接收天线接收所经过的多径衰落信道h1,1上的信噪比的估计值和通过指定的第2个发射天线发送的发送信号t2(n)到指定的第1个接收天线接收所经过的多径衰落信道h2,1上的信噪比的估计值,从图4中可以看出,利用本发明方法估计得到的盲信噪比值较之于现有的基于CP OFDM的信噪比估计方法估计得到的信噪比值更接近于实际信噪比值。
为了进行算法相对性能比较,在此定义信噪比估计值的归一化均方误差NMSE和信噪比估计值的归一化偏差NBias作为性能评估函数, 其中,1≤q≤Nnum,表示第q个符号信噪比的估计值,ρq表示第q个符号信噪比的实际值。
图5给出了利用现有的基于CP OFDM的信噪比估计方法估计得到的信噪比值和利用本发明方法估计得到的盲信噪比值各自的归一化均方误差随信噪比的变化曲线图。从图5中可以看出,利用本发明方法得到的h1,1上的信噪比的估计值和h2,1上的信噪比的估计值各自的归一化均方误差随信噪比的变大而变小,趋势是一样的,且相比于现有的基于CPOFDM的信噪比估计方法误差更小。
图6给出了利用现有的基于CP OFDM的信噪比估计方法估计得到的信噪比值和利用本发明方法估计得到的盲信噪比值各自的归一化偏差随信噪比的变化曲线图。从图6中可以看出,利用本发明方法得到的h1,1上的信噪比的估计值和h2,1上的信噪比的估计值各自的归一化偏差随着信噪比的增大而变小,趋势是相同的,相比于现有的基于CP OFDM的信噪比估计方法误差更小。
Claims (8)
1.一种基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于包括以下步骤:
步骤1:设定CP-Free MIMO-OFDM系统中的各个信道采用多径衰落模型,即CP-FreeMIMO-OFDM系统中的各个信道为多径衰落信道;设定CP-Free MIMO-OFDM系统的发送端具有MT个发射天线且接收端具有MR个接收天线;其中,MT表示发射天线的个数,MT>1,MR表示接收天线的个数,MR>1;
步骤2:在CP-Free MIMO-OFDM系统的发送端,将接收的频域数据信号变换为MT个时域数据信号,每个时域数据信号为由多个长度一致的OFDM符号构成的OFDM信号,所有时域数据信号各自包含的OFDM符号的个数相同;然后对每个OFDM信号中的各个OFDM符号进行子载波映射,得到对应的子载波映射后的OFDM信号;接着对每个子载波映射后的OFDM信号先后进行数模转换和射频模块处理,得到对应的射频模块处理后的OFDM信号;再在每个射频模块处理后的OFDM信号上叠加一个用于消除码间干扰的匹配信号,得到对应的发送信号;
步骤3:CP-Free MIMO-OFDM系统的发送端将MT个发送信号各通过一个发射天线相互平行地传输给CP-Free MIMO-OFDM系统的接收端;
步骤4:在CP-Free MIMO-OFDM系统的接收端,每个接收天线接收到的接收信号由一个发送信号通过一个发射天线发送后经过MT个多径衰落信道后形成的具有不同时延的MT路信号组成,将通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被第r个接收天线接收到的接收信号记为yr(n),其中,1≤i≤MT,ti(n)表示通过第i个发射天线发送的发送信号,1≤r≤MR,1≤m≤Lh,Lh表示多径衰落信道的多径阶数,Si,r,m表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的信号功率因子,e为自然基数,j表示复数中的虚数单位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的到达相位,τi,r,m表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m径的时延,ti(n-τi,r,m)表示通过第i个发射天线发送的发送信号ti(n)经过hi,r的第m径的时延后的信号,N表示高斯白噪声的噪声功率因子,vi,r(n)表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r上的平稳高斯白噪声;
步骤5:在CP-Free MIMO-OFDM系统的接收端,根据自相关函数的定义,获取yr(n)的自相关函数,记为 ;其中,τ表示延时变量,E()表示数学期望,yr(n+τ)表示yr(n)延时τ后的信号,表示yr(n+τ)的共扼,为yr(n)的自相关函数,1≤m1≤Lh,1≤m2≤Lh,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的信号功率因子,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的信号功率因子,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的到达相位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的到达相位,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径的时延,表示通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m2径的时延,表示通过第i个发射天线发送的发送信号ti(n)经过hi,r的第m1径的时延后的信号,表示的共扼,表示ti(n)延时τ后的信号ti(n+τ)经过hi,r的第m2径的时延后的信号,表示ti(n)的自相关函数延时后的值,δ(τ)为vi,r(n)的自相关函数,δ(τ)是延时变量τ的冲击函数;
步骤6:在CP-Free MIMO-OFDM系统的接收端,对作傅立叶级数展开,得到yr(n)的周期自相关函数,记为 其中,k为循环频率,P表示OFDM信号的循环周期,亦为一个OFDM符号的长度,P=M+2L,M表示一个OFDM符号的数据长度,即表示OFDM信号的子载波数,Λ为引入的中间符号,表示ti(n)的周期自相关函数延时后的值,δ(τ)δ(k)为vi,r(n)的周期自相关函数,δ(k)是变量k的冲击函数;
步骤7:在CP-Free MIMO-OFDM系统的接收端,根据通过第i个发射天线发送的发送信号ti(n)到第r个接收天线接收所经过的多径衰落信道hi,r的第m1径和第m2径的时间间隔相等与否这两种情况,将转化为 然后根据发送信号的周期自相关函数的能量分布,令得到其中,表示ti(n)的周期自相关函数;
步骤8:在CP-Free MIMO-OFDM系统的接收端,假设需要估计通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比,那么通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数为然后根据yb(n)是否有通过指定的第a个发射天线发送的发送信号ta(n),将转化为其中,a∈[1,MT],b∈[1,MR],Si,b,m表示通过第i个发射天线发送的发送信号ti(n)到指定的第b个接收天线接收所经过的多径衰落信道hi,b的第m径的信号功率因子,τi,b,m表示通过第i个发射天线发送的发送信号ti(n)到指定的第b个接收天线接收所经过的多径衰落信道hi,b的第m径的时延,1≤J≤MT,且J≠a,SJ,b,m表示通过第J个发射天线发送的发送信号tJ(n)到指定的第b个接收天线接收所经过的多径衰落信道hJ,b的第m径的信号功率因子,Sa,b,m表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的第m径的信号功率因子,表示通过第J个发射天线发送的发送信号tJ(n)的周期自相关函数,表示通过指定的第a个发射天线发送的发送信号ta(n)的周期自相关函数,τJ,b,m表示通过第J个发射天线发送的发送信号tJ(n)到指定的第b个接收天线接收所经过的多径衰落信道hJ,b的第m径的时延,τa,b,m表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的第m径的时延;
步骤9:在CP-Free MIMO-OFDM系统的接收端,对通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b上的盲信噪比进行估计,具体过程为:
步骤9_1:在的基础上,令ta(n)的循环频率等于ka且ka≠0,在排除ha,b的干扰信道的影响后,得到ta(n)的循环频率等于ka时yb(n)的周期自相关函数,记为 然后对的等号两边同时取绝对值,得到其中,表示循环频率k等于ka时通过指定的第a个发射天线发送的发送信号ta(n)的周期自相关函数,符号“||”为取绝对值符号,ξa,b表示通过指定的第a个发射天线发送的发送信号ta(n)到指定的第b个接收天线接收所经过的多径衰落信道ha,b的信号功率因子;
步骤9_2:分析当k=0且τ=0时,转化为然后对的等号两边同时取绝对值,得到其中,表示k=0且τ=0时通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数,表示k=0且τ=0时通过第i个发射天线发送的发送信号ti(n)的周期自相关函数,ξ1,b表示通过第1个发射天线发送的发送信号t1(n)到指定的第b个接收天线接收所经过的多径衰落信道h1,b的信号功率因子,表示k=0且τ=0时通过第1个发射天线发送的发送信号t1(n)的周期自相关函数,ξ2,b表示通过第2个发射天线发送的发送信号t2(n)到指定的第b个接收天线接收所经过的多径衰落信道h2,b的信号功率因子,表示k=0且τ=0时通过第2个发射天线发送的发送信号t2(n)的周期自相关函数,表示通过第MT个发射天线发送的发送信号到指定的第b个接收天线接收所经过的多径衰落信道的信号功率因子,表示k=0且τ=0时通过第MT个发射天线发送的发送信号的周期自相关函数;
步骤9_3:对通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线接收到的接收信号yb(n)的周期自相关函数进行统计量估计,将统计量估计后的值记为 其中,Nnum表示每个OFDM信号中包含的OFDM符号的个数,wb(n)表示通过第i个发射天线发送的发送信号ti(n)经过MT个多径衰落信道后被指定的第b个接收天线实际接收到的接收信号,wb(n+τ)表示wb(n)延时τ后的信号,wb *(n+τ)为wb(n+τ)的共轭;
2.根据权利要求1所述的基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于所述的步骤2中,将接收的频域数据信号变换为MT个时域数据信号的具体过程为:
步骤2_1:对接收的频域数据信号进行串并转化,分成MT个子符号流;
步骤2_2:对每个子符号流进行无失真压缩编码,并在无失真压缩编码过程中在该子符号流中加入用于克服该子符号流在多径衰落信道中受到的干扰和噪声影响的冗余信息,得到对应的无失真压缩编码后的频域数据信号;
步骤2_3:对每个无失真压缩编码后的频域数据信号进行空时调制,得到对应的空时调制后的频域数据信号;
步骤2_4:对每个空时调制后的频域数据信号先后进行傅里叶逆变换和正交频分复用调制处理,得到对应的时域数据信号,即共得到MT个时域数据信号。
3.根据权利要求2所述的基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于所述的步骤2_2中,采用信道编码技术对每个子符号流进行无失真压缩编码。
4.根据权利要求2或3所述的基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于所述的步骤2_3中,利用调制器对每个无失真压缩编码后的频域数据信号进行空时调制。
5.根据权利要求1或2所述的基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于所述的步骤2中,对每个OFDM信号中的各个OFDM符号进行子载波映射的具体过程为:对于任一个OFDM符号,将该OFDM符号中的第1个采样点至第L个采样点按序构成第一映射集,并将该OFDM符号中的第L+1个采样点至第2L个采样点按序构成第二映射集;然后将该OFDM符号中的第1个采样点作为第一起始点,间隔p个采样点后将第一映射集中的各个采样点按序替换该OFDM符号中的第p+1个采样点至第p+L个采样点;并将该OFDM符号中的第L+1个采样点作为第二起始点,间隔p个采样点后将第二映射集中的各个采样点按序替换该OFDM符号中的第p+L+1个采样点至第p+2L个采样点,经过两次采样点替换操作后实现了该OFDM符号的子载波映射过程;其中,L表示第一映射集和第二映射集各自中包含的采样点的个数,即为第一映射集和第二映射集各自的长度,M表示一个OFDM符号的数据长度,即表示OFDM信号的子载波数,p=M。
6.根据权利要求1所述的基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于所述的步骤2中,在每个射频模块处理后的OFDM信号上叠加一个用于消除码间干扰的匹配信号,得到对应的发送信号的具体过程为:将在第i'个射频模块处理后的OFDM信号上叠加匹配信号后得到的发送信号记为ti'(n),ti'(n)=si'(n)+ci'(n);其中,1≤i'≤MT,n表示离散时间点,ti'(n)为离散时间信号,si'(n)表示第i'个射频模块处理后的OFDM信号,ci'(n)表示在si'(n)上叠加的用于消除码间干扰的匹配信号,C表示信道矩阵,Cinter为C的干扰部分,CH为C的厄米特矩阵。
7.根据权利要求5所述的基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于所述的步骤9_1中,排除ha,b的干扰信道的影响的具体过程为:令通过指定的第a个发射天线发送的发送信号ta(n)中的OFDM符号进行子载波映射时的第一映射集和第二映射集的长度为La,并令通过第J个发射天线发送的发送信号tJ(n)中的OFDM符号进行子载波映射时的第一映射集和第二映射集的长度为LJ,其中,a∈[1,MT],1≤J≤MT,且J≠a,La∈[1,L],LJ∈[1,L];然后令再使得通过第J个发射天线发送的发送信号tJ(n)的周期自相关函数处于零点,即使得且从而使得转化为排除了ha,b的干扰信道的影响。
8.根据权利要求1所述的基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法,其特征在于所述的步骤9执行完毕后,重复执行步骤2至步骤9共Number次,得到ha,b上的盲信噪比的Number个估计值,然后将Number个估计值的平均值作为ha,b上的盲信噪比的最终估计值;其中,Number∈[500,1000]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210209866.8A CN114785645B (zh) | 2022-03-03 | 2022-03-03 | 基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210209866.8A CN114785645B (zh) | 2022-03-03 | 2022-03-03 | 基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114785645A CN114785645A (zh) | 2022-07-22 |
CN114785645B true CN114785645B (zh) | 2023-05-30 |
Family
ID=82423910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210209866.8A Active CN114785645B (zh) | 2022-03-03 | 2022-03-03 | 基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114785645B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115442199B (zh) * | 2022-08-30 | 2024-04-16 | 西安电子科技大学 | 无cp的mimo-ofdm一体化信号设计及处理方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102255836A (zh) * | 2011-07-19 | 2011-11-23 | 宁波大学 | 一种基于mimo-ofdm信号循环平稳特性的盲信噪比估计方法 |
CN102790744A (zh) * | 2012-07-18 | 2012-11-21 | 宁波大学 | 一种正交频分复用系统中的信号干扰噪声比估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009130713A1 (en) * | 2008-04-23 | 2009-10-29 | Indian Institute Of Technology | A system and method for estimation and correction of carrier frequency offset in mimo-ofdm based wireless communications systems |
US8290462B2 (en) * | 2009-02-04 | 2012-10-16 | Agere Systems Llc | Receiver and method for estimating a plurality of estimated transfer functions corresponding to wireless channels in a multiple-input system |
-
2022
- 2022-03-03 CN CN202210209866.8A patent/CN114785645B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102255836A (zh) * | 2011-07-19 | 2011-11-23 | 宁波大学 | 一种基于mimo-ofdm信号循环平稳特性的盲信噪比估计方法 |
CN102790744A (zh) * | 2012-07-18 | 2012-11-21 | 宁波大学 | 一种正交频分复用系统中的信号干扰噪声比估计方法 |
Non-Patent Citations (1)
Title |
---|
基于相关性的OFDM盲识别与参数估计方法;张孟伯;杨晓静;张欣冉;樊斌斌;;火力与指挥控制(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114785645A (zh) | 2022-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3116090B1 (ja) | 通信システム、送信装置、受信装置、送信方法、受信方法、および、情報記録媒体 | |
US20110228884A1 (en) | Method and device for filtering channel estimation result in orthogonal frequency division multiplexing | |
JP2005124125A (ja) | Ofdm伝送方式におけるキャリアの配置方法、送信装置及び受信装置 | |
CN100521665C (zh) | 一种用于固定训练序列填充调制系统的迭代分解方法 | |
JP5102305B2 (ja) | 周波数選択性で高速フェージングチャネル上でシングルキャリア空間周波数のブロック符号化した伝送を行うためのシステム及び方法 | |
Nasir et al. | Performance of coarse and fine timing synchronization in OFDM receivers | |
CN104168227A (zh) | 一种应用于正交频分复用系统的载波同步的方法 | |
CN113259291B (zh) | 利用水声连续信号动态多普勒跟踪实现的相位补偿方法 | |
CN101729479B (zh) | 一种基于ofdm信号循环平稳特性的盲信道估计方法 | |
CN114785645B (zh) | 基于CP-Free MIMO-OFDM信号循环平稳特性的盲信噪比估计方法 | |
CN108650005B (zh) | 一种用于mimo-fbmc/oqam系统中的导频结构以及信道估计方法 | |
CN1913396B (zh) | 单/多载波共融数字广播系统通信方法 | |
CN102255836B (zh) | 一种基于mimo-ofdm信号循环平稳特性的盲信噪比估计方法 | |
CN106453192A (zh) | 一种基于屏蔽互补二元序列偶的符号同步方法及系统 | |
CN109194594A (zh) | 一种基于连续载波聚合的相位噪声抑制方法 | |
CN1933464B (zh) | 多天线接收系统中降低复杂度的均衡接收装置及接收方法 | |
CN1235362C (zh) | 估计下行链路信道的方法 | |
CN1984109A (zh) | 通信系统中的信道估计器及其信道估计方法 | |
Bhoyar et al. | Leaky least mean square (LLMS) algorithm for channel estimation in BPSK-QPSK-PSK MIMO-OFDM system | |
CN102143116A (zh) | 一种基于ofdm信号循环平稳特性的盲信噪比估计方法 | |
CN101447969A (zh) | 一种多带正交频分复用超宽带系统的信道估计方法 | |
CN1933467B (zh) | 宽带单载波/多载波均衡接收装置及其接收方法 | |
CN101102114A (zh) | 基于二值全通序列保护间隔填充的频域信道估计方法 | |
CN102821078B (zh) | 基于ofdm信号循环平稳特性的盲信道估计方法 | |
Ribeiro et al. | An OFDM Symbol Design for Reduced Complexity MMSE Channel Estimation. |
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 |