CN107947899B - 单载波交织频分多址系统中的多用户信噪比盲估计方法 - Google Patents

单载波交织频分多址系统中的多用户信噪比盲估计方法 Download PDF

Info

Publication number
CN107947899B
CN107947899B CN201711144165.6A CN201711144165A CN107947899B CN 107947899 B CN107947899 B CN 107947899B CN 201711144165 A CN201711144165 A CN 201711144165A CN 107947899 B CN107947899 B CN 107947899B
Authority
CN
China
Prior art keywords
user
signal
multiple access
division multiple
frequency division
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
CN201711144165.6A
Other languages
English (en)
Other versions
CN107947899A (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN201711144165.6A priority Critical patent/CN107947899B/zh
Publication of CN107947899A publication Critical patent/CN107947899A/zh
Application granted granted Critical
Publication of CN107947899B publication Critical patent/CN107947899B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/2668—Details of algorithms
    • H04L27/2673—Details of algorithms characterised by synchronisation parameters
    • H04L27/2676—Blind, i.e. without using known symbols
    • 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/0001—Systems modifying transmission characteristics according to link quality, e.g. power backoff
    • H04L1/0006—Systems modifying transmission characteristics according to link quality, e.g. power backoff by adapting the transmission format
    • 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/0056—Systems characterized by the type of code used
    • H04L1/0071—Use of interleaving
    • 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
    • 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/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2692—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with preamble design, i.e. with negotiation of the synchronisation sequence with transmitter or sequence linked to the algorithm used at the receiver

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Power Engineering (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开一种单载波交织频分多址系统中的多用户信噪比盲估计方法，其步骤为：(1)将输入的信号变换为频域信号；(2)将频域数据变换为时域数据；(3)发送已调制信号；(4)接收含有噪声的信号；(5)多用户信号分离获得各用户时域信号；(6)计算噪声的平均功率；(7)计算单载波交织频分多址符号的平均功率；(8)计算单载波交织频分多址符号的信噪比；(9)计算每个子帧的信噪比。本发明盲估计信噪比，提升了单载波频分多址系统的传输效率；通过大量的各用户时域信号直接计算每个子帧的信噪比，使得本发明能够提高信噪比估计的准确性，尤其是低信噪比估计的准确性。

Description

单载波交织频分多址系统中的多用户信噪比盲估计方法
技术领域
本发明属于通信技术领域，更进一步涉及信号处理技术领域中的单载波交织频分多址系统中的多用户信噪比盲估计方法。本发明可用于卫星通信、移动通信、电离层散射通信、对流层散射通信等通信技术领域中，采用单载波交织频分多址系统(SC-IFDMA)，并对系统中各用户的信噪比进行盲估计的方法。
背景技术
在通信系统中，常见的表征信道传输质量的指标有：信道时延扩展、信噪比、误码率、误帧率等。其中误码率、误帧率是通信系统传输质量的直接衡量指标，往往需要统计大量的数量才能得到，耗时较长，不利于即时使用。但信噪比(Signal to Noise Ratio，SNR)可以实时测量，它与误码率、误帧率也有直接的关系，所以常采用信噪比作为通信系统传输质量的衡量指标。高精度的信噪比估计可以为功率控制、越区切换和信道速率匹配等提供准确的信道质量信息，充分发掘信道潜在的传输能力，提升通信系统的性能。
现有的信噪比估计方法总体上可以分为两种：一种是基于已知数值的导频数据估计信噪比的方法，另一种是盲估计方法。前者需要不断地在每条子载波上发送大量的导频数据，提高估计的准确性，这样会降低通信系统的传输效率；后者是通过统计有用信号和噪声干扰的特征值估计信噪比，插入的导频数据也作为部分样本数据盲估计信噪比。
中国科学院上海微系统与信息技术研究所在其申请的专利文献“一种信噪比估计方法”(申请号201410152838.2，申请日2014.04.16，公开号103916342A)中提出一种根据接收信号中插入的导频数据估计信道信噪比的方法。该方法首先从接收到的数据符号中获取到导频所在位置的导频数值、信道估计值和高斯白噪声，然后根据导频数值、信道估计值和高斯白噪声获取相应导频所在位置的导频位置信息，最后根据导频位置、信道估计值和导频位置信息获取接收信号的信噪比估计值。该方法存在的不足之处是，在估计信噪比时，需要一定的先验信息，为提高信噪比估计的准确性，需要在频域信号中插入较多的导频数据，这样会降低数据传输的效率，但过少的导频会导致信噪比估计的准确度下降，不能很好的反应通信链路质量，影响系统的性能。
上海鑫皇实业有限公司在其申请的专利文件“一种OFDM系统的信噪比盲估计方法”(申请号201510288359.8，申请日2015.05.29，公开号106302286A)中提出一种OFDM系统的信噪比盲估计方法。该方法首先对OFDM前导序列进行傅里叶变换，然后经傅里叶变换后的前导符号获取数据子载波平均功率，经傅里叶变换后的前导符号获取虚载波平均功率，最后根据信号的平均功率和虚载波平均功率确定信噪比。该方法存在的不足之处是，需要使用前导符号估计信号的功率和噪声功率，若信道变化剧烈，前导符号估计出的信噪比不再适用于数据符号部分，而且前导符号的样本数据有限，导致低信噪比估计不准确。
发明内容
本发明的目的在于针对上述现有技术存在的不足，提出一种单载波交织频分多址系统中使用各用户时域信号盲估计信噪比的方法，适用于信道具有非线性功放，信号包络波动小，低信噪比特性的通信环境，可以提升估计低信噪比的准确性，提高单载波交织频分多址系统的传输效率。
为实现上述目的，本发明采取的技术方案，包括如下步骤：
(1)将输入的信号变换为频域信号：
将各用户的的二进制信号，输入到单载波交织频分多址系统中，分别经交织和编码、数字调制、串并转换后，进行点数等于M的离散傅里叶变换DFT，得到各用户的频域信号，其中M表示单载波交织频分多址系统中各用户占用的子载波数目；
(2)将频域信号变换为时域信号：
在各用户的频域信号中，插入在接收端各用户时域信号中作为部分样本数据盲估计信噪比的块状导频，按照交织式映射公式，分别将各用户的M条子载波映射到总数为N的子载波上，进行点数与N相等的逆快速傅里叶变换IFFT，得到多用户时域信号，其中N表示单载波交织频分多址系统中多用户子载波的总数；
(3)发送已调制信号：
在多用户时域信号中加入长度为C的循环前缀，将经过并串转换后生成的多用户已调制信号以子帧为单位发送出去，其中C表示单载波交织频分多址系统中循环前缀占用的子载波数目；
(4)接收含有噪声的信号：
多用户已调制信号经过平坦信道或者频率选择性衰落信道，受到信道高斯白噪声的干扰，到达接收端，接收端接收到含有噪声的信号；
(5)多用户信号分离获得各用户时域信号：
(5a)将含有噪声的信号进行串并转换，去除长度为C的循环前缀，得到含有噪声的多用户时域信号；
(5b)将含有噪声的多用户时域信号进行点数等于N的快速傅里叶变换FFT，得到含有噪声的多用户频域信号；
(5c)将含有噪声的多用户频域信号，按照交织式解映射公式，将N条子载波分别解映射到M条子载波上，得到含有噪声的各用户频域信号；
(5d)将含有噪声的各用户频域信号进行点数等于M的逆快速傅里叶变换IFFT，再按照交织式映射公式，分别将各用户的M条子载波映射到总数为N的子载波上，进行点数与N相等的逆快速傅里叶变换IFFT，得到各用户时域信号；
(6)计算噪声的平均功率：
利用噪声平均功率公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号噪声的平均功率；
(7)计算单载波交织频分多址符号的平均功率：
利用符号平均功率公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号的平均功率；
(8)计算单载波交织频分多址符号的信噪比：
利用符号信噪比公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号的信噪比；
(9)计算每个子帧的信噪比：
利用子帧信噪比公式，计算各用户时域信号每个子帧的信噪比，所估计的信噪比指的是多用户已调制信号中各个用户受信道高斯白噪声干扰的信噪比。
本发明与现有技术相比具有以下优点：
第一，由于本发明在单载波交织频分多址系统中，在各用户的频域信号中，插入在接收端各用户时域信号中作为部分样本数据盲估计信噪比的块状导频，克服了现有技术中基于已知数值的导频数据估计信噪比的方法，需在频域信号中插入较多导频数据，并在接收端用这些已知数值的导频数据估计信道的信噪比，会降低单载波交织频分多址系统传输效率的缺点，使得本发明能够直接使用各用户时域信号对信道的信噪比进行盲估计，提升了单载波交织频分多址系统的传输效率。
第二，由于本发明单载波交织频分多址系统中，直接计算各用户时域信号每个子帧中每个单载波交织频分多址符号噪声的平均功率和单载波交织频分多址符号的平均功率，进而估计单载波交织频分多址系统中每个子帧的信噪比，克服了现有技术中样本数据有限，不能准确估计信噪比的缺点，使得本发明能够利用大量的各用户时域信号估计信噪比，提高信噪比估计的准确性，尤其是低信噪比估计的准确性。
附图说明
图1是本发明的流程框图；
图2是本发明单载波交织频分多址系统子帧结构示意图；
图3是本发明各用户时域信号估计信噪比与现有技术的基于已知数值的导频数据估计信噪比仿真对比图；
图4是本发明各用户时域信号信噪比估计结果与理论信噪比归一化均方误差及现有技术的基于已知数值的导频数据估计信噪比结果与理论信噪比归一化均方误差对比图。
具体实施方式
下面结合附图对本发明做进一步的描述。
参照附图1，对本发明实现的具体步骤描述如下。
步骤1，将输入的信号变换为频域信号。
将各用户的的二进制信号，输入到单载波交织频分多址系统中，分别经交织和编码、数字调制、串并转换后，进行点数等于M的离散傅里叶变换DFT，得到各用户的频域信号，其中M表示单载波交织频分多址系统中各用户占用的子载波数目。
步骤2，将频域信号变换为时域信号。
参照附图2所示的单载波交织频分多址系统子帧结构示意图，白色实线方框表示单载波交织频分多址系统中各用户的数据，黑色实线方框表示在各用户的频域信号中，插入在接收端各用户时域信号中作为部分样本数据盲估计信噪比的块状导频，白色虚线勾画出的方框表示单载波交织频分多址系统中单载波交织频分多址符号。
在各用户的频域信号中，插入在接收端各用户时域信号中作为部分样本数据盲估计信噪比的块状导频，按照交织式映射公式，分别将各用户的M条子载波映射到总数为N的子载波上，进行点数与N相等的逆快速傅里叶变换IFFT，得到多用户时域信号，其中N表示单载波交织频分多址系统中多用户子载波的总数。
交织式映射公式如下：
Figure BDA0001472116260000051
其中，Xk表示交织式映射后总数为N的第k条子载波，k表示多用户子载波的序号，k的取值是k＝0,1,2,...,N-1，Xp表示交织式映射前总数为M第p条子载波，p表示各用户子载波的序号，p的取值是p＝0,1,2,...,M-1，Q表示子载波映射间隔，其数值大小等于单载波交织频分多址系统中可容纳的用户数，
Figure BDA0001472116260000052
u表示不同用户映射的起始位置，其取值范围是0≤u≤Q-1。
多用户时域信号是指，子载波的总数为N，单载波交织频分多址符号为F*L二维时频结构的多用户时域信号，其中，F表示单载波交织频分多址系统中子帧的总数，*表示相乘操作，L表示每个子帧中单载波交织频分多址符号的总数。
步骤3，发送已调制信号。
在多用户时域信号中加入长度为C的循环前缀，将经过并串转换后生成的多用户已调制信号以子帧为单位发送出去，其中C表示单载波交织频分多址系统中循环前缀占用的子载波数目。
步骤4，接收含有噪声的信号。
多用户已调制信号经过平坦信道或者频率选择性衰落信道，受到信道高斯白噪声的干扰，到达接收端，接收端接收到含有噪声的信号。
步骤5，多用户信号分离获得各用户时域信号。
将含有噪声的信号进行串并转换，去除长度为C的循环前缀，得到含有噪声的多用户时域信号；
将含有噪声的多用户时域信号进行点数等于N的快速傅里叶变换FFT，得到含有噪声的多用户频域信号；
将含有噪声的多用户频域信号，按照交织式解映射公式，将N条子载波分别解映射到M条子载波上，得到含有噪声的各用户频域信号；
将含有噪声的各用户频域信号进行点数等于M的逆快速傅里叶变换IFFT，再按照交织式映射公式，分别将各用户的M条子载波映射到总数为N的子载波上，进行点数与N相等的逆快速傅里叶变换IFFT，得到各用户时域信号。
交织式解映射公式如下：
Figure BDA0001472116260000061
其中，
Figure BDA0001472116260000062
表示交织式解映射后总数为M第p条子载波，p表示各用户子载波的序号，p的取值是p＝0,1,2,...,M-1，
Figure BDA0001472116260000063
表示交织式解映射前总数为N的第k条子载波，k表示多用户子载波的序号，k的取值是k＝0,1,2,...,N-1，u表示不同用户映射的起始位置，其取值范围是0≤u≤Q-1，Q表示子载波映射间隔，其数值大小等于单载波交织频分多址系统中可容纳的用户数，
Figure BDA0001472116260000064
步骤6，计算噪声的平均功率。
根据各用户时域信号按子载波序号具有周期性重复的性质，重复形式与映射的起始位置u有关，利用噪声平均功率公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号噪声的平均功率。
噪声平均功率公式如下：
Figure BDA0001472116260000065
其中，
Figure BDA0001472116260000071
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号噪声的平均功率，v表示单载波交织频分多址系统中各用户的序号，其取值是v＝0,1,2,...,Q-1，Q表示子载波映射间隔，其数值大小等于单载波交织频分多址系统中可容纳的用户数，
Figure BDA0001472116260000072
d表示单载波交织频分多址系统中子帧的序号，其取值是d＝0,1,2,...,F-1，F表示单载波交织频分多址系统中子帧的总数，t表示单载波交织频分多址符号的序号，其取值是t＝0,1,2,...,L-1，L表示每个子帧中单载波交织频分多址符号的总数，k表示多用户子载波的序号，k的取值是k＝0,1,2,...,N-1，∑表示求和操作，||·||2表示求二阶矩操作，[·]表示求实部和虚部的绝对值操作，Yt,k表示各用户时域信号每个子帧的第t个单载波交织频分多址符号中第k条子载波，
Figure BDA0001472116260000073
表示各用户时域信号每个子帧的第t个单载波交织频分多址符号中第
Figure BDA0001472116260000074
条子载波，b表示按含有噪声的时域信号重复次数的序号，其取值是b＝0,1,2，
Figure BDA0001472116260000075
表示各用户时域信号中每个子帧的第t个单载波交织频分多址符号中第
Figure BDA0001472116260000076
条子载波，
Figure BDA0001472116260000077
表示各用户时域信号中每个子帧的第t个单载波交织频分多址符号中第
Figure BDA0001472116260000078
条子载波。
步骤7，计算单载波交织频分多址符号的平均功率。
利用符号平均功率公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号的平均功率。
符号平均功率公式是指，利用以下两种公式中的任意一种，计算各用户时域信号每个子帧的每个单载波交织频分多址符号的平均功率：
Figure BDA0001472116260000079
Figure BDA0001472116260000081
其中，
Figure BDA0001472116260000082
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的平均功率，v表示单载波交织频分多址系统中各用户的序号，其取值是v＝0,1,2,...,Q-1，Q表示子载波映射间隔，其数值大小等于单载波交织频分多址系统中可容纳的用户数，
Figure BDA0001472116260000083
d表示单载波交织频分多址系统中子帧的序号，其取值是d＝0,1,2,...,F-1，F表示单载波交织频分多址系统中子帧的总数，t表示单载波交织频分多址符号的序号，其取值是t＝0,1,2,...,L-1，L表示每个子帧中单载波交织频分多址符号的总数，k表示多用户子载波的序号，k的取值是k＝0,1,2,...,N-1，∑表示求和操作，||·||2表示求二阶矩操作，Yt,k表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的第k条子载波，
Figure BDA0001472116260000084
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号噪声的平均功率，[·]表示求实部和虚部的绝对值操作，b表示按含有噪声的时域信号重复次数的序号，其取值是b＝0,1,2，
Figure BDA0001472116260000085
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的第
Figure BDA0001472116260000086
条子载波，
Figure BDA0001472116260000087
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的第
Figure BDA0001472116260000088
条子载波。
步骤8，计算单载波交织频分多址符号的信噪比。
利用符号信噪比公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号的信噪比。
符号信噪比公式如下：
Figure BDA0001472116260000089
其中，
Figure BDA0001472116260000091
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的信噪比，v表示单载波交织频分多址系统中各用户的序号，其取值是v＝0,1,2,...,Q-1，Q表示子载波映射间隔，其数值大小等于单载波交织频分多址系统中可容纳的用户数，
Figure BDA0001472116260000092
d表示单载波交织频分多址系统中子帧的序号，其取值是d＝0,1,2,...,F-1，F表示单载波交织频分多址系统中子帧的总数，t表示单载波交织频分多址符号的序号，其取值是t＝0,1,2,...,L-1，L表示每个子帧中单载波交织频分多址符号的总数，k表示多用户子载波的序号，k的取值是k＝0,1,2,...,N-1，
Figure BDA0001472116260000093
表示表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的平均功率，
Figure BDA0001472116260000094
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号噪声的平均功率。
步骤9，计算每个子帧的信噪比。
利用子帧信噪比公式，计算各用户时域信号每个子帧的信噪比，所估计的信噪比指的是多用户已调制信号中各个用户受信道高斯白噪声干扰的信噪比。
子帧信噪比公式如下：
Figure BDA0001472116260000095
其中，
Figure BDA0001472116260000096
表示第v个用户时域信号中第d个子帧的信噪比，v表示单载波交织频分多址系统中各用户的序号，其取值是v＝0,1,2,...,Q-1，Q表示子载波映射间隔，其数值大小等于单载波交织频分多址系统中可容纳的用户数，
Figure BDA0001472116260000097
d表示单载波交织频分多址系统中子帧的序号，其取值是d＝0,1,2,...,F-1，F表示单载波交织频分多址系统中子帧的总数，L表示每个子帧中单载波交织频分多址符号的总数，表示求和操作，
Figure BDA0001472116260000098
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的信噪比，t表示单载波交织频分多址符号的序号，其取值是t＝0,1,2,...,L-1。
下面结合仿真图对本发明的效果做进一步的描述。
1.仿真条件：
运用MATLABR2013a软件对本发明和现有技术的基于已知数值的导频数据估计信噪比的方法进行对比仿真，其结果如图3和图4所示。
2.仿真内容与结果分析：
图3是本发明各用户时域信号估计信噪比与现有技术的基于已知数值的导频数据估计信噪比仿真对比图，其横轴表示信噪比，单位是dB，纵轴表示估计出的信噪比，单位是dB。本发明的仿真实验中取单载波交织频分多址系统中单用户占用的子载波数目M＝32，单载波交织频分多址系统中子载波的总数N＝256，单载波交织频分多址系统中子帧的总数F＝24，每个子帧中单载波交织频分多址符号的总数L＝15，子载波映射间隔Q＝8，映射的起始位置u＝0，使用如下的本发明符号平均功率公式，计算单载波交织频分多址符号的平均功率，进而估计各用户时域信号每个子帧的信噪比。
Figure BDA0001472116260000101
与本发明估计信噪比的方法对应的现有技术是基于已知数值的导频数据估计每个子帧的信噪比，将上述两种方法估计信噪比的结果与理论信噪比进行比较。图3中以星型标示的实线表示理论信噪比曲线，以三角形标示的实线表示基于已知数值的导频估计出的信噪比曲线，以五角星标示的实线表示本发明估计出的信噪比曲线。
由三条曲线可以看出，基于已知数值的导频数据估计出的信噪比，可以估计到-6dB，采用本发明估计出的信噪比，可以估计到-18dB。据此，本发明通过使用大量的各用户时域信号，能够有效提高信噪比估计的范围，获得更好的估计性能，同时，使用本发明技术与现有技术相比可以节省单载波交织频分多址系统的开销，提高信号的传输效率。
图4是本发明各用户时域信号信噪比估计结果与理论信噪比归一化均方误差及现有技术的基于已知数值的导频数据估计信噪比结果与理论信噪比归一化均方误差对比图，其横轴表示信噪比，单位是dB，纵轴表示归一化均方误差。本发明的仿真实验仿真条件同上，分别计算估计出的信噪比与真实信噪比的相对误差，得到归一化均方误差。
图4中以三角形标示的实线表示已知数值的导频估计出的信噪比与理论信噪比的归一化均方误差，以五角星标示的实线表示本发明估计出的信噪比与理论信噪比的归一化均方误差。由两条曲线可以看出，基于已知数值的导频数据估计出的信噪比的归一化均方误差要高于本发明，这是由于使用现有技术可采用已知数值的导频数据数量有限，估计的误差较大。使用本发明在信噪比估计正确的前提下，归一化均方误差在10-1以下。进一步证明采用本发明的方法可以提升信噪比估计的准确性。

Claims (6)

1.一种单载波交织频分多址系统中的多用户信噪比盲估计方法，其特征在于，包括如下步骤：
(1)将用户信号变换为频域信号：
将各用户的的二进制信号，输入到单载波交织频分多址系统中，分别经交织和编码、数字调制、串并转换后，进行点数等于M的离散傅里叶变换DFT，得到各用户的频域信号，其中M表示单载波交织频分多址系统中各用户占用的子载波数目；
(2)将频域信号变换为时域信号：
在各用户的频域信号中，插入在接收端各用户时域信号中作为部分样本数据盲估计信噪比的块状导频，按照交织式映射公式，分别将各用户的M条子载波映射到总数为N的子载波上，进行点数与N相等的逆快速傅里叶变换IFFT，得到多用户时域信号，其中N表示单载波交织频分多址系统中多用户子载波的总数；
(3)发送已调制信号：
在多用户时域信号中加入长度为C的循环前缀，将经过并串转换后生成的多用户已调制信号以子帧为单位发送出去，其中C表示单载波交织频分多址系统中循环前缀占用的子载波数目；
(4)接收含有噪声的信号：
多用户已调制信号经过平坦信道或者频率选择性衰落信道，受到信道高斯白噪声的干扰，到达接收端，接收端接收到含有噪声的信号；
(5)多用户信号分离获得各用户时域信号：
(5a)将含有噪声的信号进行串并转换，去除长度为C的循环前缀，得到含有噪声的多用户时域信号；
(5b)将含有噪声的多用户时域信号进行点数等于N的快速傅里叶变换FFT，得到含有噪声的多用户频域信号；
(5c)将含有噪声的多用户频域信号，按照交织式解映射公式，将N条子载波分别解映射到M条子载波上，得到含有噪声的各用户频域信号；
(5d)将含有噪声的各用户频域信号进行点数等于M的逆快速傅里叶变换IFFT，再按照交织式映射公式，分别将各用户的M条子载波映射到总数为N的子载波上，进行点数与N相等的逆快速傅里叶变换IFFT，得到各用户时域信号；
(6)计算噪声的平均功率：
利用下述噪声平均功率公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号噪声的平均功率；
Figure FDA0002352779260000021
其中，
Figure FDA0002352779260000022
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号噪声的平均功率，v表示单载波交织频分多址系统中各用户的序号，其取值是v＝0,1,2,...,Q-1，Q表示子载波映射间隔，其数值大小等于单载波交织频分多址系统中可容纳的用户数，
Figure FDA0002352779260000023
d表示单载波交织频分多址系统中子帧的序号，其取值是d＝0,1,2,...,F-1，F表示单载波交织频分多址系统中子帧的总数，t表示单载波交织频分多址符号的序号，其取值是t＝0,1,2,...,L-1，L表示每个子帧中单载波交织频分多址符号的总数，k表示多用户子载波的序号，k的取值是k＝0,1,2,...,N-1，∑表示求和操作，||·||2表示求二阶矩操作，[·]表示求实部和虚部的绝对值操作，Yt,k表示各用户时域信号每个子帧的第t个单载波交织频分多址符号中第k条子载波，
Figure FDA0002352779260000024
表示各用户时域信号每个子帧的第t个单载波交织频分多址符号中第
Figure FDA0002352779260000025
条子载波，b表示按含有噪声的时域信号重复次数的序号，其取值是b＝0,1,2，
Figure FDA0002352779260000026
表示各用户时域信号中每个子帧的第t个单载波交织频分多址符号中第
Figure FDA0002352779260000027
条子载波，
Figure FDA0002352779260000028
表示各用户时域信号中每个子帧的第t个单载波交织频分多址符号中第
Figure FDA0002352779260000031
条子载波，u表示不同用户映射的起始位置，其取值范围是0≤u≤Q-1；
(7)计算单载波交织频分多址符号的平均功率：
利用下述两种公式中的任意一种符号平均功率公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号的平均功率：
Figure FDA0002352779260000032
Figure FDA0002352779260000033
其中，
Figure FDA0002352779260000034
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的平均功率；
(8)计算单载波交织频分多址符号的信噪比：
利用符号信噪比公式，计算各用户时域信号每个子帧中每个单载波交织频分多址符号的信噪比；
(9)计算每个子帧的信噪比：
利用子帧信噪比公式，计算各用户时域信号每个子帧的信噪比，所估计的信噪比指的是多用户已调制信号中各个用户受信道高斯白噪声干扰的信噪比。
2.根据权利要求1所述的单载波交织频分多址系统中的多用户信噪比盲估计方法，其特征在于：步骤(2)中所述的交织式映射公式如下：
Figure FDA0002352779260000035
其中，Xk表示交织式映射后总数为N的第k条子载波，Xp表示交织式映射前总数为M第p条子载波，p表示各用户子载波的序号，p的取值是p＝0,1,2,...,M-1。
3.根据权利要求1所述的单载波交织频分多址系统中的多用户信噪比盲估计方法，其特征在于：步骤(2)中所述的多用户时域信号是指，子载波的总数为N，单载波交织频分多址符号为F*L二维时频结构的多用户时域信号，其中，*表示相乘操作。
4.根据权利要求2所述的单载波交织频分多址系统中的多用户信噪比盲估计方法，其特征在于：步骤(5c)中所述的交织式解映射公式如下：
Figure FDA0002352779260000041
其中，
Figure FDA0002352779260000042
表示交织式解映射后总数为M第p条子载波，
Figure FDA0002352779260000043
表示交织式解映射前总数为N的第k条子载波。
5.根据权利要求1所述的单载波交织频分多址系统中的多用户信噪比盲估计方法，其特征在于：步骤(8)中所述的符号信噪比公式如下：
Figure FDA0002352779260000044
其中，
Figure FDA0002352779260000045
表示第v个用户时域信号中第d个子帧的第t个单载波交织频分多址符号的信噪比。
6.根据权利要求5所述的单载波交织频分多址系统中的多用户信噪比盲估计方法，其特征在于：步骤(9)中所述的子帧信噪比公式如下：
Figure FDA0002352779260000046
其中，
Figure FDA0002352779260000047
表示第v个用户时域信号中第d个子帧的信噪比。
CN201711144165.6A 2017-11-17 2017-11-17 单载波交织频分多址系统中的多用户信噪比盲估计方法 Active CN107947899B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711144165.6A CN107947899B (zh) 2017-11-17 2017-11-17 单载波交织频分多址系统中的多用户信噪比盲估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711144165.6A CN107947899B (zh) 2017-11-17 2017-11-17 单载波交织频分多址系统中的多用户信噪比盲估计方法

Publications (2)

Publication Number Publication Date
CN107947899A CN107947899A (zh) 2018-04-20
CN107947899B true CN107947899B (zh) 2020-04-07

Family

ID=61932753

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711144165.6A Active CN107947899B (zh) 2017-11-17 2017-11-17 单载波交织频分多址系统中的多用户信噪比盲估计方法

Country Status (1)

Country Link
CN (1) CN107947899B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109412984B (zh) * 2018-11-15 2021-01-01 中国电子科技集团公司第三十八研究所 一种基于Aitken加速法的多天线场景下盲信噪比估算方法
CN114338328B (zh) * 2022-01-05 2023-03-17 电子科技大学 一种非协作短波突发通信信号的失步检测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101534281A (zh) * 2009-04-20 2009-09-16 北京交通大学 一种ofdm系统基于梳状导频的分集式信道估计方法
CN103634260A (zh) * 2013-12-06 2014-03-12 电子科技大学 一种基于sim调制的sc-ifdma通信方法
CN103916342A (zh) * 2014-04-16 2014-07-09 中国科学院上海微系统与信息技术研究所 一种信噪比估计方法
WO2016068666A1 (ko) * 2014-10-31 2016-05-06 엘지전자 주식회사 Bpsk 또는 qbpsk 방식을 사용하는 무선 통신 시스템에서 블라인드 방식의 cfo 추정 방법
CN106302286A (zh) * 2015-05-29 2017-01-04 上海鑫皇实业有限公司 一种ofdm系统的信噪比盲估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101534281A (zh) * 2009-04-20 2009-09-16 北京交通大学 一种ofdm系统基于梳状导频的分集式信道估计方法
CN103634260A (zh) * 2013-12-06 2014-03-12 电子科技大学 一种基于sim调制的sc-ifdma通信方法
CN103916342A (zh) * 2014-04-16 2014-07-09 中国科学院上海微系统与信息技术研究所 一种信噪比估计方法
WO2016068666A1 (ko) * 2014-10-31 2016-05-06 엘지전자 주식회사 Bpsk 또는 qbpsk 방식을 사용하는 무선 통신 시스템에서 블라인드 방식의 cfo 추정 방법
CN106302286A (zh) * 2015-05-29 2017-01-04 上海鑫皇实业有限公司 一种ofdm系统的信噪比盲估计方法

Also Published As

Publication number Publication date
CN107947899A (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN101945066B (zh) 一种ofdm/oqam系统的信道估计方法
CN102724147B (zh) 一种水声正交频分复用的信道估计方法
CN102387115B (zh) 一种ofdm导频方案设计及信道估计方法
CN111147407B (zh) 基于信道预测的tmsbl水声ofdm时变信道估计方法
CN103200139B (zh) 一种ofdm信号带宽盲估计方法
CN104735014B (zh) 一种基于前导符号差分相关的定时同步方法
CN104410590A (zh) 一种基于压缩感知的短波ofdm抑制干扰联合信道估计方法
CN101808053A (zh) 基于ofdm的信道估计方法及装置
CN108562883B (zh) 一种多载波雷达系统的最大似然距离估计算法
CN107454032A (zh) 一种基于子载波间幅值积的ofdm频偏估计方法
CN101729456A (zh) 正交频分复用通信系统的信道估计方法
CN110736977A (zh) 一种适用于声呐通探的一体化信号产生方法
CN107947899B (zh) 单载波交织频分多址系统中的多用户信噪比盲估计方法
CN105187352A (zh) 一种基于ofdm前导的整数频偏估计方法
CN103117978A (zh) 基于多普勒频移补偿的水声通信多载波匹配方法及装置
CN100493056C (zh) 带时域包络加权的正交频分复用系统的频域信道估计方法
CN102377726B (zh) Ofdm系统的定时同步方法
CN102215184B (zh) 一种上行定时偏差的估计方法及系统
CN111416782A (zh) 一种基于空载波的ofdm系统频偏估计分析方法
CN105162533B (zh) 一种发射机幅度不平衡和相位不平衡测量方法
CN101026433B (zh) 一种用于自适应调制编码的信噪比估算方法
CN103916342A (zh) 一种信噪比估计方法
CN104270328B (zh) 一种信噪比实时估计方法
CN101945070B (zh) 一种噪声测量的方法及装置
US8731109B2 (en) Methods and systems for effective channel estimation in OFDM systems

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