CN102723978A - 多天线信道特征参数联合估计的下界建立方法 - Google Patents
多天线信道特征参数联合估计的下界建立方法 Download PDFInfo
- Publication number
- CN102723978A CN102723978A CN2012102251611A CN201210225161A CN102723978A CN 102723978 A CN102723978 A CN 102723978A CN 2012102251611 A CN2012102251611 A CN 2012102251611A CN 201210225161 A CN201210225161 A CN 201210225161A CN 102723978 A CN102723978 A CN 102723978A
- Authority
- CN
- China
- Prior art keywords
- partiald
- dtri
- antenna
- channel
- log
- 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
Images
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明多天线信道特征参数联合估计的下界建立方法,涉及无线通信和参数估计技术领域。在对多天线信道采用估计方法获得包括信道衰落、时延和频率偏置这三类时变信道特征参数的估计值和估计方差时,可用本发明的下界作为评价指标。包括以下步骤:初始化待估计参数及其估计方差;根据信道特征参数模型推导一步状态转移概率密度函数和观测似然函数;然后利用递推方式得到费希尔信息矩阵,对其求逆后得到克拉美-罗下界。
Description
技术领域
本发明涉及无线通信系统的接收机设计技术领域,具体是针对目前时变的多天线信道,结合非线性滤波理论中用递推方式得到的克拉美-罗下界(CRLB)公式,提供的信道时变参数联合获取时的估计下界,能给接收机设计提供参考。
背景技术
在多天线系统即多输入多输出系统(MIMO,Multiple InputMultiple Output)的接收端获得信道状态信息(CSI)进而正确地解调出发射信号,可使接收机性能获得较大提升。特别是采用分布式天线配置后,可以克服集中式多天线系统面临的天线间隔大,较难在接收端布置的问题。由于发射机和接收机之间的传播路径非常复杂,多天线接收机在性能上很大程度还是受到无线信道的影响,如阴影衰落,多普勒频移和时延。
无线信道参数估计下界是信道估计中重要的一环,它的主要功能是给各信道估计算法提供性能参考,也为给定指标能否实现提供依据。参数估计下界主要有克拉美-罗下界、巴塔卡里亚界和巴兰金界,由于巴塔卡里亚界和巴兰金界计算过于复杂,使学术界大多研究无偏估计时参数的最小均方误差下界,即克拉美-罗下界。
目前针对集中式MIMO系统时不变特征参数获取时的估计下界研究较多,但是多天线无线信道是典型的时变信道,针对时变特征参数的估计下界还很少,还缺乏多天线信道包括衰落、时延和由多普勒频移引起的频率偏置这三类时变特征参数的估计下界。
发明内容
本发明的目的在于提供一种多天线信道特征参数联合估计的下界建立方法。它利用公知的用递推方法得到的CRLB公式,联系多天线时变信道的具体个性,将时延和信道衰落用直角坐标方式表示为复信道衰落,时延等效为复信道衰落的相位,序贯地获得复信道衰落和频率偏置的下界。该方法对高精度接收机的设计具有参考价值。
本发明是通过以下技术方案实现的:
A、记待估计的NT根发送天线和某一根接收天线间信道的特征参数为xk=[Re(hk)T,Im(hk)T,ωk T]T,其中为信道衰落的实部,为信道衰落的虚部,为频率偏置,下标k表示离散时刻,上标T表示转置,初始化参数值为x0=[Re(h0)T,Im(h0)T,γ0 T]T,根据先验概率密度函数p(x0),初始化参数的方差矩阵值为初始化费希尔信息矩阵的值为
B、由多天线信道参数的随机模型xk+1=fk(xk)+nk,式中,fk(□)为时变函数,nk为过程噪声。由于信道衰落的幅值服从瑞利分布,相位服从(0,2π]间的均匀分布,对于慢衰落,当接收机带宽B和符号周期Ts一定时,瑞利衰落的中值场强只产生比较平缓的变化,进一步考虑信道频偏缓慢变化时,可将信道时变特征参数参数化为在均值附近抖动的一阶自回归模型xk+1=Akxk+nk,式中系数矩阵Ak由多天线信道模型确定,过程噪声nk服从均值uk是0及方差矩阵是Qk的正态分布,所述的p(xk+1|xk)为其中|·|表示行列式;
C、由某一根接收天线接收到的信号 式中,Wp为接收端已知的第p根发送天线的发射功率,hp,k为第p根发送天线和某一根接收天线间在时刻k的复信道衰落,sp,k为时刻k在第p根发送天线上的发送符号,γp,k=fp,kTS为归一化后的载波频率偏置,其中,fp,k为时刻k的绝对实际频偏,0<γp,k<0.5,ηk为均值是0及实部和虚部的方差均是的复高斯观测噪声,根据ηk的分布特征推导出以xk+1为条件的观测量yk+1的观测似然函数p(yk+1|xk+1)为其中
D、由步骤B得到的p(xk+1|xk)和步骤C得到的p(yk+1|xk+1),采用公知的用递推方式获得费希尔信息矩阵的公式: 式中:
E、求取克拉美-罗下界Ck+1为随时刻k重复步骤B、C、D、E,上面所述的克拉美-罗下界Ck+1为某一根接收天线可获得的估计下界,适用于多天线系统中各个接收天线的估计下界获取。
本发明将适用于非线性滤波的具有递推形式的CRLB公式与具有时变特征参数的多天线信道模型相结合,相对于非时变特征参数的CRLB,本方法具有更广的适用性;其次,本方法中的费希尔信息矩阵为递推形式,从而具有高效的实时性;最后,本方法在观测时长较短的范围内就能达到收敛效果,具有良好的稳健性,且较易实现。因此,本发明适合在实际系统中应用,例如各类集中式MIMO系统、分布式MIMO系统(如协同通信系统,利用中继传输的通信系统等)。
附图说明
图1为具有时变特征参数的分布式多天线系统信道模型;
图2为随着观测值个数增加,信道衰落实部和信道衰落虚部的CRLB变化曲线;
图3为随着观测值个数增加,信道频率偏置的CRLB变化曲线;
图4为不同SNR下,信道衰落实部和信道衰落虚部的CRLB变化曲线
图5为不同SNR下,信道频率偏置的CRLB曲线
图6为不同过程噪声方差对信道衰落实部和信道衰落虚部CRLB的影响;
图7为不同过程噪声方差对信道频率偏置CRLB的影响。
具体实施方式
以下结合附图和实施例对本发明作进一步的描述
该实施例以本发明技术方案为前提进行实施,给出了具体的实施方式和具体的操作过程,但本发明的保护范围不限于下述实例。
实施例中,具有时变特征参数的多天线系统信道模型(如附图1所示),本发明采用2发2收的多天线系统,待传输的信号为由哈达码矩阵随机产生S=[S1,S2],采用哈达码矩阵保证了S1和S2在任意时刻的相互正交性,发射信号的总功率为W,W=10^(SNR/10),各天线的发射功率采用平均分配方案。
如图1所示,具有时变特征参数的多天线系统信道模型的特征参数联合获取时下界的建立方法,包括以下步骤:
A、记待估计的2根发送天线和某一根接收天线问信道的特征参数为xk=[Re(h1,k),Re(h2,k).Im(h1,k).Im(h2,k),γ1,k,γ2,k]T,初始化参数值为x0=[Re(h0)T,Im(h0)T,γ0 T]T=[0.1,0.3,0.6,0.4,0.25,0.25]T,令根据先验概率密度函数p(x0),初始化参数的方差矩阵值为 初始化费希尔信息矩阵的值为 xk为矩阵Xk=[x1,k,x2,k]T中任一元素,表示任一根接收天线与2根发送天线之间的信道特征参数;
B、由于信道衰落的幅值服从瑞利分布,相位服从(0,2π]间的均匀分布,对于慢衰落,当接收机带宽B和符号周期Ts一定时,瑞利衰落的中值场强只产生比较平缓的变化,进一步考虑信道频偏缓慢变化时,可将信道时变特征参数参数化为在均值附近抖动的一阶自回归模型xk+1=Akxk+nk,式中系数矩阵Ak由多天线信道模型确定,令Ak为适合维数的单位矩阵,过程噪声nk服从均值uk是0及方差矩阵是Qk=diag{10-5,10-5,10-5,10-5,10-5,10-5}的正态分布,所述的p(xk+1|xk)为 其中|·|表示行列式;
C、第一根接收天线端的观测信号式中,Wp为接收端己知的第p根发送天线的发射功率,hp,k为第p根发送天线和某一根接收天线问在时刻k的复信道衰落,sp,k为时刻k在第p根发送天线上的发送符号,γp,k=fp,kTs为归一化后的载波频率偏置,其中,fp,k为时刻k的绝对实际频偏,Ts为符号周期,0<γp,k<0.25,ηk为均值是0及实部和虚部的方差均是1/2的复高斯观测噪声,根据ηk的分布特征推导出以xk+1为条件的观测量yk+1的观测似然函数p(yk+1|xk+1) 为 其中
D、由步骤B得到的p(xk+1|xk)和步骤C得到的p(yk+1|xk+1),采用公知的用递推方式获得费希尔信息矩阵的公式: 式中:
E、按时刻k更新Jk+1和Ck+1,上面所述的克拉美-罗下界Ck+1为第一根接收天线可获得的估计下界,适用于多天线系统中各个接收天线的估计下界获取。。
本发明中,附图2和附图3所示为信噪比SNR=0dB且过程噪声的方差 时,在接收端第一根接收天线可获得的信道衰落实部、虚部和频率偏置的CRLB曲线,由附图可见,在观测值的个数N<20,就可以接近稳态。仿真还表明当信噪比较高时,在N=10就能到达稳态。
附图4和附图5是在观测噪声方差为Qk=diag{10-5,10-5,10-5,10-5,10-5,10-5}时的CRLB与信噪比的关系曲线,该图表明随着SNR增大,参数的CRLB均减小,但减小的幅度有所减缓。
附图6和附图7为SNR=20dB时,在过程噪声方差为中每个元素分别10-8时、10-6时和10-4时的CRLB曲线,可以得到随着过程噪声方差的增加,其CRLB随着观测值个数均能到达稳态。但是当噪声方差为Qk中的每个元素大于10-4时CRLB随着观测值个数迅速增加,且呈现发散状态。
尽管本发明的内容已经通过上述实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (2)
1.一种多天线信道特征参数联合估计的下界建立方法,其特征是运用递推方式求取多天线信道特征参数联合估计的费希尔信息矩阵,然后求逆得到下界,包括:
A、记待估计的NT根发送天线和某一根接收天线间信道的特征参数为xk=[Re(hk)T,Im(hk)T,γk T]T,其中为信道衰落的实部,为信道衰落的虚部,为频率偏置,下标k表示离散时刻,上标T表示转置,初始化参数值为x0=[Re(h0)T,Im(h0)T,γ0 T]T,初始化参数的方差矩阵值为初始化费希尔信息矩阵的值为J0;
B、由多天线信道参数的随机模型xk+1=fk(xk)+nk,式中,fk(□)为时变函数,nk为过程噪声,根据nk的分布特征推导出以xk为条件的一步递推量xk+1的状态转移概率密度函数p(xk+1|xk);
C、由某一根接收天线接收到的信号 式中,Wp为接收端已知的第p根发送天线的发射功率,hp,k为第p根发送天线和某一根接收天线间在时刻k的复信道衰落,Sp,k为时刻k在第p根发送天线上的发送符号,γp,k=fp,kTs为归一化后的载波频率偏置,其中,fp,k为时刻k的绝对实际频偏,Ts为符号周期,0<γp,k<0.5,ηk为均值是0及实部和虚部的方差均是的复高斯观测噪声,根据ηk的分布特征推导出以xk+1为条件的观测量yk+1的观测似然函数p(yk+1|xk+1);
D、由步骤B得到的p(xk+1|xk)和步骤C得到的p(yk+1|xk+1),采用公知的用递推方式获得费希尔信息矩阵的公式: 式中:
E、求取克拉美-罗下界Ck+1为随时刻k重复步骤B、C、D、E。
2.如权利要求1所述的多天线信道特征参数联合估计的下界建立方法,其特征是运用递推方式求取多天线信道特征参数联合估计的费希尔信息矩阵,然后求逆得到下界,还包括:
步骤C中所述观测似然函数p(yk+1|xk+1)为其中
步骤E中所述克拉美-罗下界Ck+1为某一根接收天线可获得的估计下界。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210225161.1A CN102723978B (zh) | 2012-06-29 | 2012-06-29 | 多天线信道特征参数联合估计的下界建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210225161.1A CN102723978B (zh) | 2012-06-29 | 2012-06-29 | 多天线信道特征参数联合估计的下界建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102723978A true CN102723978A (zh) | 2012-10-10 |
CN102723978B CN102723978B (zh) | 2014-12-10 |
Family
ID=46949662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210225161.1A Expired - Fee Related CN102723978B (zh) | 2012-06-29 | 2012-06-29 | 多天线信道特征参数联合估计的下界建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102723978B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103499811A (zh) * | 2013-08-28 | 2014-01-08 | 电子科技大学 | 提高雷达目标估计性能的天线个数分配方法 |
CN108667152A (zh) * | 2018-06-22 | 2018-10-16 | 南京大学 | 一种鲁棒的无线能量传输方法 |
CN108957422A (zh) * | 2018-06-01 | 2018-12-07 | 电子科技大学 | 一种基于量化数据的云mimo雷达的根克拉美罗下界计算方法 |
CN111585645A (zh) * | 2019-02-18 | 2020-08-25 | 富士通株式会社 | 光链路导致的偏振态变化监测装置及方法、光接收机 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222458A (zh) * | 2008-01-22 | 2008-07-16 | 上海师范大学 | Mimo-ofdm信道的低阶递推最小均方误差估计 |
CN101697495A (zh) * | 2009-11-05 | 2010-04-21 | 湖南大学 | 一种基于博弈论的mimo信道跟踪方法 |
WO2010140101A1 (en) * | 2009-06-04 | 2010-12-09 | Telefonaktiebolaget L M Ericsson (Publ) | Improved mobile radio channel estimation |
-
2012
- 2012-06-29 CN CN201210225161.1A patent/CN102723978B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222458A (zh) * | 2008-01-22 | 2008-07-16 | 上海师范大学 | Mimo-ofdm信道的低阶递推最小均方误差估计 |
WO2010140101A1 (en) * | 2009-06-04 | 2010-12-09 | Telefonaktiebolaget L M Ericsson (Publ) | Improved mobile radio channel estimation |
CN101697495A (zh) * | 2009-11-05 | 2010-04-21 | 湖南大学 | 一种基于博弈论的mimo信道跟踪方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103499811A (zh) * | 2013-08-28 | 2014-01-08 | 电子科技大学 | 提高雷达目标估计性能的天线个数分配方法 |
CN108957422A (zh) * | 2018-06-01 | 2018-12-07 | 电子科技大学 | 一种基于量化数据的云mimo雷达的根克拉美罗下界计算方法 |
CN108957422B (zh) * | 2018-06-01 | 2022-07-29 | 电子科技大学 | 一种基于量化数据的云mimo雷达的根克拉美罗下界计算方法 |
CN108667152A (zh) * | 2018-06-22 | 2018-10-16 | 南京大学 | 一种鲁棒的无线能量传输方法 |
CN108667152B (zh) * | 2018-06-22 | 2021-04-27 | 南京大学 | 一种鲁棒的无线能量传输方法 |
CN111585645A (zh) * | 2019-02-18 | 2020-08-25 | 富士通株式会社 | 光链路导致的偏振态变化监测装置及方法、光接收机 |
CN111585645B (zh) * | 2019-02-18 | 2023-03-31 | 富士通株式会社 | 光链路导致的偏振态变化监测装置及方法、光接收机 |
Also Published As
Publication number | Publication date |
---|---|
CN102723978B (zh) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102223327B (zh) | CoMP多用户系统中基于广义交替最大的信道估计方法 | |
CN104506259B (zh) | 宽带多天线通信系统接收通道间的时延差估计及校正方法 | |
CN102158272B (zh) | 一种射频通道的校准方法、装置及系统 | |
CN102647387B (zh) | 同频干扰的消除方法及装置 | |
CN102723978B (zh) | 多天线信道特征参数联合估计的下界建立方法 | |
CN105656819A (zh) | 一种基于压缩感知和大规模mimo的自适应信道估计方法 | |
CN101741508B (zh) | Lte发射模式7下cqi获取及使用的方法 | |
CN103888209A (zh) | 一种宽带接收阵列天线通道幅相误差时域校正方法 | |
CN101359950A (zh) | 一种基于奇异值分解迫零波束成型的多天线中继传输方法 | |
CN101621813A (zh) | 一种面向下行相干协作多点传输的上行信道估计方法 | |
CN101222458A (zh) | Mimo-ofdm信道的低阶递推最小均方误差估计 | |
CN101188448B (zh) | 一种智能天线的校准方法、装置及系统 | |
CN105450559A (zh) | 一种自干扰信道估计方法和设备 | |
CN105391664A (zh) | 基于mimo的lte全双工系统中模拟干扰抵消方法 | |
CN103929382A (zh) | 一种大规模mimo系统的大尺度衰落估计方法与装置 | |
CN101383643A (zh) | 一种mimo实时平台及其同步和信道联合估计方法 | |
Socheleau et al. | Concise derivation of scattering function from channel entropy maximization | |
CN106130938B (zh) | Tdd大规模mimo系统多用户联合信道估计方法 | |
CN101364961A (zh) | 基于mimo非实时平台的自适应信道估计方法 | |
CN102710565B (zh) | 分布式多天线移动信道特征参数的联合估计方法 | |
CN103166881B (zh) | 智能天线校准方法及系统 | |
Sim et al. | Low-complexity nonlinear self-interference cancellation for full-duplex radios | |
CN103701539A (zh) | 一种cdma用户数量测定方法 | |
CN101312597B (zh) | 一种信道估计单元和信道估计方法 | |
CN107809399A (zh) | 一种针对量化接收信号的多天线毫米波信道估计方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141210 Termination date: 20170629 |
|
CF01 | Termination of patent right due to non-payment of annual fee |