CN101662441B - 信号估计方法、装置及检测系统 - Google Patents

信号估计方法、装置及检测系统 Download PDF

Info

Publication number
CN101662441B
CN101662441B CN200910174733A CN200910174733A CN101662441B CN 101662441 B CN101662441 B CN 101662441B CN 200910174733 A CN200910174733 A CN 200910174733A CN 200910174733 A CN200910174733 A CN 200910174733A CN 101662441 B CN101662441 B CN 101662441B
Authority
CN
China
Prior art keywords
variance
signal
frequency domain
transmission signal
estimated value
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
CN200910174733A
Other languages
English (en)
Other versions
CN101662441A (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.)
Huawei Technologies Co Ltd
Original Assignee
Huawei Technologies Co Ltd
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 Huawei Technologies Co Ltd filed Critical Huawei Technologies Co Ltd
Priority to CN200910174733A priority Critical patent/CN101662441B/zh
Publication of CN101662441A publication Critical patent/CN101662441A/zh
Application granted granted Critical
Publication of CN101662441B publication Critical patent/CN101662441B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明实施例公开了一种信号估计方法、装置及检测系统,属于通信领域。所述方法包括:接收经过离散傅里叶变换的频域接收信号;根据先验信息重建发送信号的均值和方差;根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差,根据重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和等效噪声的方差对频域的接收信号进行线性滤波和干扰抵消,得到发送信号的频域估计值;对发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值。本发明实施例根据接收信号估计发送信号时,在信道估计中考虑了信道衰落系数估计误差对信号估计过程的影响,使信号的估计结果更加准确。

Description

信号估计方法、装置及检测系统
技术领域
本发明涉及通信领域,尤其涉及一种信号估计方法、装置及检测系统。
背景技术
目前,在无线通信系统中使用多天线发送和多天线接收,可以成倍地提高通信系统的传输能力。为了进一步提高数据传输速率,采用带宽传输技术,例如OFDM(Orthogonal Frequency Division Multiplexing,正交频分复用)技术实现数据传输。由于OFDM系统发射信号峰均比较高,为了降低采用OFDM技术的系统中发射信号的峰均比,现有带宽传输采用具有低峰均比的多天线DFT-S-OFDM(Discrete Fourier Transform Spread OFDM,离散傅里叶变换扩展的正交频分复用)技术,这一技术虽然具有单载波传输技术低峰均比的特性,却又引入同一天线上不同时刻的信号干扰。
针对天线间存在的干扰以及同一天线上不同时刻的信号干扰存在多种估计方法,在现有的信号估计方法中,信道参数必须已知才能够对接收端的接收信号进行检测,进而估计出发送端发送的信号。信道参数在仿真时通常都假设精确已知,但在实际系统中,信道频域衰落系数通常通过信道估计方法得到。但是这样的信道估计存在误差,包含误差的信道估计值用于根据接收端的接收信号估计发送端的发送信号时,会影响估计的精度,降低接收机的性能。
发明内容
本发明的实施例提供一种信号估计方法及装置,在根据接收信号估计发送信号时,在信道估计中考虑信道衰落系数估计误差对信号估计过程的影响,使信号的估计结果更加准确。
本发明实施例提供一种信号估计方法,包括:
接收经过离散傅里叶变换的频域接收信号;
根据先验信息重建发送信号的均值和发送信号的方差;
根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差,其中,所述等效噪声包括信道频域衰落系数的估计误差和加性高斯噪声;
根据所述重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声的方差对所述频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值;
对所述发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值。
本发明实施例提供一种信号估计装置,包括:
接收模块,用于接收经过离散傅里叶变换的频域接收信号;
重建模块,用于根据先验信息重建发送信号的均值和发送信号的方差;
等效噪声方差获取模块,用于根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差,其中,所述等效噪声包括信道频域衰落系数的估计误差和加性高斯噪声;
处理模块,用于根据所述重建模块重建的发送信号的均值的傅里叶变换、所述重建模块重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声方差获取模块获取的等效噪声的方差对所述接收模块接收的频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值;
反变换模块,用于对所述处理模块得到的发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值。
本发明实施例提供一种检测系统,包括译码器和上述的信号估计装置;所述译码器,用于接收所述信号估计装置输出的信息;
所述译码器还用于向所述信号估计装置反馈先验信息。
本发明实施例提供的信号估计方法、装置及检测系统,在无线通信系统基站侧根据接收信号估计发送信号时,在信道估计中考虑了信道衰落系数估计误差对信号估计过程的影响,使信号的估计结果更加准确。
附图说明
图1为本发明实施例提供的一种信号估计方法流程图;
图2为本发明实施例提供的又一种信号估计方法流程图;
图3为本发明实施例提供的信号估计装置结构框图;
图4为本发明实施例提供的检测系统的结构框图。
具体实施方式
下面结合附图对本发明实施例信号估计方法及装置进行详细描述。
本发明实施例提供的信号估计方法,可以应用于基站侧的接收机中。为了抵抗噪声和干扰,信号接收机通常采用差错控制编码,而高效的差错控制编码通常使用软判决译码,即,为译码器输入比特位的软信息,相应地,要求检测器为译码器提供软信息。其中,软信息可以是信号中比特位1或0出现的概率,而非比特位(1或0)本身。例如似然比信息是软信息的一种,其为比特位中1出现的概率/0出现的概率。
在基站侧的信号接收端,检测器同译码器迭代工作的迭代检测译码接收机中,检测器是软输入软输出的。其中,软输入软输出表示检测器不仅能够输出软信息至译码器,而且能够使用译码器反馈的软信息作为先验信息输入检测器。其中,译码器为检测器反馈的先验信息可以认为是原始发送信号的信息比特经过编码所得编码比特的软信息。
如图1所示,为本发明实施例提供的一种信号估计方法流程图,该方法包括以下内容。
101:接收经过离散傅里叶变换的频域接收信号。
在本发明实施例中,可以是对基站侧的接收天线上的接收信号进行去循环前缀操作后,再对进行离散傅里叶变换,可以得到频域的接收信号,而检测器接收该频域接收信号,对该频域接收信号进行进一步处理。
102:根据先验信息重建发送信号的均值和发送信号的方差。
在本发明实施例中,发送信号可以是终端的发送天线发出的信号,该信号经过信道变换后,到达基站侧的接收天线,成为接收天线的接收信号。
其中,先验信息可以为译码器向检测器反馈的软信息,可以认为是终端发送天线的发送信号经过编码所得编码比特的软信息。
其中,根据先验信息重建发送信号的均值和发送信号的方差可以根据贝叶斯原理得到,例如可以为:
假设每Mc个编码比特映射到星座点的集合χ中的某个星座点α(c)上,形成发送信号xk,xk为发送信号矢量xt中的元素。若译码器生成发送信号xk的先验信息为Lk,j α,j=1,2,...M。根据先验信息Lk,j α,可以计算发送信号xk的概率如下:
P [ x k = α ( c ) ] = Π j = 1 M c 1 1 + exp ( - c ~ j L k , j α )
其中,α(c)是星座点集合χ中的某一点, c ~ j = + 1 c j = 1 - 1 c j = 0 , cj表示映射到星座点α(c)所对应的第j个比特,exp()表示取指数,∏表示连乘。
根据上述发送信号xk的概率,可以重建发送信号xk的均值μk和方差vk
μ k = Σ α ( c ) ∈ χ α ( c ) P [ x k = α ( c ) ]
v k = ( Σ α ( c ) ∈ χ | α ( c ) | 2 P [ x k = α ( c ) ] ) - | μ k | 2 ,
其中,||2表示绝对值的平方。
103:根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差,其中,等效噪声包括信道频域衰落系数的估计误差和加性高斯噪声。
其中,等效噪声的方差可以通过以下方式获取:
σ z ~ 2 ( k ) = N T σ λ ~ 2 + σ z 2 , 其中,表示信道频域衰落系数的估计误差的方差,表示等效噪声的方差,σz 2表示加性高斯噪声的方差,NT表示终端侧的发送天线个数,取自然数。当然,还可以通过其他方式获取等效噪声的方差,只要是在等效噪声的方差中可以体现加性高斯噪声的方差和信道频域衰落系数的估计误差的方差即可。
104:根据重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和等效噪声的方差对频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值;
其中,该步骤中,可以通过以下方式得到发送信号的频域估计值:
X ^ = { diag { α } Σ - 1 Λ ^ H ( y - Λ ^ η ) + η } ,
其中,η表示重建的发送信号的均值的傅里叶变换,y表示频域的接收信号,
Figure G2009101747336D00052
表示信道频域衰落系数的估计值矩阵, Σ - 1 = diag { Σ 1 - 1 , Σ 2 - 1 , . . . , Σ M - 1 } , Σ k = Λ ^ k H Λ ^ k V ‾ + σ z ~ 2 ( k ) I N T , k = 1,2 , . . . M , M表示每个用户在频域上所占用子载波的个数,()H表示对矩阵进行转置,
Figure G2009101747336D00055
表示等效噪声的方差,
Figure G2009101747336D00056
表示表示NT维的单位矩阵, V ‾ = diag { v ‾ 1 , v ‾ 2 , . . . , v ‾ M } , 表示重建的发送信号的方差,diag{}表示取对角阵, α = [ α 1,1 , . . . , α N T , 1 , . . . , α 1 , M , . . . , α N T , M ] T , α t , 1 = α t , 2 = . . . = α t , M = α t = 1 ρ t , ρ t = 1 M Σ k = 1 M e t H Σ - 1 Λ ^ k H Λ ^ k e t , et是NT维的列向量,其第t个元素为1,其余为零。
105:对发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值。
本发明实施例提供的信号估计方法,由于在无线通信系统基站侧根据接收信号估计发送信号时,在信道估计中考虑了信道衰落系数估计误差对信号估计过程的影响,使信号的估计结果更加准确。
如图2所示,为本发明实施例提供的信号估计方法流程图,该方法可以应用于多天线DFT-S-OFDM通信系统的软输入软输出检测器中。
本发明实施例假设多天线DFT-S-OFDM通信系统中有NT个终端侧的发送天线,基站侧有NR个接收天线,此时传输信道为多天线发送和多天线接收信道,即MIMO(Multiple Input Multiple Out-put)信道。需要说明的是,本发明实施例以MIMO信道为例说明,但并一定要限定在MIMO的场景中。该信号估计方法可以包括以下内容。
201:对接收信号进行离散傅里叶变换,得到频域的接收信号。
其中,在基站侧的NR个接收天线中的第r个接收天线上的接收信号,经过离散傅里叶变换为频域的接收信号,该信号可以为:
y r = Ω T F N Y r = Σ t = 1 N T Λ r , t F x t + z r - - - ( 1 )
其中,ΩT表示子载波反映射矩阵,FN为N×N的归一化DFT矩阵,N表示频域的子载波总个数,Yr表示第r个接收天线去除循环前缀的接收信号,xt表示第t个终端的发送天线的发送信号矢量,其元素个数为M;M表示每个用户在频域上所占用子载波的个数;F为M×M的归一化DFT矩阵;Λr,t为M×M的对角阵,其对角元素为第r个接收天线与第t个用户发送天线之间信道的频域衰落系数λr,t(k),k=1,2,...M,记为Λr,t=diag{λr,t(1),λr,t(2),...,λr,t(M)};zr为第r个接收天线上的加性高斯噪声,假设该加性高斯噪声的分布参数中均值为0,方差为σz 2;以vec{A}表示由矩阵A逐列排成的列矢量,设 y = vec { [ y 1 , y 2 , . . . , y N R ] T } , x = vec { [ x 1 , x 2 , . . . , x N T ] T } , z = vec { [ z 1 , z 2 , . . . , z N R ] T } , 则公式(1)可以写成矩阵和向量的表达式:
y = Λ F ~ x + z - - - ( 2 )
其中,Λ=diag{Λ1,Λ2,...,ΛM},
Figure G2009101747336D00066
F ~ = F ⊗ I N T , 其中,
Figure G2009101747336D00068
表示Kronecker乘积(矩阵的直积),表示NT维的单位矩阵。
202:根据先验信息重建发送信号的均值和发送信号的方差;
其中,发送信号可以是终端的发送天线发出的信号,该信号经过信道变换后,到达基站侧的接收天线,成为接收天线的接收信号。
先验信息可以为译码器向检测器反馈的软信息,可以认为是终端发送天线的发送信号经过编码所得编码比特的软信息。
重建的发送信号的均值和重建的发送信号的方差可以是根据先验信息重建得到的。根据先验重建发送信号的均值和发送信号的方差可以通过步骤102中的方法实现,此处不赘述。
假设μ表示利用先验信息重建发送信号的均值矢量,V表示利用先验信息重建发送信号的协方差矩阵:
μ = E { x } = [ μ 1,1 , . . . , μ N T , 1 , . . . , μ 1 , M , . . . , μ N T , M ] T - - - ( 3 )
V = Cov { x , x } = diag { v 1,1 , . . . , v N T , 1 , . . . , v 1 , M , . . . , v N T , M } - - - ( 4 )
其中,μt,m和vt,m分别表示终端侧的第t个发送天线的第m个发送符号的均值和方差。
在本发明实施例中,将信号估计的方法应用于检测器同译码器迭代工作的迭代检测译码接收机中,可以根据译码器输出的先验信息,重建发送信号的均值和发送信号的方差作为检测器的软输入,与译码器之间迭代工作。
203:根据重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、加性高斯噪声的方差以及信道频域衰落系数的估计误差的方差获取等效噪声的方差。
在本发明实施例中,等效噪声包括信道频域衰落系数的估计误差和加性高斯噪声,具体说明如下:
首先,构建包含信道频域衰落系数的估计误差的接收信号模型;
其中,信道频域衰落系数可以通过信道估计方法得到,而信道估计一般是存在误差,在信道估计存在误差的情况下,假设信道频域衰落系数λr,t(k)的估计误差模型为:
λ r , t ( k ) = λ ^ r , t ( k ) + λ ~ r , t ( k ) - - - ( 5 )
其中,
Figure G2009101747336D00081
为信道频域衰落系数λr,t(k)的估计值,
Figure G2009101747336D00082
为信道频域衰落系数λr,t(k)的估计误差,假设该估计误差
Figure G2009101747336D00083
的均值为0,方差为
Figure G2009101747336D00084
且各个之间统计独立。
定义: Λ ^ = diag { Λ ^ 1 , Λ ^ 2 , . . . , Λ ^ M } , Λ ~ = diag { Λ ~ 1 , Λ ~ 2 , . . . , Λ ~ M } , 其中,
Figure G2009101747336D00088
Figure G2009101747336D00089
其中,k=1,2,...M;M、NT、NR表示的含义均与上述201中相同,此处不再赘述。
因此,信道频域衰落系数的估计误差模型可记为矩阵形式,如下:
Λ = Λ ^ + Λ ~ - - - ( 6 )
Λ表示信道频域衰落系数矩阵;
Figure G2009101747336D000811
为信道频域衰落系数矩阵Λ的估计矩阵;
Figure G2009101747336D000812
为信道频域衰落系数矩阵Λ的估计误差矩阵。
然后,将信道频域衰落系数的估计误差模型(6)代入公式(2)中,可得包含信道频域衰落系数的估计误差的接收信号模型为:
y = Λ ^ F ~ x + Λ ~ F ~ x + z = Λ ^ F ~ x + z ~ - - - ( 7 )
其中, z ~ = Λ ~ F ~ x + z 表示信道频域衰落系数的估计误差对接收信号的影响与加性高斯噪声对接收信号的影响之和,可以将其称之为接收信号的等效噪声。
具体地,步骤203可以包括:
1)对重建的发送信号的均值矢量μ进行傅里叶变换,得到傅里叶变换之后的矢量,记为η;同时定义 η k = [ η 1 , k , η 2 , k , . . . , η N T , k ] T ; 其中,k=1,2,...M,表示子载波序号;
2)对重建的发送信号的方差vt,m在同一天线内进行平均处理,即:
v ‾ t = 1 M Σ m = 1 M v t , m ,
3)计算每个子载波上的等效噪声方差为 σ z ~ 2 ( k ) = σ λ ~ 2 ( η k H η k + Σ t = 1 N T v ‾ t ) + σ z 2 , 其中,k=1,2,...M表示子载波序号,
Figure G2009101747336D00092
表示信道频域衰落系数的估计误差的方差; η k = [ η 1 , k , η 2 , k , . . . , η N T , k ] T , σz 2表示加性高斯噪声的方差。
204:根据重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和等效噪声的方差对接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值和发送信号的估计误差的方差。
205:对线性滤波和干扰抵消后的发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值,输出发送信号的估计值。
上述步骤204的内容和步骤205的内容可以用公式表示:
由包含信道频域衰落系数的估计误差的接收信号模型(7),利用最小均方准则,可得发送信号x的估计值:
x ^ = F ~ H { diag { α } Σ - 1 Λ ^ H ( y - Λ ^ η ) + η } - - - ( 8 )
其中, Σ - 1 = diag { Σ 1 - 1 , Σ 2 - 1 , . . . , Σ M - 1 } , Σ k = Λ ^ k H Λ ^ k V ‾ + σ z ~ 2 ( k ) I N T , V ‾ = diag { v ‾ 1 , v ‾ 2 , . . . , v ‾ M } , α = [ α 1,1 , . . . , α N T , 1 , . . . , α 1 , M , . . . , α N T , M ] T , α t , 1 = α t , 2 = . . . = α t , M = α t = 1 ρ t , ρ t = 1 M Σ k = 1 M e t H Σ - 1 Λ ^ k H Λ ^ k e t , 各个参数中,( )H表示对矩阵进行转置运算;et是NT维的列向量,其第t个元素为1,其余为零;其他各个参数的含义均与上述相同,此处不再赘述。
公式(8)中,表示对接收信号在频域进行线性滤波;
Figure G2009101747336D000912
表示对接收信号在频域进行干扰抵消;
Figure G2009101747336D000913
表示对接收信号进行反离散傅里叶变换。
在公式(8)中,残余的干扰和噪声可以用高斯分布近似,将公式(8)的结果
Figure G2009101747336D00101
等效为发送信号x经过加性高斯噪声信道的输出,即:
x ^ t , m = Δ x t , m + n t , m - - - ( 9 )
其中,
Figure G2009101747336D00103
表示发送信号x的估计值矩阵中的元素;xt,m表示发送信号x矩阵中的元素;nt,m表示发送信号的估计误差,即残余的干扰和噪声矩阵中的元素;
Figure G2009101747336D00105
表示近似等于。
进一步地,根据公式(8)还可以得到发送信号估计误差的方差为:
σ n , t , m 2 = α t - v ‾ t - - - ( 10 )
其中,各个参数的含义与上述表示相同,此处不再赘述。
206:根据发送信号的估计值和发送信号估计误差的方差,计算每个比特的对数似然比。
其中,步骤206的计算方法与现有的方法基本相同,此处不再赘述。
此外,计算后得到的对数似然比即比特的软信息,可以将该软信息作为译码器的输入,经过译码器输出后再经过解交织器得到先验信息,利用该先验信息可以用于后续过程中重建发送信号的均值和方差。
在迭代检测译码接收机中,基于信道频域衰落系数的估计误差的信号估计方法,可以将信道频域衰落系数的估计误差的方差、利用译码器反馈的先验信息重建的发送信号的均值和方差作为检测器的软输入,将
Figure G2009101747336D00107
和σn,t,m 2作为检测器的软输出,与译码器之间进行迭代。在检测器与译码器迭代工作的过程中,译码器输出的先验信息是不断更新的。对于初始的检测,由于译码器没有输出,此时,可以假设根据先验信息重建的发送信号的均值为0,根据先验信息重建的发送信号的方差为1,将该均值和方差作为检测器的软输入。
本实施例在多天线DFT-S-OFDM系统的信号估计方法中,由于在无线通信系统基站侧根据接收信号估计发送信号时,在信道估计中考虑了信道衰落系数估计误差对信号估计过程的影响,使信号的估计结果更加准确,同时满足了无线通信系统迭代接收机对高性能以及低处理复杂度的要求。
本发明实施例提供一种信号估计装置,如图3所示,该装置包括:
接收模块31,用于接收经过离散傅里叶变换的频域接收信号;
重建模块32,用于根据先验信息重建发送信号的均值和发送信号的方差;
等效噪声方差获取模块33,用于根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差,其中,等效噪声包括信道频域衰落系数的估计误差和加性高斯噪声;
处理模块34,用于根据重建模块32重建的发送信号的均值的傅里叶变换、重建模块32重建的发送信号的方差、信道频域衰落系数的估计值和等效噪声方差获取模块获取的等效噪声的方差对接收模块接收的频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值;
反变换模块35,用于对处理模块34得到的发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值。
其中,处理模块34可以包括:
线性滤波单元,用于根据重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和等效噪声的方差对频域的接收信号在频域进行线性滤波;
干扰抵消单元,用于根据重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和等效噪声的方差对频域的接收信号在频域进行干扰抵消。
进一步地,处理模块34,还可以用于根据重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和等效噪声的方差对频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号估计误差的方差。
进一步地,该信号估计装置还可以包括:
对数似然比计算模块,根据处理模块得到的发送信号的估计值和发送信号估计误差的方差,计算每个比特的对数似然比,后续的装置(例如译码器)可以根据对数似然比更新先验信息。
其中,等效噪声方差获取模块33可以通过以下方式获取等效噪声方差:
根据加性高斯噪声的方差σz 2和信道频域衰落系数的估计误差的方差
Figure G2009101747336D00121
获取等效噪声的方差 σ z ~ 2 ( k ) = N T σ λ ~ 2 + σ z 2 , 其中,NT表示终端侧的发送天线个数;或者,
根据重建的发送信号的均值的傅里叶变换η、重建的发送信号的方差
Figure G2009101747336D00123
加性高斯噪声的方差σz 2以及信道频域衰落系数的估计误差的方差
Figure G2009101747336D00124
获取等效噪声的方差 σ z ~ 2 ( k ) = σ λ ~ 2 ( η k H η k + Σ t = 1 N T v ‾ t ) + σ z 2 ,
其中,k表示子载波序号,NT表示终端侧的发送天线个数,
Figure G2009101747336D00126
表示对所有发送天线内重建的发送信号的方差求和。
上述信号估计装置中的各个模块的工作原理可以参考前述的信号估计方法实施例中的相关描述,此处不再赘述。
本实施例提供的信号估计装置可以是基于信道估计存在误差的检测器,由于在无线通信系统基站侧根据接收信号估计发送信号时,在信道估计中考虑了信道衰落系数估计误差对信号估计过程的影响,使信号的估计结果更加准确。
本发明实施例还提供一种检测系统,该系统可以包括译码器以及图3实施例所示的信号估计装置,该译码器用于接收信号估计装置输出的信息;
该译码器还用于向信号估计装置反馈先验信息。
进一步地,该检测系统可以是检测器同译码器迭代工作的迭代检测译码接收机。在该接收机中,信号估计装置可以是检测器,该检测器是软输入软输出的,即检测器不仅能够输出软判决信息至译码器,而且能够使用译码器反馈的软信息作为先验信息。如图4所示,为本发明实施例提供的一种检测系统,该系统可以是存在信道估计误差的迭代接收机,该迭代接收机中包括:对接收信号进行去循环前缀部分,OFDM解调部分,基于信道估计存在误差的检测器(即信号估计装置)、交织器、解交织器以及译码器。
其中,OFDM解调部分包括对接收信号的串并转换、傅里叶变换、子载波反映射以及并串部分;
交织器与解交织器可以把信号序列重新组织,以提高信号接收性能;
基于信道估计存在误差的检测器(即信号估计装置)向译码器方向输出信号比特的软信息,例如比特的对数似然比;由于在实际通信系统中,信道频域衰落系数是通过信道估计方法得到的,而信道估计总是存在误差,该误差会降低接收性能,因此本发明实施例在检测器内部引入信道衰落系数估计误差对接收信号检测的影响,以提高信号的接收性能;
译码器向检测器反馈软信息,该软信息为先验信息,在译码器与检测器迭代工作的过程中,该先验信息是不断更新的,检测器可以根据先验信息重建发送信号的均值和方差作为检测器的软输入。
本实施例提供的包含信号估计装置的迭代接收机,在根据接收信号估计发送信号时,在信道估计中考虑了信道衰落系数估计误差对信号估计过程的影响,使信号的估计结果更加准确,满足了无线通信系统迭代接收机对高性能以及低处理复杂度的要求。
本发明实施例可以利用软件实现,相应的软件程序可以存储在可读取的存储介质中,例如,计算机的硬盘、缓存或光盘中。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,在没有超过本申请的精神和范围内,可以通过其他的方式实现。当前的实施例只是一种示范性的例子,不应该作为限制,所给出的具体内容不应该限制本申请的目的。例如,多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。
另外,所描述系统,装置和方法以及不同实施例的示意图,在不超出本申请的范围内,可以与其它系统,模块,技术或方法结合或集成。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。

Claims (11)

1.一种信号估计方法,其特征在于,所述方法包括:
接收经过离散傅里叶变换的频域接收信号;
根据先验信息重建发送信号的均值和发送信号的方差;
根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差,其中,所述等效噪声包括信道频域衰落系数的估计误差和加性高斯噪声;
根据所述重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声的方差对所述频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值;
对所述发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值;
所述根据所述重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声的方差对所述频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值包括:
根据
Figure FDA0000130814540000011
对接收信号y在频域进行线性滤波,根据
Figure FDA0000130814540000012
对接收信号在频域进行干扰抵消,获得所述发送信号的频域估计值 X ^ = { diag { α } Σ - 1 Λ ^ H ( y - Λ ^ η ) + η } ,
其中,η表示重建的发送信号的均值的傅里叶变换,y表示频域的接收信号,
Figure FDA0000130814540000014
表示信道频域衰落系数的估计值矩阵,
Figure FDA0000130814540000015
Figure FDA0000130814540000016
k=1,2,...,M,M表示每个用户在频域上所占用子载波的个数,()H表示对矩阵进行转置,
Figure FDA0000130814540000017
表示等效噪声的方差,
Figure FDA0000130814540000018
表示表示NT维的单位矩阵,
Figure FDA0000130814540000019
表示重建的发送信号的方差,diag{ }表示取对角阵, α = [ α 1,1 , . . . , α N T , 1 , . . . , α 1 , M , . . . , α N T , M ] T , α t , 1 = α t , 2 = . . . = α t , M = α t = 1 ρ t ,
Figure FDA00001308145400000112
et是NT维的列向量,其第t个元素为1,其余为零。
2.根据权利要求1所述的信号估计方法,其特征在于,所述根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差包括:
根据加性高斯噪声的方差
Figure FDA0000130814540000021
和信道频域衰落系数的估计误差的方差
Figure FDA0000130814540000022
获取等效噪声的方差其中,NT表示终端侧的发送天线个数。
3.根据权利要求1所述的信号估计方法,其特征在于,所述根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差包括:
根据重建的发送信号的均值的傅里叶变换η、重建的发送信号的方差
Figure FDA0000130814540000024
加性高斯噪声的方差
Figure FDA0000130814540000025
以及信道频域衰落系数的估计误差的方差
Figure FDA0000130814540000026
获取等效噪声的方差
Figure FDA0000130814540000027
其中,k表示子载波序号,NT表示终端侧的发送天线个数,
Figure FDA0000130814540000028
表示对所有发送天线内重建的发送信号的方差求和。
4.根据权利要求1-3中任一项所述的信号估计方法,其特征在于,所述方法还包括:对所述频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号估计误差的方差
Figure FDA0000130814540000029
其中,
5.根据权利要求4所述的信号估计方法,其特征在于,根据所述发送信号的估计值和所述发送信号估计误差的方差,计算每个比特的对数似然比,根据所述对数似然比更新所述先验信息。
6.一种信号估计装置,其特征在于,包括:
接收模块,用于接收经过离散傅里叶变换的频域接收信号;
重建模块,用于根据先验信息重建发送信号的均值和发送信号的方差;
等效噪声方差获取模块,用于根据加性高斯噪声的方差和信道频域衰落系数的估计误差的方差获取等效噪声的方差,其中,所述等效噪声包括信道频域衰落系数的估计误差和加性高斯噪声;
处理模块,用于根据所述重建模块重建的发送信号的均值的傅里叶变换、所述重建模块重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声方差获取模块获取的等效噪声的方差对所述接收模块接收的频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号的频域估计值;
反变换模块,用于对所述处理模块得到的发送信号的频域估计值进行反离散傅里叶变换,得到发送信号的估计值;
所述处理模块具体用于:
根据
Figure FDA0000130814540000031
对接收信号y在频域进行线性滤波,根据
Figure FDA0000130814540000032
对接收信号在频域进行干扰抵消,获得所述发送信号的频域估计值 X ^ = { diag { α } Σ - 1 Λ ^ H ( y - Λ ^ η ) + η } ,
其中,η表示重建的发送信号的均值的傅里叶变换,y表示频域的接收信号,
Figure FDA0000130814540000034
表示信道频域衰落系数的估计值矩阵, Σ - 1 = diag { Σ 1 - 1 , Σ 2 - 1 , . . . Σ M - 1 } , Σ k = Λ ^ k H Λ ^ k V ‾ + σ z ~ 2 ( k ) I N T , k=1,2,...,M,M表示每个用户在频域上所占用子载波的个数,( )H表示对矩阵进行转置,表示等效噪声的方差,
Figure FDA0000130814540000038
表示表示NT维的单位矩阵,
Figure FDA0000130814540000039
表示重建的发送信号的方差,diag{ }表示取对角阵, α = [ α 1,1 , . . . , α N T , 1 , . . . , α 1 , M , . . . , α N T , M ] T , α t , 1 = α t , 2 = . . . = α t , M = α t = 1 ρ t ,
Figure FDA00001308145400000312
et是NT维的列向量,其第t个元素为1,其余为零。
7.根据权利要求6所述的信号估计装置,其特征在于,所述处理模块包括:
线性滤波单元,用于根据所述重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声的方差对所述频域的接收信号在频域进行线性滤波;
干扰抵消单元,用于根据所述重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声的方差对所述频域的接收信号在频域进行干扰抵消。
8.根据权利要求6所述的信号估计装置,其特征在于,所述处理模块,还用于根据所述重建的发送信号的均值的傅里叶变换、重建的发送信号的方差、信道频域衰落系数的估计值和所述等效噪声的方差对所述频域的接收信号在频域进行线性滤波和干扰抵消,得到发送信号估计误差的方差。
9.根据权利要求8所述的信号估计装置,其特征在于,所述信号估计装置还包括:
对数似然比计算模块,根据所述处理模块得到的发送信号的估计值和所述发送信号估计误差的方差,计算每个比特的对数似然比。
10.根据权利要求6-9任一项所述的信号估计装置,其特征在于,所述等效噪声方差获取模块通过以下方式获取所述等效噪声方差:
根据加性高斯噪声的方差
Figure FDA0000130814540000041
和信道频域衰落系数的估计误差的方差
Figure FDA0000130814540000042
获取等效噪声的方差
Figure FDA0000130814540000043
其中,NT表示终端侧的发送天线个数;或者,
根据重建的发送信号的均值的傅里叶变换η、重建的发送信号的方差、加性高斯噪声的方差
Figure FDA0000130814540000045
以及信道频域衰落系数的估计误差的方差
Figure FDA0000130814540000046
获取等效噪声的方差
Figure FDA0000130814540000047
其中,k表示子载波序号,NT表示终端侧的发送天线个数,
Figure FDA0000130814540000048
表示对所有发送天线内重建的发送信号的方差求和。
11.一种检测系统,其特征在于,所述系统包括译码器和如权利要求6-10中任意一项权利要求所述的信号估计装置;所述译码器,用于接收所述信号估计装置输出的信息;
所述译码器还用于向所述信号估计装置反馈先验信息。
CN200910174733A 2009-09-11 2009-09-11 信号估计方法、装置及检测系统 Active CN101662441B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910174733A CN101662441B (zh) 2009-09-11 2009-09-11 信号估计方法、装置及检测系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910174733A CN101662441B (zh) 2009-09-11 2009-09-11 信号估计方法、装置及检测系统

Publications (2)

Publication Number Publication Date
CN101662441A CN101662441A (zh) 2010-03-03
CN101662441B true CN101662441B (zh) 2012-10-17

Family

ID=41790236

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910174733A Active CN101662441B (zh) 2009-09-11 2009-09-11 信号估计方法、装置及检测系统

Country Status (1)

Country Link
CN (1) CN101662441B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2533159C2 (ru) * 2010-09-16 2014-11-20 ЗетТиИ Уистрон Телеком АБ Способ и система для улучшенного подавления помех посредством выбора пути
CN102546489B (zh) * 2012-01-10 2015-01-28 华为技术有限公司 一种无线通信中解调有效噪声的计算方法及装置
CN113542164B (zh) * 2021-07-13 2023-03-31 电子科技大学 一种盲信号检测中的通信参数估计方法及系统
CN114500192B (zh) * 2022-03-23 2023-06-02 广州海格通信集团股份有限公司 一种信道估计方法及终端设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101312440A (zh) * 2007-05-24 2008-11-26 中国科学院微电子研究所 一种对正交频分复用通信系统信号信噪比进行估计的方法
CN101330358A (zh) * 2007-06-22 2008-12-24 中兴通讯股份有限公司 一种干扰抑制分集合并的方法及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101312440A (zh) * 2007-05-24 2008-11-26 中国科学院微电子研究所 一种对正交频分复用通信系统信号信噪比进行估计的方法
CN101330358A (zh) * 2007-06-22 2008-12-24 中兴通讯股份有限公司 一种干扰抑制分集合并的方法及系统

Also Published As

Publication number Publication date
CN101662441A (zh) 2010-03-03

Similar Documents

Publication Publication Date Title
CN101841375B (zh) 一种多输入多输出单载波分块传输系统的检测方法及装置
CN101505290B (zh) 改进的宽带mimo中频偏估计方法
CN101350800B (zh) 载波间干扰处理装置和方法以及无线通信系统的接收机
KR100918717B1 (ko) 다입다출력 직교주파수분할다중화 이동통신 시스템에서의신호 시퀀스 추정 방법 및 장치
CN100373840C (zh) 归一化迭代软干扰抵消信号检测方法和装置
US8325588B2 (en) Iterative reception method and iterative receiver
CN102624666B (zh) 稀疏信道模型下多路收发的正交多载波水声通信循环译码方法
CN106302298A (zh) 一种消除正交频分复用水声通信系统限幅噪声的方法
CN104767587A (zh) 基于ofdm系统下联合信道编译码的压缩感知信道估计方法
CN100571098C (zh) 通信系统中低复杂度的极大似然检测方法及装置
CN101662441B (zh) 信号估计方法、装置及检测系统
CN105322991A (zh) 基于wfrft预编码的多输入多输出传输系统及传输方法
CN101589571A (zh) 用于频率选择性和快速衰落信道上的单载波空间频率块编码传输的系统和方法
CN106972875B (zh) 一种mimo系统下多维联合估计动态稀疏信道的方法
US8040967B2 (en) Noise power interpolation in a multi-carrier system
CN101267409A (zh) 一种mimo-ofdm双选择性信道的跟踪方法
CN102790746B (zh) Ofdm系统的信道估计方法
CN100499611C (zh) 无线通信系统空域最大后验概率检测方法
CN109802901B (zh) 基于到达角测量的3d mimo信道估计方法及系统
CN103326987A (zh) 基于循环星座移位的峰平比抑制方法
EP2148485B1 (en) Method and apparatus for OFDM based on 3D signal constellations
CN102916924B (zh) 可调带宽的偏移调制多载波频分多址传输方法
CN110138461A (zh) 自适应多入多出与正交频分复用相结合的水声通信方法
US9166841B2 (en) Receiving apparatus and receiving method
CN102065047B (zh) 一种双发射天线ofdm信号的发射方法及其发射装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant