[发明专利]一种适用于直接序列扩频超宽带系统的信道估计方法

权利要求书

1.一种适用于直接序列扩频超宽带系统的信道估计方法，在该系统的发射机与接收机之间采用基于数据包的传输策略，发射机发出的每个数据包以1个导频数据块开始,后面为M个净荷数据块，其中导频数据块由收发双方约定的已知数据组成，用于信道估计，净荷数据块由信源产生的随机业务数据组成；导频数据块和每个净荷数据块具有相同的长度，均由N个数据符号组成；

发射机工作流程如下：

信源产生的随机业务数据被组合成净荷数据块，每M个净荷数据块附加1个导频符号块后构成1个数据包，数据包经二进制相移键控调制、直接序列扩频和超宽带脉冲成形后，被超宽带天线发送到无线信道，一个数据块对应的发送信号s(t)表示为

s(t)=Σi=1Ndiw(t-iT)---(1)]]>

w(t)=Σj=1Nscjg(t-jTc)---(2)]]>

其中，d_i∈{+1,-1}为BPSK调制后的数据符号，i表示数据符号的序号，N表示一个数据块的长度；T表示符号周期；c_j表示扩频序列的一个码片，下标j表示码片的序号，N_s表示扩频序列长度，T_c表示扩频序列的码片周期；g(t)为超宽带脉冲波形；w(t)为扩频序列与超宽带脉冲g(t)卷积得到的脉冲序列波形；t表示时间；w(t-iT)表示平移i个符号周期后的脉冲序列波形；g(t-jT_c)表示平移j个码片周期后的超宽带脉冲波形；

超宽带无线信道的结构如下：

信道冲激响应h(t)表示为一个抽头延迟线模型

h(t)=Σl=0L-1αlδ(t-lTc)---(3)]]>

其中，l表示多径分量的序号，L代表多径分量的数目，α_l表示第l条径的增益，lT_c为第l条径的时间延迟，T_c是码片周期，t表示时间，δ(t)是Dirac函数，时间延迟的最大值(L-1)T_c大于一个符号周期T但小于一个数据块周期NT，显然，前N_s条多径即多径分量0,1,…,N_S-1的时间延迟小于一个符号周期T，后面L-N_s条径将延伸到后续一个至几个符号周期；

接收机工作流程如下：

一个数据块期间，接收机天线从无线信道接收到的信号r(t)表示为

r(t)=s(t)*h(t)+n(t)=Σi=1NdiΣl=0L-1αlw(t-iT-lTc)+n(t)---(4)]]>

其中，*表示卷积运算，n(t)是零均值、单边功率谱密度为N₀（瓦/赫兹）的高斯噪声，接收信号r(t)首先经过波形相关器，采样电路以码片速率对波形相关器输出进行采样并储存，所得采样值为

zi,j=∫0Tcr(t+iT+jTc)g(t)dt---(5)]]>

其中，z_i,j表示数据块中第i个数据符号的第j个采样值，i表示数据符号的序号，j表示采样值的序号，公式（4）中r(t+iT+jT_c)表示接收信号r(t)平移i个符号周期和j个码片周期后的结果，信道估计模块利用采样值序列z_i,j,i＝1,2,...,N;j＝1,2,...,N_S，估计出多径信道信息即各多径分量的增益；波形级Rake合并模块根据信道估计结果选出一定数量的多径分量，采用最大比合并将对应采样值z_i,j进行加权求和，对于某净荷数据块第i个符号的第j个码片，Rake合并后的输出变量为

zi,jr=Σf=1Nfαfzi,j-τf---(6)]]>

其中，f表示Rake合并支路的序号，N_f为Rake合并支路数目，α_f和τ_f分别表示第f支路的加权值和时间延迟，表示第f支路对应采样值，对于某净荷数据块第i个符号，波形级Rake合并后的输出变量序列经过解扩处理，得到解扩后变量

Zi=Σj=1Nscjzi,jr---(7)]]>

其中，Z_i表示第i个符号的解扩后变量，即发射机采用的扩频序列，c_j表示扩频序列的一个码片，j表示码片的序号，为了消除符号间干扰，采用各种均衡器对（7）式中解扩后变量Z_i进行处理，均衡器的抽头系数来自于信道估计结果，对于第i个符号，假设均衡器输出变量用表示，则通过检测Z_i的极性可恢复数据符号，即

d^i=+1,Z~i>0-1,Z~i<0---(8)]]>

最后，（8）式中恢复的数据符号经BPSK逆映射还原为二进制数据，当整个数据包接收完成后，按数据包结构解封装后输出至信源；

本信道估计方法步骤如下：

(1)利用数据包中导频数据块对应的波形相关器输出采样值序列，接收机估计出时间延迟小于一个符号周期的所有多径分量的增益；

时间延迟小于一个符号周期的多径分量即公式（3）中第0至第N_S-1条径，将导频数据块中每个数据符号对应的波形相关器采样序列表示为N_s维列向量，第i个符号对应的列向量Z_1,i为

Z1,i=zi,1zi,2···zi,NS,i=1,2,...,N---(9)]]>

利用公式（9）对导频数据块中N个符号的列向量进行加权平均，得加权平均值

Z‾1=1NΣi=1Ndi·Z1,i---(10)]]>

其中，d_i∈{+1,-1}为BPSK调制后的导频数据符号，i表示导频数据符号的序号，利用发射机扩频序列码片构成如下三角矩阵

其中，c_j,j＝1,2,…,N_S为扩频序列的码片，j表示码片的序号，基于公式（10）和（11），计算时间延迟小于一个符号周期的所有多径分量的增益如下

α^0α^1······α^NS-1=C-1Z‾1---(12)]]>

其中，为第j条径的增益估计值，C^-1表示矩阵C的逆矩阵；

（2）利用数据包中导频数据块对应的波形相关器输出采样值序列以及步骤（1）估计结果，接收机估计出时间延迟大于一个符号周期的所有多径分量的增益；

时间延迟大于一个符号周期的多径分量分布在第2个至第p个符号周期，每个符号周期内包含N_s条径，通常p小于数据块长度N；在步骤（1）估计结果的基础上，我们通过一个迭代过程，从第2个符号周期开始，逐个符号周期分别估计出对应多径分量的增益；对于第k（2≤k≤p）个符号周期，待估计的多径分量为第(k-1)N_s至第kN_s-1条径，具体算法如下：首先，对导频数据块内所有N个导频数据符号，利用已有估计结果消除公式（9）Z_1,i中前面已估计出的多径分量，即作如下计算：

zk,i=zk-1,i-di-(k-2)·Σn=11cnα^(k-2)NS+(1-n)Σn=12cnα^(k-2)NS+(2-n)···Σn=1Nscnα^(k-2)NS+(NS-n)-di-(k-1)·Σn=2Nscnα^(k-1)NS-n+1Σn=3Nscnα^(k-1)NS-n+2···Σn=NsNscnα^(k-1)NS-n+(NS-1)0,i=k,k+1,...,N---(13)]]>

其中，Z_k,i表示消除前面(k-1)N_s条径后波形相关器采样序列组成的向量，k表示待估计的符号周期的序号，i表示导频数据块中符号的序号，d_i-(k-2)表示第i-(k-2)个导频数据符号，d_i-(k-1)表示第i-(k-1)个导频数据符号，c_n，n＝1,2,…,N_S为扩频序列的码片，n表示码片的序号，N_S为扩频序列长度，α^(k-2)NS+(1-n),]]>α^(k-2)NS+(2-n),]]>α^(k-2)NS+(NS-n),]]>α^(k-1)NS-n+1,]]>α^(k-1)NS-n+1,]]>α^(k-1)NS-n+2,]]>α^(k-1)NS-n+(NS-1)]]>表示已估计出的某条径的增益，它们的下标计算结果表示对应径的具体序号；其次，在导频数据块第k个至第N个导频符号周期内，对公式（13）中Z_k,i进行加权求和计算

Z‾k=1N-(k-1)Σi=kNdi-(k-1)Zk,i---(14)]]>

最后，估计出第k个符号周期内所有多径分量的增益如下

α^(k-1)NSα^(k-1)NS+1···α^kNS-1=C-1Z‾k---(15)]]>

其中，α^j,j=(k-1)NS,(k-1)NS+1,...,KNS-1]]>为第j条径的增益估计值，C^-1表示矩阵C的逆矩阵；

（3）以步骤（1）和（2）估计结果作为信道估计初始值，对数据包中的净荷数据块进行逐块接收判决，每个净荷数据块的判决结果又被视作导频数据，并采用步骤（1）和步骤（2）方法重复估计多径分量的增益，通过迭代计算更新估计值，并将更新后的结果用于下一净荷数据块的接收判决，当达到规定的迭代次数后，停止估计与更新过程；

设用向量H₀表示步骤（1）和（2）估计结果，即信道估计初始值；为了进一步降低由噪声引起的估计误差，我们采用判决反馈机制对信道估计结果进行I次（I＜＜M）迭代修正，以提高估计精度，具体步骤如下：首先，用已有的信道估计初始值H₀进行Rake合并和均衡，经判决得到第1个净荷数据块的数据，这些已判决的数据被视作导频数据符号用于信道估计，具体估计过程如上述步骤（1）和（2）；假设估计结果为H’₁，利用H’₁修正信道估计初始值H₀如下

H₁＝(H'₁+H₀)/2         (16)

修正后的信道信息H₁用于第2个净荷数据块的接收与判决，判决恢复的数据又用于信道估计和修正，估计过程如上述步骤（1）和步骤（2），修正方法同公式（16），假设估计结果为H’₂，则修正后的信道信息为

H₂＝(H'₂+H₁)/2        (17)

上述迭代修正过程使信道估计精度将逐渐提高，经过I步迭代后，得到最终的估计值表示为H_I；剩余的净荷数据块都用H_I进行接收与判决；迭代次数I的取值根据实际需求选择，以便在计算复杂度和系统性能之间取得合理的折中，当达到规定的迭代次数后，停止估计与更新过程。