CN115276938A - 用于多天线子载波索引调制ofdm的导频传输方法 - Google Patents

用于多天线子载波索引调制ofdm的导频传输方法 Download PDF

Info

Publication number
CN115276938A
CN115276938A CN202210859170.XA CN202210859170A CN115276938A CN 115276938 A CN115276938 A CN 115276938A CN 202210859170 A CN202210859170 A CN 202210859170A CN 115276938 A CN115276938 A CN 115276938A
Authority
CN
China
Prior art keywords
pilot
signal
antenna
channel estimation
frequency
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.)
Pending
Application number
CN202210859170.XA
Other languages
English (en)
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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202210859170.XA priority Critical patent/CN115276938A/zh
Publication of CN115276938A publication Critical patent/CN115276938A/zh
Pending legal-status Critical Current

Links

Images

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/0003—Two-dimensional division
    • H04L5/0005—Time-frequency
    • H04L5/0007—Time-frequency the frequencies being orthogonal, e.g. OFDM(A), DMT
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0054—Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • 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/0202—Channel estimation
    • H04L25/024—Channel estimation channel estimation algorithms
    • H04L25/0256—Channel estimation using minimum mean square error criteria
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/3405—Modifications of the signal space to increase the efficiency of transmission, e.g. reduction of the bit error rate, bandwidth, or average power
    • H04L27/3411—Modifications of the signal space to increase the efficiency of transmission, e.g. reduction of the bit error rate, bandwidth, or average power reducing the peak to average power ratio or the mean power of the constellation; Arrangements for increasing the shape gain of a signal set
    • 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/003—Arrangements for allocating sub-channels of the transmission path
    • H04L5/0048—Allocation of pilot signals, i.e. of signals known to the receiver
    • 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)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Artificial Intelligence (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明属于通信抗干扰技术领域，涉及一种用于多天线子载波索引调制OFDM的导频传输方法。本发明根据一定的导频数据能量比在频域将导频稀疏叠加在调制信号上，与导频单独占用子载波的放置方式相比，减少了导频开销，提高了系统的频谱效率，且通过迭代更新也改善了系统的BER性能。

Description

用于多天线子载波索引调制OFDM的导频传输方法
技术领域
本发明属于通信抗干扰技术领域，具体涉及一种用于多天线下子载波索引调制正交频分复用(Orthogonal Frequency Division Multiplexing with Subcarrier IndexModulation,OFDM-SIM)系统中的导频传输方法。
背景技术
索引调制(Index Modulation,IM)技术是一种无线通信领域的新型传输技术，其本质是用额外的比特选择激活的索引资源来传输信息比特，OFDM-SIM系统利用子载波作为索引资源，选择部分子载波进行激活传输星座符号，从索引域和数字调制域两个维度进行信息的传递。一方面，OFDM-SIM系统通过子载波的激活位置来携带额外的信息，以弥补静默子载波造成的频谱效率的损失，另一方面，OFDM-SIM系统子载波稀疏映射结构具有更低的PAPR，静默部分子载波可以有效降低载波间干扰。因此，相较于传统的OFDM系统，OFDM-SIM系统具有更强的鲁棒性以及更高的频偏忍耐度，在未来可见光通信和大规模MIMO通信中具有良好的应用前景。
为了在接收端恢复发送的比特信息，接收机必须获取信道状态信息。通常，使用发射机和接收机均已知的导频信号进行信道估计。OFDM系统中，传统的导频放置采用频分复用的方式，额外占用系统子载波，因此，随着通信技术不断发展，系统传输速率不断提高，通信节点不断增长，在系统天线数增多时，导频开销成为了频谱消耗的主要部分。
发明内容
本发明基于提高系统频谱效率和信道估计性能的目的，提出一种用于多天线OFDM-SIM系统的索引调制数据与导频配置及收发机设计方案。该方案在发射端通过将稀疏导频与调制信号在频域叠加后通过天线发送，在接收端利用导频进行信道估计获取信道信息，通过最大似然检测进行解调。该系统在降低导频传输所需资源消耗的同时，通过高效信道估计与检测方法保证接收端的信号检测性能。
为了便于理解，对多天线OFDM-SIM系统进行以下定义：
对于一个NT×NR的多天线OFDM-SIM系统，NT是发送天线的数量，NR是接收天线的数量。在接收端设置信道估计迭代次数为I。对于每一根发射天线，设每一帧有Ns个OFDM符号，N为OFDM系统子载波数，将N个子载波平均分为G个子块，每个子块包含L＝N/G个子载波，选择K个子载波进行激活，将此系统记作(L,K)OFDM-SIM系统。发射天线编号用t(t＝1,2,…,NT)表示，接收天线编号用r(r＝1,2,…,NR)表示，一根天线上的一个OFDM符号编号用j(j＝1,2,…,Ns)表示，子块编号用g(g＝1,2,…,G)表示，l(l＝1,2,…,L)表示子块的子载波编号。
对于天线t上的一个OFDM符号的任意一个子块g传输的比特由索引比特和调制比特组成，从L个子载波中选择K个进行激活的组合数为
Figure BDA0003757223260000021
则需要
Figure BDA0003757223260000022
(
Figure BDA0003757223260000023
表示向下取整)位索引比特，有K个激活子载波传输M阶星座符号，所以调制比特数为b2＝K log2 M，一个块传输的总比特数为b＝b1+b2。
本发明的技术方案是：
用于多天线子载波索引调制OFDM的导频传输方法，包括：
发射端：
步骤1：产生信息比特流。以第t根发射天线的第j个OFDM符号的子块g为例，产生信息比特
Figure BDA0003757223260000024
通过索引比特
Figure BDA0003757223260000025
选择激活的K个子载波，用以传输调制比特
Figure BDA0003757223260000026
进行M阶星座调制的符号，剩余的(L-K)个子载波保持静默，不传送信息。
步骤2：星座符号映射和子载波激活选择。对于发射天线t的第j个OFDM符号的子块g，调制比特为
Figure BDA0003757223260000027
经过M阶星座映射后得到调制符号向量
Figure BDA0003757223260000028
其中
Figure BDA0003757223260000029
Φ表示M阶星座符号集合。再通过索引比特
Figure BDA00037572232600000210
进行激活子载波选择，得到激活子载波位置向量
Figure BDA00037572232600000211
最后，通过映射规则：
Figure BDA00037572232600000212
将星座调制符号放置到对应的激活子载波上，得到子块发射符号
Figure BDA00037572232600000213
步骤3：生成导频图案。在系统发送端，同一根天线上的OFDM符号使用相同的导频序列，发射天线t使用导频序列
Figure BDA0003757223260000031
(其中V为导频数量，v＝1,2,…,V)。采用一种稀疏导频放置方案，导频等间隔放置，一帧包含Ns个OFDM符号，设置导频位置集合为
Figure BDA0003757223260000032
其中
Figure BDA0003757223260000033
表示发射天线t的导频位置，其中
Figure BDA0003757223260000034
表示发射天线t的第j个OFDM符号的第v个导频位置)。则发射天线t的第j个OFDM符号的第g个子块导频信号为
Figure BDA0003757223260000035
令l'＝l+(g-1)*L，表示系统子载波编号，其中
Figure BDA0003757223260000036
表示为：
Figure BDA0003757223260000037
则发射天线t的第j个OFDM符号的导频信号为
Figure BDA0003757223260000038
发射天线t的导频信号为
Figure BDA0003757223260000039
则生成的导频信号为
Figure BDA00037572232600000310
步骤4：导频与数据叠加。将导频信号叠加在数据符号上，不单独分配频率资源。假设导频符号能量分配因子为Ep，调制符号能量分配因子Ed，为满足能量归一化，满足：G·K·Ed+V·Ep＝N。则发射天线t的第j个OFDM符号的第g个子块频域信号为
Figure BDA00037572232600000311
其中
Figure BDA00037572232600000312
则发射天线t的第j个OFDM-SIM频域信号为
Figure BDA00037572232600000313
发射天线t的频域信号为
Figure BDA00037572232600000314
则发射端的频域发射信号为
Figure BDA00037572232600000315
步骤5：频域-时域变换。将步骤4得到的发送符号向量S进行IFFT变换，得到时域信号
Figure BDA00037572232600000316
其中
Figure BDA00037572232600000317
Figure BDA00037572232600000318
表示发射天线t上第j个OFDM-SIM时域符号，
Figure BDA00037572232600000319
其中
Figure BDA00037572232600000320
Figure BDA00037572232600000321
是
Figure BDA00037572232600000322
的第n(n＝1,2,…,N)个元素。再在时域信号
Figure BDA00037572232600000323
上添加循环前缀(CP)得到最终时域发送符号
Figure BDA00037572232600000324
通过发射天线发送。
接收端：
步骤6：时域-频域转换。时域接收信号为
Figure BDA0003757223260000041
其中
Figure BDA0003757223260000042
进行去循环前缀(CP)操作，
Figure BDA0003757223260000043
再进行FFT操作得到频域待检测信号
Figure BDA0003757223260000044
Figure BDA0003757223260000045
表示接收天线r上的频域接收信号，其中
Figure BDA0003757223260000046
变换公式为：
Figure BDA0003757223260000047
Figure BDA0003757223260000048
是
Figure BDA0003757223260000049
的第m个元素。接收符号向量可以用矩阵表示为：
Y＝HS+W
Figure BDA00037572232600000410
其中，H为信道频域响应矩阵，S为频域发送符号，W为加性复高斯白噪声向量。
步骤7：初始信道估计。按照已知的导频位置A，从第r根接收天线的接收信号Yr中取出对应的检测信号
Figure BDA00037572232600000411
利用频域检测信号
Figure BDA00037572232600000412
导频Pt和导频位置At，使用某种信道估计方法得到接收天线r与发射天线t之间的初始信道估计值
Figure BDA00037572232600000413
然后得到初始信道估计
Figure BDA00037572232600000414
若迭代次数I＝0，最终信道估计值
Figure BDA00037572232600000415
转到步骤11，否则，设置i＝0，表示当前迭代次数，令
Figure BDA00037572232600000416
转到步骤8。
步骤8：信号粗检测。根据信道估计值
Figure BDA00037572232600000417
频域接收信号Y和导频P进行信号粗检测得到检测结果
Figure BDA0003757223260000051
步骤9：干扰抵消操作。为了更加准确得到接收天线r与发射天线t之间的信道估计，根据步骤8得到的检测结果
Figure BDA0003757223260000052
和导频信号P，从接收信号Yr中减去其他发射天线数据和导频的影响，更新用于检测发射天线t到接收天线r之间信道系数的接收信号
Figure BDA0003757223260000053
表示如下。
Figure BDA0003757223260000054
其中Yr表示接收天线r的接收信号。
步骤10：更新信道估计值。当前迭代次数i＝i+1，使用步骤8得到的粗检测结果
Figure BDA0003757223260000055
导频Pt和步骤9更新的频域检测信号
Figure BDA0003757223260000056
进行信道估计，得到接收天线r与发射天线t之间的信道估计更新值
Figure BDA0003757223260000057
更新所有接收天线与发射天线间的信道信息，得到
Figure BDA0003757223260000058
当i＝I，迭代结束，最终信道估计值
Figure BDA0003757223260000059
转到步骤11，否则，转到步骤8。
步骤11：数据最终检测。使用最终的信道估计值
Figure BDA00037572232600000510
和导频P，频域接收信号Y进行信号检测，得到发送的数据比特流。
本发明的有益效果为：本发明根据一定的导频数据能量比在频域将导频稀疏叠加在调制信号上，与导频单独占用子载波的放置方式相比，减少了导频开销，提高了系统的频谱效率，且通过迭代更新也改善了系统的BER性能。
附图说明
图1为实施例1的叠加稀疏梳状导频原理图示意。
图2为实施例2的叠加稀疏菱形导频原理图示意。
图3为应用叠加稀疏导频方法的OFDM-SIM系统的发射端框图。
图4为应用叠加稀疏导频方法的OFDM-SIM系统的接收端框图。
具体实施方式
下面结合附图和实施例，详细描述本发明的技术方案：
实施例1：
下面以发送天线NT为2，接收天线NR为1，循环前缀CP＝64，调制符号采用BPSK，每个子块子载波L＝2，激活子载波K＝1，总载波数N＝1024的(2,1)OFDM-SIM系统，一帧包含Ns＝4个OFDM符号，设插入导频间隔d＝4，则一个OFDM符号含有导频数量V＝256，设发射天线1使用的导频序列为
Figure BDA0003757223260000061
发射天线2使用的导频序列为
Figure BDA0003757223260000062
假设两根发射天线的导频位置相同，均为A，设置叠加梳状导频，每个OFDM符号的导频位置均为Aj＝1:4:1024(j＝1,2,3,4)，则
Figure BDA0003757223260000063
信道估计采用LS信道估计，信号检测采用ML检测方法，迭代次数I＝5，为例介绍本发明的具体实施方式。
系统子块数为G＝N/L＝1024/2＝512，一个子块内索引比特长度
Figure BDA0003757223260000064
Figure BDA0003757223260000065
表示向下取整，BPSK调制符号的调制阶数M＝2，调制比特长度b2＝K log2(M)＝1×log2(2)＝1。
用t(t＝1,2)表示发射天线编号，r(r＝1)表示接收天线编号，j(j＝1,2,3,4)表示OFDM符号编号，g(t＝1,2,…,512)表示子块编号，l(l＝1,2)表示子块内子载波编号。
步骤1：产生信息比特流。以第t根发射天线的第j个OFDM符号的子块g为例，产生信息比特
Figure BDA0003757223260000066
步骤2：星座符号映射和子载波激活选择。对于发射天线t第j个OFDM符号的子块g，它的发送比特
Figure BDA0003757223260000067
将调制比特
Figure BDA0003757223260000068
进行BPSK星座映射得到调制符号
Figure BDA0003757223260000069
子载波激活规则为：索引比特
Figure BDA00037572232600000610
为‘0’，表示激活该子块的第一个子载波，索引比特
Figure BDA00037572232600000611
为‘1’，表示激活该子块的第二个子载波，则通过索引比特
Figure BDA00037572232600000612
获得激活子载波位置
Figure BDA00037572232600000613
通过索引映射将调制符号放置到激活子载波上得到频域信号
Figure BDA00037572232600000614
其中
Figure BDA00037572232600000615
步骤3：生成导频图案。利用导频位置A，生成对应发射天线t第j个OFDM符号的子块g的导频信号
Figure BDA00037572232600000616
令l′＝l+(g-1)*2，其中l＝1,2，映射规则为：
Figure BDA0003757223260000071
则发射天线t的第j个OFDM符号的导频信号为
Figure BDA0003757223260000072
发射天线t的导频信号为
Figure BDA0003757223260000073
则生成的导频信号为P＝[P1；P2]。
步骤4：导频与数据叠加。对应现有参数设置的系统下，对应发射天线t，一个OFDM符号的数据和导频功率总和为1024，导频功率Ep与数据功率Ed满足512×1×Ed+256×Ep＝1024，设置Ep＝0.6，Ed＝1.7。将导频Pt(t＝1,2)与经过BPSK调制和索引调制的数据叠加放置，对应发射天线t的第j个OFDM符号子块g的频域符号表示为：
Figure BDA0003757223260000074
其中
Figure BDA0003757223260000075
则发射天线t的第j个OFDM-SIM频域信号为
Figure BDA0003757223260000076
发射天线t的频域符号表示为：
Figure BDA0003757223260000077
则最终发送的符号为S＝[S1；S2]。
步骤5：频域-时域变换。将步骤4得到的发送符号向量S进行IFFT变换，得到时域发送信号s＝[s1；s2]，其中
Figure BDA0003757223260000078
而
Figure BDA0003757223260000079
IFFT变换公式为：
Figure BDA00037572232600000710
Figure BDA00037572232600000711
是
Figure BDA00037572232600000712
的第n(n＝1,2,…,1024)个元素。再在时域信号
Figure BDA00037572232600000713
上添加循环前缀(CP)得到最终时域发送符号
Figure BDA00037572232600000714
通过发射天线发送。
接收端：
步骤6：时域-频域转换。将时域接收信号r＝[r1]，其中r＝[r1,r2,r3,r4]，进行去循环前缀(CP)操作，yj＝rj[64+(1:1024)](j＝1,2,3,4)，再进行FFT操作得到频域待检测信号Y＝[Y1]＝[Y1,Y2,Y3,Y4]，用矩阵的形式可以表示为：
Figure BDA00037572232600000715
其中，Yj＝[Y1,j,Y2,j,…,Y1024,j]T，FFT计算公式为：
Figure BDA00037572232600000716
步骤7：初始信道估计。根据插入导频位置A取出对应的检测信号Yp＝Y[A]和有效导频
Figure BDA00037572232600000717
使用频域检测信号Yp，导频
Figure BDA00037572232600000718
和导频位置A，采用LS信道估计方法得到导频位置的第一次信道估计值
Figure BDA0003757223260000081
其中
Figure BDA0003757223260000082
再通过频域线性插值得到整个频带的初始信道估计值
Figure BDA0003757223260000083
迭代次数I＝5，设置i＝0，表示当前迭代次数，令
Figure BDA0003757223260000084
转到步骤8。
步骤8：信号粗检测。根据信道估计值
Figure BDA0003757223260000085
和频域检测信号Y，导频P，以子块为单位，使用ML检测算法：
Figure BDA0003757223260000086
Ω表示一个子块的所有可能发送符号集合。检测完所有子块，最终得到信号粗检测结果
Figure BDA0003757223260000087
其中
Figure BDA0003757223260000088
步骤9：干扰抵消操作。为了检测第t根发射天线的数据，用步骤8得到的检测结果
Figure BDA0003757223260000089
和导频信号P，从检测信号Y中减去其他天线发射数据的影响。
Figure BDA00037572232600000810
表示第t(t＝1,2)根发射天线到接收天线之间消除其他发射天线数据干扰的更新接收信号。
Figure BDA00037572232600000811
步骤10：更新信道估计值。当前迭代次数i＝i+1，使用步骤8得到的粗检测结果
Figure BDA00037572232600000812
导频Pt和步骤9更新的频域检测信号
Figure BDA00037572232600000813
使用LS信道估计算法更新所有接收天线与发射天线间的信道信息，得到信道估计更新值
Figure BDA00037572232600000814
其中
Figure BDA00037572232600000815
当i＝5，迭代结束，最终信道估计值
Figure BDA00037572232600000816
转到步骤11，否则，转到步骤8。
步骤11：数据最终检测。使用最终的信道估计值
Figure BDA00037572232600000820
导频P和频域检测信号Y以子块为单位进行ML检测：
Figure BDA00037572232600000817
Ω表示一个子块的所有可能发送符号集合。检测完所有子块，最终得到信号检测结果
Figure BDA00037572232600000818
其中
Figure BDA00037572232600000819
恢复出发送的数据比特流。
实施例2：
下面以发送天线NT为2，接收天线NR为1，循环前缀CP＝64，调制符号采用BPSK，每个子块子载波L＝2，激活子载波K＝1，总载波数N＝1024的(2,1)OFDM-SIM系统，一帧包含Ns＝4个OFDM符号，设插入导频间隔d＝4，则一个OFDM符号含有导频数量V＝256，设发射天线1使用的导频为
Figure BDA0003757223260000091
发射天线2使用的导频为
Figure BDA0003757223260000092
假设两根发射天线的导频位置相同，均为A，设置导频形状为菱形导频，相邻OFDM符号的导频交错放置，导频位置为Aj＝1:4:1024(j＝1,3)和Aj＝3:4:1024(j＝2,4)，则
Figure BDA0003757223260000093
信道估计采用MMSE和LS信道估计，信号检测采用ML检测方法，迭代次数I＝5，为例介绍本发明的具体实施方式。
系统子块数为G＝N/L＝1024/2＝512，一个子块内索引比特长度
Figure BDA0003757223260000094
Figure BDA0003757223260000095
表示向下取整，BPSK调制符号的调制阶数M＝2，调制比特长度b2＝K log2(M)＝1×log2(2)＝1。
用t(t＝1,2)表示发射天线编号，r(r＝1)表示接收天线编号，j(j＝1,2,3,4)表示OFDM符号编号，g(t＝1,2,…,512)表示子块编号，l(l＝1,2)表示子块内子载波编号。
步骤1：产生信息比特流。以发射天线t的第j个OFDM符号的子块g为例，产生信息比特
Figure BDA0003757223260000096
步骤2：星座符号映射和子载波激活选择。对于发射天线t第j个OFDM符号的子块g，它的发送比特
Figure BDA0003757223260000097
将调制比特
Figure BDA0003757223260000098
进行BPSK星座映射得到调制符号
Figure BDA0003757223260000099
子载波激活规则为：索引比特
Figure BDA00037572232600000910
为‘0’，表示激活该子块的第一个子载波，索引比特
Figure BDA00037572232600000911
为‘1’，表示激活该子块的第二个子载波，则通过索引比特
Figure BDA00037572232600000912
获得激活子载波位置
Figure BDA00037572232600000913
通过索引映射将调制符号放置到激活子载波上得到频域信号
Figure BDA00037572232600000914
其中
Figure BDA00037572232600000915
步骤3：生成导频图案。利用导频位置A，生成对应发射天线t第j个OFDM符号的子块g的导频信号
Figure BDA00037572232600000916
令l′＝l+(g-1)*2，其中l＝1,2，映射规则为：
Figure BDA0003757223260000101
则发射天线t的第j个OFDM符号的导频信号为
Figure BDA0003757223260000102
发射天线t的导频信号为
Figure BDA0003757223260000103
则生成的导频信号为P＝[P1；P2]。
步骤4：导频与数据叠加。对应现有参数设置的系统下，对应发射天线t，一个OFDM符号的数据和导频功率总和为1024，导频功率Ep与数据功率Ed满足512×1×Ed+256×Ep＝1024，设置Ep＝0.6，Ed＝1.7。将导频Pt(t＝1,2)与经过BPSK调制和索引调制的数据叠加放置，对应发射天线t的第j个OFDM符号子块g的频域符号表示为：
Figure BDA0003757223260000104
其中
Figure BDA0003757223260000105
则发射天线t的第j个OFDM-SIM频域信号为
Figure BDA0003757223260000106
发射天线t的频域符号表示为：
Figure BDA0003757223260000107
则发送符号为S＝[S1；S2]。
步骤5：频域-时域变换。将步骤4得到的发送符号向量S进行IFFT变换，得到时域发送信号s＝[s1；s2]，其中
Figure BDA0003757223260000108
而
Figure BDA0003757223260000109
IFFT变换公式为：
Figure BDA00037572232600001010
Figure BDA00037572232600001011
是
Figure BDA00037572232600001012
的第n(n＝1,2,…,1024)个元素。再在时域信号
Figure BDA00037572232600001013
上添加循环前缀(CP)得到最终时域发送符号
Figure BDA00037572232600001014
通过发射天线发送。
接收端：
步骤6：时域-频域转换。将时域接收信号r＝[r1]，其中r＝[r1,r2,r3,r4]，进行去循环前缀(CP)操作，yj＝rj[64+(1:1024)](j＝1,2,3,4)，再进行FFT操作得到频域待检测信号Y＝[Y1]＝[Y1,Y2,Y3,Y4]，用矩阵的形式可以表示为：
Figure BDA00037572232600001015
其中，
Figure BDA00037572232600001016
FFT计算公式为：
Figure BDA00037572232600001017
步骤7：初始信道估计。根据插入导频位置A取出对应的检测信号Yp＝Y[A]和有效导频
Figure BDA00037572232600001018
使用频域检测信号Yp，导频
Figure BDA00037572232600001019
和导频位置A，先采用LS信道估计方法得到导频位置的初始信道估计值
Figure BDA0003757223260000111
其中
Figure BDA0003757223260000112
再采用MMSE信道估计得到整个频带的信道估计值
Figure BDA0003757223260000113
MMSE信道估计方法为：
Figure BDA0003757223260000114
Figure BDA0003757223260000115
是频域上真实信道向量和临时信道估计向量之间的互相关矩阵，RHH是信道自相关矩阵，E是维度与RHH相同的单位阵，
Figure BDA0003757223260000116
表示信噪比。其中
Figure BDA0003757223260000117
和RHH中的元素满足：
Figure BDA0003757223260000118
其中，n(n＝1,2,…,1024)和j(j＝1,2,3,4)分别表示子载波(频率)和OFDM符号(时间)的编号。在一个呈指数衰减的多径功率时延谱(PDP)中，频域相关
Figure BDA0003757223260000119
△f是子载波间隔，τrms是信道均方时延。同时，对于具有最大多普勒频率fmax和Jakes功率谱的衰落信道，时域相关rt[j]为：rt[j]＝J0(2πfmaxjTsym)，Tsym为OFDM符号周期，J0(·)为第一类零阶贝塞尔函数。
迭代次数I＝5，设置i＝0，表示当前迭代次数，令
Figure BDA00037572232600001110
转到步骤8。
步骤8：信号粗检测。根据信道估计值
Figure BDA00037572232600001111
和频域检测信号Y，导频P，以子块为单位，使用ML检测算法：
Figure BDA00037572232600001112
Ω表示一个子块的所有可能发送符号集合。检测完所有子块，最终得到信号粗检测结果
Figure BDA00037572232600001113
其中
Figure BDA00037572232600001114
步骤9：干扰抵消操作。为了检测第t根发射天线的数据，用步骤8得到的检测结果
Figure BDA00037572232600001115
和导频信号P，从检测信号Y中减去其他天线发射数据的影响。
Figure BDA00037572232600001116
表示第t(t＝1,2)根发射天线到接收天线之间消除其他发射天线数据干扰的更新接收信号。
Figure BDA00037572232600001117
步骤10：更新信道估计值。当前迭代次数i＝i+1，使用步骤8得到的粗检测结果
Figure BDA00037572232600001118
导频Pt和步骤9更新的频域检测信号
Figure BDA00037572232600001119
使用LS信道估计算法更新所有接收天线与发射天线间的信道信息，得到信道估计更新值
Figure BDA00037572232600001120
其中
Figure BDA00037572232600001121
当i＝5，迭代结束，最终信道估计值
Figure BDA0003757223260000121
转到步骤11，否则，转到步骤8。
步骤11：数据最终检测。使用最终的信道估计值
Figure BDA0003757223260000122
导频P和频域检测信号Y以子块为单位进行ML检测：
Figure BDA0003757223260000123
Ω表示一个子块的所有可能发送符号集合。检测完所有子块，最终得到信号检测结果
Figure BDA0003757223260000124
其中
Figure BDA0003757223260000125
恢复出发送的数据比特流。

Claims (1)

1.用于多天线子载波索引调制OFDM的导频传输方法，定义多天线OFDM-SIM系统发射天线数量为NT，接收天线数量为NR，系统总子载波数为N，每个子块载波数为L，其中K个激活载波用于传输M阶星座调制符号，其余(L-K)个子载波静默，因此含有G＝N/L个子块；设置V为有效导频长度，Ed为数据的归一化发送功率，Ep为导频的归一化发送功率，且满足G·K·Ed+V·Ep＝N，接收机信道估计迭代次数设为I；其特征在于，所述方法包括以下步骤：
发射端：
步骤1：产生信息比特流：对于任意一个子块g，信息比特由索引比特和调制比特组成，其中
Figure FDA0003757223250000011
个索引比特用于指示激活子载波位置，b2＝Klog2(M)个调制比特进行星座符号映射，因此，一个OFDM-SIM符号传输有m＝G(b1+b2)位比特；
步骤2：星座符号映射和子载波激活选择：将调制比特进行M阶星座符号映射得到调制符号向量Z，利用索引比特选择激活的子载波，得到激活子载波位置向量Γ，根据Γ将星座符号放置到激活子载波上，得到频域信号
Figure FDA0003757223250000012
步骤3：生成导频信号：通过发射天线t的导频位置At将长度为V的导频序列pt放置到对应位置，t＝1,2,…,NT，得到导频信号
Figure FDA0003757223250000013
步骤4：在系统发送端，通过导频与数据叠加得到发射天线t的系统发送符号：
Figure FDA0003757223250000014
从而得到发送端频域符号
Figure FDA0003757223250000015
步骤5：频域-时域变换：将步骤4得到的发送符号向量S依次进行IFFT、和加循环前缀得到OFDM-SIM时域符号，通过发射天线发送；
接收端：
步骤6：时域-频域转换：将接收信号进行去循环前缀，FFT操作后获得频域待检测信号Y；
步骤7：初始信道估计：按照已知的有效导频位置A，取出对应的检测信号Yp和有效导频Pp，使用频域检测信号Yp和导频Pp，依次进行信道估计和频域插值得到初始信道估计值
Figure FDA0003757223250000021
若迭代次数I＝0，最终信道估计值
Figure FDA0003757223250000022
转到步骤11，否则，设置i＝0，表示当前迭代次数，令
Figure FDA0003757223250000023
转到步骤8；
步骤8：信号粗检测：根据信道估计值
Figure FDA0003757223250000024
频域接收信号Y和导频P进行信号粗检测得到检测结果
Figure FDA0003757223250000025
步骤9：干扰抵消操作：根据步骤8得到的检测结果
Figure FDA0003757223250000026
和导频信号P，从检测信号Yr中减去其他天线数据和导频的影响，更新用于检测发射天线t到接收天线r之间信道系数的接收信号
Figure FDA0003757223250000027
表示为：
Figure FDA0003757223250000028
步骤10：更新信道估计值：当前迭代次数i＝i+1，使用步骤8得到的粗检测结果
Figure FDA0003757223250000029
导频Pt和步骤9更新的频域检测信号
Figure FDA00037572232500000210
进行信道估计，得到接收天线r与发射天线t之间的信道估计更新值
Figure FDA00037572232500000211
更新所有接收天线与发射天线间的信道信息，得到
Figure FDA00037572232500000212
当i＝I，迭代结束，最终信道估计值
Figure FDA00037572232500000213
转到步骤11，否则，转到步骤8；
步骤11：使用最终的信道估计值
Figure FDA00037572232500000214
和导频P，频域检测信号Y进行检测，得到发送的数据比特流。
CN202210859170.XA 2022-07-21 2022-07-21 用于多天线子载波索引调制ofdm的导频传输方法 Pending CN115276938A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210859170.XA CN115276938A (zh) 2022-07-21 2022-07-21 用于多天线子载波索引调制ofdm的导频传输方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210859170.XA CN115276938A (zh) 2022-07-21 2022-07-21 用于多天线子载波索引调制ofdm的导频传输方法

Publications (1)

Publication Number Publication Date
CN115276938A true CN115276938A (zh) 2022-11-01

Family

ID=83767809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210859170.XA Pending CN115276938A (zh) 2022-07-21 2022-07-21 用于多天线子载波索引调制ofdm的导频传输方法

Country Status (1)

Country Link
CN (1) CN115276938A (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW201210226A (en) * 2010-02-23 2012-03-01 Qualcomm Inc Multi-term demapping for multi-channel wireless communication
CN106656875A (zh) * 2017-02-21 2017-05-10 电子科技大学 用于子载波索引调制ofdm系统的信道估计方法
CN107147483A (zh) * 2017-06-21 2017-09-08 电子科技大学 一种基于导频图案调制的sim‑ofdm系统通信方法
CN108847917A (zh) * 2018-06-28 2018-11-20 电子科技大学 一种导频图案调制正交频分复用传输方法
CN109412998A (zh) * 2018-12-25 2019-03-01 电子科技大学 导频图案调制系统中的位置图案设计方法
CN111163026A (zh) * 2019-12-26 2020-05-15 电子科技大学 用于空间调制多天线系统的导频放置与信道估计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW201210226A (en) * 2010-02-23 2012-03-01 Qualcomm Inc Multi-term demapping for multi-channel wireless communication
CN106656875A (zh) * 2017-02-21 2017-05-10 电子科技大学 用于子载波索引调制ofdm系统的信道估计方法
CN107147483A (zh) * 2017-06-21 2017-09-08 电子科技大学 一种基于导频图案调制的sim‑ofdm系统通信方法
CN108847917A (zh) * 2018-06-28 2018-11-20 电子科技大学 一种导频图案调制正交频分复用传输方法
CN109412998A (zh) * 2018-12-25 2019-03-01 电子科技大学 导频图案调制系统中的位置图案设计方法
CN111163026A (zh) * 2019-12-26 2020-05-15 电子科技大学 用于空间调制多天线系统的导频放置与信道估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵矗: "高效导频传输与检测技术研究与仿真", 《中国优秀硕士学位论文全文数据库信息科技辑》, no. 1 *

Similar Documents

Publication Publication Date Title
CN101951353B (zh) 一种干扰环境下的正交频分复用系统信道估计方法
CN101005475A (zh) 正交频分复用通信中时间和频率同步的方法及系统
CN101222470B (zh) 双天线离散傅立叶扩频广义多载波系统的信道估计方法
CN101815042B (zh) 一种正交频分复用系统信道估计方法和装置
JP4650104B2 (ja) 通信装置及び通信方法
CN101204030A (zh) 一种正交频分时分发射机、接收机及其方法
WO2009107146A1 (en) Optimal training sequence and channel estimation method and system for superimposed training based ofdm systems
US8391427B2 (en) Channel estimation methods and apparatus utilizing the same
CN109600327B (zh) 一种基于虚部干扰利用的信道估计方法
JP2005304040A (ja) 直交周波数分割多重接続システムにおける高速周波数ホッピングのための送受信装置
CN111585688B (zh) 一种基于索引调制的ocdm水声通信方法
TWI410090B (zh) 正交分頻多工系統之訊號發送、接收方法及裝置
CN101083515A (zh) 一种发射分集的正交频分复用的信道估计方法及装置
CN102752253A (zh) 时频域联合处理抑制正交频分复用系统子载波间干扰的方法
CN102143101A (zh) 镜像扩展的频域加窗正交频分多址信道估计方法
CN108847917B (zh) 一种导频图案调制的正交频分复用传输方法
CN107634922A (zh) 基于小波变换调制系统的压缩感知时域信道估计方法
CN108234102A (zh) 一种具有低复杂度检测算法的sm-gfdm系统
CN101807954A (zh) 上行多用户时域同步频分多址接入方法
JP4830613B2 (ja) マルチユーザ通信システム、通信装置及びそれらを用いるマルチパス伝送路推定方法
US8098749B2 (en) CFR estimation method for multi-band OFDM-based UWB systems
KR20090074591A (ko) 사이클릭 프리픽스가 없는 신호를 처리하는 중계기 기반의통신 시스템 및 그 방법
CN102487364B (zh) 一种信道估计方法及装置
WO2011009399A1 (zh) 导频插入方法、导频插入模块和发射机
CN101888359B (zh) 基于tds-ofdm的多址接入系统及其信号帧结构、接收方法

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20221101