CN101626254A - 一种高速联合检测的方法 - Google Patents

Abstract

本发明公开了一种高速联合检测的方法，包括：确定每个用户设备UE的消除了频偏的信道冲击响应；根据每个UE的消除了频偏的信道冲击响应检测发送符号，得到每个UE发送的消除了频偏的符号序列估计值。本发明通过对每个UE的频率偏移进行估计，并在此基础上计算出消除了频偏后的信道响应估计后，基于所述消除了频偏后的信道响应估计进行联合检测，并对联合检测的结果进行纠正频偏，使多普勒频移在本发明的联合检测中进行了补偿。因此，本发明能克服多普勒频移对联合检测结果的影响，具有很好的联合检测性能，不仅适用于UE处于静止或慢速移动的情况，还适用于UE处于高速移动的情况。

Description

一种高速联合检测的方法

技术领域

本发明涉及时分同步码分多址接入(TD-SCDMA)系统中的信号联合检测技术，特别是涉及一种适用于用户设备(UE)处于高速移动状态时的信号联合检测(Joint Detection，JD)方法。

背景技术

在TD-SCDMA系统中，每个时隙的接收信号包括864个码片。其中，开始的352个码片是第一个数据域的接收信号，接着的144个码片是训练序列(MIDMABLE)域的接收信号，再接着的352个码片是第二个数据域的接收信号，最后16个码片是GP(保护周期)。当信号接收端接收到用户设备(UE)发送的信号后需要根据接收信号进行联合检测，从而从接收信号中将每个UE在第一个和第二个数据域发送的符号序列求解出来。

通常，接收端在对接收信号进行联合检测时，首先需要根据MIDMABLE域的接收信号估计出每个UE的信道冲击响应。这里，每个用户在发送信号时，根据系统事先为其分配的一个训练序列偏移(MIDAMBLE SHIFT)值确定它发送的MIDAMBLE域信号。MIDAMBLE域接收信号是每个UE发送的MIDAMBLE域信号经过信道以后的信号的迭加。由于每个UE的MIDAMBLE SHIFT值在接收端是已知的，因此，每个UE发送的MIDMABLE域信号在接收端也是已知的。接收端利用其接收的MIDAMBLE域接收信号和已知的每个UE发送的MIDAMBLE域信号，通过信道估计可以得到每个UE的信道冲击响应。然后，接收端利用每个UE的信道冲击响应构造出用于求解每个UE发送的符号序列的系统矩阵。最后，接收端根据所述系统矩阵求解相关的失量方程即可将每个UE发送的符号序列求解出来。

在N根接收天线的TD-SCDMA系统中，假设有K个MIDMABLESHIFT，每个UE分配一个MIDMABLE SHIFT，第k个UE采用的MIAMBLESHIFT值为k，1≤k≤K，第k个UE的信道冲击响应为，

$h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},...,h_{k,l}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}]---\left(1\right)$

其中，W是信道冲击响应的窗长，n为天线编号，1≤n≤N，h_k，l ⁽ⁿ⁾为第k个UE在第n根天线的第l个路径上的信道冲击响应，l为路径编号，1≤l≤W。

由所有K个UE在第n根天线上的信道冲击响应构成的第n根天线的信道冲击响应矢量为：

$h^{\left(n\right)}={[h_1^{\left(n\right)},...,h_k^{\left(n\right)},...,h_K^{\left(n\right)}]}^T---\left(2\right)$

在上述TD-SCDMA系统中，假设e_n表示第n根天线的接收信号，e_n为包括864个元素的列矢量，其中，每个元素对应一个码片的接收信号；e_n，mid表示第n根天线的MIDAMBLE域接收信号，它是由e_n中间的MIDAMBLE域的144个元素(码片)构成的列矢量。e_n，i表示第i个数据域接收信号，i＝1，2。e_n，1是由e_n的前(352+W-1)个元素构成的列矢量；e_n，2是由e_n的最后352+16个元素中的前(352+W-1)个元素构成的列矢量。由于e_n，mid的最开始的W-1个元素受到第一个数据域的干扰，W的最大值可能是16，因此，在进行信道冲击响应估计的时候，只选取e_n，mid最后的128个元素，构成列矢量e′_n，mid。e′_n，mid可以由如下方程表示：

e′_n，mid＝Gh⁽ⁿ⁾+v⁽ⁿ⁾         (3)

其中，v⁽ⁿ⁾为第n根天线上的MIDMABLE域接收信号的噪声失量，G是维数为128×KW矩阵，G矩阵由该时隙采用的基本MIDAMBLE码构成，G矩阵可以由下式描述，

G＝[G⁽¹⁾，G⁽²⁾，...，G^(K)]        (4)

这里，G^(k)具有如下结构

$G^{\left(k\right)}=\left[\begin{array}{ccccc}{\underset{\‾}{m}}_{W+(K-k)W}& {\underset{\‾}{m}}_{W-1+(K-k)W}& ...& {\underset{\‾}{m}}_{2+(K-k)W}& {\underset{\‾}{m}}_{1+(K-k)W}\\ {\underset{\‾}{m}}_{W+1+(K-k)W}& {\underset{\‾}{m}}_{W+(K-k)W}& ...& {\underset{\‾}{m}}_{3+(K-k)W}& {\underset{\‾}{m}}_{2+(K-k)W}\\ {\underset{\‾}{m}}_{W+2+(K-k)W}& {\underset{\‾}{m}}_{W+1+(K-k)W}& ...& {\underset{\‾}{m}}_{4+(K-k)W}& {\underset{\‾}{m}}_{3+(K-k)W}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ {\underset{\‾}{m}}_{W+127+(K-k)W}& {\underset{\‾}{m}}_{W+126+(K-k)W}& ...& {\underset{\‾}{m}}_{129+(K-k)W}& {\underset{\‾}{m}}_{128+(K-k)W}\end{array}\right]$

其中，m＝[m ₁，m ₂，m ₃，…，m ₁₂₈]是小区的基本MIDAMBLE码经过QPSK调制而成的复值符号序列，m _i＝m _i-P，i＝P+1，P+2，…，L_m+(K-1)W，P＝128，L_m＝144。

具体地，现有的联合检测方法包括：

步骤一、估计每个UE在每根天线上的信道冲击响应。

根据公式(3)可以推出

h⁽ⁿ⁾＝G^-1e′_n，mid               (5)

由公式(4)可以计算G^-1，然后，根据G^-1和e′_n，mid求解公式(5)即可得到第n根天线的信道冲击响应矢量h⁽ⁿ⁾，根据公式(2)进而可以确定第k个UE在第n根天线上的信道冲击响应h_k ⁽ⁿ⁾。

步骤二、根据步骤一中得到的h_k ⁽ⁿ⁾，对接收信号的数据域部分信号e_n，i进行联合检测，求解出每个UE发送的符号序列。

这里，需要首先将MIDAMBLE域的接收信号对第一个数据域e_n，1的最后W-1个元素的干扰消除掉，同时将MIDAMBLE域的接收信号对第二个数据域e_n，2的最开始的W-1的元素的干扰消除掉，经过干扰消除的第一、第二数据域信号仍旧用e_n，1和e_n，2表示，具体的干扰消除方法参见相关文献，此处不再赘述。具体地，根据h⁽ⁿ⁾求解出每个UE发送的符号序列的方法包括：

步骤a、构造系统矩阵A。

系统矩阵A的结构为，

$A=\left[\begin{array}{c}{A}^{\left(1\right)}\\ A^{\left(2\right)}\\ .\\ .\\ .\\ A^{\left(N\right)}\end{array}\right]$

其中，A⁽ⁿ⁾具有如下结构：

其中，构成系统矩阵A⁽ⁿ⁾的b_k ⁽ⁿ⁾包括W+SF_k-1个元素，SF_k为第k个UE的扩频因子，b_k ⁽ⁿ⁾具有如下结构：

$b_k^{\left(n\right)}=[b_{k,1}^{\left(n\right)},b_{k,2}^{\left(n\right)},...b_{k,W+{\mathrm{SF}}_k-1}^{\left(n\right)}],$

由公式 $b_k^{\left(n\right)}=h_k^{\left(n\right)}\⊗u_k$ 可以计算b_k ⁽ⁿ⁾，这里，

表示卷积， $u_k=[u_{k,1},u_{k,2},...,u_{k,{\mathrm{SF}}_k}],$ u_k为第k个UE采用的扩频码，它是该用户分配的OVSF信道码和小区扰码之积。

A⁽ⁿ⁾的前N₁列对应第1个UE；仅接着的N₂列对应第2个UE；从第

列到第

列对应第k个UE；最后第N_K列对应第K个UE。前N₁列中，第一列的前SF₁+W-1行由矢量b₁ ⁽ⁿ⁾的元素构成，其余元素都为0；第二列是由第一列向下平移SF₁行得到。第m列由第m-1列元素向下平移SF₁行得到。从第

列到第

列，第

列的前SF_k+W-1是由b_k ⁽ⁿ⁾的元素构成，其余元素都为0；第

列由其前一列向下平移SF_k行得到；第列由其前一列向下平移SF_k行得到，这里m∈{2，...，N_k}，N_k表示第k个UE发送的符号总数。

步骤b、根据系统矩阵A估计每个UE发送的符号序列，得到每个UE发送的符号序列估计值。

假设第k个UE在第i(i＝1，2)个数据域发送的符号序列为： $d_{k,i}=[d_{k,i,l},d_{k,i,2},...,d_{k,i,N_k}],$ 其中，d_k，i，m表示发送的第m个符号。所有K个UE发送的符号序列为d_i＝[d_1，i，d_2，i，…，d_K，i]，i＝1，2。可以建立如下方程：

e_i＝Ad_i+q_i    (6)

其中，e_i表示所有N根天线上的第i个数据域的接收信号，i＝1，2，e_i具有如下结构，

$e_i={[e_{1,i}^T,e_{2,i}^T,...,e_{N,i}^T]}^T$

q_i表示所有N根天线上的第i个数据域的接收信号中的噪声矢量，q_i具有如下结构，

$q_i=\left[\begin{array}{c}{q}_i^{\left(1\right)}\\ q_i^{\left(2\right)}\\ .\\ .\\ .\\ q_i^{\left(N\right)}\end{array}\right]$

其中，q_i ⁽ⁿ⁾是第n根天线上的第i个数据域的接收信号中的噪声矢量。基于MMSE算法求解方程(6)得到下述公式，

${\hat{d}}_{\mathrm{MMSE}-i}={(A^H{R}_n^{-1}A+R_d^{-1})}^{-1}{A}^H{R}_n^{-1}{e}_i---\left(7\right)$

其中，

表示第i个数据域的检测符号估计失量，其结构为：

${\hat{d}}_{\mathrm{MMSE}-i}=[{\hat{d}}_{1,i},...,{\hat{d}}_{k,i},...,{\hat{d}}_{K,i}],$ 其中，

是d_k，i的估计。

R_n表示第i个数据域的噪声协方差矩阵，通常，第一个数据域的噪声协方差矩阵和第二个数据域的噪声协方差矩阵近似相同，因此这里不再区别标记，统一用R_n表示，R_n由 $R_n=E\left[q_i{q}_i^T\right]$ 计算得到；R_d表示在第i个数据域所有K个UE发送的符号序列的协方差矩阵， $R_d=E\left[d_i{d}_i^T\right]=I.$

利用公式(7)，根据系统矩阵A、R_n、R_d和e_i可以得到每个数据流域的发送符号序列，从而可以得到每个UE在每个数据域的发送符号序列。当每根天线的噪声不相关时，R_n是对角矩阵，对角矩阵的元素由每根天线的平均噪声功率组成。具体地求R_n的方法参见有关文献，这里不再赘述。

当UE高速移动时，它的多普勒频移(Doppler frequency deviation)非常大，这将对现有的联合检测方法中计算的信道冲击响应的准确性有很大影响。因此，在UE高速移动的情况下，如果采用现有的方法进行联合检测，则联合检测的性能将会非常差，并且随着UE移动的速度越大，多普勒频移越严重，联合检测的性能也越差。

由此可见，现有的联合检测方法不适用于UE处于高速移动的情况。

发明内容

有鉴于此，本发明的主要目的在于提供一种联合检测的方法，该方法不仅适用于UE处于静止或慢速移动的情况，还适用于UE处于高速移动的情况。

为了达到上述目的，本发明提出的技术方案为：

一种高速联合检测的方法，包括以下步骤：

a、确定每个用户设备UE的消除了频偏的信道冲击响应；

b、根据每个UE的消除了频偏的信道冲击响应检测发送符号，得到每个UE发送的消除了频偏的符号序列估计值。

较佳地，所述步骤a包括：

a1、对每根天线上的接收信号进行信道冲击响应的估计，得到每根天线的信道冲击响应 $h^{\left(n\right)}={[h_1^{\left(n\right)},...,h_k^{\left(n\right)},...,h_K^{\left(n\right)}]}^T,$ 其中， $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},...,h_{k,l}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}]$ 为第k个

UE在第n根天线上的信道冲击响应，n为天线编号，k为UE编号，1≤k≤K，

l为路径编号，1≤l≤W；

a2、对每根天线的信道冲击响应h⁽ⁿ⁾进行频偏补偿，得到每个UE的消除了频偏的信道冲击响应h′⁽ⁿ⁾。

较佳地，所述步骤a2包括：

根据每根天线的信道冲击响应h⁽ⁿ⁾，确定每个UE在每条路径上的频率偏差Δf_k，l；

根据每个UE在每条路径上的频率偏差Δf_k，l，对G矩阵进行纠正频偏得到矩阵G′＝[G′⁽¹⁾，…，G′^(k)，...，G′^(K)]，其中，

$G^{\left(k\right)}=\left[\begin{array}{ccccc}{e}^{{\mathrm{ja}}_1}{\underset{\‾}{m}}_{W+(K-k)W}& e^{{\mathrm{ja}}_2}{\underset{\‾}{m}}_{W-1+(K-k)W}& ...& e^{{\mathrm{ja}}_{W-1}}{\underset{\‾}{m}}_{2+(K-k)W}& e^{{\mathrm{ja}}_W}{\underset{\‾}{m}}_{1+(K-k)W}\\ e^{j(a_1+{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+1+(K-k)W}& e^{j(a_2+{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+(K-k)W}& ...& e^{j(a_{W-1}+{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{3+(K-k)W}& e^{j(a_W+{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{2+(K-k)W}\\ e^{j(a_1+2{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+2+(K-k)W}& e^{j(a_2+2{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+1+(K-k)W}& ...& e^{j(a_{W-1}+2{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{4+(K-k)W}& e^{j(a_W+{2\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{3+(K-k)W}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ e^{j(a_1+127{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+127+(K-k)W}& e^{j(a_2+127{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+126+(K-k)W}& ...& e^{j(a_{W-1}+127{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{129+(K-k)W}& e^{j(a_W+127{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{128+(K-k)W}\end{array}\right]$

m＝[m ₁，m ₂，m ₃，…，m ₁₂₈]是由小区的基本MIDAMBLE码经过QPSK调制而成的复值符号序列，Δ_l＝2π·Δf_k，l·T_c，1≤l≤W， $a_l=e^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_{k,l}\·(352+16)\·T_c}$ 为第l条路径上MIDAMBLE域的第17个码片的初始相位，T_c为每个码片的持续时间；

根据G′和每根天线的信道冲击响应h⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

较佳地，所述步骤a2包括：

根据每根天线的信道冲击响应h⁽ⁿ⁾，确定每个UE在每条路径上的频率偏差Δf_k，l，其中，k为UE编号，l为路径编号；

根据每个UE在每条路径上的频率偏差，对G矩阵进行纠正频偏得到矩阵G′＝[G′⁽¹⁾，G′⁽²⁾，...，G′^K]，其中，

$G^{\left(k\right)}=\left[\begin{array}{ccccc}{e}^{{\mathrm{ja}}_1}{\underset{\‾}{m}}_{W+(K-k)W}& e^{{\mathrm{ja}}_2}{\underset{\‾}{m}}_{W-1+(K-k)W}& ...& e^{{\mathrm{ja}}_{W-1}}{\underset{\‾}{m}}_{2+(K-k)W}& e^{{\mathrm{ja}}_W}{\underset{\‾}{m}}_{1+(K-k)W}\\ e^{j(a_1+{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+1+(K-k)W}& e^{j(a_2+{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+(K-k)W}& ...& e^{j(a_{W-1}+{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{3+(K-k)W}& e^{j(a_W+{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{2+(K-k)W}\\ e^{j(a_1+2{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+2+(K-k)W}& e^{j(a_2+2{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+1+(K-k)W}& ...& e^{j(a_{W-1}+2{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{4+(K-k)W}& e^{j(a_W+{2\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{3+(K-k)W}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ e^{j(a_1+127{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+127+(K-k)W}& e^{j(a_2+127{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+126+(K-k)W}& ...& e^{j(a_{W-1}+127{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{129+(K-k)W}& e^{j(a_W+127{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{128+(K-k)W}\end{array}\right]$

其中，m＝[m ₁，m ₂，m ₃，…，m ₁₂₈]是小区的基本MIDAMBLE码经过QPSK调制而成的复值符号序列，Δ_l＝2π·Δf_k，l·T_c， $a_l=e^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_{k,l}\·(352+16)T_c}$ 为第l条路径上MIDAMBLE域的第17个码片的初始相位，T_c为每个码片持续的时间；

从每个UE的路径中选取M条强径l_k，1…l_k，M；

在h_k ⁽ⁿ⁾中仅保留第k个UE的M条强径所对应的信道冲击响应值，得到新的矢量h′_k ⁽ⁿ⁾；根据h′_k ⁽ⁿ⁾得到 $h^{\′\left(n\right)}={[h_1^{\′\left(n\right)},...,h_k^{\′\left(n\right)},...,h_K^{\′\left(n\right)}]}^T;$

在G′^(k)中仅保留第k个UE的M条强径的路径编号所对应的M列元素得到矩阵G″^(k)；根据G″^(k)得到G″＝[G″⁽¹⁾，…，G″^(k)，…，G″^(K)]；

根据G″和h′⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

较佳地，所述步骤a2中根据G′和每根天线的信道冲击响应h⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应包括：

根据G′计算(G′)^H；

根据R_h，n＝E(v⁽ⁿ⁾(v⁽ⁿ⁾)^H)＝σ²I计算R_h，n ^-1，其中，v⁽ⁿ⁾为第n根天线的训练序列MIDMABLE域接收信号的噪声失量；

计算R_h＝E(h⁽ⁿ⁾(h⁽ⁿ⁾)^H)；

根据(G′)^H、R_h，n ^-1、R_h和第n根天线的MIDAMBLE域接收信号e′_n，mid，基于MMSE算法计算 $h_{\mathrm{MMSE}}^{\left(n\right)}={({\left(G^{\′}\right)}^H{R}_{h,n}^{-1}{G}^{\′}+R_h^{-1})}^{-1}{\left(G^{\′}\right)}^H{R}_{h,n}^{-1}{e}_{n,\mathrm{mid}}^{\′},$ 得到每根天线消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

较佳地，所述步骤a2中根据G″和h′⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应包括：

根据G″计算(G″)^H；

根据R_h，n＝E(v⁽ⁿ⁾(v⁽ⁿ⁾)^H)＝σ²I计算R_h，n ^-1，其中，v⁽ⁿ⁾为第n根天线的MIDMABLE域接收信号的噪声失量；

计算 $R_{h^{\′}}=E\left(h^{\′\left(n\right)}{\left(h^{\′\left(n\right)}\right)}^H\right)=\left[\begin{array}{cccc}{p}_{1,l_{\mathrm{1,1}}}^{\left(n\right)}& 0& ...& 0\\ 0& p_{1,l_{1,2}}^{\left(n\right)}& ...& 0\\ .& .& & .\\ .& .& ...& .\\ .& .& & .\\ 0& 0& ...& p_{K{,l}_{K,M}}^{\left(n\right)}\end{array}\right],$ 其中， $p_{k.l}^{\left(n\right)}={\left|\right|h_{k,l}^{\left(n\right)}\left|\right|}^2$ 为第k个UE的第l条路径的功率值，l∈{l_k，1，…，l_k，M}；

根据(G″)^H、R_h，n ^-1、R_h′和第n根天线的MIDAMBLE域接收信号e′_n，mid，基于MMSE算法计算 $h_{\mathrm{MMSE}}^{\′\left(n\right)}={({\left(G^{\′\′}\right)}^H{R}_{h,n}^{-1}{G}^{\′\′}+R_{h^{\′}}^{-1})}^{-1}{\left(G^{\′\′}\right)}^H{R}_{h,n}^{-1}{e}_{n,\mathrm{mid}}^{\′},$ 得到h′_MMSE ⁽ⁿ⁾，h′_MMSE ⁽ⁿ⁾由每个UE的M条强径的消除了频偏的信道冲击响应值构成；用h′_MMSE ⁽ⁿ⁾更新h_(n)得到h_MMSE ⁽ⁿ⁾，所述更新h⁽ⁿ⁾的方法为：用h′_MMSE ⁽ⁿ⁾中每个UE的M条强径的消除了频偏的信道冲击响应值替代h_k ⁽ⁿ⁾中对应路径的信道冲击响应值，得到h″_k ⁽ⁿ⁾；根据h″_k ⁽ⁿ⁾得到每根天线消除了频偏的信道冲击响应 $h_{\mathrm{MMSE}}^{\left(n\right)}={[h_1^{\′\′\left(n\right)},...,h_k^{\′\′\left(n\right)},...,h_K^{\′\′\left(n\right)}]}^T.$

较佳地，所述步骤b包括：

b1x、根据每个UE的消除了频偏的信道冲击响应，构造系统矩阵A；

b2x、根据系统矩阵A，得到每个UE发送的符号序列估计值

其中，k为UE编号，i＝1，2为数据域编号，m∈{1，...，N_k}为符号编号，N_k为第k个UE发送的符号总数；

b3x、利用公式 $d_{k,i,m}^{\′}={\hat{d}}_{k,i,m}{e}^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_k\·\left((m-0.5){\mathrm{SF}}_k\right)\·T_c}$ 对

进行纠正频偏，得到每个UE发送的消除了频偏的符号序列估计值d′_k，i，m；

其中，SF_k为第k个UE的扩频因子，Δf_k为第k个UE的路径频率偏差平均值，Tc为每个码片的持续时间。

较佳地，所述步骤b3x中 $\mathrm{\Δ}{f}_k=\frac{1}{W}\underset{i=1}{\overset{W}{\mathrm{\Σ}}}\mathrm{\Δ}{f}_{k,l},$ 其中，W为信道冲击响应的窗长，Δf_k，l为第k个UE的第l个路径的频率偏差，Tc为每个码片的持续时间。

较佳地，所述步骤b3x中 $\mathrm{\Δ}{f}_k=\frac{1}{M}\underset{i=1}{\overset{M}{\mathrm{\Σ}}}\mathrm{\Δ}{f}_{k,l},$ M为每个UE的强径数量，Δf_k，l为第k个UE的路径编号为l的强径的频率偏差，Tc为每个码片的持续时间。

较佳地，所述步骤b包括：

b1y、根据每个UE的消除了频偏的信道冲击响应 $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},...,h_{k,l}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}],$ 确定b_k，m ⁽ⁿ⁾(i，l)，

其中，b_k，m ⁽ⁿ⁾(i，l)是由W+SF_k-1个元素构成的列矢量，b_k，m ⁽ⁿ⁾(i，l)的前l-1个元素为0，后W-l个元素是0，n为天线编号，k为UE编号，m为符号编号，m∈{1，...，N_k}，N_k表示第k个发送的符号总数，i为数据域编号，i＝1，2，l为路径编号，1≤l≤W，SF_k为第k个UE的扩频因子，W为信道冲击响应的窗长，b_k，m ⁽ⁿ⁾(1，l)的具体结构为：

$b_{k,m}^{\left(n\right)}\left(\mathrm{1,1}\right)=\left[\begin{array}{c}{u}_{k,1}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·((m-1)S{F}_k+1)T_c}\\ u_{k,2}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·((m-1)S{F}_k+2)T_c}\\ .\\ .\\ .\\ u_{k,S{F}_k}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·\left(\mathrm{mS}{F}_k\right)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

$b_{k,m}^{\left(n\right)}\left(\mathrm{1,2}\right)=\left[\begin{array}{c}0\\ u_{k,1}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·((m-1)S{F}_k+2)T_c}\\ u_{k,2}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·((m-1){\mathrm{SF}}_k+3)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·(\mathrm{mS}{F}_k+1)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

……

$b_{k,m}^{\left(n\right)}(1,l)=\left[\begin{array}{c}0\\ .\\ .\\ .\\ 0\\ u_{k,1}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,l}\·((m-1){\mathrm{SF}}_k+l)T_c}\\ u_{k,2}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,l}\·((m-1){\mathrm{SF}}_k+l+1)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,l}\·(\mathrm{mS}{F}_k+l-1)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

……

$b_{k,m}^{\left(n\right)}(1,W)=\left[\begin{array}{c}0\\ .\\ .\\ .\\ 0\\ u_{k,1}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}\·((m-1){\mathrm{SF}}_k+W)T_c}\\ u_{k,2}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}\·((m-1){\mathrm{SF}}_k+W+1)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}({\mathrm{mSF}}_k+W-1)T_c}\end{array}\right],$

其中， $u_k=[u_{k,1},u_{k,2},...,u_{k,{\mathrm{SF}}_k}],$ u_k为第k个UE的扩频码，Δf_k，l为第k个UE的第l条路径的频率偏差；

所述b_k，m ⁽ⁿ⁾(2，l)根据公式 $b_{k,m}^{\left(n\right)}(2,l)=b_{k,m}^{\left(n\right)}(1,l)\·e^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_{k,l}\·(352+144)\·T_c}$ 确定；

b2y、按照 $b_{k,m}^{\left(n\right)}\left(i\right)=\underset{l=1}{\overset{W}{\mathrm{\Σ}}}{b}_{k,m}^{\left(n\right)}(i,l),$ 确定b_k，m ⁽ⁿ⁾(i)；

b3y、利用b_k，m ⁽ⁿ⁾(i)确定A⁽ⁿ⁾′(i)，并按照 $A^{\′}\left(i\right)=\left[\begin{array}{c}{A}^{\left(1\right)\′}\left(i\right)\\ .\\ .\\ .\\ A^{\left(n\right)\′}\left(i\right)\\ .\\ .\\ .\\ A^{\left(N\right)\′}\left(i\right)\end{array}\right]$ 确定包含了路径频偏的新的系统矩阵A′(i)；其中，所述利用b_k，m ⁽ⁿ⁾(i)确定A⁽ⁿ⁾′(i)的方法为：

将A⁽ⁿ⁾′(i)的前N₁列对应第1个UE，紧接着的N₂列对应第2个UE，从第

列到第

列对应第k个UE，1≤k≤K，最后第N_K列对应第K个UE；

其中，在第k个UE对应的N_k列中，将列矢量b_k，n ⁽ⁿ⁾(i)的SF_k+W-1个元素依次放入第n列的从第(n-1)SF_k+1行开始的W+SF_k-1行，将第n列其余行的元素均设为0，n为所述第k个UE对应的N_k列中的第n列，1≤n≤N_k；

b4y、根据新的系统矩阵A′(i)，得到每个UE发送的消除了频偏的符号序列估计值d′_k，i，m。

综上所述，本发明提出的高速联合检测的方法通过对每个UE的频率偏移进行估计，并在此基础上计算出消除了频偏后的信道响应估计后，基于所述消除了频偏后的信道响应估计进行联合检测，并对联合检测的结果进行纠正频偏，使多普勒频移在本发明的联合检测中进行了补偿，因此，本发明能克服多普勒频移对联合检测结果的影响，具有很好的联合检测性能。由此可见，本发明不仅适用于UE处于静止或慢速移动的情况，还适用于UE处于高速移动的情况。

附图说明

图1为本发明方法的流程图。

图2为本发明方法实施例一的流程图。

图3为图2中步骤202的一种实施流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图及具体实施例对本发明作进一步地详细描述。

本发明的主要思想是，通过基于消除了频偏后的信道响应估计进行联合检测，并对联合检测的结果进行纠正频偏，使本发明能够克服多普勒频移对联合检测性能的影响。

图1为本发明的高速联合检测方法流程图。如图1所示，本发明主要包括：

步骤101、确定每个UE的消除了频偏的信道冲击响应。

本步骤中，首先，按照现有的信道估计方法对每根天线上的接收信号进行信道冲击响应的估计，然后，对每根天线的信道冲击响应进行频偏补偿，得到每个UE的消除了频偏的信道冲击响应。

步骤102、根据每个UE的消除了频偏的信道冲击响应检测发送符号，得到每个UE发送的消除了频偏的符号序列估计值。

本步骤，通过基于步骤101中得到的每个UE的消除了频偏的信道冲击响应对每个UE的发送符号进行检测，减少频移对联合检测性能的影响。

上述为本发明的高速联合检测方法的总体描述，下面通过本发明的较佳实施例对本发明进行详细说明。

图2为本发明实施例一的流程图。如图2所示，实施例一用于实现本发明的步骤101和步骤102，具体包括：

步骤201、对每根天线上的接收信号进行信道冲击响应的估计，得到每根天线的信道冲击响应。

本步骤中，可以按照现有的信道估计方法对每根天线上的接收信号进行信道冲击响应的估计，得到每根天线的信道冲击响应为：

$h^{\left(n\right)}={[h_1^{\left(n\right)},...,h_2^{\left(n\right)},...,h_K^{\left(n\right)}]}^T,$

其中， $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},h_{k,2}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}]$ 为第k个UE在第n根天线上的信道冲击响应，n为天线编号，k为UE编号，1≤k≤K，l为路径编号，1≤l≤W；

具体的信道估计方法为现有技术，此处不再赘述。

步骤202、对每根天线的信道冲击响应h⁽ⁿ⁾进行频偏补偿，得到每个UE的消除了频偏的信道冲击响应。

在实际应用中，本步骤的一种实现方法是：先根据h⁽ⁿ⁾确定每个UE在每条路径上的频率偏差Δf_k，l；然后再根据Δf_k，l对G矩阵进行纠正频偏得到矩阵G′，最后根据G′和h⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应。具体的操作流程如图3所示，包括

步骤2021、根据h⁽ⁿ⁾确定每个UE在每条路径上的频率偏差Δf_k，l。

本步骤中，可以按照现有方法得到Δf_k，l，此处不再赘述。

步骤2022、根据Δf_k，l对G矩阵进行纠正频偏得到矩阵G′。

本步骤中，通过对现有联合检测方法中G矩阵的每个元素进行纠正频偏得到G′＝[G′⁽¹⁾，G′⁽²⁾，...，G′^K]，其中，

$G^{\left(k\right)}=\left[\begin{array}{ccccc}{e}^{{\mathrm{ja}}_1}{\underset{\‾}{m}}_{W+(K-k)W}& e^{{\mathrm{ja}}_2}{\underset{\‾}{m}}_{W-1+(K-k)W}& ...& e^{{\mathrm{ja}}_{W-1}}{\underset{\‾}{m}}_{2+(K-k)W}& e^{{\mathrm{ja}}_W}{\underset{\‾}{m}}_{1+(K-k)W}\\ e^{j(a_1+{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+1+(K-k)W}& e^{j(a_2+{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+(K-k)W}& ...& e^{j(a_{W-1}+{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{3+(K-k)W}& e^{j(a_W+{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{2+(K-k)W}\\ e^{j(a_1+2{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+2+(K-k)W}& e^{j(a_2+2{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+1+(K-k)W}& ...& e^{j(a_{W-1}+2{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{4+(K-k)W}& e^{j(a_W+{2\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{3+(K-k)W}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ e^{j(a_1+127{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+127+(K-k)W}& e^{j(a_2+127{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+126+(K-k)W}& ...& e^{j(a_{W-1}+127{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{129+(K-k)W}& e^{j(a_W+127{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{128+(K-k)W}\end{array}\right]$

这里，G′^(k)中的第l列与第k个UE的第l条路径相对应，m＝[m ₁，m ₂，m ₃，…，m ₁₂₈]是小区的基本MIDAMBLE码经过QPSK调制而成的复值符号序列，Δ_l＝2π·Δf_k，l·T_c， $a_l=e^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_{k,l}\·(352+16)\·T_c}$ 为第l条路径上MIDAMBLE域的第17个码片的初始相位，Δf_k，l为第k个UE的第l条路径的频率偏差，T_c为每个码片的持续时间。

步骤2023、根据G′和h⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应。

本步骤中，较佳地，采用MMSE算法确定每个UE的消除了频偏的信道冲击响应，具体操作流程包括：

根据G′计算(G′)^H；

根据R_h，n＝E(v⁽ⁿ⁾(v⁽ⁿ⁾)^H)计算R_h，n ^-1，其中，v⁽ⁿ⁾为第n根天线的训练序列MIDMABLE域接收信号的噪声失量；R_h，n＝E(v⁽ⁿ⁾(v⁽ⁿ⁾)^H)的计算方法参见现有文献。

$R_h=E\left(h^{\left(n\right)}{\left(h^{\left(n\right)}\right)}^H\right)=\left[\begin{array}{cccc}{P}_{\mathrm{1,1}}& 0& ...& 0\\ 0& p_{\mathrm{1,2}}& ...& 0\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ 0& 0& ...& P_{K,W}\end{array}\right]$ 为对角矩阵。

其中，h⁽ⁿ⁾为步骤201中得到的第n根天线的信道冲击响应矢量， $p_{k,l}^{\left(n\right)}={\left|\right|h_{k,l}^{\left(n\right)}\left|\right|}^2,$ h_k，l ⁽ⁿ⁾为第k个UE的第n根天线的第l条路径的信道冲击响应；

根据(G′)^H、R_h，n ^-1、R_h和第n根天线的MIDAMBLE域接收信号e′_n，mid计算 $h_{\mathrm{MMSE}}^{\left(n\right)}={({\left(G^{\′}\right)}^H{R}_{h,n}^{-1}{G}^{\′}+R_h^{-1})}^{-1}{\left(G^{\′}\right)}^H{R}_{h,n}^{-1}{e}_{n,\mathrm{mid}}^{\′}$ 得到每根天线消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

需要说明的是步骤2023也可以利用公式h⁽ⁿ⁾＝(G′)^-1e′_n，mid求解，但是相比上述方法，该方法的性能差，此处不再赘述。

另外，在实际应用中，为减少算法的复杂度，步骤202中还可以只对h⁽ⁿ⁾中一部分元素进行上述运算，得到其中一部分元素的消除了频偏的信道冲击响应。较佳的具体方法如下：

x1、从每个UE的路径中选取M条强径l_k，1…l_k，M；

这里，根据每个UE的每条路径的功率 $p_{k,l}^{\left(n\right)}={\left|\right|h_{k,l}^{\left(n\right)}\left|\right|}^2,$ 从每个UE的W条路径中选出功率最大的M条路径，其中，k为UE编号，l为路径编号；所述M条强径的编号为l_k，m，m＝1，2，…，M。在实际应用中，M可以直接预先设定，也可以通过预先设定的功率门限值确定。

x2、在h_k ⁽ⁿ⁾中仅保留第k个UE的M条强径的信道冲击响应值，得到新的矢量h′_K ⁽ⁿ⁾；根据h′_k ⁽ⁿ⁾得到 $h^{\′\left(n\right)}={[h_1^{\′\left(n\right)},...,h_k^{\′\left(n\right)},...,h_K^{\′\left(n\right)}]}^T;$

具体的， $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},h_{k,2}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}]$ 的第l个元素是第k个UE的第l条路径的信道冲击响应值；将 $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},h_{k,2}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}]$ 中与第l_k，m条被选择的路径所对应的第l_k，m个元素保存，其余元素删除，得到一个包含M个元素的列矢量h′_k ⁽ⁿ⁾。

x3、在G′^(k)矩阵中仅保留第k个UE的M条强径的路径编号所对应的M列元素得到矩阵G″^(k)；根据G″^(k)得到G″＝[G″⁽¹⁾，…，G″^(k)，…，G″^(K)]；

具体的，将G^(k)中第l_k，m列保存，m＝1，2，…，M，其余列全部删除，得到的矩阵G″^(k)。

x4、根据G″和h′⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

较佳地，这里采用MMSE算法确定每个UE的消除了频偏的信道冲击响应，具体操作流程包括：

根据G″计算(G″)^H；

根据R_h，n＝E(v⁽ⁿ⁾(v⁽ⁿ⁾)^H)＝σ²I计算R_h，n ^-1，其中，v⁽ⁿ⁾为第n根天线的MIDMABLE域接收信号的噪声失量；

计算 $R_{h^{\′}}=E\left(h^{\′\left(n\right)}{\left(h^{\′\left(n\right)}\right)}^H\right)=\left[\begin{array}{cccc}{p}_{1,l_{\mathrm{1,1}}}^{\left(n\right)}& 0& ...& 0\\ 0& p_{1,l_{1,2}}^{\left(n\right)}& ...& 0\\ .& .& & .\\ .& .& ...& .\\ .& .& & .\\ 0& 0& ...& p_{K{,l}_{K,M}}^{\left(n\right)}\end{array}\right],$ 其中，

$p_{k,l}^{\left(n\right)}={\left|\right|h_{k,l}^{\left(n\right)}\left|\right|}^2$ 为第k个UE的第l条路径的功率值，l∈{l_k，1，…，l_k，M}；

根据(G″)^H、R_h，n ^-1、R_h′和第n根天线的MIDAMBLE域接收信号e′_n，mid，基于MMSE算法计算 $h_{\mathrm{MMSE}}^{\′\left(n\right)}={({\left(G^{\′\′}\right)}^H{R}_{h,n}^{-1}{G}^{\′\′}+R_{h^{\′}}^{-1})}^{-1}{\left(G^{\′\′}\right)}^H{R}_{h,n}^{-1}{e}_{n,\mathrm{mid}}^{\′},$ 得到h′_MMSE ⁽ⁿ⁾。

h′_MMSE ⁽ⁿ⁾由每个UE的M条强径的消除了频偏的信道冲击响应值构成；用h′_MMSE ⁽ⁿ⁾更新h⁽ⁿ⁾，所述更新h⁽ⁿ⁾的方法为：用h′_MMSE ⁽ⁿ⁾中每个UE的M条强径的消除了频偏的信道冲击响应值h′_k ⁽ⁿ⁾替代h_k ⁽ⁿ⁾中对应路径的信道冲击响应值，得到h″_k ⁽ⁿ⁾；例如：假设第k个UE的M条强径的编号为：l_k，m，m＝1，2，…，M，从h′_MMSE ⁽ⁿ⁾中h′_k ⁽ⁿ⁾可以得到第l_k，m条路径的消除了频偏的信道冲击响应值，m＝1，2，…，M；用h′_MMSE ⁽ⁿ⁾中h′_k ⁽ⁿ⁾的第l_k，m条路径的信道冲击响应值替代 $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},h_{k,2}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}]$ 中第l_k，m条路径的信道冲击响应值。根据h″_k ⁽ⁿ⁾即可得到每根天线消除了频偏的信道冲击响应 $h_{\mathrm{MMSE}}^{\left(n\right)}={[h_1^{\′\′\left(n\right)},...,h_k^{\′\′\left(n\right)},...,h_K^{\′\′\left(n\right)}]}^T.$

这里，由于G″和h′⁽ⁿ⁾仅为G′和h⁽ⁿ⁾中的部分元素，因此，减少了计算h′_MMSE ⁽ⁿ⁾时算法复杂度。

步骤203、根据每个UE的消除了频偏的信道冲击响应，构造系统矩阵A。

本步骤中，构造系统矩阵A的方法与现有技术相同，此处不再赘述。

步骤204、根据系统矩阵A，得到每个UE发送的符号序列估计值

其中，k为UE编号，i为数据域编号，m∈{1，...，N_k}为符号编号，N_k为第k个UE发送的符号总数。

本步骤的实现方法为现有技术，此处不再赘述。

步骤205、对所述每个UE发送的符号序列估计值进行纠正频偏，得到每个UE发送的消除了频偏的符号序列估计值d′_k，i，m。

本步骤中，利用公式 $d_{k,i,m}^{\′}={\hat{d}}_{k,i,m}{e}^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_k\·\left((m-0.5)S{F}_k\right){\·T}_c}$ 对步骤204中得到

进行纠正频偏，其中，d′_k，i，m为第k个UE在第i个数据域发送的第m个符号的消除了频偏的估计值，SF_k为第k个UE的扩频因子，Δf_k为第k个UE的路径频率偏差平均值， $\mathrm{\Δ}{f}_k=\frac{1}{W}\underset{l=1}{\overset{W}{\mathrm{\Σ}}}\mathrm{\Δ}{f}_{k,l},$ W为信道冲击响应的窗长，Δf_k，l为第k个UE的第l条路径的频率偏差，Tc为一个码片持续的时间。

这里需要说明的是，为提高本发明的算法性能，还可以仅根据每个UE的M条强径确定该UE的路径频率偏差平均值Δf_k，即 $\mathrm{\Δ}{f}_k=\frac{1}{M}\underset{l=1}{\overset{M}{\mathrm{\Σ}}}\mathrm{\Δ}{f}_{k,l},$ 其中，Δf_k，l为第k个UE的编号为l的强径的频率偏差。

上述方案中，步骤201、202实现了本发明的步骤101；步骤203、204、205实现了本发明的步骤102。

在实际应用中，为提高本发明的算法性能，步骤102还可以采用新的联合检测算法直接得到每个UE在数据域i(i＝1，2)发送的符号序列的消除了频偏的估计值。新的联合检测算法表现在：采用了新的系统矩阵A′，新的系统矩阵A′考虑了每个UE的每条路径的频偏。每个数据域的系统矩阵A⁽ⁿ⁾′(i)并不相同。具体利用新的系统矩阵A′实现步骤102的方法如下：

b1y、根据每个UE的消除了频偏的信道冲击响应 $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},...,h_{k,l}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}],$ 确定b_k，m ⁽ⁿ⁾(i，l)，

其中，b_k，m ⁽ⁿ⁾(i，l)是由W+SF_k-1个元素构成的列矢量，b_k，m ⁽ⁿ⁾(i，l)的前l-1个元素为0，后W-l个元素是0，n为天线编号，k为UE编号，m为符号编号，m∈{1，...，N_k}，N_k表示第k个发送的符号总数，i为数据域编号，i＝1，2，l为路径编号，1≤l≤W，SF_k为第k个UE的扩频因子，W为信道冲击响应的窗长，b_k，m ⁽ⁿ⁾(1，l)的具体结构为：

$b_{k,m}^{\left(n\right)}\left(\mathrm{1,1}\right)=\left[\begin{array}{c}{u}_{k,1}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·((m-1)S{F}_k+1)T_c}\\ u_{k,2}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·((m-1)S{F}_k+2)T_c}\\ .\\ .\\ .\\ u_{k,S{F}_k}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·\left(\mathrm{mS}{F}_k\right)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

$b_{k,m}^{\left(n\right)}\left(\mathrm{1,2}\right)=\left[\begin{array}{c}0\\ u_{k,1}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·((m-1)S{F}_k+2)T_c}\\ u_{k,2}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·((m-1){\mathrm{SF}}_k+3)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·(\mathrm{mS}{F}_k+1)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

……

$b_{k,m}^{\left(n\right)}(1,l)=\left[\begin{array}{c}0\\ .\\ .\\ .\\ 0\\ u_{k,1}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,l}\·((m-1){\mathrm{SF}}_k+l)T_c}\\ u_{k,2}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,l}\·((m-1){\mathrm{SF}}_k+l+1)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,l}\·(\mathrm{mS}{F}_k+l-1)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

……

$b_{k,m}^{\left(n\right)}(1,W)=\left[\begin{array}{c}0\\ .\\ .\\ .\\ 0\\ u_{k,1}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}\·((m-1){\mathrm{SF}}_k+W)T_c}\\ u_{k,2}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}\·((m-1){\mathrm{SF}}_k+W+1)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}({\mathrm{mSF}}_k+W-1)T_c}\end{array}\right]$

其中， $u_k=[u_{k,1},u_{k,2},...,u_{k,{\mathrm{SF}}_k}],$ u_k为第k个UE的扩频码，是该UE分配的OVSF信道码和小区扰码之积，Δf_k，l为第k个UE的第l条路径的频率偏差；

所述b_k，m ⁽ⁿ⁾(2，l)根据公式 $b_{k,m}^{\left(n\right)}(2,l)=b_{k,m}^{\left(n\right)}(1,l)\·e^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_{k,l}\·(352+144)\·T_c}$ 确定；

按照 $b_{k,m}^{\left(n\right)}\left(i\right)=\underset{l=1}{\overset{W}{\mathrm{\Σ}}}{b}_{k,m}^{\left(n\right)}(i,l),$ 确定b_k，m ⁽ⁿ⁾(i)；

利用b_k，m ⁽ⁿ⁾(i)确定A⁽ⁿ⁾′(i)，并按照 $A^{\′}\left(i\right)=\left[\begin{array}{c}{A}^{\left(1\right)\′}\left(i\right)\\ .\\ .\\ .\\ A^{\left(n\right)\′}\left(i\right)\\ .\\ .\\ .\\ A^{\left(N\right)\′}\left(i\right)\end{array}\right]$ 确定包含了路径频偏的新的

系统矩阵A′(i)；其中，所述利用b_k，m ⁽ⁿ⁾(i)确定A⁽ⁿ⁾′(i)的方法为：

将A⁽ⁿ⁾′(i)的前N₁列对应第1个UE，紧接着的N₂列对应第2个UE，从第

列到第

列对应第k个UE，1≤k≤K，最后第N_k列对应第K个UE；

其中，在第k个UE对应的N_k列中，将列矢量b_k，n ⁽ⁿ⁾(i)的SF_k+W-1个元素依次放入第n列的从第(n-1)SF_k+1行开始的W+SF_k-1行，将第n列其余行的元素均设为0，n为所述第k个UE对应的N_k列中的第n列，1≤n≤N_k；

这里，具体利用b_k，m ⁽ⁿ⁾(i)构建A⁽ⁿ⁾′(i)的流程为：

A⁽ⁿ⁾′(i)的前N₁列中，第一列的前SF₁+W-1行由矢量b_1，1 ⁽ⁿ⁾(i)的元素构成，其余元素都为0；第二列是由矢量b_1，2 ⁽ⁿ⁾(i)构成，b_1，2 ⁽ⁿ⁾(i)放入第二列的起始位置相对b_1，1 ⁽ⁿ⁾(i)放入第一列的起始位置(第一列的第一行)向下平移SF₁行得到，即：b_1，2 ⁽ⁿ⁾(i)的元素从第二列的第SF_k+1行开始放入第二列，第二列的其余元素为0。第n列由b_1，n ⁽ⁿ⁾(i)构成，相对于第n-1列第n列向下平移SF₁行，即：b_1，n ⁽ⁿ⁾(i)从第(n-1)SF_k+1行开始放入第n列。从第

列到第

列，第

列的前SF_k+W-1是由b_k，1 ⁽ⁿ⁾(i)的元素构成，其余元素都为0；第

列由b_k，2 ⁽ⁿ⁾(i)的元素构成，b_k，2 ⁽ⁿ⁾(i)放入第

列的起始位置相对于前一列下移SF_k行；第

列由b_k，m ⁽ⁿ⁾(i)的元素构成，b_k，m ⁽ⁿ⁾(i)放入第列的起始位置相对于b_k，m-1 ⁽ⁿ⁾(i)放入前一列的起始位置向下平移SF_k行，这里m∈{2，...，N_k}，N_k表示UE发送的符号总数。

根据新的系统矩阵A′(i)，得到每个UE发送的消除了频偏的符号序列估计值d′_k，i，m。

这里，基于MMSE算法，利用新的系统矩阵A′(i)分别求解方程 $d_{\mathrm{MMSE}-i}^{\′}={(A^{\text{\′H}}\left(i\right)R_n^{-1}{A}^{\′}\left(i\right)+R_d^{-1})}^{-1}{A}^{\′H}\left(i\right)R_n^{-1}{e}_i,$ 得到每个UE在每个数据域发送的符号序列的消除了频偏的估计值d′_k，i，m。具体求解方程 $d_{\mathrm{MMSE}-i}^{\′}={(A^{\text{\′H}}\left(i\right)R_n^{-1}{A}^{\′}\left(i\right)+R_d^{-1})}^{-1}{A}^{\′H}\left(i\right)R_n^{-1}{e}_i$ 的方法与现有技术相同，此处不再赘述。

综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1、一种高速联合检测的方法，其特征在于，该方法包括以下步骤：

a、确定每个用户设备UE的消除了频偏的信道冲击响应；

b、根据每个UE的消除了频偏的信道冲击响应检测发送符号，得到每个UE发送的消除了频偏的符号序列估计值。

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

a1、对每根天线上的接收信号进行信道冲击响应的估计，得到每根天线的信道冲击响应 $h^{\left(n\right)}={[h_1^{\left(n\right)},...,h_k^{\left(n\right)},...,h_K^{\left(n\right)}]}^T,$ 其中， $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},...,h_{k,l}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}]$ 为第k个UE在第n根天线上的信道冲击响应，n为天线编号，k为UE编号，1≤k≤K，l为路径编号，1≤l≤W；

a2、对每根天线的信道冲击响应h⁽ⁿ⁾进行频偏补偿，得到每个UE的消除了频偏的信道冲击响应h′⁽ⁿ⁾。

3、根据权利要求2所述的方法，其特征在于，所述步骤a2包括：

根据每根天线的信道冲击响应h⁽ⁿ⁾，确定每个UE在每条路径上的频率偏差Δf_k，l；

根据每个UE在每条路径上的频率偏差Δf_k，l，对G矩阵进行纠正频偏得到矩阵G′＝[G′⁽¹⁾，…，G′^(k)，...，G′^(K)]，其中，

$G^{\left(k\right)}=\left[\begin{array}{ccccc}{e}^{j{a}_1}{\underset{\‾}{m}}_{W+(K-k)W}& e^{j{a}_2}{\underset{\‾}{m}}_{W-1+(K-k)W}& ...& e^{j{a}_{W-1}}{\underset{\‾}{m}}_{2+(K-k)W}& e^{j{a}_W}{\underset{\‾}{m}}_{1+(K-k)W}\\ e^{j(a_1+{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+1+(K-k)W}& e^{j(a_2+{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+(K-k)W}& ...& e^{j(a_{W-1}+{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{3+(K-k)W}& e^{j(a_W+{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{2+(K-k)W}\\ e^{j(a_1+2{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+2+(K-k)W}& e^{j(a_2+2{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+1+(K-k)W}& ...& e^{j(a_{W-1}+2{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{4+(K-k)W}& e^{j(a_W+2{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{3+(K-k)W}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ e^{j(a_1+127{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+127+(K-k)W}& e^{j(a_2+127{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+126+(K-k)W}& ...& e^{j(a_{W-1}+127{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{129+(K-k)W}& e^{j(a_W+127{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{128+(K-k)W}\end{array}\right]$

m＝[m ₁，m ₂，m ₃，…，m ₁₂₈]是由小区的基本MIDAMBLE码经过QPSK调制而成的复值符号序列，Δ_l＝2π·Δf_k，l·T_c，1≤l≤W， $a_l=e^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_{k,l}\·(352+16)\·T_c}$ 为第l条路径上MIDAMBLE域的第17个码片的初始相位，T_c为每个码片的持续时间；

根据G′和每根天线的信道冲击响应h⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

4、根据权利要求2所述的方法，其特征在于，所述步骤a2包括：

根据每根天线的信道冲击响应h⁽ⁿ⁾，确定每个UE在每条路径上的频率偏差Δf_k，l，其中，k为UE编号，l为路径编号；

根据每个UE在每条路径上的频率偏差，对G矩阵进行纠正频偏得到矩阵G′＝[G′⁽¹⁾，G′⁽²⁾，...，G′^K]，其中，

$G^{\left(k\right)}=\left[\begin{array}{ccccc}{e}^{j{a}_1}{\underset{\‾}{m}}_{W+(K-k)W}& e^{j{a}_2}{\underset{\‾}{m}}_{W-1+(K-k)W}& ...& e^{j{a}_{W-1}}{\underset{\‾}{m}}_{2+(K-k)W}& e^{j{a}_W}{\underset{\‾}{m}}_{1+(K-k)W}\\ e^{j(a_1+{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+1+(K-k)W}& e^{j(a_2+{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+(K-k)W}& ...& e^{j(a_{W-1}+{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{3+(K-k)W}& e^{j(a_W+{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{2+(K-k)W}\\ e^{j(a_1+2{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+2+(K-k)W}& e^{j(a_2+2{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+1+(K-k)W}& ...& e^{j(a_{W-1}+2{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{4+(K-k)W}& e^{j(a_W+2{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{3+(K-k)W}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ e^{j(a_1+127{\mathrm{\Δ}}_1)}{\underset{\‾}{m}}_{W+127+(K-k)W}& e^{j(a_2+127{\mathrm{\Δ}}_2)}{\underset{\‾}{m}}_{W+126+(K-k)W}& ...& e^{j(a_{W-1}+127{\mathrm{\Δ}}_{W-1})}{\underset{\‾}{m}}_{129+(K-k)W}& e^{j(a_W+127{\mathrm{\Δ}}_W)}{\underset{\‾}{m}}_{128+(K-k)W}\end{array}\right]$

其中，m＝[m ₁，m ₂，m ₃，…，m ₁₂₈]是小区的基本MIDAMBLE码经过QPSK调制而成的复值符号序列，Δ_l＝2π·Δf_k，l·T_c， $a_l=e^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_{k,l}\·(352+16)T_c}$ 为第l条路径上MIDAMBLE域的第17个码片的初始相位，T_c为每个码片持续的时间；

从每个UE的路径中选取M条强径l_k，1…l_k，M；

在h_k ⁽ⁿ⁾中仅保留第k个UE的M条强径所对应的信道冲击响应值，得到新的矢量h′_k ⁽ⁿ⁾；根据h′_k ⁽ⁿ⁾得到 $h^{\′\left(n\right)}={[h_1^{\′\left(n\right)},...,h_k^{\′\left(n\right)},...,h_K^{\′\left(n\right)}]}^T;$

在G′^(k)中仅保留第k个UE的M条强径的路径编号所对应的M列元素得到矩阵G″^(k)；根据G″^(k)得到G″＝[G″⁽¹⁾，…，G″^(k)，…，G″^(K)]；

根据G″和h′⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

5、根据权利要求3所述的方法，其特征在于，所述步骤a2中根据G′和每根天线的信道冲击响应h⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应包括：

根据G′计算(G′)^H；

根据R_h，n＝E(v⁽ⁿ⁾(v⁽ⁿ⁾)^H)＝σ²I计算R_h，n ^-1，其中，v⁽ⁿ⁾为第n根天线的训练序列MIDMABLE域接收信号的噪声失量；

计算R_h＝E(h⁽ⁿ⁾(h⁽ⁿ⁾)^H)；

根据(G′)^H、R_h，n ^-1、R_h和第n根天线的MIDAMBLE域接收信号e′_n，mid，基于MMSE算法计算 $h_{\mathrm{MMSE}}^{\left(n\right)}={({\left(G^{\′}\right)}^H{R}_{h,n}^{-1}{G}^{\′}+R_h^{-1})}^{-1}{\left(G^{\′}\right)}^H{R}_{h,n}^{-1}{e}_{n,\mathrm{mid}}^{\′},$ 得到每根天线消除了频偏的信道冲击响应h_MMSE ⁽ⁿ⁾。

6、根据权利要求4所述的方法，其特征在于，所述步骤a2中根据G″和h′⁽ⁿ⁾确定每个UE的消除了频偏的信道冲击响应包括：

根据G″计算(G″)^H；

根据R_h，n＝E(v⁽ⁿ⁾(v⁽ⁿ⁾)^H)＝σ²I计算R_h，n ^-1，其中，v⁽ⁿ⁾为第n根天线的MIDMABLE域接收信号的噪声失量；

计算 $R_{h^{\′}}=E\left(h^{\′\left(n\right)}{\left(h^{\′\left(n\right)}\right)}^H\right)=\left[\begin{array}{cccc}{p}_{1,l_{1,1}}^{\left(n\right)}& 0& ...& 0\\ 0& p_{1,l_{\mathrm{1,2}}}^{\left(n\right)}& ...& 0\\ .& .& & .\\ .& .& ...& .\\ .& .& & .\\ 0& 0& ...& p_{K,l_{K,M}}^{\left(n\right)}\end{array}\right],$ 其中， $p_{k,l}^{\left(n\right)}={\left|\right|h_{k,l}^{\left(n\right)}\left|\right|}^2$ 为第k个UE的第l条路径的功率值，l∈{l_k，1，…，l_k，M}；

根据(G″)^H、R_h，n ^-1、R_h′和第n根天线的MIDAMBLE域接收信号e′_n，mid，基于MMSE算法计算 $h_{\mathrm{MMSE}}^{\′\left(n\right)}={({\left(G^{\′\′}\right)}^H{R}_{h,n}^{-1}{G}^{\′\′}+R_{h^{\′}}^{-1})}^{-1}{\left(G^{\′\′}\right)}^H{R}_{h,n}^{-1}{e}_{n,\mathrm{mid}}^{\′},$ 得到h′_MMSE ⁽ⁿ⁾，h′_MMSE ⁽ⁿ⁾由每个UE的M条强径的消除了频偏的信道冲击响应值构成；用h′_MMSE ⁽ⁿ)更新h⁽ⁿ⁾得到h_MMSE ⁽ⁿ⁾，所述更新h⁽ⁿ⁾的方法为：用h′_MMSE ⁽ⁿ⁾中每个UE的M条强径的消除了频偏的信道冲击响应值替代h_k ⁽ⁿ⁾中对应路径的信道冲击响应值，得到h″_k ⁽ⁿ⁾；根据h″_k ⁽ⁿ⁾得到每根天线消除了频偏的信道冲击响应 $h_{\mathrm{MMSE}}^{\left(n\right)}={[h_1^{\′\′\left(n\right)},...,h_k^{\′\′\left(n\right)},...,h_K^{\′\′\left(n\right)}]}^T.$

7、根据权利要求1至6所述的任一方法，其特征在于，所述步骤b包括：

b1x、根据每个UE的消除了频偏的信道冲击响应，构造系统矩阵A；

b2x、根据系统矩阵A，得到每个UE发送的符号序列估计值其中，k为UE编号，i＝1，2为数据域编号，m∈{1，...，N_k}为符号编号，N_k为第k个UE发送的符号总数；

b3x、利用公式 $d_{k,i,m}^{\′}={\hat{d}}_{k,i,m}{e}^{j2\mathrm{\π}\·\mathrm{\Δ}{f}_k\·\left((m-0.5)S{F}_k\right)\·T_c}$ 对

进行纠正频偏，得到每个UE发送的消除了频偏的符号序列估计值d′_k，i，m；

其中，SF_k为第k个UE的扩频因子，Δf_k为第k个UE的路径频率偏差平均值，Tc为每个码片的持续时间。

8、根据权利要求7所述的方法，其特征在于，所述步骤b3x中 $\mathrm{\Δ}{f}_k=\frac{1}{W}\underset{i=1}{\overset{W}{\mathrm{\Σ}}}\mathrm{\Δ}{f}_{k,l},$ 其中，W为信道冲击响应的窗长，Δf_k，l为第k个UE的第l个路径的频率偏差，Tc为每个码片的持续时间。

9、根据权利要求7所述的方法，其特征在于，所述步骤b3x中 $\mathrm{\Δ}{f}_k=\frac{1}{M}\underset{i=1}{\overset{M}{\mathrm{\Σ}}}\mathrm{\Δ}{f}_{k,l},$ M为每个UE的强径数量，Δf_k，l为第k个UE的路径编号为l的强径的频率偏差，Tc为每个码片的持续时间。

10、根据权利要求1至6所述的任一方法，其特征在于，所述步骤b包括：

b1y、根据每个UE的消除了频偏的信道冲击响应 $h_k^{\left(n\right)}=[h_{k,1}^{\left(n\right)},...,h_{k,l}^{\left(n\right)},...,h_{k,W}^{\left(n\right)}],$ 确定b_k，m ⁽ⁿ⁾(i，l)，

其中，b_k，m ⁽ⁿ⁾(i，l)是由W+SF_k-1个元素构成的列矢量，b_k，m ⁽ⁿ⁾(i，l)的前l-1个元素为0，后W-l个元素是0，n为天线编号，k为UE编号，m为符号编号，m∈{1，...，N_k}，N_k表示第k个发送的符号总数，i为数据域编号，i＝1，2，l为路径编号，1≤l≤W，SF_k为第k个UE的扩频因子，W为信道冲击响应的窗长，b_k，m ⁽ⁿ⁾(1，l)的具体结构为：

$b_{k,m}^{\left(n\right)}\left(\mathrm{1,1}\right)=\left[\begin{array}{c}{u}_{k,1}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·((m-1)S{F}_k+1)T_c}\\ u_{k,2}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,1}\·((m-1)S{F}_k+2)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,1}^{\left(n\right)}\·e^{j2\mathrm{\π}{\mathrm{\Δf}}_{k,1}\·\left(m{\mathrm{SF}}_k\right)T_c}\\ 0\\ .\\ .\end{array}\right]$

$b_{k,m}^{\left(n\right)}\left(\mathrm{1,2}\right)=\left[\begin{array}{c}0\\ u_{k,1}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·((m-1){\mathrm{SF}}_k+2)T_c}\\ u_{k,2}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·((m-1){\mathrm{SF}}_k+3)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,2}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,2}\·(m{\mathrm{SF}}_k+1)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

......

$b_{k,m}^{\left(n\right)}(1,l)=\left[\begin{array}{c}0\\ .\\ .\\ .\\ 0\\ u_{k,1}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π}{\mathrm{\Δf}}_{k,l}\·((m-1){\mathrm{SF}}_k+l)T_c}\\ u_{k,2}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π}{\mathrm{\Δf}}_{k,l}\·((m-1){\mathrm{SF}}_k+l+1)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,l}^{\left(n\right)}\·e^{j2\mathrm{\π}{\mathrm{\Δf}}_{k,l}(m{\mathrm{SF}}_k+l-1)T_c}\\ 0\\ .\\ .\\ .\\ 0\end{array}\right]$

......

$b_{k,m}^{\left(n\right)}(1,W)=\left[\begin{array}{c}0\\ .\\ .\\ .\\ 0\\ u_{k,1}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}\·((m-1){\mathrm{SF}}_k+W)T_c}\\ u_{k,2}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π\Δ}{f}_{k,W}((m-1){\mathrm{SF}}_k+W+1)T_c}\\ .\\ .\\ .\\ u_{k,{\mathrm{SF}}_k}{h}_{k,W}^{\left(n\right)}\·e^{j2\mathrm{\π}{\mathrm{\Δf}}_{k,W}\·(m{\mathrm{SF}}_k+W-1)T_c}\end{array}\right],$

其中， $u_k=[u_{k,1},u_{k,2},...,u_{k,{\mathrm{SF}}_k}],$ u_k为第k个UE的扩频码，Δf_k，l为第k个UE的第l条路径的频率偏差；

所述b_k，m ⁽ⁿ⁾(2，l)根据公式 $b_{k,m}^{\left(n\right)}(2,l)=b_{k,m}^{\left(n\right)}(1,l)\·e^{j2\mathrm{\π}\·{\mathrm{\Δf}}_{k,l}\·(352+144)\·T_c}$ 确定；

b2y、按照 $b_{k,m}^{\left(n\right)\left(i\right)}=\underset{l=1}{\overset{W}{\mathrm{\Σ}}}{b}_{k,m}^{\left(n\right)}(i,l),$ 确定b_k，m ⁽ⁿ⁾(i)；

b3y、利用b_k，m ⁽ⁿ⁾(i)确定A⁽ⁿ⁾′(i)，并按照 $A^{\′}\left(i\right)=\left[\begin{array}{c}{A}^{\left(1\right)\′}\left(i\right)\\ .\\ .\\ .\\ A^{\left(n\right)\′}\left(i\right)\\ .\\ .\\ .\\ A^{\left(N\right)\′}\left(i\right)\end{array}\right]$ 确定包含了路径频偏的新的系统矩阵A′(i)；其中，所述利用b_k，m ⁽ⁿ⁾(i)确定A⁽ⁿ⁾′(i)的方法为：

将A⁽ⁿ⁾′(i)的前N₁列对应第1个UE，紧接着的N₂列对应第2个UE，从第

列到第

列对应第k个UE，1≤k≤K，最后第N_K列对应第K个UE；

其中，在第k个UE对应的N_k列中，将列矢量b_k，n ⁽ⁿ⁾(i)的SF_k+W-1个元素依次放入第n列的从第(n-1)SF_k+1行开始的W+SF_k-1行，将第n列其余行的元素均设为0，n为所述第k个UE对应的N_k列中的第n列，1≤n≤N_k；

b4y、根据新的系统矩阵A′(i)，得到每个UE发送的消除了频偏的符号序列估计值d′_k，i，m。

Images