CN103825857A - 一种lte上行频域ls信道估计定点实现方法 - Google Patents

Abstract

本发明公开一种LTE上行频域LS信道估计定点实现方法，该方法包括：S1.对接收信号中的导频数据进行信道估计，获得初始信道响应I₀；S2.根据初始信道响应I₀，生成信道估计预处理表；S3.根据信道估计预处理表对初始信道响应I₀进行预处理，得到信道响应I。

Description

一种LTE上行频域LS信道估计定点实现方法

技术领域

本发明涉及移动通信领域，具体涉及一种LTE上行频域LS信道估计定点实现方法。

背景技术

LTE无线通信系统的上行链路中，由于传输信号受到信道干扰的影响，导致接收到的信号在相位和幅度上发生畸变，所述信道干扰主要是无线信道中的多径衰落和多普勒频移。因此，接收机需要通过信道均衡技术消除信道干扰。

传统的LTE上行频域信道估计定点实现方法，如图1所示，该方法包括：

对接收到的数据去除循环前缀(CP)后，经过快速傅里叶变换(FFT)变换到频域，从频域数据中提取导频的频域数据Y_RS；

将获得的Y_RS和用户的本地导频数据X_RS进行匹配，获得上行频域信道估计I，所述匹配方法为：

$I=\frac{Y_{\mathrm{RS}}}{X_{\mathrm{RS}}};$

在得到LTE上行频域信道估计I之后，对于接收数据中不是导频数据的用户数据Y_R，满足Y_R=I×X_T，其中，X_T为上行UE(UserEquipment)侧发送的数据，则

$X_T=\frac{Y_R}{I}.$

现有的最小二乘法(LS)信道估计方法在定点数字信号处理(DSP)平台上实现时，存在的问题是：一般来说，定点数字信号处理器内没有硬件除法器，大量的除法运算会严重影响接收机处理效率。

发明内容

本发明所要解决的技术问题是现有的信道估计方法中大量的除法运算影响到定点运算平台处理效率。

为此目的，本发明提出一种LTE上行频域LS信道估计定点实现方法，该方法包括：

S1.对接收信号中的导频数据进行信道估计，获得初始信道响应I₀；

S2.根据初始信道响应I₀，生成信道估计预处理表；

S3.根据信道估计预处理表对初始信道响应I₀进行预处理，得到信道响应I。

其中，所述步骤S1包括：

对接收到的数据去除循环前缀后，经过快速傅里叶变换将时域接收数据变换为频域接收数据，从频域接收数据中提取导频的频域数据Y_RS；

将获得的Y_RS和用户的本地导频数据X_RS进行匹配，获得信道响应I₀，所述匹配为：

$\begin{array}{c}{I}_0=\{\frac{Y_{\mathrm{RS}}\left(0\right)}{X_{\mathrm{RS}}\left(0\right)},\frac{Y_{\mathrm{RS}}\left(1\right)}{X_{\mathrm{RS}}\left(1\right)},...,\frac{Y_{\mathrm{RS}}\left(i\right)}{X_{\mathrm{RS}}\left(i\right)},...,\frac{Y_{\mathrm{RS}}(n-1)}{X_{\mathrm{RS}}(n-1)}\}\\ =\{\frac{Y_R\left(0\right)\×X_{\mathrm{RS}}{\left(0\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(0\right)\right|}^2},\frac{Y_R\left(1\right)\×X_{\mathrm{RS}}{\left(1\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(1\right)\right|}^2},...,\frac{Y_R\left(i\right)\×X_{\mathrm{RS}}{\left(i\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(i\right)\right|}^2},...,\frac{Y_R(n-1)\×X_{\mathrm{RS}}{(n-1)}^{*}}{{\left|X_{\mathrm{RS}}(n-1)\right|}^2}\};\end{array}$

其中，X_RS(i)^*与X_RS(i)互为共轭复数，|X_RS(i)|²为本地导频数据X_RS(i)的幅度的平方，所述|X_RS(i)|²=2^R，其中，i=1,2,3…n，n为导频数据的数量，R为正整数；

所述匹配等价于：I₀(i)=[Y_R(i)×X_RS(i)^*]>>R，其中“>>”为右移操作。

其中，在步骤S1和步骤S2之间，该方法还包括：

S11.根据计算出的初始信道响应I₀，获得信道功率的变化范围[P_L，P_H]；

S12.对所述信道功率的变化范围[P_L，P_H]进行修正，得到修正后的信道功率变化范围[P_min，P_max]，其中P_min=P_L-w×(P_H-P_L)，P_max=P_H+w×(P_H-P_L)，所述w根据LTE系统接收机的精确度确定，所述P_min>0；

S13.判断接收信号是否为第一个接收信号，若是，则执行步骤S2，否则，执行步骤S13；

S14.判断初始信道响应I₀中是否存在某个信道估计值所对应的信道功率在[P_min，P_max]外，若成立，执行S2，否则执行S3。

其中，所述步骤S2包括：

将[P_min，P_max]平均分为N份，每一等分点的功率值为P_j，j=1,2,3,…,N-1，P₀=P_min，其中，所述N的取值由定点实现平台所需精确度确定；

根据定点实现平台表示的最大正整数λ，通过计算P_j×M_j=λ，得到M_j，j=0,1,2,3,…,N-1，所述M_j，j=0,1,2,3,…,N-1为预处理表。

其中，所述步骤S3包括：

根据初始信道响应I₀中的信道估计值H_i，计算H_i×(H_i)^*，其中，所述H_i为I₀中的任意一个信道估计值，(H_i)^*与H_i互为共轭复数；

通过不等式P_j≤H_i×(H_i)^*，j=0,1,2,3,…,N-1，选择最大的P_j，得到P_j对应的M_j；

根据所述信道估计值H_i计算对应的H，得到I₀中每个元素对应的H，所述H的计算公式如下：

H=(H_i×M_j)>>R₁

其中，“>>”为右移操作，R₁为满足不等式

的最大整数，j=0,1,2,…,N-1；

所有H构成信道响应I。

较佳的，该方法进一步包括：

S4.根据得到信道响应I，计算上行UE侧发送的数据X_T。

其中，所述X_T中第k个符号X_T(k)的计算公式如下：

$X_T\left(k\right)=\frac{Y_R\left(k\right)}{I\left(k\right)}=\frac{Y_R\left(k\right)\×I{\left(k\right)}^{*}}{{\left|I\left(k\right)\right|}^2}$

其中，Y_R(k)为接收端接收到的数据，I(k)^*与I(k)互为共轭复数，k≥0，k为正整数。

其中，所述步骤S4包括：

通过对Y_R(k)×I(k)^*右移(R₂+1)位，得到X_T中第k个符号X_T(k)，其中，R₂为满足

的最大整数。

相比于现有技术，本发明提供的方法的有益效果是：

1.借助导频信号获得初始信道估计矩阵I₀，由于传输信道的影响，不同频点时信道估计H_i的模值的平方|H_i|²可能会有微小偏差，本发明提供的方法通过对|H_i|²进行补偿调整，保证每个频点上|H_i|²近似相等，有效地补偿定点运算带来的频率选择性衰落对信道估计I的影响，又避免进行除法运算，节省上行接收端处理时间。

2.本发明提供的方法在定点实现时避免进行除法运算，节省上行接收端处理时间，仅生成动态更新信道估计预处理表涉及少量除法运算，待预处理表稳定后能保持较长时间不更新，提高运算效率。

3.本发明提供的方法在定点实现时，能根据接收机实际字长及接收到空口数据的功率动态调整信道估计的补偿值，具有可操作性。当接收到数据的幅度发生变化时，能快速调整至合适的补偿值，具有自适应处理能力。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1示出了传统的LTE上行频域信道估计方法流程图；

图2示出了上行链路一个子帧的结构图；

图3示出了LTE系统上行子帧的时频二维结构图；

图4示出了一种LTE上行频域LS信道估计定点实现方法流程图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

对于LTE系统上行链路，传输信号采用SC-FDMA(Single-carrierFrequency-Division Multiple Access)技术和块状导频结构，导频RS(Reference Signals)符号和用户数据符号时分复用。根据LTE协议TS36.211第5.5.2节，对于Normal循环前缀，上行链路一个子帧结构如图2所示，其中导频符号与用户数据符号共有14个，记作0-13，其中第4个符号(Symbol)和第11个符号为导频符号，一个时隙(Slot)包含7个符号。

LTE系统上行子帧的时频二维结构，如图3所示，示意图假设一个子帧占用了12个子载波（Subcarrier），一个时隙的12个子载波为一个资源栅格(RB)。

假设导频发送信号的系统模型为：

Y_RS=X_RS×I+N

其中，Y_RS为接收到的参考信号矩阵，X_RS为本地的参考信号矩阵，I为信道响应，N为加性高斯白噪声(AWGN)矢量。

本发明实施例中采用最小二乘(Least Squares，LS)算法得到信道估计I的估计，具体过程如下：

令J=(Y_RS-X_RS×I_LS)^H×(Y_RS-X_RS×I_LS)，其中，(Y_RS-X_RS×I_LS)^H为(Y_RS-X_RS×I_LS)的共轭转置矩阵；

计算 $\frac{\∂J}{{\∂H}_{\mathrm{LS}}}=0;$

得到LS算法的信道估计I_LS， $I_{\mathrm{LS}}=\{I_{\mathrm{LS}}\left(0\right),I_{\mathrm{LS}}\left(1\right),...,I_{\mathrm{LS}}(n-1)\}$ $=\{\frac{Y_{\mathrm{RS}}\left(0\right)}{X_{\mathrm{RS}}\left(0\right)},\frac{Y_{\mathrm{RS}}\left(1\right)}{X_{\mathrm{RS}}\left(1\right)},......,\frac{Y_{\mathrm{RS}}(n-1)}{X_{\mathrm{RS}}(n-1)}\},$ 其中n为导频信号的数量。

对于时不变信道，同一频点不同符号间信道响应仅存在微小偏差，因而在定点实现时同一时隙内各符号可以使用同一组信道估计值。

随着通信系统大规模商业，尤其是LTE/LTE-A产业化不断推进，定点DSP凭借优秀的数字处理能力在通信设备中被广泛使用。在定点运算平台上，往往没有专门的除法运算器，导致除法计算速度远不及乘法运算。经实测，在TI TMS320C6474芯片上二者运算速度相差15倍以上。因而，当需要处理的数据量较大时，会影响系统的性能。所以工程实践中，在使用定点运算平台时需要尽量避免大规模的除法运算。还需要说明的是数字实现信号处理时，数字量化误差是时刻存在的，但是需要控制在不影响系统性能的范围内。

一种LTE上行频域LS信道估计定点实现方法，如图4所示，该方法包括：

S1.对接收信号中的导频数据进行信道估计，获得初始信道响应I₀；

S2.根据初始信道响应I₀，生成信道估计预处理表；

S3.根据信道估计预处理表对初始信道响应I₀进行预处理，得到信道响应I。

所述步骤S1包括：

对接收到的数据去除循环前缀后，经过快速傅里叶变换将时域接收数据变换为频域接收数据，从频域接收数据中提取导频的频域数据Y_RS；

将获得的Y_RS和用户的本地导频数据X_RS进行匹配，获得信道响应I₀，所述匹配为：

$\begin{array}{c}{I}_0=\{\frac{Y_{\mathrm{RS}}\left(0\right)}{X_{\mathrm{RS}}\left(0\right)},\frac{Y_{\mathrm{RS}}\left(1\right)}{X_{\mathrm{RS}}\left(1\right)},...,\frac{Y_{\mathrm{RS}}\left(i\right)}{X_{\mathrm{RS}}\left(i\right)},...,\frac{Y_{\mathrm{RS}}(n-1)}{X_{\mathrm{RS}}(n-1)}\}\\ =\{\frac{Y_R\left(0\right)\×X_{\mathrm{RS}}{\left(0\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(0\right)\right|}^2},\frac{Y_R\left(1\right)\×X_{\mathrm{RS}}{\left(1\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(1\right)\right|}^2},...,\frac{Y_R\left(i\right)\×X_{\mathrm{RS}}{\left(i\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(i\right)\right|}^2},...,\frac{Y_R(n-1)\×X_{\mathrm{RS}}{(n-1)}^{*}}{{\left|X_{\mathrm{RS}}(n-1)\right|}^2}\};\end{array}$

其中，X_RS(i)^*与X_RS(i)互为共轭复数，|X_RS(i)|²为本地导频数据X_RS(i)的幅度的平方，所述|X_RS(i)|²=2^R，其中，i=1,2,3…n，n为导频数据的数量，R为正整数；

所述匹配等价于：I₀(i)=[Y_R(i)×X_RS(i)^*]>>R，其中“>>”为右移操作。

在步骤S1和步骤S2之间，该方法还包括：

S11.根据计算出的初始信道响应I₀，获得信道功率的变化范围[P_L，P_H]；

S12.对所述信道功率的变化范围[P_L，P_H]进行修正，得到修正后的信道功率变化范围[P_min，P_max]，其中P_min=P_L-w×(P_H-P_L)，P_max=P_H+w×(P_H-P_L)，所述w根据LTE系统接收机的精确度确定，所述P_min>0；

S13.判断接收信号是否为第一个接收信号，若是，则执行步骤S2，否则，执行步骤S13；

S14.判断初始信道响应I₀中是否存在某个信道估计值所对应的信道功率在[P_min，P_max]外，若成立，执行S2，否则执行S3。

所述步骤S2包括：

将[P_min，P_max]平均分为N份，每一等分点的功率值为P_j，j=1,2,3,…,N-1，P₀=P_min，其中，所述N的取值由定点实现平台所需精确度确定；

根据定点实现平台表示的最大正整数λ，通过计算P_j×M_j=λ，得到M_j，j=0,1,2,3,…,N-1，所述M_j，j=0,1,2,3,…,N-1为预处理表。

所述步骤S3包括：

根据初始信道响应I₀中的信道估计值H_i，计算H_i×(H_i)^*，其中，所述H_i为I₀中的任意一个信道估计值，(H_i)^*与H_i互为共轭复数；

通过不等式P_j≤H_i×(H_i)^*，j=0,1,2,3,…,N-1，选择最大的P_j，得到P_j对应的M_j；

根据所述信道估计值H_i计算对应的H，得到I₀中每个元素对应的H，所述H的计算公式如下：

H=(H_i×M_j)>>R₁

其中，“>>”为右移操作，R₁为满足不等式

的最大整数，j=0,1,2,…,N-1；

所有H构成信道响应I。

较佳的，该方法进一步包括：

S4.根据得到信道响应I，计算上行UE（User Equipment）侧发送的数据X_T。

所述X_T中第k个符号X_T(k)的计算公式如下：

$X_T\left(k\right)=\frac{Y_R\left(k\right)}{I\left(k\right)}=\frac{Y_R\left(k\right)\×I{\left(k\right)}^{*}}{{\left|I\left(k\right)\right|}^2}$

其中，Y_R(k)为接收端接收到的数据，I(k)^*与I(k)互为共轭复数，k≥0，k为正整数。

所述步骤S4包括：

通过对Y_R(k)×I(k)^*右移(R₂+1)位，得到X_T中第k个符号X_T(k)，其中，R₂为满足

的最大整数。

实施例1：

本实施例中接收机运行后执行第一次信道估计运算或者信道估计矩阵中的某个信道估计所对应的功率在[P_min，P_max]外，需重新生成信道估计预处理表。所述接收机为带符号的定点实现平台，其参数为：运算字长为32位，定标为Q₀格式。则该定点实现平台中能表示的最大正整数λ=2³¹-1。

接收到的数据为1个资源栅格(RB)的数据，12位导频数据计算出的H_i分别为：

序号i	1	2	3	4	5	6
							Hi	229+284j	306+228j	353+126j	372+23j	377-71j	328-159j
序号i	7	8	9	10	11	12
							Hi	284-237j	227-297j	139-350j	50-385j	-56-378j	-158-332j

根据定点实现平台的精确度要求，w取值0.2，N取值32。

统计得到P_L=132865，P_H=150725，则P_min=129293，P_max=154297。将[P_min，P_max]均分为32份，P_j及其对应的M_j为：

j	1	2	3	4	5	6	7	8
									Pj	129293	130074	130856	131637	132419	133200	133981	134763
Mj	16609	16510	16411	16314	16217	16122	16028	15935
									j	9	10	11	12	13	14	15	16
Pj	135544	136325	137107	137888	138670	139451	140232	141014
									Mj	15843	15753	15663	15574	15486	15400	15314	15229
j	17	18	19	20	21	22	23	24
									Pj	141795	142576	143358	144139	144921	145702	146483	147265
Mj	15145	15062	14980	14899	14818	14739	14660	14582
									j	25	26	27	28	29	30	31	32
Pj	148046	148827	149609	150390	151172	151953	152734	153516
									Mj	14506	14429	14354	14279	14206	14133	14060	13989

满足

的R₁=13，满足的

的R₂=8。对于非导频位置的信号，接收到的有效数据为：

序号i	1	2	3	4	5	6
							数据Yi	277-117j	-163-423j	155+33j	76-222j	-249-259j	-170+188j
序号i	7	8	9	10	11	12
							数据Yi	-285-211j	288+10j	-71+59j	74-142j	165-303j	438+67j

对H_i进行(H_i×M_j)>>R₁处理后得到H，所有H构成信道估计矩阵I，对Y_R×I^*右移(R₂+1)位，恢复原发送序列，则原始发送的数据为：

序号i	1	2	3	4	5	6
							Xi	117-408i	-517-326i	215-29i	86-311i	-264-403i	-331+134i
序号i	7	8	9	10	11	12
							Xi	-116-479i	229+322i	-110-60i	199+73i	370+279i	-347+512i

实施例2：

在实施例1的基础上，本实施例中接收机执行第q次信道估计运算，q>1，且信道估计矩阵中的所有信道估计所对应的功率均在[P_min，P_max]范围内。

接收到的数据为1个资源栅格(RB)的数据，12位导频数据计算出的H_i分别为：

序号i	1	2	3	4	5	6
							Hi	229+284j	306+228j	353+126j	372+23j	377-71j	328-159j
序号i	7	8	9	10	11	12
							Hi	284-237j	227-297j	139-350j	50-385j	-56-378j	-158-332j

根据实施例1中得到的P_min=129293，P_max=154297，P_j及其对应的M_j为：

j	1	2	3	4	5	6	7	8
									Pj	129293	130074	130856	131637	132419	133200	133981	134763
Mj	16609	16510	16411	16314	16217	16122	16028	15935
									j	9	10	11	12	13	14	15	16
Pj	135544	136325	137107	137888	138670	139451	140232	141014
									Mj	15843	15753	15663	15574	15486	15400	15314	15229
j	17	18	19	20	21	22	23	24
									Pj	141795	142576	143358	144139	144921	145702	146483	147265
Mj	15145	15062	14980	14899	14818	14739	14660	14582
									j	25	26	27	28	29	30	31	32
Pj	148046	148827	149609	150390	151172	151953	152734	153516
									Mj	14506	14429	14354	14279	14206	14133	14060	13989

满足的R₁=13，满足的

的R₂=8。对于非导频位置的信号，接收到的有效数据为：

序号i	1	2	3	4	5	6
							数据Yi	277-117j	-163-423j	155+33j	76-222j	-249-259j	-170+188j
序号i	7	8	9	10	11	12
							数据Yi	-285-211j	288+10j	-71+59j	74-142j	165-303j	438+67j

对H_i进行(H_i×M_j)>>R₁处理后得到H，所有H构成信道估计矩阵I，对Y_R×I^*右移(R₂+1)位，恢复原发送序列，则原始发送的数据为：

序号i	1	2	3	4	5	6
							Xi	117-408i	-517-326i	215-29i	86-311i	-264-403i	-331+134i
序号i	7	8	9	10	11	12
							Xi	-116-479i	229+322i	-110-60i	199+73i	370+279i	-347+512i

虽然结合附图描述了本发明的实施方式，但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。

Claims (8)

1.一种LTE上行频域LS信道估计定点实现方法，其特征在于，该方法包括：

S1.对接收信号中的导频数据进行信道估计，获得初始信道响应I₀；

S2.根据初始信道响应I₀，生成信道估计预处理表；

S3.根据信道估计预处理表对初始信道响应I₀进行预处理，得到信道响应I。

2.根据权利要求1所述的方法，其特征在于，所述步骤S1包括：

对接收到的数据去除循环前缀后，经过快速傅里叶变换将时域接收数据变换为频域接收数据，从频域接收数据中提取导频的频域数据Y_RS；

将获得的Y_RS和用户的本地导频数据X_RS进行匹配，获得信道响应I₀，所述匹配为：

$\begin{array}{c}{I}_0=\{\frac{Y_{\mathrm{RS}}\left(0\right)}{X_{\mathrm{RS}}\left(0\right)},\frac{Y_{\mathrm{RS}}\left(1\right)}{X_{\mathrm{RS}}\left(1\right)},...,\frac{Y_{\mathrm{RS}}\left(i\right)}{X_{\mathrm{RS}}\left(i\right)},...,\frac{Y_{\mathrm{RS}}(n-1)}{X_{\mathrm{RS}}(n-1)}\}\\ =\{\frac{Y_R\left(0\right)\×X_{\mathrm{RS}}{\left(0\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(0\right)\right|}^2},\frac{Y_R\left(1\right)\×X_{\mathrm{RS}}{\left(1\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(1\right)\right|}^2},...,\frac{Y_R\left(i\right)\×X_{\mathrm{RS}}{\left(i\right)}^{*}}{{\left|X_{\mathrm{RS}}\left(i\right)\right|}^2},...,\frac{Y_R(n-1)\×X_{\mathrm{RS}}{(n-1)}^{*}}{{\left|X_{\mathrm{RS}}(n-1)\right|}^2}\};\end{array}$

其中，X_RS(i)^*与X_RS(i)互为共轭复数，|X_RS(i)|²为本地导频数据X_RS(i)的幅度的平方，所述|X_RS(i)|²=2^R，其中，i=1,2,3…n，n为导频数据的数量，R为正整数；

所述匹配等价于：I₀(i)=[Y_R(i)×X_RS(i)^*]>>R，其中“>>”为右移操作。

3.根据权利要求1所述的方法，其特征还在于，在步骤S1和步骤S2之间，该方法还包括：

S11.根据计算出的初始信道响应I₀，获得信道功率的变化范围[P_L，P_H]；

S12.对所述信道功率的变化范围[P_L，P_H]进行修正，得到修正后的信道功率变化范围[P_min，P_max]，其中P_min=P_L-w×(P_H-P_L)，P_max=P_H+w×(P_H-P_L)，所述w根据LTE系统接收机的精确度确定，所述P_min>0；

S13.判断接收信号是否为第一个接收信号，若是，则执行步骤S2，否则，执行步骤S13；

S14.判断初始信道响应I₀中是否存在某个信道估计值所对应的信道功率在[P_min，P_max]外，若成立，执行S2，否则执行S3。

4.根据权利要求3所述的方法，其特征在于，所述步骤S2包括：

将[P_min，P_max]平均分为N份，每一等分点的功率值为P_j，j=1,2,3,…,N-1，P₀=P_min，其中，所述N的取值由定点实现平台所需精确度确定；

根据定点实现平台表示的最大正整数λ，通过计算P_j×M_j=λ，得到M_j，j=0,1,2,3,…,N-1，所述M_j，j=0,1,2,3,…,N-1为预处理表。

5.根据权利要求4所述的方法，其特征在于，所述步骤S3包括：

根据初始信道响应I₀中的信道估计值H_i，计算H_i×(H_i)^*，其中，所述H_i为I₀中的任意一个信道估计值，(H_i)^*与H_i互为共轭复数；

通过不等式P_j≤H_i×(H_i)^*，j=0,1,2,3,…,N-1，选择最大的P_j，得到P_j对应的M_j；

根据所述信道估计值H_i计算对应的H，得到I₀中每个元素对应的H，所述H的计算公式如下：

H=(H_i×M_j)>>R₁

其中，“>>”为右移操作，R₁为满足不等式

的最大整数，j=0,1,2,…,N-1；

所有H构成信道响应I。

6.根据权利要求1所述的方法，其特征还在于，该方法进一步包括：

S4.根据得到信道响应I，计算上行UE侧发送的数据X_T。

7.根据权利要求6所述的方法，其特征在于，所述X_T中第k个符号X_T(k)的计算公式如下：

$X_T\left(k\right)=\frac{Y_R\left(k\right)}{I\left(k\right)}=\frac{Y_R\left(k\right)\×I{\left(k\right)}^{*}}{{\left|I\left(k\right)\right|}^2}$

其中，Y_R(k)为接收端接收到的数据，I(k)^*与I(k)互为共轭复数，k≥0，k为正整数。

8.根据权利要求7所述的方法，其特征在于，所述步骤S4包括：

通过对Y_R(k)×I(k)^*右移(R₂+1)位，得到X_T中第k个符号X_T(k)，其中，R₂为满足

的最大整数。

Images