[go: up one dir, main page]

CN104836652B - 一种低信噪比下的空时分组编码mimo‑ofdm系统时频同步新方法 - Google Patents

一种低信噪比下的空时分组编码mimo‑ofdm系统时频同步新方法 Download PDF

Info

Publication number
CN104836652B
CN104836652B CN201510234000.2A CN201510234000A CN104836652B CN 104836652 B CN104836652 B CN 104836652B CN 201510234000 A CN201510234000 A CN 201510234000A CN 104836652 B CN104836652 B CN 104836652B
Authority
CN
China
Prior art keywords
mrow
msub
munderover
mfrac
sequence
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.)
Expired - Fee Related
Application number
CN201510234000.2A
Other languages
English (en)
Other versions
CN104836652A (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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum 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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201510234000.2A priority Critical patent/CN104836652B/zh
Publication of CN104836652A publication Critical patent/CN104836652A/zh
Application granted granted Critical
Publication of CN104836652B publication Critical patent/CN104836652B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00—Arrangements affording multiple use of the transmission path
    • H04L5/0001—Arrangements for dividing the transmission path
    • H04L5/0014—Three-dimensional division
    • H04L5/0023—Time-frequency-space
    • 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
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2657—Carrier synchronisation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2662—Symbol synchronisation
    • H04L27/2663—Coarse synchronisation, e.g. by correlation
    • 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/0204—Channel estimation of multiple channels

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radio Transmission System (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明是一种低信噪比下的STBC MIMO‑OFDM系统时频同步新方法，其特征在于：在低信噪比条件下，针对STBC MIMO‑OFDM系统很难获得精确的时频同步问题，本发明提出一种低信噪比条件下的STBC MIMO‑OFDM系统时频同步新方法。该方法首先利用一组复互补序列集和正交矩阵构造一个复正交互补序列，通过正交矩阵变换扩展该序列，从而获得合适长度的复正交互补序列；为提高该序列的自相关性和互相关性，训练序列是由该序列的共轭取反特性构造，并叠加于数据信号上。在接收端通过训练序列与数据信号相关获取良好的同步信号，并且整数频偏估计范围可以达到ε∈(‑N/4,N/4)。在低信噪比条件下，本发明提出的同步方法与常规方法比较，在‑15dB时可以获得更加精确的时间同步和较大的频偏范围，且计算复杂度更低。

Description

一种低信噪比下的空时分组编码MIMO-OFDM系统时频同步新 方法
技术领域
本发明涉及MIMO-OFDM技术领域，特别是一种低信噪比下的STBC MIMO-OFDM系统时频同步新方法。
背景技术
随着无线通信技术的快速发展，人们对无线通信技术的传输速率和稳定性以及可靠性具体越来越高的要求。现有的正交频分复用(OFDM)技术是一种多载波传输技术，其具有频谱利用率较高、抗多径衰落等特点。多输入多输出(MIMO)技术则是充分利用空间资源，实现多发射多接收的功能，在不增加频谱资源和发射天线功率的情况下，提高信道的容量。将MIMO与OFDM结合的系统具有信号传输稳定、频谱利用率高、高速传输率等特点，能够很好地满足下一代无线传输网络的要求，并且也是后期无线通信技术发展的必然趋势。其主要应用的领域包括：无线局域网(WLAN)、无线宽带接入(Wimax)、下一代移动无线通信系统(B3G/4G)等。
为了提高无线通信技术的传输速率和稳定性以及可靠性，在低信噪比条件下的空时分组编码(STBC)多输入多输出正交频分复用(MIMO-OFDM)系统中，MIMO-OFDM技术具有信号传输稳定、频谱利用率高、高速传输率等特点，并且采用STBC编码可以实现多发天线联合编码技术，天线分集可抵抗信道的衰落，提高无线通信链路的可靠性；通过STBC编码增加传输的空时冗余信息，提高无线通信传输的稳定性与可靠性。
在低信噪比条件下，现有的MIMO-OFDM同步算法时频同步算法性能均受到影响。因此，在低信噪比条件下，要想获得性能良好的时频同步性能，该时频同步算法是值得深入研究。Ali Rachini等人提出STBC MIMO OFDM系统时间同步算法，采用正交性良好的CAZAC序列和Walsh-Hadamard序列作为训练序列，利用本地训练序列与数据信号相关获取良好的时间同步性能。仿真结果表明在低信噪比条件下，时间同步性能较优。但没有考虑频偏估计对系统性能的影响。FAN Hui-li等人提出MIMO-OFDM系统时频同步算法，采用共轭结构的CAZAC序列作为训练序列，利用接收端数据信息自相关获取同步信息。仿真结果表明该方法在低信噪比条件下时频同步较好，但该方法不足是频偏估计范围较小。对比文献：[1]R.N.Ali,B.D.Ali,N.Fabienne and B.Bilal.Timing synchronization method forMIMO-OFDM system using orthogonal preamble[C].19th International Conferenceon Telecommunications(ICT 2012),2012:1-6.[2]H.L.Fan,J.F.Sun,P.Yang,D.S.Li.ARobust Timing and Frequency Synchronization Algorithm for HF MIMO OFDMSystems[C].IEEE Conference on Global Mobile Congress(GMC),2010:1-4.这些同步方法都是将训练序列插入到OFDM数据符号上的同步方法，其不足处是降低系统的频谱利用率和降低发射机的发射功率。
因此，为了克服现有MIMO-OFDM系统同步方法存在的不足，本发明提出一种低信噪比下的STBC MIMO-OFDM系统时频同步新算法，该同步方法的训练序列是叠加到OFDM数据符号上，其可以有效的提高频带利用率和发射机的发射功率。
发明内容
本发明的目的是在于提出一种低信噪下的STBC MIMO-OFDM系统时频同步新方法，从而克服现有的MIMO-OFDM系统同步方法的不足。本发明提出同步新方法，首先是设计一种正交性良好的复正交互补序列，其次利用该序列的共轭取反特性构造出新的训练序列，该训练序列具有良好的自相关性和较弱的互相关性，并且有优良的正交性。将训练序列叠加到一个完整的OFDM数据符号上成为发射信号，在接收端，利用接收端数据信号和本地训练序列相关可以获得精确的同步信息，利用接收端数据信号与本地训练序列相关运算求相位来获取整数频偏估计值。精确的时间同步和频率同步可以保证传输的数据信息正确解调。
为了实现上述目的，本发明提出的一种低信噪比下的STBC MIMO-OFDM系统时频同步新方法的创新点：
(1)本发明设计一种基于复正交互补序列的训练序列，该训练序列采用共轭取反特性提高训练序列的自相关性。
(2)本发明提出的新时频同步算法，时间同步算法是采用接收端数据信息和本地训练序列获取同步信息，用最大值确定精确的定时同步位置，在实现系统定时同步基础上，进行频率同步，且在时域可估计较大的频偏范围，从而保证传输的数据信息的正确解调。
本发明提出一种低信噪比下的STBC MIMO-OFDM系统时频同步新方法，该系统的天线设置NT个发射天线和NR个接收天线，其具体步骤如下：
步骤1：复正交互补序列设计：
(1)设计一个复互补序列集{A1,B1,C1,D1},其中A1＝[1,1,j,-j]，B1＝[1,1,-j,j]，和C2＝B1,D2＝-A1，其中是A1的共轭取反运算，-A1表示A1的负值，j表示虚数符号，j2＝-1；
(2)利用上述复互补序列集构造行向量为4，列向量为8的矩阵E4×8，其中式中(·)T表示矩阵转置；
(3)利用正交矩阵重新构造行向量为16，列向量为16新矩阵E16×16，其矩阵式中(·)T表示矩阵转置；
(4)利用正交矩阵重新构造行向量为32，列向量为32新矩阵E32×32，获得矩阵
(5)重复上述第(4)步，可以扩展矩阵从而获得更长的正交序列；
步骤2：利用上述复正交互补序列设计步骤，则可获取正交矩阵通过矩阵构造一组序列ai(n)，其中n的取值范围n∈[0,N/4-1]，i的取值范围i∈[1,N/4]，N表示子载波数；
步骤3：当i＝1时，构造长度为N/4的序列a1(n)，将序列a1(n)取共轭取反运算获取新的序列b1(n)，有：
式中：表示序列a1()序列取共轭运算；
将序列a1(n)和b1(n)构成长度为N/2的序列t1(n)，有：
步骤4：将序列t1(n)重复一次构造长度为N的序列c1(n)，有：
步骤5：将训练序列叠加到一个完整的OFDM符号上信号表示为：
式中：P表示发射机的总功率，β是功率分配因子，其取值范围0<β<1，定义为其中xi(n)为第i帧传输信号，x(n)表示传输信号，ci(n)为第i帧传输信号叠加的训练序列，c(n)为传输信号叠加的训练序列，E(●)表示对括号中序列做均值计算；
步骤6：利用接收端接收到的本地训练序列和数据信息相关获取时间同步，定时度量函数表达式为：
其中：
式中：Nt为发射天线的数目，Nr为接收天线的数目，Ng为循环前缀的长度，表示接收端数据信息取共轭运算，n是采样点的数目，ci(n)表示发射天线上的训练序列，d为整数值，其表示本地训练序列和接收端数据信息的相对滑动位置，m表示循环次数；
为了方便数学分析，假设发射机的总功率P为1，且训练序列与数据信息及信道噪声无关。
将(8)式带入到(5)式得到：
式中：表示训练序列的发射功率。当接收端数据信号的叠加训练序列与本地训练序列数据块对齐时，B{P(d)}达到最大值。由于公式(5)中定时度量函数R(d)是归一化作用，则M(d)达到最大值时，即可判断此刻的d是OFDM符号开始位置，同步算法表示为：
式中：表示定时同步估计值；
步骤6：由(9)式可以看出同步精确位置和训练序列功率分配因子有关，当功率分配因子越大，系统定时同步性能越好；
步骤7：采用接收端数据信息自相关性，获取小数频偏估计，其数学公式表示为：
式中：real(·)表示取实部运算，imag(·)表示取虚部运算，表示小数频偏估计值，其估计范围：
步骤8：是将估计出的小数频率偏移补偿后进行频率的整数频偏估计，通过本地训练序列与接收端数据信号相关运算求相位来获取整数频偏估计值，表达式：
式中：表示整数频偏估计，其估计范围是：
附图说明
为了更加详细说明一种低信噪比下的STBC MIMO-OFDM系统时频同步新方法，针对本发明所涉及到的图表进行详细说明。
图1是本发明STBC MIMO-OFDM系统框图。图中，STBC MIMO-OFDM系统框图主要由两部分组成，在发送端主要有：输入的数据信息经过MPSK调制102、空时编码103、IFFT调制104等模块；在接收端主要有：同步和信道估计112、FFT解调116、空时解码117、MPSK解调118等模块。
图2是本发明训练序列叠加到OFDM符号上的结构框图。图中，训练序列是叠加到各发送天线OFDM数据符号之上。
图3是本发明复正交互补序列的构造示意图。图中，复正交互补序列是利用复互补序列集构造矩阵，并通过正交矩阵和矩阵变换扩展合适长度的正交序列。
图4是本发明基于复正交互补序列的训练序列构造示意图。图中，采用合适长度的复正交互补序列，并通过共轭取反和重复特性重新构造基于复正交互补序列的训练序列。
图5是本发明STBC MIMO-OFDM系统时频同步算法流程图。图中，利用接收端的本地训练序列和数据信息相关获取定时同步信息，用最大值判断定时同步位置，随后进行小数频偏估计和补偿，最后采用本地训练序列与接收端数据信号相关运算求相位获取整数频偏估计值。
图6是本发明不同功率分配因子的误码率仿真图。图中，横坐标表示信噪比(SNR)，纵坐标表示误码率性能(BER)，β表示功率分配因子，由图中看出，随着功率分配因子β的增大，系统误码率性能降低。
图7是本发明不同功率分配因子的同步正确概率图。图中，横坐标表示信噪比(SNR)，纵坐标表示时间同步正确概率，β表示功率分配因子，由图中看出，随着功率分配因子β的增大，系统时间同步正确概率性能提高。
图8是本发明时间同步正确概率性能比较图。图中，横坐标表示信噪比(SNR)，纵坐标表示时间同步正确概率，在文献[1]中是Ali Rachini等人采用正交性良好的Walsh-Hadamard序列和CAZAC序列分别作为训练序列的两种时间同步算法，在文献[2]中是FanHui-li等人采用移位正交性的CAZAC序列时间同步算法，由图看出，改进算法的时间同步正确概率的性能明显优于文献[1]和[2]。对比文献：[1]R.N.Ali,B.D.Ali,N.Fabienne andB.Bilal.Timing synchronization method for MIMO-OFDM system using orthogonalpreamble[C].19th International Conference on Telecommunications(ICT 2012),2012:1-6.[2]H.L.Fan,J.F.Sun,P.Yang,D.S.Li.A Robust Timing and FrequencySynchronization Algorithm for HF MIMO OFDM Systems[C].IEEE Conference onGlobal Mobile Congress(GMC),2010:1-4.
图9是本发明频偏估计均方误差性能比较图。图中，横坐标表示信噪比(SNR)，纵坐标表示频偏估计均方误差(MSE)，在文献[2]中是Fan Hui-li等人提出采用接收信号自相关获取小数频偏估计算法，由图看出，改进算法的频偏估计均方误差性能更优于文献[2]。对比文献：[2]H.L.Fan,J.F.Sun,P.Yang,D.S.Li.A Robust Timing and FrequencySynchronization Algorithm for HF MIMO OFDM Systems[C].IEEE Conference onGlobal Mobile Congress(GMC),2010:1-4.
具体实施方式
下面结合附图对本发明专利做进一步详细说明。
图1是本发明具有NT根发射天线和NR根接收天线的STBC MIMO-OFDM系统基本框图。该系统框图主要是由发射机和接收机组成，在发射机端，主要包括的模块有数据比特流模块101，调制模块(MPSK)102、空时编码模块103、傅里叶逆变换(IFFT)模块104、并串变换模块105、并行训练序列模块106、功率分配模块107、叠加模块108、插入保护间隔模块109和数模转换(D/A)模块。在接收机端，主要包括的模块有模数转换(A/D)模块111、同步和信道估计模块112、训练序列与数据分离模块113、去保护间隔模块114、串并变换模块115、傅里叶变换(FFT)模块116，空时解码模块117，解调(MPSK)模块118，数据比特流模块119。本发明主要是在发射机端采用将并行训练序列模块106叠加到OFDM数据符号上，这样可以提高系统的频谱利用率和传输效率，节省带宽资源。在之前采用空时编码模块103可以实现多发天线编码技术，天线分集还可以抵抗信道的衰落，从而提高通信链路的可靠性，此外，还增加传输的空时冗余信息，提高无线通信传输的稳定性与可靠性。在接收机端，同步和信道估计模块主要是确定FFT的窗口的起始位置和保证各子载波的正交性，随后进行空时解码模块117和解调(MPSK)模块118实现对OFDM数据模块的正确解调。
图2是本发明训练序列叠加到OFDM符号上的结构框图，将基于复正交互补的训练序列叠加到各发射天线上一个完整的OFDM数据符号上，201是训练序列c1(n)，202是训练序列c2(n)。
图3是本发明复正交互补序列的构造示意图，设计复正交互补序列步骤：301首先是设计一个复互补序列集，302利用上述复互补序列集构造矩阵E4×8，303利用正交矩阵H2×4构造矩阵E16×16，304利用正交矩阵F2×2重新构造矩阵E32×32，305重复上述步骤第四步可以获得更长的正交序列。
图4是本发明基于复正交互补序列的训练序列构造示意图，复正交互补序列的训练序列的构造方法：401利用复正交互补序列的设计步骤获得正交矩阵并利用该矩阵构造一组序列ai(n),i＝0,…,N/4，402当i＝1时举例，将长度为N/4的序列a1(n)作共轭取反运算获取新的序列b1(n)，403将序列a1(n)和b1(n)组合构造成长度为N/2的序列t1(n)，404将序列t1(n)重复一次构造长度为N的训练序列c1(n)。
图5是本发明STBC MIMO-OFDM系统时频同步算法流程图，501是接收信号，502利用接收端的本地训练序列和数据信息相关获取定时同步，503用最大值判断定时同步位置，504小数频偏估计，505频偏补偿，506利用本地训练序列与接收端数据信号相关运算求相位来获取整数频偏估计值。
图6是本发明不同功率分配因子的误码率仿真图，图7是不同功率分配因子的同步正确概率仿真图，本发明的主要仿真参数设置：发射天线数目NT＝2，接收天线数目NR＝2，蒙卡罗仿真次数是10000次，调制方式采用MPSK，子载波数N＝1024，循环前缀的长度为Ng＝128，选取信道环境是多径衰落信道，频偏设置为ε＝10.35。通过系统误码率性能和同步正确概率性能权衡叠加训练序列功率分配因子β的最优功率分配因子。当在误码率较小及同步正确概率较高时，确定此时的功率分配因子是最优值。由图6和图7可以看出，在相同信噪比下，随着功率分配因子增大β，误码率性能降低，但时间同步正确概率却增强。为了算法性能折中，本文选取功率分配因子与理论值相近的值β＝0.3作为后续的仿真分析将在此基础上进行。
图8是本发明时间同步正确概率性能仿真图，本发明的主要仿真参数设置：发射天线数目NT＝2，接收天线数目NR＝2，蒙卡罗仿真次数是10000次，调制方式采用MPSK，子载波数N＝1024，循环前缀的长度为Ng＝128，选取信道环境是多径衰落信道，频偏设置为ε＝10.35，功率分配因子设置为β＝0.3。由图可以看出，在SNR＝-19dB时，本文正确概率可达到90％以上。与对比文献[1]和[2]比较，本文提出的时间同步方法性能较好。
图9是本发明频偏估计均方误差性能仿真图，本发明的主要仿真参数设置：发射天线数目NT＝2，接收天线数目NR＝2，蒙卡罗仿真次数是10000次，调制方式采用MPSK，子载波数N＝1024，循环前缀的长度为Ng＝128，选取信道环境是多径衰落信道，频偏设置为ε＝10.35，功率分配因子设置为β＝0.3。由图可看出，本文小数频偏估计算法的频偏范围为|εf|≤1具有较低的均方误差，克服了常规算法频偏范围是|εf|≤0.5的缺陷。在相同信噪比条件下，本文提出的频偏估计算法与对比文献[2]的算法性能比较。具有较低的均方误差。

Claims (1)

1.一种低信噪比下的空时分组编码MIMO-OFDM系统时频同步新方法，其特征在于：
步骤1：复正交互补序列设计：
(1)设计一个复互补序列集{A1,B1,C1,D1},其中A1＝[1,1,j,-j]，B1＝[1,1,-j,j]，和C2＝B1,D2＝-A1，其中是A1的共轭取反运算，-A1表示A1的负值，j表示虚数符号，j2＝-1；
(2)利用上述复互补序列集构造行向量为4，列向量为8的矩阵E4×8，其中式中(·)T表示矩阵转置；
(3)利用正交矩阵重新构造行向量为16，列向量为16新矩阵E16×16，其矩阵
(4)利用正交矩阵重新构造行向量为32，列向量为32新矩阵E32×32，获得矩阵
(5)重复上述第(4)步方法，改变行向量大小，可以获得更长的正交序列；
步骤2：利用上述复正交互补序列设计步骤，则可获取正交矩阵通过矩阵构造一组序列ai(n)，其中n的取值范围n∈[0,N/4-1]，i的取值范围i∈[1,N/4]，N表示子载波数；
步骤3：当i＝1时，构造长度为N/4的序列a1(n)，将序列a1(n)取共轭取反运算获取新的序列b1(n)，有：
<mrow> <msub> <mi>b</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msubsup> <mi>a</mi> <mn>1</mn> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>n</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>4</mn> <mo>-</mo> <mn>1</mn> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式中：表示序列a1()序列取共轭运算；
将序列a1(n)和b1(n)构成长度为N/2的序列t1(n)，有：
<mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>n</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>4</mn> <mo>-</mo> <mn>1</mn> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>b</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>n</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mi>N</mi> <mo>/</mo> <mn>4</mn> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>-</mo> <mn>1</mn> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
步骤4：将序列t1(n)重复一次构造长度为N的序列c1(n)，有：
<mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>n</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>-</mo> <mn>1</mn> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>n</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>,</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
步骤5：将训练序列叠加到一个完整的OFDM符号上信号表示为：
<mrow> <msub> <mi>r</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;beta;</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <msub> <mi>x</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msqrt> <mrow> <mi>P</mi> <mi>&amp;beta;</mi> </mrow> </msqrt> <msub> <mi>c</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>n</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mo>...</mo> <mi>N</mi> <mo>+</mo> <mi>N</mi> <mi>g</mi> <mo>-</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中：P表示发射机的总功率，β是功率分配因子，其取值范围0<β<1，定义为其中xi(n)为第i帧传输信号，x(n)表示传输信号，ci(n)为第i帧传输信号叠加的训练序列，c(n)为传输信号叠加的训练序列，E(●)表示对括号中序列做均值计算；
步骤6：利用接收端接收到的本地训练序列和数据信息相关获取时间同步，定时度量函数表达式为：
<mrow> <mi>M</mi> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> <mrow> <mo>|</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
其中：
<mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>T</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>R</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>1</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>/</mo> <mn>4</mn> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mo>&amp;lsqb;</mo> <msubsup> <mi>r</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>d</mi> <mo>+</mo> <mfrac> <mi>N</mi> <mn>4</mn> </mfrac> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mi>m</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mi>N</mi> <mn>4</mn> </mfrac> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mi>m</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>r</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>d</mi> <mo>+</mo> <mi>m</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>T</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>R</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mo>|</mo> <msub> <mi>r</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式中：NT为发射天线的数目，NR为接收天线的数目，Ng为循环前缀的长度，表示接收端数据信息rj()取共轭运算，n是采样点的数目，ci(n)表示发射天线上的训练序列，d为整数值，其表示本地训练序列和接收端数据信息的相对滑动位置，m表示循环次数；
步骤6：由于公式(5)中定时度量函数R(d)是归一化作用，则M(d)达到最大值时，即可判断此刻的d是OFDM符号开始位置，同步算法表示为：
<mrow> <mover> <mi>&amp;tau;</mi> <mo>^</mo> </mover> <mo>=</mo> <munder> <mi>argmax</mi> <mi>d</mi> </munder> <mo>{</mo> <mo>|</mo> <mi>M</mi> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>|</mo> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
式中：表示定时同步估计值；
步骤7：采用接收端数据信息自相关性，获取小数频偏估计，其数学公式表示为：
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <mover> <mi>&amp;tau;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>T</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>R</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>r</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mover> <mi>&amp;tau;</mi> <mo>^</mo> </mover> <mo>+</mo> <mi>n</mi> <mo>+</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> <msub> <mi>r</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mi>l</mi> <mi>e</mi> <mo>{</mo> <mfrac> <mrow> <mi>r</mi> <mi>e</mi> <mi>a</mi> <mi>l</mi> <mo>&amp;lsqb;</mo> <mi>F</mi> <mrow> <mo>(</mo> <mover> <mi>&amp;tau;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mi>i</mi> <mi>m</mi> <mi>a</mi> <mi>g</mi> <mo>&amp;lsqb;</mo> <mi>F</mi> <mrow> <mo>(</mo> <mover> <mi>&amp;tau;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
式中：real(·)表示取实部运算，imag(·)表示取虚部运算，表示小数频偏估计值，其估计范围：
步骤8：将估计出的小数频率偏移补偿后进行频率的整数频偏估计，通过本地训练序列与接收端数据信号相关运算求相位来获取整数频偏估计值，表达式：
<mrow> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mi>N</mi> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mi>l</mi> <mi>e</mi> <mo>&amp;lsqb;</mo> <mi>Q</mi> <msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <mrow> <mi>d</mi> <mo>=</mo> <mover> <mi>&amp;tau;</mi> <mo>^</mo> </mover> </mrow> </msub> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>Q</mi> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>T</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>R</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>1</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>/</mo> <mn>2</mn> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>r</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>d</mi> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mi>m</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>c</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mi>m</mi> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>r</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>d</mi> <mo>+</mo> <mi>n</mi> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mi>m</mi> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>c</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mfrac> <mi>N</mi> <mn>2</mn> </mfrac> <mi>m</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
式中：表示整数频偏估计，其估计范围是：
CN201510234000.2A 2015-05-11 2015-05-11 一种低信噪比下的空时分组编码mimo‑ofdm系统时频同步新方法 Expired - Fee Related CN104836652B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510234000.2A CN104836652B (zh) 2015-05-11 2015-05-11 一种低信噪比下的空时分组编码mimo‑ofdm系统时频同步新方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510234000.2A CN104836652B (zh) 2015-05-11 2015-05-11 一种低信噪比下的空时分组编码mimo‑ofdm系统时频同步新方法

Publications (2)

Publication Number Publication Date
CN104836652A CN104836652A (zh) 2015-08-12
CN104836652B true CN104836652B (zh) 2017-12-08

Family

ID=53814310

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510234000.2A Expired - Fee Related CN104836652B (zh) 2015-05-11 2015-05-11 一种低信噪比下的空时分组编码mimo‑ofdm系统时频同步新方法

Country Status (1)

Country Link
CN (1) CN104836652B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11218236B2 (en) 2016-06-01 2022-01-04 Qualcomm Incorporated Time division multiplexing of synchronization channels
US11563505B2 (en) 2016-06-01 2023-01-24 Qualcomm Incorporated Time division multiplexing of synchronization channels
CN106385308B (zh) * 2016-09-19 2019-05-24 中国人民解放军海军航空大学 基于k-s检测的stbc-ofdm信号盲识别方法
CN107070840B (zh) * 2017-06-07 2019-10-22 北京科技大学 一种符号定时同步方法及其系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101494487A (zh) * 2009-02-23 2009-07-29 电子科技大学 一种虚拟miso系统中的频率偏移补偿方法
CN103227769A (zh) * 2013-05-06 2013-07-31 西南石油大学 一种新的降低stbc mimo-ofdm系统峰均比的方法
CN104363192A (zh) * 2014-10-21 2015-02-18 江苏中兴微通信息科技有限公司 一种兼容多种帧格式的mimo通信系统的接收方法和装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070127360A1 (en) * 2005-12-05 2007-06-07 Song Hyung-Kyu Method of adaptive transmission in an orthogonal frequency division multiplexing system with multiple antennas
US20090080558A1 (en) * 2007-09-25 2009-03-26 An John F Space-time mimo wireless system based on feedback optimum weight design

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101494487A (zh) * 2009-02-23 2009-07-29 电子科技大学 一种虚拟miso系统中的频率偏移补偿方法
CN103227769A (zh) * 2013-05-06 2013-07-31 西南石油大学 一种新的降低stbc mimo-ofdm系统峰均比的方法
CN104363192A (zh) * 2014-10-21 2015-02-18 江苏中兴微通信息科技有限公司 一种兼容多种帧格式的mimo通信系统的接收方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MIMO－OFDM 系统同步技术研究进展;罗仁泽等;《电讯技术》;20140828;1171-1176页 *
叠加共轭对称训练序列的OFDM系统时频同步方法;罗仁泽等;《电讯技术》;20140928;1231-1236页 *

Also Published As

Publication number Publication date
CN104836652A (zh) 2015-08-12

Similar Documents

Publication Publication Date Title
CN101494528B (zh) 发射分集块传输系统的训练序列设计及其信道估计方法
US8548101B2 (en) Demodulation method and demodulator for orthogonal frequency multiplexing—multiple input multiple output system
WO2017219389A1 (zh) 大规模mimo系统中实现完美全向预编码的同步信号和信号的发送与接收方法
CN101005475A (zh) 正交频分复用通信中时间和频率同步的方法及系统
CN104836652B (zh) 一种低信噪比下的空时分组编码mimo‑ofdm系统时频同步新方法
CN111245766B (zh) 基于频域双分量扩展加权傅里叶变换的计算分集方法
CN1964341B (zh) 多入多出-正交频分复用系统的接收端的频偏估计方法
CN102098262B (zh) 无线通信中数据发射方法及系统、接收方法及系统
CN101355543A (zh) 基于正交训练序列的mimo-scfde系统信道估计方法
CN107276726B (zh) 一种Massive MIMO FBMC波束空时编码下行传输方法
CN102970271A (zh) 一种联合估计载波频偏的频率同步方法
CN101291311B (zh) 多输入多输出正交频分复用系统的同步实现方法及装置
CN104796372A (zh) 一种低实现复杂度的lte系统下行辅同步信号检测方法
CN103259757B (zh) 一种有效的mimo-ofdm系统的时间与频率同步新方法
CN107332606A (zh) 基于双采样的leo系统差分空时正交频分复用编码方法
CN102377699B (zh) 一种mu-mimo系统的信道估计方法和装置
Bhoyar et al. Leaky least mean square (LLMS) algorithm for channel estimation in BPSK-QPSK-PSK MIMO-OFDM system
Li et al. A generalized analytical solution to channel estimation with intersymbol interference cancelation and co-channel interference cancelation for single input single output/multiple input single output digital terrestrial multimedia broadcasting systems
CN105049095A (zh) Sc-mimo系统双流三天线或者四天线的分集收发方法及装置
CN103929395A (zh) 基于恒定包络调制的ofdm系统频偏同步方法
CN102014101B (zh) 基于多发送天线技术的ofdm同步方法
CN105187186B (zh) 低信噪比zcz序列及互模糊函数的分布式mimo-ofdm系统同步方法
CN101494490B (zh) 基于反馈的随机分布式自组织通信方法及系统
CN101945068A (zh) 一种低发射功率的单载波频域均衡系统接收机检测方法
CN104468053A (zh) 一种基于ngb-w的多天线复用传输系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171208

Termination date: 20180511

CF01 Termination of patent right due to non-payment of annual fee