一种快变OFDM系统的渐进迭代时变信道估计和ICI消除方法
专利汇可以提供一种快变OFDM系统的渐进迭代时变信道估计和ICI消除方法专利检索，专利查询，专利分析的服务。并且一种快变OFDM系统的渐进 迭代 时变信道估计和ICI消除方法，用于无线通信信道估计领域，用于快变信道产生ICI严重影响情况下，基于导频的信道估计。其特征在于：将载波间干扰ICI与信道噪声之和(Sum of ICI and channel noise，SIN)作为Kalman 滤波器 的去噪对象，消除Kalman估计时ICI的影响。此外，用于迭代的数据采用渐进增长的方式，用于测量的子载 波数 量，在迭代的过程中沿各个导频点两侧缓慢增长，从而抑制由ICI产生的影响。在SNR较小的情况下，该方案的性能比起现有 算法 有显著的提升。，下面是一种快变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，其中包含Np个导频符号 和Nd个数据符号，Nd+Np＝N，np＝1,2,…,Np，导频在频域上的位置矩阵表示成： 其中
且保证Np≥L，L是信道多径数l的最大值，即 Np个
导频被平均的插入到N个载波之中并且在传输过程中保持不变，在N个载波中导频点符号表示为
步骤(2)，数据通过OFDM系统发送至接收端，在接收端去掉循环前缀后，按以下步骤用多项式基扩展模型P-BEM对信道进行建模：
步骤(2.1)，利用多项式基扩展模型P-BEM来描述具有时频双选特性的多径传播信道，S (s)
则第 个OFDM符号的第n个子载波第l径的信道冲激响应h (n,l)表示为：
(s) (s) (s)
h (n,l)＝QCl +ξl (n),0≤n≤N-1，
(s) -3
其中，ξl 表示建模时每一个OFDM符号的第l径的模型误差，其值小于10 ，在计算(s) (s) (s)
时忽略，即认为h (n,l)＝QCl ，Q是一个N×B的正交基函数矩阵，Cl 则是由基函数对应的B个系数 组成的向量 fmax是信道
的最高频率，Ts是采样时间，
步骤(2.2)，将在接收端的接收信号表示成以下形式：
(s) (s) (s) (s)
y ＝H x +W ，
(s) (s) (s) (s) T (s) (s) (s) (s) T
其中，x ＝[x1 ,x2 …xN ]、y ＝[y1 ,y2 ,…,yN ]分别表示频域上第s个去(s) (s)
掉循环前缀后的发送信号和接收信号，W 是其频域上的白噪声，H 是N×N的信道矩阵：
其中，矩阵的每一个元素为多径信道的信道冲击响应的和，计算方式如下：
m，k表示上述矩阵H(s)的第m行k列的值，m＝1,2,…,m,…,M，M≤N，k＝1,2,…,k,…,K，K≤N，τl是第l径的时延，Gl(s)(M,K)为信道冲击相应的频域表达矩阵，其每一个元素计算如下：
步骤(2.3)，根据P-BEM模型将接收信号进行重新建模，表示成带有P-BEM系数的表达式如下：
y(s)＝Φ(s)g(s)+W(s)，
其中，
g(s)＝[C1(s)T,C2(s)T…CL(s)T]T，表示P-BEM算法中的系数矩阵，
表示重新建模后，与发送数据相关的系数矩阵，其计算方法
如下：
其中，
是第l径的傅里叶变换，
Γ＝[Γ1,Γ2,…,ΓL]，表示L径的总傅里叶变换矩阵，
diag(x(s))表示以向量x(s)为对角元素的矩阵，
步骤(3)，利用AR模型对信道BEM系数进行建模：
步骤(3.1)，按下式计算Cl(s)相关矩阵：
其中，j表示相关的阶数，即进行相关运算的OFDM符号的符号间隔，j的取值为(s) (s-1)
[-1,0,1]， 分别表示当前OFDM信号的Cl 和前一个符号的Cl 的相关矩
(s) (s) (s+1)
阵，当前OFDM信号的Cl 的自相关矩阵，当前OFDM信号的Cl 和后一个符号的Cl 的相关矩阵(·)H表示Hermitian运算， 其中E[·]
表示均值，J0(·)表示第一类的零阶贝塞尔函数，fd＝vfc/c是终端的移动速度为v时的最大多普勒频移，fc是载波频率，c是光速， 代表第l径的信道冲击响应的方差，并假设步骤(3.2)，根据Yule-Walker方程得到信道P-BEM参数的状态转移方程：
g(s)＝Ag(s-1)+U(s)，
(s) (s)
将OFDM系统发送符号的时间顺序g 看做控制系统中状态转移过程g(s)，即g ＝g(s)，状态转移方程系数A＝diag(a1，a2，…，al，…aL)， diag(x)表示以向量x为对角元素的矩阵，U(s)代表第s个OFDM符号的AR模型的建模误差；
步骤(4)，对Kalman滤波器进行初始化并计算初始更新方程：
步骤(4.1)，按照下式对Kalman滤波器进行初始化：
形式如 和 中下标前一个s均表示当前的状态为g(s)，后一个s表示第s个OFDM符号， P(O|O)是用于计算的初始值， 表示OFDM符号的g(s)的初始值，P(O|O)表示对应的误差相关矩阵，0LB，1是LB×1的零矩阵，
步骤(4.2)，按下式计算Kalman的初始时间更新方程：
i＝1，s＝1，
H
P(s)＝AP(0|0)(A)+V[U(s)]，
i表示迭代次数， 表示Kalman方程中状态估计g(s)的中间变量，P(s)表示中间变量 对应的误差相关矩阵；用V[·]表示协方差矩阵，V[U(s)]＝diag(u1，u2…uL)，步骤(5)，进行第一次信道估计迭代运算，此时迭代次数i＝1，本次迭代中仅使用导频点所在的子载波处接收到数据做信道估计，除了导频点所在子载波外其他子载波上的数据视为ICI，用SIN方法消除未知数据对导频处信道估计的影响，实现无ICI干扰的导频辅助Kalman信道估计，具体步骤如下：
步骤(5.1)，仅将接收信号中各个导频点对应的载波位置的接收信号用于计算，并将接收信号分成各个导频点所在子载波上的数据，除导频点子载波外的其它子载波上的数据对各个导频点所在子载波的干扰和噪声三部分，如下式所示：
其中，
表示相邻的各个导频点之间的距离， 是一
2 (s)
个Np×Np的单位矩阵，σ 高斯白噪声W 的方差，上式中第二项为除导频点子载波外的其它子载波上的数据对导频点所在子载波的干扰ICI，
(s)
步骤(5.2)，将数据ICI干扰项考虑成信道噪声W (Ps)的一部分作为滤波器的去噪对象，将步骤(2)中的算法按照SIN估计的方法改写，令 则
SIN估计的Kalman观测方程表示为：
其中：
步骤(5.3)，计算 的协方差矩阵
计算中假设ICI为高斯白噪声，令 因为噪声与ICI二者相互
独立，所以
UICI矩阵中的每一个元素的计算式为：
其中，m，k表示矩阵的第m行k列，Es是发送数据的功率，
是功率为Pv时的多普勒功率普函数，f是传输频率，α表示用于迭
代的最边缘数据距对应的各个导频点的距离，第一次计算时为0，rag为计算的精度，rag＝[0,1,2,3]，并且:
ρ(α,rag,N)＝ρ(0,rag,N)-ρ1(α,rag,N)，
步骤(5.4)，分别按下述三式计算kalman增益K(s)，第s个OFDM符号转移到状态的状态估计矩阵 和与 对应的协防差矩阵P(s|s)，构成观测更新方程组，其中，Φ＝SIN
Φ ，
K(s)＝P(s)(Φ(s))H(Φ(s)P(s)(Φ(s))H+V[W(s)])-1，
P(s|s)＝P(s)-K(s)Φ(s)P(s)，
(s)
步骤(5.5)，根据下式计算出信道矩阵H 的估计值：
(s)
步骤(5.6)，利用下式对信道矩阵进行QR分解，得到矩阵R ：
(s) (s)
H ＝IR
(s)
其中I是一个单位矩阵，R 是一个上三角矩阵，
步骤(5.7)，按下式对数据进行QR数据检测：
′(s) H (s)
其中y ＝(I) y ， 和 分别是数据的检测值和检测值星座图量化后的结果，[·]m,k代表矩阵的第m行k列，[·]m是向量的第m个元素，[·]k是向量的第k个元素，(s)
O(·)表示解调运算，m，k表示矩阵H 的第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干扰，
SIN
步骤(6.2)，步骤(5.2)中所述的方法计算Φ (s)，
步骤(6.3)，随迭代次数增加，用于迭代的数据增加，此时更新的 的协方差矩阵UICI计算如下：
其中，A+m表示导频点右侧的用于计算的数据到对应的导频点的距离，A-m表示导频点左侧的用于计算的数据到对应的导频点的距离，
步骤(6.4)，按步骤(5.5)所述方法计算kalman增益K(s)，第s个OFDM符号转移到状态 的状态估计矩阵 和与 对应的协防差矩阵P(s|s)，构成观测更新方程组，其中，SIN
Φ＝Φ ，
H H -1
K(s)＝P(s)(Φ(s))(Φ(s)P(s)(Φ(s))+V[W(s)]) ，
P(s|s)＝P(s)-K(s)Φ(s)P(s)，
(s)
步骤(6.5)，按步骤(5.5)所述方法计算出信道矩阵的估计值H ：
(s)
步骤(6.6)，按步骤(5.6)所述方法对信道矩阵进行QR分解得到R ：
(s) (s)
H ＝IR ，
步骤(6.7)，按步骤(5.7)所述方法对数据进行QR数据检测：
步骤(7)，判断是否所有的数据是否已经都用于迭代，如果是，则算法结束，如果不是，则继续，
步骤(8)，通过判断迭代的次数，决定是否需要增加迭代算法的输入数据，判断算法如下：
设定步长 Δ， 比 较 迭 代 次 数 i 和 Δ·μ+2，其 中，
如果i＝Δ·μ+2，则选择已用于计算的数据所在的子载波两侧
的子载波处的数据作为下次迭代新加入的用于计算的数据，与已用于计算的数据一起带入SIN算法，带入到步骤(6)中重新计；
如果i≠Δ·μ+2，则不增加用于计算的数据，返回步骤(6)进行迭代；
结束。
一种快变OFDM系统的渐进迭代时变信道估计和ICI消除方
法
技术领域
发明内容
Np个导频被平均的插入到N个载波之中并且在传输过程中保持不变，在N个载波中导频点符号表示为
的最高频率，Ts是采样时间，
1,2,…,k,…,K，K≤N，τl是第l径的时延，Gl (M,K)为信道冲击相应的频域表达矩阵，其每一个元素计算如下：
矩阵。(·)表示 Hermitian运算， 其中E[·]
表示均值,J0(·)表示第一类的零阶贝塞尔函数，fd=vfc/c是终端的移动速度为v时的最大多普勒频移，fc是载波频率，c是光速， 代表第l径的信道冲击响应的方差，并假设[0034] 步骤(3.2)，根据Yule-Walker方程得到信道P-BEM参数的状态转移方程：
2
示相邻的各个导频点之间的距离， 是一个Np×Np的单位矩阵，σ 高
(s)
斯白噪声W 的方差，上式中第二项为除导频点子载波外的其它子载波上的数据对导频点所在子载波的干扰ICI，
则SIN估计的Kalman观测方程表示为：
迭代的最边缘数据距对应的各个导频点的距离，第一次计算时为0，rag为计算的精度，rag=[0,1,2,3]，并且:
H H -1
的子载波处的数据作为下次迭代新加入的用于计算的数据，与已用于计算的数据一起带入SIN算法，带入到步骤(6)中重新计；
3（b） 表示第二次计算时的接收信号，其中①表示导频点所在子载波 ②表示本次迭代新加入计算的数据所在的子载波 和 ③表示表示ICI所在的子载
波，④表示α，即用于迭代的最边缘的数据与导频点之间的距离，此时α＝1，⑤表示导频点之间的距离 图3（c）表示迭代次数i＝Δ·μ+2，且i＞2时的接收信号，其中
①表示上次用于计算的数据所在子载波
示本次迭代新加入计算的数据所在的子载波 和 ④表示α，此时，α
＝Δ+1；
具体实施方式
Np个导频被平均的插入到N个载波之中并且在传输过程中保持不变，在N个载波中导频点符号表示为
(s) (s) (s)
的最高频率，Ts是采样时间，
当前OFDM信号的Cl(s)的自相关矩阵，当前OFDM信号的Cl(s)和后一个符号的Cl(s+1)的相关矩阵。(·)H表示Hermitian运算， 其中E[·]
表示均值,J0(·)表示第一类的零阶贝塞尔函数，fd=vfc/c是终端的移动速度为v时的最大多普勒频移，fc是载波频率，c是光速， 代表第l径的信道冲击响应的方差，并假设[0144] 步骤(3.2)，根据YuleWalker方程得到信道P-BEM参数的状态转移方程：
2
示相邻的各个导频点之间的距离， 是一个Np×Np的单位矩阵，σ 高
(s)
斯白噪声W 的方差，上式中第二项为除导频点子载波外的其它子载波上的数据对导频点所在子载波的干扰ICI，
(s)
则SIN估计的Kalman观测方程表示为：
迭代的最边缘数据距对应的各个导频点的距离，第一次计算时为0，rag为计算的精度，rag=[0,1,2,3]，并且:
H H -1
的子载波处的数据作为下次迭代新加入的用于计算的数据，与已用于计算的数据一起带入SIN算法，带入到步骤(6)中重新计；
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 信息处理装置、信息处理方法及程序 | 2020-05-11 | 275 |
| 一种基于稀疏感知的自适应空时均衡方法 | 2020-05-12 | 905 |
| 一种基于心率变异性的团队协同能力的评估模型 | 2020-05-12 | 968 |
| 融合目标检测与跟踪的视频半自动目标标注方法 | 2020-05-13 | 604 |
| 一种通过建立QSAR模型预测化学品以斑马鱼胚胎为受体的毒性的方法 | 2020-05-13 | 350 |
| 一种针对抵近飞行过程的目标跟踪方法 | 2020-05-13 | 865 |
| 基于谱聚类的音频数据聚类方法 | 2020-05-12 | 865 |
| 一种基于空间误差的城市土地利用模拟元胞自动机方法 | 2020-05-08 | 54 |
| 移动通信系统中针对LTE-A系统消息进行解析处理的方法 | 2020-05-11 | 257 |
| NB_IoT系统中针对空口系统消息进行解析处理的方法 | 2020-05-08 | 99 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
