背景技术
在基于正交频分复用(OFDM)多载波技术的新一代无线通信系统中,信道估计影响到整个系统的性能。
在最小均方误差(Minimum Mean Square Error,MMSE)信道估计算法能够充分利用OFDM时域以及频域的相关性信息,获得优良的性能。
由于功率时延谱(PDP)与频域相关系数互为傅里叶变换对,如果能够准确的获取PDP就能够准确获取频域相关系数,从而可以根据部分频点上的信道冲激响应估计与该频点频域相关性较高的频点上的信道冲激响应,用于频域滤波和频域信道估计,可见,PDP是进行信道估计的重要工具。
PDP与频域相关系数的关系说明如下:
假设频域相关函数可表示为:
(公式1)
其中,L表示多径个数,τl表示第l径的传输时延,σl 2表示第l径的幅度方差,假设有用符号时间表示为Tuse,则有△f=1/Tuse。
假设多径时延为采样时间的整数倍(对于时延不是采样时间整数倍的情况,因为能量主要分散到其临近的整数倍采样时间抽头上,所以仍然可以采用该假设),即N表示有用符号时间总采样点数。
代入公式1可以得到:
(公式2)
其中,il表示第l径的抽头位置。
进一步地,若定义信道的功率时延谱向量p,满足:
(公式3)
其中,m=0,1,2,...,N-1,将公式3代入公式2中,则有:
(公式4)。
实际应用中,由于PDP无法直接获得,通过采用预设谱型的方式进行频域滤波和频域信道估计。
下面以矩形谱为例,基于预设谱型生成rf[k]的过程如下:
第一步,估计最大时延值τmax;
第二步,生成矩形功率时延谱向量p,满足:
(公式5);
第三步,根据功率时延谱向量p生成信道的频域相关系数,表示为:
(公式6)。
对于其他预设谱型仅在于第二步中生成的功率时延谱的公式不同,生成频域相关系数的过程相同,在此不做累述。
虽然根据预设谱型和最大时延确定PDP谱型,实现简单,但实际的PDP谱型包括信道冲击长度、包络形状、抽头间隔、相对功率值等几个方面的信息,并不是理想的预设谱型。预设谱型实际上仅能够满足实际PDP谱型的信道冲击长度信息,对于实际PDP谱型的其它几个方面的信息则直接假设为预设谱型包含的相应信息。
可见,直接采用预设谱型为实际PDP谱型,会使得计算的频域相关系数与实际差异较大。并且,由于受到噪声和干扰的影响,通过估计PDP谱型包含的各方面的信息得到实际PDP谱型的方式,实现难度较大,且准确性也不高。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本发明实施例中,如图1所示,进行PDP估计的详细方法流程如下:
步骤101:获取各导频的时域信道估计。
优选地,获取各天线、各端口、各导频的频域信道估计,将该频域信道估计补零为预定点数后进行离散傅里叶逆变换(IFFT),获得时域信道估计。
具体地,如果预定点数大于导频点数,将频域信道估计补零为预定点数后进行IFFT;如预定点数小于导频点数,将频域信道估计等间隔抽样为预定点数后进行IFFT。
优选地,该预定点数大于导频的个数且为2的整数次幂。
步骤102:根据各导频的时域信道估计获得时域信道估计的平均值,根据该时域信道估计的平均值以及预先确定的门限值确定有效径,并确定该有效径在OFDM的快速傅里叶变换(FFT)的点数中的位置编号。
具体地,时域信道估计的平均值为通过计算所述各天线、各端口、各导频的所述时域信道估计的平均值获得;或者通过计算所述各天线、各端口、各导频的部分所述时域信道估计的平均值获得。
优选地,计算部分导频的时域信道估计的平均值,具体为:
根据获取的各导频的时域信道估计确定功率最大径,计算该功率最大径所在的天线以及天线端口的各列导频的时域信道估计的平均值。
本发明实施例中,通过计算时域信道估计的平均值达到降低噪声影响的目的。
优选地,门限值包括功率门限值和时限门限值。
优选地,确定有效径的具体过程如下:
确定功率门限值和时限门限值;
判断是否存在时域信道估计的平均值大于功率门限值、且时延小于时延门限值的径;
若存在,确定时域信道估计的平均值大于功率门限值、且时延小于时延门限值的径,根据确定的径确定有效径;
若不存在,获取各时域信道估计的功率值中的最大值所对应的径,根据获取的径确定有效径。
具体地,将确定的径以及与确定的径的距离在预设范围内的径作为有效径;或者,将获取的径以及与获取的径的距离在预设范围内的径作为有效径。
优选地,确定功率门限值的具体过程如下:
根据测量的导频的噪声功率以及第一预设系数确定第一功率门限值;
计算各导频的时域信道估计的功率,根据各功率中的最大值以及第二预设系数确定第二功率门限值;
将第一功率门限值和第二功率门限值中的较大值确定为功率门限值。
实际应用中,噪声功率可以有多种获取方式,包括但不限于:将测量获得的各天线、各端口、各导频的噪声功率的平均值作为导频的噪声功率;或者,为测量获得的某一个导频的噪声功率;或者,将测量获得的一个以上的导频的噪声功率的平均值作为导频的噪声功率。
优选地,根据循环前缀(Cyclic Prefix,CP)个数、终端提前接收的点数以及进行OFDM的FFT的点数,确定时延门限值。
具体地,时延门限值的计算公式为:
其中,NCP为CP的个数,NTao为终端侧由定时提前量决定的提前接收点数,NFFT为进行OFDM的FFT的点数。
优选地,确定有效径在OFDM的FFT的点数中的位置编号的具体过程为:
假设获取时域信道估计时所采用的IFFT的点数为预定点数,确定有效径后获取该有效径在预定点数中的位置,根据该有效径在预定点数中的位置、预定点数以及进行OFDM的FFT的点数,确定该有效径在进行OFDM的FFT的点数中的位置编号。
具体地,有效径在进行OFDM的FFT的点数中的位置编号,计算公式如下:
其中,NFFT1表示获取时域信道估计时所采用的IFFT的点数,即预定点数,NFFT表示进行OFDM的FFT的点数,ndelay1表示有效径在预定点数中的位置,ndelay表示有效径在进行OFDM的FFT的点数中的位置编号。
步骤103:根据有效径以及有效径在OFDM的FFT点数中的位置编号,确定PDP。
其中,根据有效径以及有效径在OFDM的FFT点数中的位置编号,确定PDP有两种实现方式,具体如下:
第一种实现方式,具体为:
确定有效径的归一化功率;
计算在进行OFDM的FFT点数中的位置编号相同的各有效径的归一化功率之和,并作为包含有进行OFDM的FFT的点数的PDP中、与该位置编号对应的位置的功率。
其中,确定有效径的归一化功率,具体为:获取有效径的功率,并计算各有效径的功率之和;将有效径的功率与各有效径的功率之和的比值,作为该有效径的归一化功率。其中,有效径的功率为该有效径所对应的时域信道估计的功率。
第二种具体实现,具体为:
将PDP的FFT点的功率初始化为零;
确定包含各有效径在OFDM的FFT点数中的位置编号的位置编号集合,该位置编号集合中的元素唯一,并确定位置编号集合中包含的元素个数N;
将PDP各FFT点中、位置编号集合中的各元素所确定的位置的功率设置为N分之一。
优选地,确定PDP后,将该PDP进行快速傅里叶变换(FFT),获取频域相关系数。
较佳地,针对连续的多个子帧分别计算得到频域相关系数后,计算各子帧的频域相关系数的平均值,将该平均值作为最终的频域相关系数,以提高得到的频域相关系数的准确性;或者针对连续多个子帧分别计算得到PDP后,对各PDP进行合并归一化,将合并后的结果作为最终的PDP,以提高得到的PDP的准确性,进而提高根据PDP获得的频域相关系数的准确性。
以下通过两个具体实施例对本发明中PDP估计的过程进行完整、详细的说明。
第一具体实施例中,如图2所示,进行PDP估计的具体过程如下:
步骤1:根据各天线各端口各列导频的频域信道估计获取时域信道估计。
具体为:对各符号上的导频点的频域信道估计Hka,p补零后进行NFFT1点的IFFT变换,其中,NFFT1大于频域上的导频点数、且为2的整数次幂,用公式表示为:hka,p=IFFTn-IFFT(Hka,p),其中,n_IFFT=NFFT1。
步骤2:计算各天线各端口的各列导频对应的时域信道估计的平均h_mean,用公式表示为:
其中,k表示径的编号,且k={0,…,NFFT1-1},Nka表示天线数,Np表示天线端口数,Npilot表示导频列数。
步骤3:根据测量获得的导频的噪声功率确定第一功率门限值,以及根据步骤1获得的各时域信道估计中的功率最大值确定第二功率门限值,并选择第一功率门限值和第二功率门限值中的较大值作为功率门限值。
具体地,第一功率门限值表示为:G1=α·σ2,其中,σ2表示测量获得的导频的噪声功率,α为通过仿真预先确定的系数。
实际应用中,噪声功率可以有多种获取方式,包括但不限于:将测量获得的各天线、各端口、各导频的噪声功率的平均值作为导频的噪声功率;或者,为测量获得的某一个导频的噪声功率;或者,将测量获得的一个以上的导频的噪声功率的平均值作为导频的噪声功率。
具体地,第二功率门限值表示为:G2=β·|hmax|2,其中,|hmax|2为步骤1获得的各时域信道估计中的功率最大值,即功率最大的径的功率,β为通过仿真预先确定的系数。
步骤4:计算时延门限值,具体计算公式如下:
其中,NCP为CP的个数,NTao为终端侧由定时提前量决定的提前接收点数,NFFT为进行OFDM的FFT的点数。
步骤5:选择h_mean中满足功率大于功率门限值G且时延值在τ内的各径为有效径,并记录所有有效径在NFFT1样点内的位置ndelay1,也就是记录有效径的编号k,以及记录有效径的功率P1。
较佳地,为了减小功率泄露带来的损失,将有效径两侧在预定距离范围内(△内)的径同时作为有效径。
步骤6:如果不存在满足h_mean>G的径,则将步骤1获得的各时域信道估计中的功率最大值对应的时域信道估计hmax及其两侧距离在预定范围内(△内)的径作为有效径,记录NFFT1样点内有效径的位置ndelay1,也就是记录有效径的编号k,以及记录有效径的功率P1。
步骤7:换算为NFFT样点内的ndelay和归一化后的功率P。
具体地,有效径在OFDM的NFFT样点内的位置编号计算公式如下:
具体地,有效经的归一化后的功率为:该有效径的功率与各有效径的功率之和的比值,计算公式如下:
步骤8:获取时延功率谱PDP,对PDP功率谱各径进行归零再计算有效径功率。
具体为:对PDP功率谱各径功率初始化为零,具体公式如下:
PDP(i)=0,其中,i=1,2,…,NFFT;
根据有效径在OFDM的NFFT样点内的位置编号以及有效径的归一化后的功率计算PDP,具体公式如下:
步骤9:进行NFFT点的FFT变换,将PDP转到频域,获取频域相关系数。用公式表示为:rf=FFTn(PDP),n=NFFT。
第二具体实施例中,进行PDP估计时步骤1~步骤4与第一具体实施例中的步骤1~步骤4相同,此处不再重复,以下仅描述该具体实施例中的后续步骤。
步骤5:选择h_mean中满足功率大于功率门限值G且时延值在τ内的各径为有效径,并记录所有有效径在NFFT1样点内的位置ndelay1,也就是记录有效径的编号k。
较佳地,为了减小功率泄露带来的损失,将有效径两侧在预定距离范围内(△内)的径同时作为有效径。
步骤6:如果不存在满足h_mean>G的径,则将步骤1获得的各时域信道估计中的功率最大值对应的时域信道估计hmax及其两侧距离在预定范围内(△内)的径作为有效径,记录NFFT1样点内有效径的位置ndelay1,也就是记录有效径的编号k。
步骤7:换算为NFFT样点内的位置编号ndelay。
具体地,有效径在OFDM的NFFT样点内的位置编号计算公式如下:
步骤8:获取时延功率谱PDP。
具体地,对PDP功率谱各径功率初始化为零,具体公式如下:
PDP(i)=0,其中,i=1,2,…,NFFT;
确定包含各有效径在OFDM的FFT点数中的位置编号的位置编号集合I,该位置编号集合中的元素唯一,并确定位置编号集合中包含的元素个数N,该位置编号集合可表示为:I={i|i=ndelay},I中元素唯一,如果有多个ndelay等于i,则只记录一次,I中元素个数为N。
根据位置编号集合中包含的元素以及元素个数获取PDP,用公式表示为:
其中,,i∈I。
步骤9:进行NFFT点的FFT变换,把时延功率谱PDP转到频域,获取频域相关系数。用公式表示为:rf=FFTn(PDP),其中,n=NFFT。
基于同一发明构思,本发明实施例还提供了一种PDP估计装置,该装置的具体实施可参见上述方法部分的描述,重复之处不再赘述,如图3所示,该装置主要包括:
获取模块301,用于获取各导频的时域信道估计;
处理模块302,用于根据所述各导频的时域信道估计获得时域信道估计的平均值,根据所述时域信道估计的平均值以及预先确定的门限值确定有效径,并确定所述有效径在进行正交频分复用OFDM的快速傅里叶变换FFT的点数中的位置编号;
确定模块303,用于根据所述有效径以及所述位置编号,确定功率时延谱PDP。
优选地,所述门限值包括功率门限值和时限门限值;
所述处理模块具体用于:
确定所述功率门限值和所述时限门限值;
判断是否存在所述时域信道估计的平均值大于所述功率门限值、且时延小于所述时延门限值的径;
若存在,确定所述时域信道估计的平均值大于所述功率门限值、且时延小于所述时延门限值的径,根据确定的所述径确定有效径;
若不存在,获取各所述时域信道估计的功率值中的最大值所对应的径,根据获取的所述径确定有效径。
优选地,所述处理模块具体用于:
根据测量的导频的噪声功率以及第一预设系数确定第一功率门限值;
计算各导频的所述时域信道估计的功率,根据各功率中的最大值以及第二预设系数确定第二功率门限值;
将所述第一功率门限和所述第二功率门限值确定为所述功率门限值。
优选地,所述处理模块具体用于:
根据CP个数、终端提前接收的点数以及所述进行OFDM的FFT的点数,确定所述时延门限值。
优选地,所述处理模块具体用于:
将确定的所述径以及与确定的所述径的距离在预设范围内的径作为所述有效径;
或者,
将获取的所述径以及与获取的所述径的距离在预设范围内的径作为所述有效径。
优选地,所述获取模块具体用于:获取各天线、各端口、各导频的频域信道估计;将所述频域信道估计补零为预定点数后进行离散傅里叶逆变换IFFT,获得所述时域信道估计。
具体地,获取模块在预定点数大于导频点数时,将频域信道估计补零为预定点数后进行IFFT;在预定点数小于导频点数时,将频域信道估计等间隔抽样为预定点数后进行IFFT。
优选地,所述预定点数大于所述导频的个数且为2的整数次幂。
优选地,所述处理模块具体用于:
计算所述各天线、各端口、各导频的所述时域信道估计的平均值;或者,
计算所述各天线、各端口、各导频的部分所述时域信道估计的平均值。
优选地,处理模块具体用于:
获取有效径在预定点数中的位置,根据有效径在预定点数中的位置、预定点数以及进行OFDM的FFT的点数,确定该有效径在进行OFDM的FFT的点数中的位置编号。
其中,确定模块确定PDP有以下两种实现方式:
第一种,所述确定模块具体用于:
确定所述有效径的归一化功率;
计算所述位置编号相同的各所述有效径的归一化功率之和,作为包含有所述进行OFDM的FFT的点数的所述PDP中、与所述位置编号对应的位置的功率。
其中,所述确定模块获取所述有效径的功率,并计算各所述有效径的功率之和;将所述有效径的功率与各所述有效径的功率之和的比值,作为所述有效径的归一化功率。
第二种,所述确定模块具体用于:
将所述PDP各FFT点的功率初始化为零;
确定包含各所述有效径的所述位置编号的位置编号集合,所述位置编号集合中的元素唯一,并确定所述位置编号集合中包含的元素个数N;
将所述PDP各FFT点中、所述位置编号集合中的各元素所确定的位置的功率设置为N分之一。
基于上述技术方案,本发明实施例中,通过实际信道估计确定有效径以及有效径在进行OFDM的FFT的点数中的位置编号,即PDP谱的抽头间隔,根据有效径以及有效径在进行OFDM的FFT的点数中的位置编号得到PDP,使得能够结合PDP谱的抽头间隔确定的PDP谱型,提高了PDP估计的准确性,进而能够获得较为准确的频域相关系数,提高了吞吐量。并且,对于一些非理想的信道场景能够有效提高信道估计精度,提高信道估计算法的整体性能。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。