[发明专利]一种快变OFDM系统的渐进迭代时变信道估计和ICI消除方法

权利要求书

1.快变OFDM系统的渐进迭代时变信道估计和ICI消除方法，其特征在于，在OFDM系统的信道估计中，将ICI与噪声之和SIN作为Kalman滤波器的去噪对象，从仅使用导频点信息开始，渐进地增加用于迭代计算的信息，依以下步骤实现：

步骤(1)，发射端产生发送数据，将导频数据按照梳状导频方式插入到发送数据中：

发送端设定如下：s表示第s个OFDM符号，s＝1,2,…,s,…S，每个OFDM符号包含N个子载波，n＝1,2,…,n,…,N，其中包含N_p个导频符号和N_d个数据符号，N_d+N_p＝N，n_p＝1,2,…,N_p，导频在频域上的位置矩阵表示成：其中且保证N_p≥L，L是信道多径数l的最大值，即l＝1，2，…,,l…,L，N_p个导频被平均的插入到N个载波之中并且在传输过程中保持不变，在N个载波中导频点符号表示为xp(s)=x(s)(Ps)=[xp1(s),xp2(s),···,xpNp(s)]T,]]>

步骤(2)，数据通过OFDM系统发送至接收端，在接收端去掉循环前缀后，按以下步骤用多项式基扩展模型P-BEM对信道进行建模：

步骤(2.1)，利用多项式基扩展模型P-BEM来描述具有时频双选特性的多径传播信道，则第S个OFDM符号的第n个子载波第l径的信道冲激响应h^(s)(n,l)表示为：

h^(s)(n,l)＝QC_l^(s)+ξ_l^(s)(n),0≤n≤N-1，

其中，ξ_l^(s)表示建模时每一个OFDM符号的第l径的模型误差，其值小于10^-3，在计算时忽略，即认为h^(s)(n,l)＝QC_l^(s)，Q是一个N×B的正交基函数矩阵，C_l^(s)则是由基函数对应的B个系数组成的向量Cl(s)=[c1,l(s),c2,l(s),···,cb,l(s),···,cB,l(s)]T,]]>f_max是信道的最高频率，T_s是采样时间，

步骤(2.2)，将在接收端的接收信号表示成以下形式：

y^(s)＝H^(s)x^(s)+W^(s)，

其中，x^(s)＝[x₁^(s),x₂^(s)…x_N^(s)]^T、y^(s)＝[y₁^(s),y₂^(s),…,y_N^(s)]^T分别表示频域上第s个去掉循环前缀后的发送信号和接收信号，W^(s)是其频域上的白噪声，H^(s)是N×N的信道矩阵：

其中，矩阵的每一个元素为多径信道的信道冲击响应的和，计算方式如下：

H(s)(m,k)=Σl=0L-1Gl(s)(M,K)e-j2π(k-1N-12)τl,]]>

m，k表示上述矩阵H^(s)的第m行k列的值，m＝1,2,…,m,…,M，M≤N，k＝1,2,…,k,…,K，K≤N，τ_l是第l径的时延，G_l^(s)(M,K)为信道冲击相应的频域表达矩阵，其每一个元素计算如下：

Gl(s)(m,k)=1NΣn=0N-1h(s)(n,l)e-j2π(m-k)n/N,]]>

步骤(2.3)，根据P-BEM模型将接收信号进行重新建模，表示成带有P-BEM系数的表达式如下：

y^(s)＝Φ^(s)g^(s)+W^(s)，

其中，

g^(s)=[C₁^(s)T,C₂^(s)T…C_L^(s)T]^T，表示P-BEM算法中的系数矩阵，

表示重新建模后，与发送数据相关的系数矩阵，其计算方法如下：

Zl(s)=1N[D1diag(x(s))Γl,···,Dbdiag(x(s))Γl,···,DBdiag(x(s))Γl],]]>其中，

Γl=e-j2π(p1N-12)τle-j2π(p2-1N-12)τl···e-j2π(pNp-1N-12)τlT,]]>是第l径的傅里叶变换，

Γ＝[Γ₁,Γ₂,…,Γ_L]，表示L径的总傅里叶变换矩阵，

diag(x^(s))表示以向量x^(s)为对角元素的矩阵，

步骤(3)，利用AR模型对信道BEM系数进行建模：

步骤(3.1)，按下式计算C_l^(s)相关矩阵：

RCl(j)=((Q)HQ)-1(Q)HRh(n,l)(j)Q((Q)HQ)-1,]]>

其中，j表示相关的阶数，即进行相关运算的OFDM符号的符号间隔，j的取值为[-1,0,1]，分别表示当前OFDM信号的C_l^(s)和前一个符号的C_l^(s-1)的相关矩阵，当前OFDM信号的C_l^(s)的自相关矩阵，当前OFDM信号的C_l^(s)和后一个符号的C_l^(s+1)的相关矩阵。(·)^H表示Hermitian运算，Rh(n,l)(j)=E[h(n,l)h*(n+j,l)]=σh(n,l)2J0(2πfdTsj),]]>其中E[·]表示均值,J₀(·)表示第一类的零阶贝塞尔函数，f_d=vf_c/c是终端的移动速度为v时的最大多普勒频移，f_c是载波频率，c是光速，代表第l径的信道冲击响应的方差，并假设

步骤(3.2)，根据Yule-Walker方程得到信道P-BEM参数的状态转移方程：

g_(s)＝Ag_(s-1)+U_(s)，

将OFDM系统发送符号的时间顺序g^(s)看做控制系统中状态转移过程g_(s)，即g^(s)＝g_(s)，状态转移方程系数A＝diag(a₁,a₂,…,a_l,…a_L)，diag(x)表示以向量x为对角元素的矩阵，U_(s)代表第s个OFDM符号的AR模型的建模误差；

步骤(4)，对Kalman滤波器进行初始化并计算初始更新方程：

步骤(4.1)，按照下式对Kalman滤波器进行初始化：

g^(0|0)=OLB,1P(0|0)=diag(RC1(0),RC2(0)···RCL(0)),]]>

形式如和P_(s|s)中下标前一个s均表示当前的状态为g_(s)，后一个s表示第s个OFDM符号，P_(0|0)是用于计算的初始值，表示OFDM符号的g_(s)的初始值，P_(0|0)表示对应的误差相关矩阵，O_LB,1是LB×1的零矩阵，

步骤(4.2)，按下式计算Kalman的初始时间更新方程：

i＝1，s＝1，

g^(s)=Ag^(0|0),]]>

P_(s)＝AP_(0|0)(A)^H+V[U_(s)]，

i表示迭代次数，表示Kalman方程中状态估计g_(s)的中间变量，P_(s)表示中间变量对应的误差相关矩阵；用V[·]表示协方差矩阵，V[U_(s)]=diag(u₁,u₂…u_L)

步骤(5)，进行第一次信道估计迭代运算，此时迭代次数i＝1，本次迭代中仅使用导频点所在的子载波处接收到数据做信道估计，除了导频点所在子载波外其他子载波上的数据视为ICI，用SIN方法消除未知数据对导频处信道估计的影响，实现无ICI干扰的导频辅助Kalman信道估计，具体步骤如下：

步骤(5.1)，仅将接收信号中各个导频点对应的载波位置的接收信号用于计算，并将接收信号分成各个导频点所在子载波上的数据，除导频点子载波外的其它子载波上的数据对各个导频点所在子载波的干扰和噪声三部分，如下式所示：

yp(s)=y(s)(Ps)=H(s)[Ps,Ps]xpnp(s)+H(s)[Ps,dnpn′]xdnpn′(s)+W(s)(Ps)]]>

其中，Ps=[p1,p2,···,pNp],]]>xdnpn′(s)=[x(s)(d11),x(s)(d12),···,x(s)(d1N′),x(s)(d21),···,x(s)(dNpN′)]T,]]>n'＝1,2,…,N'表示相邻的各个导频点之间的距离，是一个N_p×N_p的单位矩阵，σ²高斯白噪声W^(s)的方差，上式中第二项为除导频点子载波外的其它子载波上的数据对导频点所在子载波的干扰ICI，

步骤(5.2)，将数据ICI干扰项考虑成信道噪声W^(s)(P_s)的一部分作为滤波器的去噪对象，将步骤（2）中的算法按照SIN估计的方法改写，令则SIN估计的Kalman观测方程表示为：

yp(s)=ΦSIN(s)g(s)+W(s)SIN,]]>

其中：

ΦSIN(s)=1N[z1(s)SIN,Z2(s)SIN···ZL(s)SIN],]]>

Zl(s)SIN=1N[D1SINdiag(xp(s))ΓlSIN···DBSINdiag(xp(s))ΓlSIN],]]>

ΓlSIN=e-j2π(p1N-12)τle-j2π(p2-1N-12)τl···e-j2π(pNp-1N-12)τlT,]]>

Γ^SIN＝[Γ₁^SIN,Γ₂,^SIN…,Γ_l^SIN]

步骤(5.3)，计算的协方差矩阵

计算中假设ICI为高斯白噪声，令因为噪声与ICI二者相互独立，所以V[W(s)SIN]=UICI+V[W(s)(Ps)];]]>

U_ICI矩阵中的每一个元素的计算式为：

UICI(m,k)=RICI(j)≈4π2Ts2Es(Σl=0L-1σh(n,l)2σDl)ρ(α,rag,N),]]>

其中，m，k表示矩阵的第m行k列，E_s是发送数据的功率，是功率为P_v时的多普勒功率普函数，f是传输频率，α表示用于迭代的最边缘数据距对应的各个导频点的距离，第一次计算时为0，rag为计算的精度，rag=[0,1,2,3]，并且:

ρ(α，rag，N)=ρ(0,rag,N)-ρ₁(α，rag,N)

步骤(5.4)，分别按下述三式计算kalman增益K_(s)，第s个OFDM符号转移到状态的状态估计矩阵和与对应的协防差矩阵P_(s|s)，构成观测更新方程组，其中，Φ=Φ^SINW(s)=W(s)SIN,]]>

K_(s)＝P_(s)(Φ_(s))^H(Φ_(s)P_(s)(Φ_(s))^H+V[W_(s)])^-1，

g^(s|s)=g^(s)+K(s)(y(s)-Φ(s)g^(s)),]]>

P_(s|s)＝P_(s)-K_(s)Φ_(s)P_(s)，

步骤(5.5)，根据下式计算出信道矩阵H^(s)的估计值：

H(s)=1NΣb=1BDbdiag(Γg^(s|s)),]]>

步骤(5.6)，利用下式对信道矩阵进行QR分解，得到矩阵R^(s)：

H^(s)＝IR^(s)

其中I是一个单位矩阵，R^(s)是一个上三角矩阵，

步骤(5.7)，按下式对数据进行QR数据检测：

其中y′^(s)＝(I)^Hy^(s)，和分别是数据的检测值和检测值星座图量化后的结果，[·]_m，k代表矩阵的第m行k列，[·]_m是向量的第m个元素，[·]_k是向量的第k个元素，O(·)表示解调运算，m，k表示矩阵H^(s)的第m行k列的值，m＝1,2,…,m,…，M，M≤N，k＝1,2,…,k,…,K，K≤N，

步骤(6)，迭代次数i=i+1，迭代计算次数i>1时，第二次迭代用导频点所在处的子载波和其两侧相邻的一个子载波处接收到的数据做信道估计，其余子载波上的数据视为ICI，其后迭代每次增加的用于计算的数据，均为上次迭代用于计算的数据所在的子载波两侧的子载波处接收到的数据，其余子载波上的数据视为ICI，SIN方法进行信道估计迭代运算如下所述：

步骤(6.1)，将接收信号分成导频点所在子载波数据，用于计算的数据所在的子载波数据，未用于计算的数据对导频点和用于计算的数据所在子载波的干扰和噪声四部分，

如下式所示：

上式中第三项为更新后的ICI干扰，

步骤(6.2)，步骤(5.2)中所述的方法计算Φ^SIN_(s)，

步骤(6.3)，随迭代次数增加，用于迭代的数据增加，此时更新的的协方差矩阵U_ICI计算如下：

UICI≈4π2Ts2Es(Σl=0L-1σh(n,l)2σDl)[12ρ(A+m,rag,N)+12ρ(A-m,rag,N)]]]>

其中，A+m表示导频点右侧的用于计算的数据到对应的导频点的距离，A-m表示导频点左侧的用于计算的数据到对应的导频点的距离，

步骤(6.4)，按步骤（5.5）所述方法计算kalman增益K_(s)，第s个OFDM符号转移到状态的状态估计矩阵和与对应的协防差矩阵P_(s|s)，构成观测更新方程组，其中，Φ=Φ^SIN，W(s)=W(s)SIN,]]>

K_(s)＝P_(s)(Φ_(s))^H(Φ_(s)P_(s)(Φ_(s))^H+V[W_(s)])^-1，

g^(s|s)=g^(s)+K(s)(y(s)-Φ(s)g^(s)),]]>

P_(s|s)＝P_(s)-K_(s)Φ_(s)P_(s)，

步骤(6.5)，按步骤（5.5）所述方法计算出信道矩阵的估计值H^(s)：

H(s)=1NΣb=1BDbdiag(Γg^(s|s)),]]>

步骤(6.6)，按步骤（5.6）所述方法对信道矩阵进行QR分解得到R^(s)：

H^(s)＝IR^(s)，

步骤(6.7)，按步骤（5.7）所述方法对数据进行QR数据检测：

步骤(7)，判断是否所有的数据是否已经都用于迭代，如果是，则算法结束，如果不是，则继续，

步骤(8)，通过判断迭代的次数，决定是否需要增加迭代算法的输入数据，判断算法

如下：

设定步长Δ，比较迭代次数i和Δ·μ+2，其中，μ=1,2,···,Δ·μ+2<NNp,]]>如果i＝Δ·μ+2，则选择已用于计算的数据所在的子载波两侧的子载波处的数据作为下次迭代新加入的用于计算的数据，与已用于计算的数据一起带入SIN算法，带入到步骤(6)中重新计；

如果i≠Δ·μ+2，则不增加用于计算的数据，返回步骤(6)进行迭代；

结束。