具体实施方式
下面参照附图来对本发明进行详细描述。
图3a-b示出MIMO-OFDM系统中发射机及接收机的基带功能框图。本发明的内容主要涉及发射机的导频插入部分和接收机中的信道估计部分。
OFDM系统是一种多载波系统,为了描述方便起见,我们将OFDM系统中的子载波分成有效子载波和虚拟子载波两类,前者用来传输数据和导频,后者是作为保护带宽不作传输用(即传输0)。根据与直流子载波的关系,有效子载波可进一步分成正的有效子载波,直流子载波和负的有效子载波。
图4示出OFDM系统中的各个子载波的相互关系。
系统参数说明
为了方便描述,首先给出下面将要用到的系统参数及其说明,如下表1所示。另外,需要指出的是,单天线OFDM系统实际上可认为是多天线OFDM系统的一个特例(也即N=1),因此下面将仅针对多天线MIMO-OFDM系统阐述本发明的内容。
表1:系统参数说明
参数 |
说明 |
Nfft |
OFDM系统的FFT变换的尺度大小 |
2Nh+1 |
OFDM系统中有效子载波的数目(其中包含直流子载波,不含虚拟子载波) |
Df |
频域导频间隔(以子载波数计算),一般为2的整数次幂,如4,8等。 |
NT |
发送天线数 |
NR |
接收天线数 |
2Mp+1 |
有效子载波中导频子载波的个数 |
系统模型描述:
考虑一个NT个发射天线,NR个接收天线的MIMO-OFDM系统,子载波总数(即FFT变换矩阵的大小)假设为Nfft,其中有效子载波(即实际用于数据和导频传输的子载波)个数假设为2Nh+1,其他的子载波作为保护带宽被称为虚拟子载波。假设所有子载波的序号取值范围为0到Nfft-1,其中0号子载波对应直流子载波。为了描述方便,这里首先定义几个特定的子载波序号集合,如下所示。
这里,
和
分别表示正有效子载波(包括直流子载波)的序号集合,负有效子载波序号集合以及虚拟子载波序号的集合。那么
表示所有有效子载波的序号集合。假设在接收机OFDM解调后接收天线υ的第l个OFDM符号的第i个子载波上的符号表示为Y
l,i (υ),则该符号可以表示为
这里,Hl,i (μ,υ)表示发射天线μ与接收天线υ之间信道的频率响应(CFR),Nl,i (υ)表示零均值,方差为σn 2的加性高斯白噪声(AWGN)。在此,我们仅仅描述一个接收天线在一个OFDM符号上的信道估计方法,其他的接收天线和其他的OFDM符号上的信道估计方法与之完全相同。因此,以下描述将略掉接收天线序号υ和OFDM符号序号l。
相位偏移的多天线导频序列
假设导频序列以Df为间隔均匀分布于OFDM符号的频域子载波中,并假设有效子载波中包含的有效导频子载波的个数为2Mp+1,序号集合表示为
其中
和
分别表示正的和负的导频子载波的序号集合,可以分别表示为公式(6)和(7)
同时,处于虚拟子载波中的导频(以下称作虚拟导频)序号集合可由公式(8)来表示:
那么NT个发射天线上的相位偏移的导频序列可表示为公式(9)
这里,
表示发送天线0上的导频序列,该导频序列的平均功率表示为
增强的MIMO-OFDM系统信道估计算法
接收机在OFDM解调后,有效导频子载波位置处的接收符号可以用矩阵表示为
这里,
和
分别表示发射天线μ发送的导频序列,信道频率响应CFR和AWGN.diag(.)表示以输入的序列作为主对角元素而构成的对角方阵。上标T表示矩阵转置操作。
图5示出根据本发明一个具体实施方式的在OFDM接收机中的增强的信道估计方法的流程图。下面参照图5并结合图3-4以及上文的OFDM系统参数和模型描述来对该方法进行描述。
如图3a所示,在OFDM发射机中有N个(其中N=1,2,3...)发射天线,用于发射多路OFDM符号序列。而如图3b所示,在OFDM接收机中,有M个(其中M=1,2,3...)个接收天线,其中每个接收天线都接收到来自上述N个发射天线发来的N路OFDM符号序列并进行合并处理,得到一个N天线OFDM符号序列,在对所述N天线OFDM符号进行时频域同步、循环前缀去除和逆离散傅立叶变换(可成为OFDM解调)后,获得一路经OFDM解调的OFDM符号序列。根据本发明的在OFDM接收机中的信道估计方法即从对所述经OFDM解调的符号序列进行处理的步骤S1处开始。应该理解,以下处理过程是基于相位偏移的多天线导频序列,其仅针对第1个接收天线与所有发射天线(第1-第N个发射天线)之间的所有信道的信道估计,其他接收天线的信道估计过程与此相同,另外,单发射天线OFDM系统仅仅是多发射天线OFDM系统的一个特例,可以由下面的描述很容易导出在单发射天线OFDM系统中的信道估计方法。
在步骤S11中,首先由经OFDM解调后的符号序列获得其中的有效导频子载波出的信号值Yep。
随后,在步骤S12中,通过将所述有效导频子载波的信号值除以发射端第一个天线所发射的OFDM信号的有效导频子载波的信号值,以生成一个有效导频子载波的信道频率响应估计值,如公式(11)所示:
对于多发射天线OFDM系统,该信道频率响应估计值应是该接收天线与所有发射天线之间的所有信道的信道频率响应的相位偏移叠加的估计值(最小二乘估计),对于单发射天线OFDM系统,信道频率响应估计算即为该接收天线与发射天线之间的信道频率响应的估计值(最小二乘估计)。
接着,在步骤S13中,利用有效导频子载波与虚拟导频子载波的相关性,由步骤S12中所获得的有效导频子载波的信道频率响应估计值来估算虚拟导频子载波处的信道频率响应估计值。
具体的,可采用公式(16)所示的Robust MMSE算法将虚拟子载波中没有传输的导频点上的信道频率响应(CFR)进行估计和预测为
也即,将有效导频子载波的信道频率响应估计值(对于多发射天线OFDM系统,为所有信道的信道频率响应的相位偏移叠加的估计值)左乘一个相关预测矩阵,得到虚拟导频子载波上的信道频率响应估计值。
这里,
其中,公式(17)和(18)表示该接收天线与所有发射天线之间的相位偏移叠加的所有信道的相关矩阵,对于单发射天线OFDM系统来说,
公式(17,18)中,上标H表示矩阵的共扼转置,Λvp (μ)和Λep (μ)是对角阵,定义如下
注意,对于单发射天线OFDM来说,上述公式(19)和(20)都变成了单位矩阵,因为μ=0。
RHvp(μ)Hep(μ)和RHep(μ)Hep(μ)表示信道的稳健的相关矩阵,其计算方法参见稳健的相关矩阵的计算部分。γ表示导频功率与噪声功率之比。
接着,在步骤S14中,对在步骤S12所获得的在有效导频子载波处的信道频率响应估计值与在步骤S13中所获得的在虚拟导频子载波处的信道频率响应估计值进行频时变换,以获得时域的信道冲击响应(CIR)估计,具体如下列公式(21)来描述:
求增强的信道估计(算法A)所得的信道冲击响应(CIR)
也即,对由步骤S12获得的正负有效导频子载波处的信道频率响应估计值(最小二乘估计)和步骤S13获得的虚拟导频子载波处的信道频率响应估计值(MMSE估计)所构成的向量进行离散傅立叶反变换(IDFT变换)。
这里,FIDFT表示Nfft/Df点的IDFT变换矩阵,如下所示
对于多发射天线OFDM系统,步骤S14所获得的时域的信道冲击响应估计应为该接收天线与所有发射天线之间的所有信道的信道冲击响应估计的时间偏移的叠加,如图7所示。而对于单发射天线OFDM系统,步骤S14所获得的时域的信道冲击响应估计即是该该接收天线与所有发射天线之间的所有信道的信道冲击响应估计。
最后,在步骤S15中,按时间偏移由所述时域的信道冲击响应估计提取出各个信道的信道冲击响应(对于单发射天线,所述时域的的信道冲击响应估计只包含该接收天线与发射天线之间的一个信道的信道冲击响应),并将所述各个信道的冲击响应分别进行时频变换,以生成各个信道的所有子载波的信道频率响应估计。
具体的,通过下面的公式(23)得到接收天线υ与所有发射天线所构成的NT个信道的CFR估计为
这里,M=N
fft/D
f/N
T,
表示列向量
的序号为μ·M到(μ+1)·M-1的元素组成的子向量,F
DFT表示N
fft点的DFT变换矩阵,表示为
公式(23)的含义是取信道冲击向量
中相应部分
然后后面添加N
fft-M个零,最后进行N
fft长度的离散傅立叶变换(DFT),即得到一个收发天线对之间的信道频率响应估计。
在另一个优选实施方式中(在此,将前述实施方式称为方案1,而本实施方式称为方案2,以示区别),在步骤S14之前,还可对步骤S12所获得的有效导频子载波处的信道频率响应进行滤波处理,以生成经过滤波的有效子载波的信道频率响应估计值。随后,在步骤S14中,对所述经过滤波的有效导频子载波和虚拟导频子载波的信道频率响应估计值进行频时变换,以获得精度更高的时域信道冲击响应估计。从而最后获得精度更高的信道频率响应估计。具体实现方式如下所述:
在步骤S13中,可采用稳健的MMSE算法将虚拟子载波中没有传输的导频点上的信道信息(CFR)进行估计和预测,同时将有效子载波位置上的最小二乘估计进行滤波,如下面公式(25)所示
(这个公式的物理意义公式前面已经说明,即对由有效导频子载波上的最小二乘估计构成的列向量左乘一个预测和滤波矩阵,将得到预测的虚拟子载波中导频位置上的信道频率响应和经滤波的有效导频子载波处的信道频率响应。)
这里,
对于单发OFDM来说,公式(26)变为
其中,
是有效导频和虚拟导频的子载波序号集合。RHp(μ)Hep(μ)表示信道的稳健的相关矩阵,其计算方法参见稳健的相关矩阵的计算部分。
在步骤S14中,求增强的信道估计所得的信道冲击响应(CIR)
该公式的物理意义是对列向量
进行离散傅立叶反变换(IDFT),得到由要估计的多个信道的冲击响应组成的信道冲击向量
图6示出根据本发明一个具体实施方式的在OFDM接收机中的增强的信道估计装置的框图。所述信道估计装置1包括一个第一计算装置(在一个优选实施例中,第一计算装置可由一个获取装置1l与一个除法装置12来实现)、一个第二计算装置13、一个频时变换装置14和一个时频变换装置15。下面参照图6并结合图3-4以及上文的OFDM系统参数和模型描述来对该信道估计装置1进行描述。
如图3a所示,在OFDM发射机中有N个(其中N=1,2,3...)发射天线,用于发射多路OFDM符号序列。而如图3b所示,在OFDM接收机中,有M个(其中M=1,2,3...)个接收天线,其中每个接收天线都接收到来自上述N个发射天线发来的N路OFDM符号序列并进行合并处理,得到一个N天线OFDM符号序列,在对所述N天线OFDM符号进行时频域同步、循环前缀去除和逆离散傅立叶变换(可成为OFDM解调)后,获得一路经OFDM解调的OFDM符号序列。根据本发明的在OFDM接收机中的信道估计方法即从对所述经OFDM解调的符号序列进行处理的步骤S1处开始。应该理解,以下处理过程是基于相位偏移的多天线导频序列,其仅针对第1个接收天线与所有发射天线(第1-第N个发射天线)之间的所有信道的信道估计,其他接收天线的信道估计过程与此相同,另外,单发射天线OFDM系统仅仅是多发射天线OFDM系统的一个特例,可以由下面的描述很容易导出在单发射天线OFDM系统中的信道估计方法。
首先,获取装置11由经OFDM解调后的符号序列获得其中的有效导频子载波出的信号值Yep。
随后,除法装置12通过将所述有效导频子载波的信号值除以发射端第一个天线所发射的OFDM信号的有效导频子载波的信号值,以生成一个有效导频子载波的信道频率响应估计值,如公式(11)所示:
对于多发射天线OFDM系统,该信道频率响应估计值应是该接收天线与所有发射天线之间的所有信道的信道频率响应的相位偏移叠加的估计值(最小二乘估计),对于单发射天线OFDM系统,信道频率响应估计算即为该接收天线与发射天线之间的信道频率响应的估计值(最小二乘估计)。
第二计算装置利用有效导频子载波与虚拟导频子载波的相关性,由第一计算装置所得的有效导频子载波的信道频率响应估计值来估算虚拟导频子载波处的信道频率响应估计值。
具体的,可采用公式(16)所示的Robust MMSE算法将虚拟子载波中没有传输的导频点上的信道频率响应(CFR)进行估计和预测为
也即,将有效导频子载波的信道频率响应估计值(对于多发射天线OFDM系统,为所有信道的信道频率响应的相位偏移叠加的估计值)左乘一个相关预测矩阵,得到虚拟导频子载波上的信道频率响应估计值。
这里,
其中,公式(17)和(18)表示该接收天线与所有发射天线之间的相位偏移叠加的所有信道的相关矩阵,对于单发射天线OFDM系统来说,
公式(17,18)中,上标H表示矩阵的共扼转置,Λvp (μ)和Λep (μ)是对角阵,定义如下
注意,对于单发射天线OFDM来说,上述公式(19)和(20)都变成了单位矩阵,因为μ=0。
RHvp(μ)Hep(μ)和RHep(μ)Hep(μ)表示信道的稳健的相关矩阵,其计算方法参见稳健的相关矩阵的计算部分。γ表示导频功率与噪声功率之比。
频时变换装置14对在除法装置12所得的在有效导频子载波处的信道频率响应估计值与第二计算装置13所得的在虚拟导频子载波处的信道频率响应估计值进行频时变换,以获得时域的信道冲击响应估计,具体如下列公式(21)来描述:
也即,对由除法装置12所得的正负有效导频子载波处的信道频率响应估计值(最小二乘估计)和第二计算装置13所得的虚拟导频子载波处的信道频率响应估计值(MMSE估计)所构成的向量进行离散傅立叶反变换(IDFT变换)。
这里,FIDFT表示Nfft/Df点的IDFT变换矩阵,如下所示
对于多发射天线OFDM系统,步骤S14所获得的时域的信道冲击响应估计应为该接收天线与所有发射天线之间的所有信道的信道冲击响应估计的时间偏移的叠加,如图7所示。而对于单发射天线OFDM系统,步骤S14所获得的时域的信道冲击响应估计即是该该接收天线与所有发射天线之间的所有信道的信道冲击响应估计。
时频变换装置15按时间偏移由所述时域的信道冲击响应估计提取出各个信道的信道冲击响应(对于单发射天线,所述时域的的信道冲击响应估计只包含该接收天线与发射天线之间的一个信道的信道冲击响应),并将所述各个信道的冲击响应分别进行时频变换,以生成各个信道的所有子载波的信道频率响应估计。
具体的,时频变换装置15可通过下面的公式(23)得到接收天线υ与所有发射天线所构成的NT个信道的CFR估计为
这里,M=N
fft/D
f/N
T,
表示列向量
的序号为μ·M到(μ+1)·M-1的元素组成的子向量,F
DFT表示N
fft点的DFT变换矩阵,表示为
公式(23)的含义是取信道冲击向量
中相应部分
然后后面添加N
fft-M个零,最后进行N
fft长度的离散傅立叶变换(DFT),即得到一个收发天线对之间的信道频率响应估计。
在另一个优选实施方式中(在此,将前述实施方式称为方案1,而本实施方式称为方案2,以示区别),第二计算装置13还可对步骤S12所获得的有效导频子载波处的信道频率响应进行滤波处理,以生成经过滤波的有效子载波的信道频率响应估计值。从而,在频时变换装置14将对所述经过滤波的有效导频子载波和虚拟导频子载波的信道频率响应估计值进行频时变换,以获得精度更高的时域信道冲击响应估计。从而最后获得精度更高的信道频率响应估计。具体实现方式如下所述:
具体的,第二计算装置14可采用稳健的MMSE算法将虚拟子载波中没有传输的导频点上的信道信息(CFR)进行估计和预测,同时将有效子载波位置上的最小二乘估计进行滤波,如下面公式(25)所示
(这个公式的物理意义公式前面已经说明,即对由有效导频子载波上的最小二乘估计构成的列向量左乘一个预测和滤波矩阵,将得到预测的虚拟子载波中导频位置上的信道频率响应和经滤波的有效导频子载波处的信道频率响应。)
这里,
对于单发OFDM来说,公式(26)变为
其中,
是有效导频和虚拟导频的子载波序号集合。RHp(μ)Hep(μ)表示信道的稳健的相关矩阵,其计算方法参见稳健的相关矩阵的计算部分。
具体的,频时变换装置14通过下列公式(21)来求增强的信道估计所得的信道冲击响应(CIR)
该公式的物理意义是对列向量
进行离散傅立叶反变换(IDFT),得到由要估计的多个信道的冲击响应组成的信道冲击向量
稳健的相关矩阵的计算
对比文献[6]给出了稳健的相关矩阵的计算方法,在此也引入本发明中作为参考。根据所述计算方法,为了得到信道的相关矩阵,不需要测量信道的实际的信道(这在一般条件下是不容易实现的),而只需要估计出信道的最大延迟。而信道的最大延迟一般比较容易估计,如根据引导序列与本地序列进行相关运算等方法。只要实际信道的最大延迟小于本方案中设定的最大延迟(K0),那么无论实际信道具有什么样的功率延迟谱(power delay profile,PDP),计算出的相关矩阵都可适用,所以称之为“稳健的(Robust)”。
信道的Nfft个子载波的稳健的相关矩阵可表示为
其中,F为归一化的Nfft点IDFT变换矩阵,如下
这里,K0表示用采样间隔归一化的信道最大延迟,即估计的信道最大延迟除以采样间隔。
那么,相关矩阵R
Hvp(μ)Hep(μ)就是矩阵R
H(μ)H(μ)的某子矩阵,该子矩阵可以通过抽取R
H(μ)H(μ)的
标识的行和
标识的列所构成。同理,可以得到相关矩阵R
Hep(μ)Hep(μ)和R
Hp(μ)Hep(μ)。
该稳健的相关矩阵的计算方法的步骤在上述参照图5和6描述的本发明内容中也进行了描述。
技术效果
图7中直观的给出了现有技术中的传统信道估计方法和根据本发明的信道估计方法所获得多个信道的信道频率响应估计值的计算机仿真结果。各个参数取值为:Nfft=1024,Nh=448,Df=4,NT=4,和Mp=112。图7中(a)部分的图是指理想的信道频率响应,(b)部分是传统信道估计方案所得到的信道频率响应估计,(c)部分是根据本发明的信道估计方案所得到的信道频率响应估计。由图7很明显的可以看出,传统信道估计方案造成了比较明显的能量泄漏,而根据本发明的信道估计方案基本上去除了能量泄漏,从而提高了信道估计的性能。
此外,我们还对根据本发明的信道估计方案进行了信道估计均方误差(MSE)的仿真,如图8和9所示。仿真条件是:NT=4,NR=4,Nfft=1024,2Mp+1=897,K0=27,基带采样频率为10MHz,仿真中所用信道的抽头参数如表1所示,并假设已获得理想的时间和频率同步。Df=4和8时得到的MSE性能分别如图8和9所示。在图8和9中,PSP/DFT-CE表示传统的信道估计方案,E-PSP/DFT-CE-A表示根据本发明中提到的方案1,E-PSP/DFT-CE-B表示本文中提出的方案2。由图中可以看出,当导频在频域插的较密时,如Df=4,方案2的性能要明显比方案1要好,而方案1又明显的去除了传统算法中的MSE平台。但当导频在频域插的较少时,如Df=8,方案2与方案1的性能已很接近,但由于去除了传统信道估计方案中的MSE平台,提高了信道估计的性能。
以上对本发明的具体实施例进行了描述。需要理解对是,本发明并不局限于上述特定对实施方式,本领域技术人员可以在所附权利要求的范围内做出各种变形或修改。