CN109510672A - 一种时变信道参数估计方法 - Google Patents
一种时变信道参数估计方法 Download PDFInfo
- Publication number
- CN109510672A CN109510672A CN201811590965.5A CN201811590965A CN109510672A CN 109510672 A CN109510672 A CN 109510672A CN 201811590965 A CN201811590965 A CN 201811590965A CN 109510672 A CN109510672 A CN 109510672A
- Authority
- CN
- China
- Prior art keywords
- orthogonal
- matrix
- orthogonal subspaces
- impulse response
- subspaces
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/382—Monitoring; Testing of propagation channels for resource allocation, admission control or handover
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Power Engineering (AREA)
- Radio Transmission System (AREA)
Abstract
本发明公开了提出了一种时变信道参数估计方法,包括以下步骤:步骤1:接收端采集多天线上接收到的信号;步骤2:构造正交子空间投影矩阵;步骤3:利用正交子空间投影矩阵,将接收信号投影到相互正交的子空间上;步骤4:在每个正交子空间内进行信道估计,得到每个正交子空间内的信道估计参数;步骤5:对正交子空间内的信道冲激响应估计值进行多普勒补偿;步骤6:将正交子空间中的信道冲激响应估计值进行合并,得到最终信道估计结果。本发明的基于正交子空间投影的时变信道估计方法可获得更小的信道估计均方误差,提升了性能。
Description
技术领域
本发明属于多天线环境下接收领域,具体涉及种基于正交子空间投影的时变信道参数估计方法。
背景技术
随着高速铁路的发展,人们出行越来越方便快捷。目前,铁路建设目标之一为80%以上的大城市都能被高速铁路网所覆盖,将北京到大部分省会城市之间的乘车时间控制在2-8小时之内。这也就意味着,在这个快节奏的社会里,人们将更愿意选择高铁出行。但随着中国进入4G时代,铁路用户对通信的需求日益增加。现今高速铁路所使用的窄带通信系统GSM-R(GSM for Railways)已不能满足用户的宽带通信需求,如视频通话、电话会议、直播等。因此,如何优化高速铁路的无线通信系统,满足客户旅途中的通信体验,成为业内广泛关注的问题。
在OFDM(Orthogonal Frequency Division Multiplexing)系统中,带宽划分为很多频域重叠非频率选择性子信道,对应的载波成为子载波。各个子信道经过调制之后,组成一个多路复用信号。然后在时域上对该信号进行截断,相当于在频域上与sinc函数卷积。由于sinc函数在周围的频载点值为零,使得该系统的子载波之间是正交的,相互独立,互不干扰,这种并行传输的特征使得OFDM 的传输速率得到极大提升。并且由于OFDM符号时域长度远大于多径实验扩展, OFDM具有很好的抗符号间干扰特性,同时具有接收机复杂度低、带宽分配灵活、方便与多种技术结合等优点,因此在无线通信中,OFDM技术是一种广泛使用的宽带空口传输技术,并且是第四代移动通信技术的物理层核心模块。
信道估计是OFDM系统信道均衡、检测、和解码的基础,信道估计性能的好坏直接影响着接收机能否根据当前信道的频域响应来准确地恢复出发射端的发送信号,是衡量OFDM系统性能优劣的主要标准。随着高铁技术的快速发展,快时变环境下的OFDM系统信道估计技术也面临了新的巨大挑战。由于接收端和发送端的相对高速移动,导致多径增益在一个OFDM符号周期内快速变化,信道变为快时变信道。同时,在快时变环境下,多普勒频移增大,破坏了子载波间的正交性,引起子载波间干扰(Inter-Carrier Interference,ICI),使得OFDM 系统性能恶化,降低信道估计精度。因此,需要找到合适的信道估计参数算法,更加准确地恢复出发射信号。
发明内容
本发明针对现有技术的不足,提出了一种时变信道参数估计方法。本发明的基于正交子空间投影的时变信道估计方法可获得更小的信道估计均方误差,提升了性能。
为达到上述目的,本发明采用以下技术方案:
一种时变信道参数估计方法,其特征在于,包括以下步骤:
步骤1:接收端采集多天线上接收到的信号;
步骤2:根据接收到的信号构造正交子空间投影矩阵;
步骤3:利用正交子空间投影矩阵,将接收到的信号投影到相互正交的子空间上;
步骤4:在每个正交子空间内进行信道估计,得到每个正交子空间内的信道冲激响应估计值;
步骤5:对正交子空间内的信道冲激响应估计值进行多普勒补偿;
步骤6:将正交子空间中的信道冲激响应估计值进行合并,得到最终信道估计结果。
进一步的,所述步骤1还包括:接收端采集排列成均匀圆阵的(UCA:UniformCircle Array)的多天线上接收到的信号,接收信号可以表示为:
式中,αl为第l条径的散射系数,al为第l条径的接收阵列响应矩阵,fd为多普勒频率,Ts为样本周期,θl为接收信号入射方向与水平方向的夹角,为时域发送信号,dl为第l条径的相关时延,z(n)为高斯白噪声。
进一步的,所述步骤2还包括:构造正交子空间投影矩阵,首先定义空间滤波器矩阵,
式中,βk为第k个正交子空间的平均角度,nr为正交子空间个数,即接收天线数,ω为载波角频率,r为UCA的半径,c为光速,为接收信号与垂直方向的夹角。接着定义正交子空间投影矩阵,即,
式中,为第k个空间滤波器矩阵。
进一步的,所述步骤3还包括:利用正交子空间投影矩阵将接收信号投影到相互正交子空间内,可以得到投影到正交子空间上的接收信号,
式中,为正交子空间投影矩阵的转置矩阵,y(n)为接收端采集到的多天线上接收的信号。
进一步的,所述步骤4还包括:在每个正交子空间内进行信道估计,首先采用基扩展模型(BEM:Basis Expansion Models),那么第k个正交子空间的接收信号可以表示为:
其中Q为BEM的阶数,δq(n)为BEM基函数,为待估计的BEM系数,为时域发送信号,dl为第l条径的相关时延,z'(n)为正交子空间中的高斯白噪声。通过利用BEM,将大大减少高速移动环境下的信道待估计参数的数量,将估计大量的信道参数转化为对少量的BEM系数进行估计,大大降低算法的复杂度。
进一步的,所述步骤3还包括将投影到每个正交子空间的接收信号表示成矩阵形式,即,
式中,z'为高斯白噪声矩阵,Β为Βq和组成的增广矩阵,即
其中,
Βq=diag[δq(0)δq(1)δq(Nc-1)],
式中,Βq为由δq(n)组成的对角矩阵,δq(n)为BEM基函数;
式中,为由圆周移位组成的L列矩阵,为时域发送信号;
式中,为矩阵的转置,
为待估计的BEM系数。
采用最小二乘法(Least Squares,LS)对BEM系数进行估计,即,
式中,为第k个正交子空间中的BEM系数估计值,Β+为Β的伪逆矩阵,为第k个正交子空间的接收信号。
那么,可以得到每个正交子空间的信道冲激响应估计值,即,
其中,为第k个正交子空间中的BEM系数估计值的转置矩阵,为BEM基函数。
进一步的,所述步骤5还包括:对正交子空间内的信道冲激响应估计值进行多普勒补偿,首先多普勒补偿矩阵可以表示为;
式中,βk为第k个正交子空间的平均角度,fd为多普勒频率,Ts为样本周期。
然后对每个正交子空间内的信道冲激响应估计值进行多普勒补偿,每个子空间内补偿后的信道冲激响应值为:
其中为两个矩阵对应元素相乘(Hadamard product),即两个矩阵对应元素相乘,为第k个正交子空间的信道冲激响应估计值,为第k个正交子空间的多普勒补偿矩阵。
进一步的,步骤6:将正交子空间中的信道冲激响应估计值进行合并,定义一个大小为nr×1的合并矩阵,即
将各个子空间中的信道冲激响应进行合并,得到最终的信道冲激响应估计值,即
式中,为第k个子空间内补偿后的信道冲激响应估计值,Ρ为合并矩阵。
由于采用上述技术方案,本发明的有益效果为:本发明利用正交子空间投影进行信道估计,对子空间中的投影接收信号进行估计,减小了多普勒频移,从而有效的抑制了载波间干扰,提升了系统误码率和均方误差性能。
附图说明
图1是本发明提供的一种时变信道参数估计方法流程图。
具体实施方案
下面结合附图对本发明优选实施例作详细说明。
步骤1:接收端采集排列成UCA(Uniform Circle Array)的多天线上接收到的信号,接收信号可以表示为
式中,αl为第l条径的散射系数,al为第l条径的接收阵列响应矩阵,fd为多普勒频率,Ts为样本周期,θl为接收信号入射方向与水平方向的夹角,为时域发送信号,dl为第l条径的相关时延,z(n)为高斯白噪声。
步骤2:构造正交子空间投影矩阵,首先定义空间滤波器矩阵,
式中,βk为第k个正交子空间的平均角度,nr为正交子空间个数,即接收天线数,ω为载波角频率,r为UCA的半径,c为光速,为接收信号与垂直方向的夹角。接着定义正交子空间投影矩阵,即,
式中,为第k个空间滤波器矩阵。
步骤3:利用正交子空间投影矩阵将接收信号投影到相互正交子空间内,可以得到投影到正交子空间上的接收信号:
步骤4:在每个正交子空间内进行信道估计,首先采用BEM(Basis ExpansionModels),那么第k个正交子空间的接收信号可以表示为:
其中Q为BEM的阶数,δq(n)为BEM基函数,为待估计的BEM系数,为时域发送信号,dl为第l条径的相关时延,z'(n)为正交子空间中的高斯白噪声。通过利用BEM,将大大减少高速移动环境下的信道待估计参数的数量,将估计大量的信道参数转化为对少量的BEM系数进行估计,大大降低算法的复杂度。
将投影到每个正交子空间的接收信号表示成矩阵形式,即,
式中,z'为高斯白噪声矩阵,Β为Βq和组成的增广矩阵,即
其中,
Βq=diag[δq(0)δq(1)δq(Nc-1)],
式中,Βq为由δq(n)组成的对角矩阵,δq(n)为BEM基函数;
式中,为由圆周移位组成的L列矩阵,为时域发送信号;
式中,为矩阵的转置,
为待估计的BEM系数。
采用最小二乘法(Least Squares,LS)对BEM系数进行估计,即,
式中,为第k个正交子空间中的BEM系数估计值,Β+为Β的伪逆矩阵,为第k个正交子空间的接收信号。
那么,可以得到每个正交子空间的信道冲激响应估计值,即,
其中,为第k个正交子空间中的BEM系数估计值的转置矩阵,为BEM基函数。
步骤5:对正交子空间内的信道冲激响应估计值进行多普勒补偿。首先多普勒补偿矩阵可以表示为,
式中,βk为第k个正交子空间的平均角度,fd为多普勒频率,Ts为样本周期。
然后对每个正交子空间内的信道冲激响应估计值进行多普勒补偿,那么每个子空间内补偿后的信道冲激响应值为,
其中为Hadamard product,即两个矩阵对应元素相乘,为第k个正交子空间的信道冲激响应估计值,为第k个正交子空间的多普勒补偿矩阵。
步骤6:将正交子空间中的信道冲激响应估计值进行合并,定义一个大小为 nr×1的合并矩阵,即,
将各个子空间中的信道冲激响应进行合并,得到最终的信道冲激响应估计值,即,
式中,为第k个子空间内补偿后的信道冲激响应估计值,Ρ为合并矩阵。
以上对本发明的优选实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。
Claims (8)
1.一种时变信道参数估计方法,其特征在于,包括以下步骤:
步骤1:接收端采集多天线上接收到的信号;
步骤2:根据接收到的信号构造正交子空间投影矩阵;
步骤3:利用正交子空间投影矩阵,将接收到的信号投影到相互正交的子空间上;
步骤4:在每个正交子空间内进行信道估计,得到每个正交子空间内的信道冲激响应估计值;
步骤5:对正交子空间内的信道冲激响应估计值进行多普勒补偿;
步骤6:将正交子空间中的信道冲激响应估计值进行合并,得到最终信道估计结果。
2.根据权利要求1所述的一种时变信道参数估计方法,其特征在于,所述步骤1还包括:接收端采集排列成均匀圆阵的多天线上接收到的信号,接收信号可以表示为:
式中,αl为第l条径的散射系数,al为第l条径的接收阵列响应矩阵,fd为多普勒频率,Ts为样本周期,θl为接收信号入射方向与水平方向的夹角,为时域发送信号,dl为第l条径的相关时延,z(n)为高斯白噪声。
3.根据权利要求1所述的一种时变信道参数估计方法,其特征在于,所述步骤2还包括:构造正交子空间投影矩阵,首先定义空间滤波器矩阵,
式中,βk为第k个正交子空间的平均角度,nr为正交子空间个数,即接收天线数,ω为载波角频率,r为均匀圆阵的半径,c为光速,为接收信号与垂直方向的夹角。接着定义正交子空间投影矩阵,即,
式中,为第k个空间滤波器矩阵。
4.根据权利要求1所述的一种时变信道参数估计方法,其特征在于,所述步骤3还包括:利用正交子空间投影矩阵将接收信号投影到相互正交子空间内,可以得到投影到正交子空间上的接收信号,
式中,为正交子空间投影矩阵Ur的转置矩阵,y(n)为接收端采集到的多天线上接收的信号。
5.根据权利要求1所述的一种时变信道参数估计方法,其特征在于,所述步骤4还包括:在每个正交子空间内进行信道估计,首先采用基扩展模型,那么第k个正交子空间的接收信号可以表示为:
其中,Q为基扩展模型的阶数,δq(n)为基扩展模型基函数,为待估计的基扩展模型系数,为时域发送信号,dl为第l条径的相关时延,z'(n)为正交子空间中的高斯白噪声。
6.根据权利要求1或4所述的一种时变信道参数估计方法,其特征在于,所述步骤3还包括将投影到每个正交子空间的接收信号表示成矩阵形式,即,
式中,z'为高斯白噪声矩阵,B为Bq和组成的增广矩阵,即,
其中,
式中,Bq为由δq(n)组成的对角矩阵,δq(n)为基扩展模型基函数;
式中,为由圆周移位组成的L列矩阵,为时域发送信号;
式中,为矩阵的转置,
为待估计的基扩展模型系数。
采用最小二乘法对基扩展模型系数进行估计,即,
式中,为第k个正交子空间中的基扩展模型系数估计值,B+为B的伪逆矩阵,为第k个正交子空间的接收信号。
那么,可以得到每个正交子空间的信道冲激响应估计值,即,
其中,为第k个正交子空间中的基扩展模型系数估计值的转置矩阵, 为基扩展模型的基函数。
7.根据权利要求1所述的一种时变信道参数估计方法,其特征在于,所述步骤5还包括:对正交子空间内的信道冲激响应估计值进行多普勒补偿,首先多普勒补偿矩阵可以表示为;
式中,βk为第k个正交子空间的平均角度,fd为多普勒频率,Ts为样本周期;
然后对每个正交子空间内的信道冲激响应估计值进行多普勒补偿,每个子空间内补偿后的信道冲激响应估计值为:
其中为两个矩阵对应元素相乘,为第k个正交子空间的信道冲激响应估计值,为第k个正交子空间的多普勒补偿矩阵。
8.根据权利要求1所述的一种时变信道参数估计方法,其特征在于,步骤6:将正交子空间中的信道冲激响应估计值进行合并,定义一个大小为nr×1的合并矩阵,即
将各个子空间中的信道冲激响应进行合并,得到最终的信道冲激响应估计值,即,
式中,为第k个子空间内补偿后的信道冲激响应估计值,P为合并矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811590965.5A CN109510672B (zh) | 2018-12-21 | 2018-12-21 | 一种时变信道参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811590965.5A CN109510672B (zh) | 2018-12-21 | 2018-12-21 | 一种时变信道参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109510672A true CN109510672A (zh) | 2019-03-22 |
CN109510672B CN109510672B (zh) | 2020-12-15 |
Family
ID=65754578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811590965.5A Active CN109510672B (zh) | 2018-12-21 | 2018-12-21 | 一种时变信道参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109510672B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111600665A (zh) * | 2020-05-12 | 2020-08-28 | 长安大学 | 基于arima滤波器的sage信道参数估计方法 |
CN113014525A (zh) * | 2021-03-15 | 2021-06-22 | 之江实验室 | 一种智能反射表面系统中干扰抑制方法 |
CN115706684A (zh) * | 2021-08-03 | 2023-02-17 | 大唐移动通信设备有限公司 | 信道参数的估计方法、装置及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102938669A (zh) * | 2012-11-20 | 2013-02-20 | 杭州电子科技大学 | 一种高速移动环境下抑制多普勒扩展的方法 |
CN105262526A (zh) * | 2015-10-30 | 2016-01-20 | 哈尔滨工程大学 | 一种基于最小化干扰投影矩阵核范数的干扰对齐方法 |
US20170168988A1 (en) * | 2015-12-09 | 2017-06-15 | The Aerospace Corporation | Signal/noise separation using frft rotational parameter obtained in relation to wigner distribution |
-
2018
- 2018-12-21 CN CN201811590965.5A patent/CN109510672B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102938669A (zh) * | 2012-11-20 | 2013-02-20 | 杭州电子科技大学 | 一种高速移动环境下抑制多普勒扩展的方法 |
CN105262526A (zh) * | 2015-10-30 | 2016-01-20 | 哈尔滨工程大学 | 一种基于最小化干扰投影矩阵核范数的干扰对齐方法 |
US20170168988A1 (en) * | 2015-12-09 | 2017-06-15 | The Aerospace Corporation | Signal/noise separation using frft rotational parameter obtained in relation to wigner distribution |
Non-Patent Citations (4)
Title |
---|
RONG ZENG 等: "Joint Estimation of Frequency Offset and Doppler Shift in High Mobility Environments Based on Orthogonal Angle Domain Subspace Projection", 《IEEE TRANSACTION ON VEHICULAR TECHNOLOGY》 * |
RONG ZENG 等: "Orthogonal Angle Domain Subspace Projection-Based Receiver Algorithm for Underwater Acoustic Communication", 《IEEE COMMUNICATIONS LETTERS》 * |
朱昂昂: "OFDM技术在水声通信中的应用研究", 《中国优秀硕士学位论文全文数据库》 * |
贾慧: "基于投影的MIMO_OFDM系统信道盲估计算法研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111600665A (zh) * | 2020-05-12 | 2020-08-28 | 长安大学 | 基于arima滤波器的sage信道参数估计方法 |
CN113014525A (zh) * | 2021-03-15 | 2021-06-22 | 之江实验室 | 一种智能反射表面系统中干扰抑制方法 |
CN113014525B (zh) * | 2021-03-15 | 2022-06-24 | 之江实验室 | 一种智能反射表面系统中干扰抑制方法 |
CN115706684A (zh) * | 2021-08-03 | 2023-02-17 | 大唐移动通信设备有限公司 | 信道参数的估计方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109510672B (zh) | 2020-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100865469B1 (ko) | 공간 분할 다중 접속 방식을 지원하는 직교주파수 다중 분할 방식 시스템에서 채널 추정 장치 및 방법 | |
JP5755679B2 (ja) | 無線マルチキャリヤ通信システムに関するパイロット伝送方法 | |
CN101778069B (zh) | 一种新型ofdm信道估计联合ici自消除方法 | |
CN102291351B (zh) | 一种ofdm无线通信系统中接收机的定时同步方法 | |
KR100922980B1 (ko) | 다중 안테나를 사용하는 직교주파수분할다중 시스템에서 채널 추정 장치 및 방법 | |
CN102158459B (zh) | 基于时频二维训练信息的ofdm块传输方法 | |
CN106789773B (zh) | 高速移动场景下的频偏和多普勒扩展的确定方法及装置 | |
CN109510672A (zh) | 一种时变信道参数估计方法 | |
CN102291363B (zh) | 一种用于ofdm系统的信道估计及数据检测方法 | |
WO2012093333A1 (en) | Method of channel estimation, method of selecting pilot information, user equipment, and base station | |
CN102291342B (zh) | 一种基于多径解析的ofdm信道估计方法 | |
CN109600327A (zh) | 一种基于虚部干扰利用的信道估计方法 | |
CN103763227A (zh) | 一种基于压缩感知的窄带干扰估计和消除方法及装置 | |
Rodriguez-Pineiro et al. | Experimental validation of ICI-aware OFDM receivers under time-varying conditions | |
CN101340406B (zh) | 多输入多输出正交频分复用系统的信道估计方法 | |
CN108449291B (zh) | 一种成对载波系统中参数的联合估计方法 | |
CN104333527A (zh) | 一种飞机用cofdm调制解调方法 | |
CN103607231B (zh) | 高速移动环境下利用多天线的快速波束切换方法 | |
CN105024956A (zh) | 高速移动下基于统计信息的信道估计和ici消除方法 | |
CN101394390B (zh) | 一种扩频型pdh微波通信系统和方法 | |
CN110213183A (zh) | 一种基于Chu序列的全双工中继系统的信道估计方法 | |
CN106534035A (zh) | 低轨卫星通信系统中的定时同步方法 | |
CN102447658B (zh) | 干扰抑制合并的方法及装置 | |
WO2011112128A1 (en) | Method and apparatus for estimating a doppler frequency | |
Semmar et al. | Digital broadcasting television channel measurements and characterization for SIMO mobile reception |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |