CN103475602B - 一种考虑同频干扰的mimo‑ofdm信道估计方法 - Google Patents
一种考虑同频干扰的mimo‑ofdm信道估计方法 Download PDFInfo
- Publication number
- CN103475602B CN103475602B CN201310418937.6A CN201310418937A CN103475602B CN 103475602 B CN103475602 B CN 103475602B CN 201310418937 A CN201310418937 A CN 201310418937A CN 103475602 B CN103475602 B CN 103475602B
- Authority
- CN
- China
- Prior art keywords
- channel
- sigma
- interference
- mimo
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明公开了一种考虑同频干扰的MIMO‑OFDM信道估计方法,包括以下步骤:1)利用MIMO‑OFDM导频信号的观测样本,基于最小二乘法的估计准则,获得信道冲激响应的初始估计;2)利用信道冲激响应的初始估计,基于残差估计方法得到干扰协方差矩阵的初始估计值;3)利用干扰协方差矩阵的初始估计值,基于Cholesky分解和时域“低通滤波”,获得干扰协方差矩阵的改善估计值;4)利用干扰协方差矩阵的改善估计值,基于最大似然估计准则获得信道冲激响应和信道频率响应的改善估计值。本发明方法计算复杂度低,估计误差性能好,且随着天线数的增多,性能增益也越大。
Description
技术领域
本发明涉及MIMO-OFDM无线通信系统,特别涉及考虑同频干扰的MIMO-OFDM系统中的信道估计方法。
背景技术
在频谱资源日趋紧张的情况下,进一步提升频谱效率是无线通信系统设计的重要目标之一。采用多天线发送和多天线接收(MIMO)的无线传输技术,可以在不额外增加时间与频率资源的情况下,通过挖掘空间维度资源,提高无线传输的有效性与可靠性,最终达到提升频谱效率的目的。正交频分复用(OFDM)技术将高速数据流转换成一组并行传输的数据流,并使子信道的频带交叠,具有较强的抗多径干扰能力和较高的频带利用率。MIMO-OFDM技术在过去十余年来一直是无线通信领域研究的主流技术之一。
在移动通信系统中,对信道进行准确的估计是实现相干解调的前提,是决定移动通信系统可靠性和有效性的重要因素之一。在未来移动通信系统中,不同小区间进行频率复用以提高频谱效率,但同时也带来了同信道干扰。传统的信道估计方法是将同信道干扰当作白噪声处理,这在同信道干扰较为严重的场景(如小区边缘)中会带来严重的性能瓶颈。因此在存在同信道干扰的情况下,研究MIMO-OFDM系统中信道估计技术是非常有意义的。
发明内容
技术问题:本发明提供一种低复杂度、高精度的考虑同频干扰的MIMO-OFDM信道估计方法。
技术方案:本发明的一种考虑同频干扰的MIMO-OFDM信道估计方法,包括以下步骤:
1)利用MIMO-OFDM导频信号的观测样本,基于最小二乘法估计准则,获得导频位置的信道冲激响应初始估计值;
2)利用步骤1)获得的信道冲激响应的初始估计,基于残差估计方法计算干扰协方差矩阵的初始估计值;
3)利用步骤2)计算得到的干扰协方差矩阵的初始估计值,基于Cholesky分解和时域“低通滤波”,获得干扰协方差矩阵的改善估计值;
4)利用步骤3)获得的干扰协方差矩阵的改善估计值,基于最大似然估计准则获得信道冲激响应的改善估计值和信道频率响应的改善估计值。
本发明的优选方案中,步骤1)的详细流程为:
首先在接收端获取MIMO-OFDM信号在Np个导频符号的Nc个子载波上的NcNp点观测样本{yk,n},0≤k≤Nc-1,0≤n≤Np-1,其中每个观测样本yk,n都是Nr×1的列向量,表示在第n个导频符号的第k个子载波,Nr根接收天线上的接收信号;由收发双方协商可知导频位置发送的导频信号为{xk,n},0≤k≤Nc-1,0≤n≤Np-1,其中每个导频信号xk,n都是Nt×1的列向量,表示在第n个导频符号的第k个子载波,Nt根发送天线上发送的导频信号;其中参数:导频符号个数Np、子载波个数Nc、发送天线根数Nt、接收天线根数Nr在系统中均为确定已知的;
然后按照下式计算导频信号xk,n的扩展形式Gk,n:
其中表示Nr×Nr单位阵,表示Kronecker乘积,
最后按照下式计算MIMO-OFDM系统中导频位置的信道冲激响应的初始估计值为:
其中(·)H表示共轭转置,则表示对导频信号xk,n的扩展形式Gk,n的共轭转置;(·)-1表示求逆。
本发明的优选方案中,步骤2)的具体流程为:首先利用步骤1)中得到的导频位置的信道冲激响应的初始估计值进行求卷积运算,得到导频信号通过导频位置的信道冲激响应初始估计值所描述的信道到达接收端的信号,然后将到达接收端的信号,减去接收天线收到的观测样本{yk,n},0≤k≤Nc-1,0≤n≤Np-1,即得到干扰信号的初始估计值;最后将干扰信号的初始估计值与其共轭转置的乘积在Np个导频位置求平均,从而获得干扰协方差矩阵的初始估计值。
本发明的优选方案中,步骤3)的具体流程为:
31)分别对Nc个子载波上干扰协方差矩阵的初始估计值进行Cholesky分解,得到Nr×Nr上三角矩阵0≤k≤Nc-1
32)取出上三角矩阵中第i行第j列元素联立向量 然后对向量进行“低通滤波”,得到除去噪声影响的向量vij;
33)将除去噪声影响的向量vij的第k个元素作为第i行第j列元素置换掉上三角矩阵中原来的元素,得到重构的上三角矩阵1≤i≤j≤Nr,0≤k≤Nc-1;
34)利用步骤33)中重构的上三角矩阵和公式得到第k子载波上干扰协方差矩阵的改善估计值
本发明的优选方案中,步骤4)中,利用下式获得信道冲激响应的改善估计值:
同时,根据与的关系 利用得到对应的然后根据
计算出第k个子载波上的信道频率响应的改善估计值
本发明利用MIMO-OFDM导频信号的观测样本,估计出导频位置的信道冲激响应。估计过程主要通过改善后的干扰协方差矩阵估计值来获得信道冲激响应和信道频率响应的改善估计值。
有益效果:本发明与现有技术相比,具有如下优点:
1、基于最大似然估计方法估计信道冲激响应和信道频率响应,传统估计方法需要迭代计算信道冲激响应和干扰的协方差矩阵,计算量大。本方法通过估计信道冲激响应和干扰的协方差矩阵改善值,不采用迭代的方法,使得计算复杂度大幅降低;
2、本方法将改善后的协方差矩阵估计值再返回估计信道冲激响应,使得信道冲激响应的估计得到改善,使得这种方法的估计误差相对传统信道估计方法性能好;
3、实验结果说明本估计方法随着天线数的增多,性能增益也越大。
附图说明
图1为本发明方法的流程图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
本发明的考虑同频干扰的MIMO-OFDM信道估计方法主要指在MIMO-OFDM系统的接收端,利用接收导频位置的信号样本估计出MIMO-OFDM信道冲激响应和信道频率响应。下面以如下场景为例对本发明所阐述的方法进行说明:MIMO-OFDM系统中可用于导频估计的子载波数为64,采用QPSK调制方式,Turbo编码,码率为1/2,导频个数为2,发送天线和接收天线个数均为2,天线间信道相关矩阵均为单位阵,干扰个数为1,期望信号的传输径数为4,干扰的传输径数为6。该方法具体如下:
1)首先在接收端获取MIMO-OFDM信号Np=2个导频符号的Nc=64个子载波上的NcNp=128点观测样本{yk,n},0≤k≤Nc-1=63,0≤n≤Np-1=1,其中每个观测样本yk,n都是Nr×1=2×1的列向量,表示在第n个导频符号的第k个子载波,Nr=2根接收天线上的接收信号,观测样本为128个2×1维列向量;由收发双方协商可知导频位置发送的导频信号为{xk,n),0≤k≤Nc-1=63,0≤n≤Np-1=1,其中每个导频信号xk,n都是Nt×1=2×1的列向量,表示在第n个导频符号的第k个子载波,Nt=2根发送天线上的发送的导频信号,因此已知的导频数据也有128个2×1维列向量。
然后按照下式计算导频信号分别计算0≤k≤Nc-1=63,0≤n≤Np-1=1范围内,xk,n的扩展形式Gk,n:
其中表示Nr×Nr=2×2单位阵,表示Kronecker乘积,Kronecker乘积的定义如下:
如果A是一个m×n的矩阵,而B是一个p×q的矩阵,克罗内克积则是一个mp×nq的分块矩阵
最后按照下式计算MIMO-OFDM系统中导频位置的信道冲激响应的初始估计值为:
其中(·)H表示共轭转置,则表示对Gk,n的共轭转置;(·)-1表示求逆。
2)利用步骤1)中得到的导频位置的信道冲激响应的初始估计值,基于残差估计方法计算干扰协方差矩阵的初始估计值,具体步骤为:
首先利用步骤1)中得到的导频位置的信道冲激响应的初始估计值进行求卷积运算,得到导频信号通过导频位置的信道冲激响应初始估计值所描述的信道到达接收端的信号,然后将到达接收端的信号,减去接收天线收到的观测样本{yk,n},0≤k≤Nc-1=63,0≤n≤Np-1=1,即得到干扰信号的初始估计值;最后将干扰信号的初始估计值与其共轭转置的乘积在Np=2个导频位置求平均,从而获得干扰协方差矩阵的初始估计值
从而得到在64个子载波上分别对应的干扰协方差矩阵的初始估计值。
3)利用步骤2)中得到的干扰协方差矩阵的初始估计值基于Cholesky分解和时域“低通滤波”,获得干扰协方差矩阵的改善估计值,具体流程如下:
31)分别对Nc=64个子载波上干扰协方差矩阵的初始估计值 进行Cholesky分解: 其中是对角元素为正数的Nr×Nr=2×2上三角矩阵,从而得到64个子载波上分别对应的上三角矩阵0≤k≤NC-1=63;
32)取出上三角矩阵中第i行第j列元素联立向量 然后对向量进行“低通滤波” 即将向量中第6个到60个元素置0,得到除去噪声影响的向量 其中,P=Fdiag(e)FH,e的第k个兀系 F表示Nc=64阶归一化傅里叶变换矩阵,Lu=6为干扰信号传输路径的最大长度;
33)将除去噪声影响的向量vij的第k个元素作为第i行第j列元素置换掉上三角矩阵中原来的元素,得到重构的上三角矩阵
vij,k表示vij的第k个元素。
本实施例中,即为重构上三角矩阵64个子载波上分别对应的
34)利用步骤33)中重构的上三角矩阵和公式得到第k子载波上干扰协方差矩阵的改善估计值从而得到64个子载波上分别对应的干扰协方差矩阵的改善估计值0≤k≤Nc-1=63。
4)利用步骤3)中获得的干扰协方差矩阵的改善估计值,基于最大似然估计准则获得信道冲激响应的改善估计值,表达式为:
根据与的关系 利用得到对应的即每取出的2行元素生成对应的然后根据
计算出第k个子载波上的信道频率响应的改善估计值从而得到到64个子载波上分别对应的信道频率响应的改善估计值0≤k≤Nc-1=63。
Claims (3)
1.一种考虑同频干扰的MIMO-OFDM信道估计方法,其特征在于,该方法包括以下步骤:
1)利用MIMO-OFDM导频信号的观测样本,基于最小二乘法估计准则,获得导频位置的信道冲激响应初始估计值,具体流程为:
首先在接收端获取MIMO-OFDM信号在Np个导频符号的Nc个子载波上的NcNp点观测样本{yk,n},0≤k≤Nc-1,0≤n≤Np-1,其中每个观测样本yk,n都是Nr×1的列向量,表示在第n个导频符号的第k个子载波,Nr根接收天线上的接收信号;由收发双方协商可知导频位置发送的导频信号为{xk,n},0≤k≤Nc-1,0≤n≤Np-1,其中每个导频信号xk,n都是Nt×1的列向量,表示在第n个导频符号的第k个子载波,Nt根发送天线上发送的导频信号;其中参数:导频符号个数Np、子载波个数Nc、发送天线根数Nt、接收天线根数Nr在系统中均为确定已知的;
然后按照下式计算导频信号xk,n的扩展形式Gk,n:
其中表示Nr×Nr单位阵,表示Kronecker乘积,表示傅里叶变换系数,(·)T表示转置,L为有用信号传输路径的最大长度;
最后按照下式计算MIMO-OFDM系统中导频位置的信道冲激响应初始估计值为:
其中(·)H表示共轭转置,则表示对导频信号xk,n的扩展形式Gk,n的共轭转置,(·)-1表示求逆;
2)利用所述步骤1)获得的信道冲激响应的初始估计,基于残差估计方法计算干扰协方差矩阵的初始估计值;
3)利用所述步骤2)计算得到的干扰协方差矩阵的初始估计值,基于Cholesky分解和时域“低通滤波”,获得干扰协方差矩阵的改善估计值,具体流程为:
31)分别对Nc个子载波上干扰协方差矩阵的初始估计值进行Cholesky分解,得到Nr×Nr上三角矩阵其中Nr为接收天线根数;
32)取出上三角矩阵中第i行第j列元素联立向量然后对向量进行“低通滤波”,得到除去噪声影响的向量vij;
33)将所述除去噪声影响的向量vij的第k个元素作为第i行第j列元素置换掉上三角矩阵中原来的元素,得到重构的上三角矩阵1≤i≤j≤Nr,0≤k≤Nc-1;
34)利用所述步骤33)中重构的上三角矩阵和公式得到第k子载波上干扰协方差矩阵的改善估计值
4)利用所述步骤3)获得的干扰协方差矩阵的改善估计值,基于最大似然估计准则获得信道冲激响应的改善估计值和信道频率响应的改善估计值。
2.根据权利要求1所述的考虑同频干扰的MIMO-OFDM信道估计方法,其特征在于,所述步骤2)的具体流程为:首先利用步骤1)中得到的导频位置的信道冲激响应的初始估计值进行求卷积运算,得到导频信号通过导频位置的信道冲激响应初始估计值所描述的信道到达接收端的信号,然后将所述到达接收端的信号,减去接收天线收到的观测样本{yk,n},0≤k≤Nc-1,0≤n≤Np-1,即得到干扰信号的初始估计值;最后将干扰信号的初始估计值与其共轭转置的乘积在Np个导频位置求平均,从而获得干扰协方差矩阵的初始估计值,其中Nc为子载波个数。
3.根据权利要求1所述的考虑同频干扰的MIMO-OFDM信道估计方法,其特征在于,所述步骤4)中,利用下式获得信道冲激响应的改善估计值:
同时,根据与的关系利用得到对应的然后根据
计算出第k个子载波上的信道频率响应的改善估计值其中Nc为子载波个数,Np为导频符号个数,为第k子载波上干扰协方差矩阵的改善估计值,Gk,n为导频信号xk,n的扩展形式,yk,n为MIMO-OFDM信号在Np个导频符号的Nc个子载波上的NcNp点观测样本,L为有用信号传输路径的最大长度,Fk表示傅里叶变换系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310418937.6A CN103475602B (zh) | 2013-09-16 | 2013-09-16 | 一种考虑同频干扰的mimo‑ofdm信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310418937.6A CN103475602B (zh) | 2013-09-16 | 2013-09-16 | 一种考虑同频干扰的mimo‑ofdm信道估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103475602A CN103475602A (zh) | 2013-12-25 |
CN103475602B true CN103475602B (zh) | 2017-01-25 |
Family
ID=49800305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310418937.6A Expired - Fee Related CN103475602B (zh) | 2013-09-16 | 2013-09-16 | 一种考虑同频干扰的mimo‑ofdm信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103475602B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104022977B (zh) * | 2014-06-17 | 2017-03-08 | 东南大学 | 一种适用于mimo‑ofdm系统的信道矩阵与干扰协方差矩阵估计方法 |
CN104022979B (zh) * | 2014-06-23 | 2017-07-18 | 东南大学 | 一种联合稀疏信道估计方法、装置及系统 |
CN105337907B (zh) * | 2014-07-25 | 2019-02-26 | 华为技术有限公司 | 一种获取信道状态信息的方法和装置 |
CN105847197B (zh) * | 2015-01-16 | 2020-03-17 | 中兴通讯股份有限公司 | 一种获取发送信号的方法和装置 |
CN106330793B (zh) * | 2015-06-26 | 2019-12-06 | 华为技术有限公司 | 信号处理方法及装置 |
CN106788628B (zh) * | 2016-06-30 | 2020-07-28 | 北京紫光展锐通信技术有限公司 | Mimo系统的信道估计方法及装置、基站 |
CN108400947B (zh) * | 2017-02-08 | 2020-08-14 | 华为技术有限公司 | 干扰噪声协方差矩阵估计方法、装置及系统 |
CN108599817B (zh) * | 2018-03-26 | 2020-09-15 | 佛山市顺德区中山大学研究院 | 一种多小区mimo-ofdm同频干扰信道估计方法 |
CN108881217B (zh) * | 2018-06-14 | 2020-10-27 | 西安交通大学 | 一种基于分层二维特征编码的安全多用户导频鉴权方法 |
CN115865109B (zh) * | 2022-11-11 | 2024-05-14 | 北京智芯微电子科技有限公司 | 多接收天线的干扰抑制合并方法、装置及介质、接收终端 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1454019A (zh) * | 2002-04-23 | 2003-11-05 | 西门子移动通讯公司 | 数字无线通信系统中信道估计的方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080298480A1 (en) * | 2007-05-31 | 2008-12-04 | Qiang Li | Asynchronous MIMO-OFDM spatial covariance estimation |
-
2013
- 2013-09-16 CN CN201310418937.6A patent/CN103475602B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1454019A (zh) * | 2002-04-23 | 2003-11-05 | 西门子移动通讯公司 | 数字无线通信系统中信道估计的方法和装置 |
Non-Patent Citations (3)
Title |
---|
Efficient Spatial Covariance Estimation for Asynchronous Co-channel Interference Suppression in MIMO-OFDM Systems;Qiang Li, Jing Zhu, Qinghua Li, and C. N. Georghiades;《IEEE transaction on wireless communication》;20081231;第7卷(第12期);第4849-4853页 * |
OFDM系统的信道估计和信号均衡技术的研究;樊同亮;《中国博士学位论文全文数据库 信息科技辑》;20130215(第2期);第56页 * |
TD-SCDMA中信道估计实现及其性能分析;赵洪昌,宋常建,张水莲;《现代电子技术》;20051231(第5期);第68-69页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103475602A (zh) | 2013-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103475602B (zh) | 一种考虑同频干扰的mimo‑ofdm信道估计方法 | |
CN101494528B (zh) | 发射分集块传输系统的训练序列设计及其信道估计方法 | |
CN101692665B (zh) | 正交频分复用-多输入多输出系统的解调方法及解调器 | |
CN104618061B (zh) | 一种大规模多天线系统中多用户信号的检测方法 | |
CN101702696B (zh) | 信道估计的实现方法和装置 | |
CN104393963B (zh) | 基于循环平稳特性的空时分组码mc‑cdma信号盲识别方法 | |
CN103685096B (zh) | 一种基于最优导频的mimo-ofdm系统信道估计方法 | |
CN104410590A (zh) | 一种基于压缩感知的短波ofdm抑制干扰联合信道估计方法 | |
CN109600327B (zh) | 一种基于虚部干扰利用的信道估计方法 | |
CN103812807B (zh) | 增强余弦调制滤波器组多载波通信系统抗干扰能力的方法 | |
CN103873406B (zh) | 水声正交频分复用通信系统帧间干扰消除方法 | |
CN104767587A (zh) | 基于ofdm系统下联合信道编译码的压缩感知信道估计方法 | |
CN104320369A (zh) | 一种基于信道估计误差和数据检测误差的迭代方法 | |
CN110213185A (zh) | 一种基于原子范数最小化的三维信道参数估计方法 | |
CN101309243A (zh) | 一种新的ofdm参数化信道估计器 | |
CN104539562A (zh) | 基于多输入多输出正交频分复用的宽带短波信道估计方法 | |
CN106209716A (zh) | 一种降低大规模mu‑mimo‑ofdm系统峰均功率比的方法 | |
CN101848183B (zh) | 多入多出正交频分复用系统中信道估计方法及装置 | |
CN101155164A (zh) | 一种dft扩频的广义多载波系统的sinr估计方法 | |
CN107276726A (zh) | 一种Massive MIMO FBMC波束空时编码下行传输方法 | |
CN102025459B (zh) | 基于非参量估计ica的mimo-ofdm系统盲去卷积方法 | |
CN105812299A (zh) | 基于联合块稀疏重构的无线传感网信道估计算法及系统 | |
CN105610484A (zh) | 大规模mimo低复杂度迭代接收方法 | |
CN104702540A (zh) | 一种用于小区边缘终端的信号接收优化方法 | |
CN102685060B (zh) | 一种正交频分复用系统中多用户mimo接收方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 210093 Nanjing University Science Park, 22 Hankou Road, Gulou District, Nanjing City, Jiangsu Province Patentee after: SOUTHEAST University Address before: 211103 No. 5 Runfa Road, Jiangning District, Nanjing City, Jiangsu Province Patentee before: Southeast University |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170125 Termination date: 20210916 |