一种低复杂度的信道估计方法与流程

1.本发明涉及无线通信技术领域，特别涉及一种低复杂度的信道估计方法。


背景技术：

2.ofdm信号在无线衰落信道中传输，其幅度会发生衰减，相位也会发生偏移。为了在接收端解调时能够准确恢复出原始数据信息，需要在接收端进行信道估计。
3.根据处理域的不同，信道估计有时域和频域之分，在接收端fft变换之后，提取目标用户的导频符号和本地导频共轭点乘后得到初始信道估计结果，时域信道估计方法需要将初始信道估计结果先进行ifft变换到时域加窗进行降噪，然后将降噪后的结果做fft转回频域，在频域完成均衡。而频域信道估计结果，对初始信道估计结果直接在频域进行滤波降噪。
4.对于ofdm系统，信道估计的一种简便高效的方法是频域信道估计方法，少一次ifft和fft，计算复杂度低，便于工程实现。对于频域信道估计方法，如果存在ta(timing advance，定时提前量)时，尤其是ta较大时，现有技术一是将初始信道估计结果直接进行固定系数的频域滤波，用固定系数的频域滤波方法可能导致信号被滤除一部分，导致降低了信道估计结果的精度，导致信号功率偏低，噪声功率偏高，影响均衡结果，进而影响通信距离。现有技术二是首先将初始信道估计结果先进行ta补偿，然后再送入频域滤波器进行降噪，最后将降噪的结果把ta再添加上，增加了一次对初始信道估计结果的ta补偿和一次对降噪结果添加ta的过程。
5.基于此现有的问题，本技术提供了一种低复杂度的信道估计方法，通过估计定时提前量ta，根据ta量进行滤波器系数修正，修正后的滤波器系数对导频的初始信道估计结果进行频域滤波，得到了信道估计的噪声功率值，此方法能够在不降低信道估计精度的同时，降低信道估计的计算量。


技术实现要素：

6.本发明的目的在于提供一种低复杂度的信道估计方法，通过估计定时提前量ta，根据ta量进行滤波器系数修正，修正后的滤波器系数对导频的初始信道估计结果进行频域滤波，得到了信道估计的噪声功率值，此方法能够在不降低信道估计精度的同时，降低信道估计的计算量。
7.本发明提供了一种低复杂度的信道估计方法，包括以下步骤：
8.采用快速傅里叶变换，将基站接收到的时域数据变换到频域，获得频域ofdm信号；
9.提取目标用户的频域ofdm信号，对目标用户的频域ofdm信号进行导频的初始信道估计；
10.根据导频的初始信道估计结果，获得同一ofdm符号上不同子载波之间的相位差，根据相位差计算得到定时提前量ta；
11.根据定时提前量ta计算得到修正系数，按照修正系数对滤波器系数进行修正；
12.根据修正后的滤波器系数，对导频的初始信道估计结果进行频域滤波，获得滤波后导频的信道估计结果；
13.获取滤波的噪声功率值，根据修正后的滤波器系数，进行噪声功率值的修正，得到噪声功率的估计值。
14.进一步地，所述采用快速傅里叶变换，将基站接收到的时域数据变换到频域，获得频域ofdm信号的步骤，包括：
15.基站接收时域数据；
16.对基站接收的时域数据执行去cp操作；
17.采用快速傅里叶变换，对去cp的时域数据进行fft变换到频域，其中fft点数记作n
fft
，取2的幂次方。
18.进一步地，还包括：
19.对初始信道估计结果h_ls(k)中前kf/2个数据和后kf/2个数据进行功率守恒操作，得到导频的初始信道估计结果
[0020][0021]
其中，kf为滤波器阶数。
[0022]
进一步地，所述提取目标用户的频域ofdm信号，对目标用户的频域ofdm信号进行导频的初始信道估计的步骤，包括：
[0023]
选定目标用户；
[0024]
设置目标用户的参数，所述目标参数包括：起始资源块rb和rb的数目；
[0025]
根据目标用户的参数，提取目标用户对应的频域ofdm信号；
[0026]
计算导频符号的初始信道估计结果h_ls(k)，初始信道估计结果表示为：
[0027]
h_ls(k)＝rxpilot(k).*conj(localpilot(k)),k＝0,1,2,...,k-1
[0028]
其中，rxpilot(k)表示接收到的导频频域ofdm信号，localpilot(k)表示本地导频频域ofdm信号，conj(x)表示对数据x取共轭。
[0029]
进一步地，所述根据导频的初始信道估计结果，获得同一ofdm符号上不同子载波之间的相位差，根据相位差计算得到定时提前量ta的步骤，包括：
[0030]
获取导频的初始信道估计结果；
[0031]
将初始信道估计结果中的信道分成若干正交子信道；
[0032]
将高速数据信号转换成并行的低速子数据流，调制到在每个子信道上进行传输，获取同一ofdm符号上不同子载波之间的相位差相位差表示为：
[0033][0034]
其中，m表示估计相位差时的子载波间隔数目；conj(x)表示对x执行共轭操作，angle(x)表示对数据x执行求角度的操作，单位为弧度；
[0035]
根据相位差φ求取定时提前量ta，其中ta表示为：
[0036]
[0037]
进一步地，所述根据定时提前量ta计算得到修正系数，按照修正系数对滤波器系数进行修正的步骤，包括：
[0038]
根据定时提前量ta计算修正系数：
[0039][0040]
其中，kf表示滤波器阶数，取偶数；kf表示修正系数索引，取值kf＝0,1,...,kf，n
fft
表示fft点数；
[0041]
根据修正系数对滤波器的原始系数进行修正，得出新的滤波器系数：
[0042]
coeff
new
(kf)＝coeff
init
(kf)
·
coeff
modified
(kf),kf＝0,1,...,kf[0043]
其中，coeff
init
(kf)为原始滤波器系数，coeff
new
(kf)为修正后的滤波器系数。
[0044]
进一步地，所述获取滤波的噪声功率值，根据修正后的滤波器系数，进行噪声功率值的修正，得到噪声功率的估计值的步骤，包括：
[0045]
根据滤波后导频的信道估计结果和导频的初始信道估计结果的差值，得到噪声序列，并计算噪声序列的噪声功率值；
[0046]
根据修正后的滤波器系数，计算滤波器系数的时域结果；
[0047]
根据滤波器系数的时域结果，获得噪声功率值的修正量；
[0048]
根据噪声功率值的修正量修正噪声功率值，得到噪声功率的估计值。
[0049]
与现有技术相比，本发明具有如下显著优点：
[0050]
本发明提供了一种低复杂度的信道估计方法，通过估计定时提前量ta，根据ta量进行滤波器系数修正，修正过程的计算量较小，不会造成信道估计过程中的计算负担，修正后的滤波器系数对导频的初始信道估计结果进行频域滤波，得到了信道估计的噪声功率值；本发明提供的一种低复杂度的信道估计方法，能够在不降低信道估计精度的同时，降低信道估计的计算量。
附图说明
[0051]
图1为本发明实施例提供的一种低复杂度的信道估计方法的估计流程图；
[0052]
图2为本发明实施例提供的频域滤波器系数转到时域的功率图；
[0053]
图3为本发明实施例提供的ls信道估计结果直接经过固定系数滤波器后转到时域的功率图；
[0054]
图4为本发明实施例提供的ls信道估计结果直接经过系数修正滤波器后转到时域的功率图。
具体实施方式
[0055]
下面结合本发明中的附图，对本发明实施例的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应属于本发明保护的范围。
[0056]
参照图1-图4，本发明提供了一种低复杂度的信道估计方法，包括以下步骤：
[0057]
s1：采用快速傅里叶变换，将基站接收到的时域数据变换到频域，获得频域ofdm信号。
[0058]
获得频域ofdm信号的步骤，包括：
[0059]
基站接收时域数据；
[0060]
对基站接收的时域数据执行去cp操作(cyclic prefix：循环移位)；
[0061]
采用快速傅里叶变换，对去cp的时域数据进行fft变换(fast fourier transform：快速傅里叶变换)到频域，其中fft点数记作n
fft
，取2的幂次方。
[0062]
s2：提取目标用户的频域ofdm信号，对目标用户的频域ofdm信号进行导频的初始信道估计。
[0063]
选定目标用户；
[0064]
设置目标用户的参数，所述目标参数包括：起始资源块rb(resource block)和rb的数目；
[0065]
根据目标用户的参数，提取目标用户对应的频域ofdm信号；
[0066]
计算初始信道估计结果h_ls(k)，h_ls(k)为行向量，初始信道估计结果表示为：
[0067]
h_ls(k)＝rxpilot(k).*conj(localpilot(k)),k＝0,1,2,...,k-1
[0068]
其中，rxpilot(k)表示接收到的导频频域ofdm信号，localpilot(k)表示本地导频频域ofdm信号，conj(x)表示对数据x取共轭。
[0069]
s3：根据导频的初始信道估计结果，获得同一ofdm符号上不同子载波之间的相位差，根据相位差计算得到定时提前量ta(timing advance：定时提前量)。得到定时提前量ta的步骤包括：
[0070]
获取导频的初始信道估计结果；
[0071]
将初始信道估计结果中的信道分成若干正交子信道；
[0072]
将高速数据信号转换成并行的低速子数据流，调制到在每个子信道上进行传输，获取同一ofdm符号上不同子载波之间的相位差相位差表示为：
[0073][0074]
其中，m表示估计相位差时的子载波间隔数目；conj(x)表示对x执行共轭操作，angle(x)表示对数据x执行求角度的操作，单位为弧度；
[0075]
根据相位差φ求取定时提前量ta，其中ta表示为：
[0076][0077]
s4：根据定时提前量ta计算得到修正系数，按照修正系数对滤波器系数进行修正，修正步骤包括：
[0078]
根据定时提前量ta计算修正系数：
[0079][0080]
其中，kf表示滤波器阶数，取偶数；kf表示修正系数索引，取值kf＝0,1,...,kf，n
fft
表示fft点数；
[0081]
根据修正系数对滤波器的原始系数进行修正，得出新的滤波器系数：
[0082]
coeff
new
(kf)＝coeff
init
(kf)
·
coeff
modified
(kf),kf＝0,1,...,kf[0083]
其中，coeff
init
(kf)为原始滤波器系数，coeff
new
(kf)为修正后的滤波器系数。
[0084]
s5：根据修正后的滤波器系数，对导频的初始信道估计结果进行频域滤波，获得滤波后导频的信道估计结果，具体滤波过程见如下公式：
[0085][0086]
s6：获取滤波的噪声功率值，根据修正后的滤波器系数，进行噪声功率值的修正，得到噪声功率的估计值，得到噪声功率的估计值的步骤包括：
[0087]
根据滤波后导频的信道估计结果和导频的初始信道估计结果的差值，得到噪声序列，并计算噪声序列的噪声功率值，其中，噪声序列为：
[0088]
noiseseq(k)＝h_ls(k)-h
filter
(k),k＝0,1,...,k-1
[0089]
噪声功率值为：
[0090][0091]
根据修正后的滤波器系数，计算滤波器系数的时域结果：
[0092]
coeff
t
(n)＝ifft(coeff,n
fft
),n＝0,1,...,n
fft-1
[0093]
其中ifft(x,n
fft
)表示将x执行n
fft
点的ifft操作；
[0094]
根据滤波器系数的时域结果，获得噪声功率值的修正量：
[0095][0096]
根据噪声功率值的修正量修正噪声功率值，得到噪声功率的估计值pn＝p
n,init
δpn。
[0097]
实施例1
[0098]
s1：采用快速傅里叶变换，将基站接收到的时域数据变换到频域，获得频域ofdm信号。
[0099]
s2：提取目标用户的频域ofdm信号，对目标用户的频域ofdm信号进行导频的初始信道估计，因为h_ls(k)开头和结束的kf/2个数据滤波时会出现功率不守恒，为了保证功率守恒，对初始信道估计结果h_ls(k)中前kf/2个数据和后kf/2个数据进行功率守恒操作，得到导频的初始信道估计结果
[0100][0101]
其中，kf为滤波器阶数；
[0102]
s3：根据导频的初始信道估计结果，获得同一ofdm符号上不同子载波之间的相位差，根据相位差计算得到定时提前量ta。
[0103]
s4：根据定时提前量ta计算得到修正系数，按照修正系数对滤波器系数进行修正。
[0104]
s5：根据修正后的滤波器系数，对导频的初始信道估计结果进行频域滤波，获得滤波后导频的信道估计结果。
[0105]
s6：获取滤波的噪声功率值，根据修正后的滤波器系数，进行噪声功率值的修正，得到噪声功率的估计值。
[0106]
实施例2
[0107]
以n
fft
＝512，目标用户的起始rb为10，rb数目为10为例，此时，k＝120。
[0108]
步骤一：将接收到的时域数据，去cp(cyclic prefix：循环移位)，进行fft变换(fast fourier transform：快速傅里叶变换)到频域，根据起始rb(resource block：资源块)和rb数目提取目标用户的频域ofdm信号，计算导频符号的初始信道估计结果h_ls(k)，h_ls(k)为行向量；图1为该滤波器系数转到时域的功率图；
[0109]
h_ls(k)＝rxpilot(k).*conj(localpilot(k)),k＝0,1,2,...,k-1
[0110]
其中，rxpilot(k)表示接收到的导频频域ofdm信号，localpilot(k)表示本地导频频域ofdm信号，conj(x)表示对数据x执行取共轭的操作。
[0111]
步骤二：利用导频的初始信道估计结果估计ta；
[0112]
phasediff＝h_ls(k' m).*conj(h_ls(k')),k'＝0,1,2,...,n-m-1
[0113][0114]
其中，m表示估计相位差时的子载波间隔数目，建议取值为6；angle(x)表示对数据x执行求角度的操作，单位为弧度；
[0115]
步骤三：利用ta值对滤波器系数进行更新；
[0116][0117]
其中，coeff
init
(kf)为原始滤波器系数，取值为[0.0025 0.0148 0.0414 0.0802 0.1223 0.1550 0.1674 0.1550 0.1223 0.0802 0.0414 0.0148 0.0025]；
[0118]
其中，kf表示滤波器阶数，为偶数，取值为12；coeff
new
(kf)为更新后的滤波器系数。
[0119]
步骤四：用更新后的滤波器系数对初始信道估计结果进行频域滤波，得到滤波后的信道估计结果h
filter
(k)：
[0120][0121]
步骤五：用滤波前和滤波后的信道估计结果估计噪声功率：
[0122]
首先，计算噪声序列：
[0123]
noiseseq(k)＝h_ls(k)-h
filter
(k),k＝0,1,...,k-1
[0124]
其次，计算初始噪声功率：
[0125][0126]
然后，计算频域滤波器系数的时域结果：
[0127]
coeff
t
(n)＝ifft(coeff,n
fft
),n＝0,1,...,n
fft-1
[0128]
接着，噪声功率修正量：
[0129][0130]
最后，得到修正后的噪声功率：
[0131]
pn＝p
n,init
δpn。
[0132]
其中，参照图2，ta＝30ts对应的衰减是-3.0088db，
[0133]
参照图3，当ta＝30ts时，ls(最小二乘法)信道估计结果直接经过固定系数的滤波器后转到时域的功率图，其中ta＝30ts时，滤波前和滤波后的功率相差2.9880db和图2中的-3.0088db基本一致。蓝色线表示ls信道估计结果在滤波前的时域功率图；红色线表示ls信道估计结果在滤波后的时域功率图。
[0134]
参照图4，当ta＝30ts时，ls信道估计结果直接经过系数修正的滤波器后转到时域的功率图。蓝色线表示ls信道估计结果在滤波前的时域功率图；红色线表示ls信道估计结果在滤波后的时域功率图。
[0135]
以上公开的仅为本发明的几个具体实施例，但是，本发明实施例并非局限于此，任何本领域的技术人员能思之的变化都应落入本发明的保护范围。