CN107332797A - 一种电力线ofdm通信***中的信道估计方法 - Google Patents
一种电力线ofdm通信***中的信道估计方法 Download PDFInfo
- Publication number
- CN107332797A CN107332797A CN201710460772.7A CN201710460772A CN107332797A CN 107332797 A CN107332797 A CN 107332797A CN 201710460772 A CN201710460772 A CN 201710460772A CN 107332797 A CN107332797 A CN 107332797A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- beta
- 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
- 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
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0256—Channel estimation using minimum mean square error criteria
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B3/00—Line transmission systems
- H04B3/54—Systems for transmission via power distribution lines
-
- 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
- H04L25/024—Channel estimation channel estimation algorithms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2626—Arrangements specific to the transmitter only
- H04L27/2627—Modulators
- H04L27/2628—Inverse Fourier transform modulators, e.g. inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
- Radio Transmission System (AREA)
Abstract
一种电力线OFDM通信***中的信道估计方法,其包括:(1)在发送端,将信息调制为频域信号,并转换为并行信号;(2)***导频,通过傅里叶逆变换IFFT之后,把频域信息变换为时域信息;(3)在经过IFFT后的信号中添加循环前缀,消除OFDM***中的ISI和ICI;经过并串转换后进入到低压电力线信道中,到达信号接收端;(4)在接收端移除循环前缀后再进行傅里叶变换FFT得到频域信号;(5)提取信号中导频位置处的信息,并对其进行信道估计,采用线性最小均方误差估计法LMMSE,利用P阶多项式展开来代替LMMSE方法中的求逆运算,并对多项式的系数进行优化,得到导频位置处的信道响应估计值。本发明能够有效估计OFDM通信***下的信道响应值,提高估计精度的同时减小了计算复杂度。
Description
技术领域
本发明涉及通信技术领域,具体涉及一种电力线OFDM通信***中的信道估计方法。
背景技术
电力线通信是一种以电力线为媒介传输数据信号的通信方式,其无需铺设额外的线路,大幅降低了通信成本。但由于电网的分支结构和节点处阻抗不匹配,使得电力线通信信道具有显著的多径特性,而正交频分复用(Orthogonal Frequency DivisionMultiplexing OFDM)技术可有效克服信道多径时延,且频谱利用率高,适用于电力线通信***。电力线信道不但存在多径效应,且信道是时变***,信道的传输函数随时间变化而变化,因而电力线信道存在频率弥散性和时间选择性衰落。OFDM通信***为了能够准确地恢复出发送的信号,需要在接收端进行分集接收、最大似然检测、相干解调,并且这些技术都需要利用信道估计信息,因此,作为OFDM技术中的关键技术,信道估计对建立可靠有效的电力线通信***有着非常重要的影响。
基于OFDM***的信道估计方法众多,一般可采用传统的最小二乘LS(LeastSquares)算法、最小均方误差MMSE(Minimum Mean Square Error)算法,或基于DFT变换的估计算法。LS算法简单,但在低信噪比时性能不理想;MMSE算法具有良好的性能,但统计信道相关性和矩阵求逆的复杂度较高;基于变换域的DFT算法将LS频域估计值变换到时域,通过时域置零实现降噪,该算法复杂度低于MMSE算法,但其性能对阈值的选取较为敏感。
发明内容
为了解决电力线OFDM通信***中的信道估计问题,本发明提供了一种低复杂度的信道估计方法,叫做多项式展开信道估计,即在已有的LMMSE估计的基础上利用P阶多项式展开来代替LMMSE方法中的求逆运算,此方法在保证估计准确度的情况下大大降低了运算复杂度。并进一步通过优化多项式的系数来使得在任何P阶多项式的情况下都能保证估计的准确性。
为了实现上述目的,本发明采用如下技术方案:
一种电力线OFDM通信***中的信道估计方法,其特征在于,包括以下步骤:(1)在发送端,将输入的信息调制为频域信号,通过串并转换将信号转换为并行信号;(2)***适当的导频,得到新的信息序列,通过傅里叶逆变换IFFT操作之后,把频域信息变换为时域信息;(3)在经过IFFT后的信号中添加循环前缀,消除OFDM***中的符号间干扰ISI和载波间干扰ICI;然后经过并串转换后,进入到低压电力线信道中,到达信号接收端;(4)在接收端移除循环前缀后再进行傅里叶变换FFT操作得到信号的频域形式;(5)提取信号中导频位置处的信息,并对其进行信道估计,采用线性最小均方误差估计法LMMSE,利用P阶多项式展开来代替LMMSE方法中的求逆运算,即P阶多项式展开信道估计方法PLMMSE,P表示自然数,得到导频位置处的信道响应估计值;再利用信道响应估计值对接收的信息进行相干解调以及符号检测还原出发送的信息。
优选地,步骤(5)中的PLMMSE算法,将展开后的多项式的每一项前面加一个系数,通过选择多项式的系数来使得在有限的P时降低信道估计误差,系数选取的原则是最小化估计的均方误差MSE。
优选地,根据步骤(5)中的LMMSE算法,信道估计值为:
这里的H=[H(0),H(1),…,H(N-1)]T,H(k)指第k个子载波的信道参数,
k∈[0,1,…N-1],N指OFDM子载波的个数;RHH=E[HHH]指子信道的自相关矩阵,HH指信道矩阵H的共轭转置;是通过最小二乘算法LS估计的信道值,X(k)则表示OFDM符号的第k个子载波传递的信息,Y(k)表示OFDM符号的第k个子载波上的接收信息;
Y=[Y(0),Y(1),…Y(N-1)]T;I指N阶单位矩阵;σ2指频域噪声方差;β是所有调制方式的指数,而且有
P阶多项式展开信道估计值为:
当参数α满足条件时,且当P→∞时,上述公式的约等号为等号,其中为矩阵的任意特征值,n∈[0,1,…,N]。
优选地,进一步对多项式的系数进行优化,得到加权的PLMMSE估计即W-PLMMSE估计为:
这里w=[w0,…,wP]T,是多项式的加权系数;同时,为了避免当P增大时,
的值变得无穷大,设置
选取加权值w=[w0,…,wP]T的原则是最小化估计的MSE,即令最小时计算w:得到最优系数wopt的计算为wopt=[w0,…,wP]T=A-1b,这里A是一个P+1行、P+1列的矩阵,b是一个P+1列矢量,A的第i行第j列元素以及b的第i个元素是:
并且矩阵A和b中的元素可从下式计算得到:
Yt指在t时刻N个子载波的接收信号,假设信道参数在时间窗口T内保持不变;通过采样协方差矩阵来求得;因为只需要计算矩阵的迹,所以T无需太大,对任意的T≥1,i≥1,可以得到
因为公式(5)的矩阵A和b中的元素是公式(7)和(8)的形式,所以通过公式(7)和(8)来计算,最终得到W-PLMMSE信道估计值从而降低计算复杂度。
附图说明
图1为本发明信道估计的实现框图。
具体实施方式
低压电力线OFDM通信***基于导频的信道估计原理框图如图1所示:
在发送端,首先将输入的二进制比特信息调制为频域信号,之后通过串并转换将信号转换为并行信号,然后根据电力线通信***的特性***适当的导频,得到新的信息序列Xm(k),设调制的子载波个数为N,Xm(k)则表示第m个OFDM符号上第k个子载波传递的信息,通过IFFT操作之后,把频域信息变换为时域信息,得到第m个OFDM符号对应的输出序列xm(n),则xm(n)表示为:
为了消除OFDM***中的符号间干扰ISI和载波间干扰ICI,在经过IFFT后的信号xm(n)中添加循环前缀,所以xm(n)变为了xmg(n),如下式:
xmg(n)=xm(n+N-NP),n=0,1,…N+NP-1 (2)
式中NP是指添加的循环前缀的符号个数。
信号xmg(n)经过并串转换后,进入到低压电力线信道中,低压电力线信道冲激响应可以描述为:
式中L表示信道的路径数;cl指第l条路径的衰减系数;τl指第l条路径的时间延迟,是指冲激函数。
信号经过低压电力线信道后到达接收端,则接收端接收到的信号ym(n)可以表示为:
式中hm(n)指第m个OFDM符号传输时的信道冲激响应;wm(n)是加性噪声干扰;是循环卷积。
在接收端移除循环前缀后再进行FFT操作得到了信号的频域形式,可以表示为:
Ym(k)=FFT(ym(n))=Xm(k)×Hm(k)+Wm(n) (5)
式中Hm(k)是信道的频率响应;Wm(n)是加性噪声的频域形式。
提取信号Ym(k)中导频位置处的信息YP(k),并对其进行信道估计,得到导频位置处的信道响应估计值HP(k)。再利用HP(k)对接收的信息进行相干解调以及符号检测还原出发送的二进制比特信息。
基于MMSE准则的信道估计方法为
这里的H=[H(0),H(1),…,H(N-1)]T,RHH=E[HHH]指子信道的自相关矩阵,X是X(0),X(1),…X(N-1)为主对角线的对角矩阵,Y=[Y(0),Y(1),…Y(N-1)]T,σ2指频域噪声方差,是通过LS估计的信道,计算简单,可表示成
MMSE估计准则将信道的先验统计特性和信道中的噪声干扰同时考虑在内,使得信道估计的准确度得到了很大的改善,但是也可以看出(XXH)-1会随着输入信号序列X的变化而变化,并且每次变化时都需要对整个信号矩阵重新求逆,即在估计每一个位置的信道响应时,都要进行求逆操作,当OFDM的子信道数增大时,矩阵的求逆运算的计算的复杂度将增加,每次要消耗大量的时间和硬件资源,不利于运用到实践中,为了解决MMSE准则的这个问题,通常用(XXH)-1的期望E[(XXH)-1]来代替(XXH)-1,因此提出了LMMSE方法:
这里I指N阶单位矩阵,SNR为β是所有调制方式的指数,而且有
但是从公式(7)中也可以看出该算法还是存在一次求逆的操作,若发送端发送大量的信息数据,矩阵的阶数较大,在对矩阵进行求逆时,运算量依旧较大,为了能够广泛应用LMMSE还是要继续降低LMMSE算法的运算复杂度,因此,本发明提出了一种多项式展开信道估计方法,并通过优化多项式的系数来使得在任意阶数的多项式展开都能得到较准确的信道估计。
本发明中的这种多项式展开信道估计方法利用了下面这条定理:
对于任意的厄米特矩阵它所有特征值满足条件λn(X)<1,那么矩阵X满足
从上述公式看出,随着p的增大,Xp的值将减小,所以,多项式的个数不需要无穷个,一个P阶的多项式便可以近似代替矩阵的求逆运算,并且可以通过选择P的大小来平衡计算复杂度和近似估计误差。
将上述定理应用到没有特征值分布限制的任意矩阵中,得到下列命题:
对于任意的厄米特矩阵X满足
若参数α满足条件当P→∞时,上述公式的约等号为等号。
将上述命题的近似性质应用于公式(7)的LMMSE估计算法中,便得到了P阶多项式展开信道估计方法
从上式可以看出,矩阵的求逆运算通过多项式展开给代替了,从而大大降低了计算复杂度。
还需要选择参数α来满足公式(9)中的性质,从简化复杂度的角度考虑,本发明选取α为另一方面,若是从加快多项式收敛的角度来看,可以选取最优参数αopt为:
此时将最大化收敛速度,因为的最大和最小特征值是关于原点对称的。
当P→∞时,PLMMSE估计的性能将逼近于MMSE估计性能,但是通常情况下,当P是一个有限值时,它估计的性能并不是最好的。因此我们想出一种方法,将展开后的每一项前面加一个参数,再求和,通过最优化多项式的系数来使得在有限的P时性能最优。因此我们得到加权的PLMMSE估计(W-PLMMSE估计)为
这里w=[w0,…,wP]T,是多项式的加权系数。同时,为了避免当P增大时,的值变得无穷大,我们同样设置
我们选取加权值的原则是最小化估计的MSE。即令最小时计算w。
因为W-PLMMSE估计的MSE为:
这里给定多项式的阶数P时,可以通过优化参数w0,…,wP来最小化MSE。即求得这个问题,可令MSE函数的偏导为零求得,也就是:
这是一个包含P+1个未知参数的P+1个方程,可以写成Aw=b的形式,这里且它的第i行第j列元素以及的第i个元素是
注意到,我们这里做了一个变量代换,矩阵A中的i=p1+1,j=p2+1,矩阵b中的i=p+1。所以对于任意P阶多项式,得到最优参数的计算为wopt=[w0,…,wP]T=A-1b。
这里,在通过A-1b计算最优参数时,复杂度较高,但是我们发现
Yt指在t时刻n个子载波的接收信号。也就是说可以通过采样协方差矩阵来求得。因为我们只需要计算矩阵的迹,所以T无需太大,对任意的T≥1,i≥1,可以得到
因为公式(16)的矩阵A和b中的元素是公式(18)和(19)的形式,所以我们可以通过公式(18)和(19)来计算,最终得到W-PLMMSE信道估计值从而降低复杂度。在平滑时间窗口区间长度T之间,定义为通过接收信号Yt,…Yt-T+1得到的,从公式(18)、(19)可看出,可以从 得出,因此进一步大大降低了复杂度。
Claims (4)
1.一种电力线OFDM通信***中的信道估计方法,其特征在于,包括以下步骤:
(1)在发送端,将输入的信息调制为频域信号,通过串并转换将信号转换为并行信号;
(2)***适当的导频,得到新的信息序列,通过傅里叶逆变换IFFT操作之后,把频域信息变换为时域信息;
(3)在经过IFFT后的信号中添加循环前缀,消除OFDM***中的符号间干扰ISI和载波间干扰ICI;然后经过并串转换后,进入到低压电力线信道中,到达信号接收端;
(4)在接收端移除循环前缀后再进行傅里叶变换FFT操作得到信号的频域形式;
(5)提取信号中导频位置处的信息,并对其进行信道估计,采用线性最小均方误差估计法LMMSE,利用P阶多项式展开来代替LMMSE方法中的求逆运算,即P阶多项式展开信道估计方法PLMMSE,P表示自然数,得到导频位置处的信道响应估计值;再利用信道响应估计值对接收的信息进行相干解调以及符号检测还原出发送的信息。
2.根据权利要求1所述的方法,其特征在于:步骤(5)中的PLMMSE算法,将展开后的多项式的每一项前面加一个系数,通过选择多项式的系数来使得在有限的P时降低信道估计误差,系数选取的原则是最小化估计的均方误差MSE。
3.根据权利要求2所述的方法,其特征在于:根据步骤(5)中的LMMSE算法,信道估计值为:
<mrow>
<msub>
<mover>
<mi>H</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>L</mi>
<mi>M</mi>
<mi>M</mi>
<mi>S</mi>
<mi>E</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mover>
<mi>H</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>L</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
这里的H=[H(0),H(1),…,H(N-1)]T,H(k)指第k个子载波的信道参数,k∈[0,1,…N-1],N指OFDM子载波的个数;RHH=E[HHH]指子信道的自相关矩阵,HH指信道矩阵H的共轭转置;是通过最小二乘算法LS估计的信道值,X(k)则表示OFDM符号的第k个子载波传递的信息,Y(k)表示OFDM符号的第k个子载波上的接收信息;
Y=[Y(0),Y(1),…Y(N-1)]T;I指N阶单位矩阵;σ2指频域噪声方差;β是所有调制方式的指数,而且有
P阶多项式展开信道估计值为:
<mrow>
<msub>
<mover>
<mi>H</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>P</mi>
<mi>L</mi>
<mi>M</mi>
<mi>M</mi>
<mi>S</mi>
<mi>E</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>P</mi>
</munderover>
<mi>&alpha;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>p</mi>
</msup>
<msub>
<mover>
<mi>H</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>L</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
当参数α满足条件时,且当P→∞时,上述公式的约等号为等号,其中为矩阵的任意特征值,n∈[0,1,…,N]。
4.根据权利要求3所述的方法,其特征在于:进一步对多项式的系数进行优化,得到加权的PLMMSE估计即W-PLMMSE估计为:
<mrow>
<msub>
<mover>
<mi>H</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>W</mi>
<mo>-</mo>
<mi>P</mi>
<mi>L</mi>
<mi>M</mi>
<mi>M</mi>
<mi>S</mi>
<mi>E</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>P</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>p</mi>
</msub>
<msup>
<mi>&alpha;</mi>
<mrow>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mi>p</mi>
</msup>
<msub>
<mover>
<mi>H</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>L</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
这里w=[w0,…,wP]T,是多项式的加权系数;同时,为了避免当p增大时,
的值变得无穷大,设置
<mrow>
<mi>&alpha;</mi>
<mo>&le;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>max</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&lambda;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
选取加权值w=[w0,…,wP]T的原则是最小化估计的MSE,即令最小时计算w:
得到最优系数wopt的计算为wopt=[w0,…,wP]T=A-1b,这里A是一个P+1行、P+1列的矩阵,b是一个P+1列矢量,A的第i行第j列元素以及b的第i个元素是:
<mrow>
<msub>
<mrow>
<mo>&lsqb;</mo>
<mi>A</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mi>j</mi>
</mrow>
</msup>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mrow>
<mo>&lsqb;</mo>
<mi>b</mi>
<mo>&rsqb;</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<mi>&alpha;</mi>
<mi>i</mi>
</msup>
<mi>t</mi>
<mi>r</mi>
<mo>(</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
并且矩阵A和b中的元素可从下式计算得到:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>X</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>Y</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>X</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>Y</mi>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>=</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>E</mi>
<mo>&lsqb;</mo>
<msup>
<mi>YY</mi>
<mi>H</mi>
</msup>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<munder>
<mi>lim</mi>
<mrow>
<mi>T</mi>
<mo>&RightArrow;</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<msup>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<mi>H</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Yt指在t时刻N个子载波的接收信号,假设信道参数在时间窗口T内保持不变;通过采样协方差矩阵来求得;因为只需要计算矩阵的迹,所以T无需太大,对任意的T≥1,i≥1,可以得到
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&ap;</mo>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<msup>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<mi>H</mi>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msup>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
2
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&ap;</mo>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<msup>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<mi>H</mi>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msup>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>R</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
因为公式(5)的矩阵A和b中的元素是公式(7)和(8)的形式,所以通过公式(7)和(8)来计算,最终得到W-PLMMSE信道估计值从而降低计算复杂度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710460772.7A CN107332797B (zh) | 2017-06-18 | 2017-06-18 | 一种电力线ofdm通信***中的信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710460772.7A CN107332797B (zh) | 2017-06-18 | 2017-06-18 | 一种电力线ofdm通信***中的信道估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107332797A true CN107332797A (zh) | 2017-11-07 |
CN107332797B CN107332797B (zh) | 2020-02-18 |
Family
ID=60195234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710460772.7A Active CN107332797B (zh) | 2017-06-18 | 2017-06-18 | 一种电力线ofdm通信***中的信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107332797B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108649997A (zh) * | 2018-04-19 | 2018-10-12 | 国网重庆市电力公司电力科学研究院 | 一种多入多出电力线通信窄带噪声的自适应消除方法、*** |
CN108737299A (zh) * | 2018-05-11 | 2018-11-02 | 重庆邮电大学 | 一种低复杂度的lmmse信道估计方法 |
CN108737301A (zh) * | 2018-05-23 | 2018-11-02 | 南通大学 | 一种基于b样条神经网络的宽带通信发射机指纹估计方法 |
CN109617847A (zh) * | 2018-11-26 | 2019-04-12 | 东南大学 | 一种基于模型驱动深度学习的无循环前缀ofdm接收方法 |
CN111490951A (zh) * | 2020-03-05 | 2020-08-04 | 熊军 | 一种基于多相fft的信号估计和信道内插方法及装置 |
CN112565141A (zh) * | 2020-10-22 | 2021-03-26 | 浙江工商大学 | 一种采用时域滤波的ofdm***抗多径时变信道混合干扰方法 |
CN113973031A (zh) * | 2021-10-28 | 2022-01-25 | 重庆邮电大学 | 一种ofdm***的信道均衡方法 |
CN114143146A (zh) * | 2021-10-28 | 2022-03-04 | 西安交通大学 | 一种基于图信号方法的ofdm***信道估计***和方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101056293A (zh) * | 2006-04-13 | 2007-10-17 | 中兴通讯股份有限公司 | 一种正交频分复用***的半盲信道估计方法 |
US20100054314A1 (en) * | 2006-12-27 | 2010-03-04 | Abb Technology Ag | Initialization of and modem for an ofdm data transmission |
CN102364902A (zh) * | 2011-10-08 | 2012-02-29 | 中国科学院上海微***与信息技术研究所 | 一种用于限幅ofdm***的信道估计方法 |
CN102821078A (zh) * | 2009-12-31 | 2012-12-12 | 宁波大学 | 基于ofdm信号循环平稳特性的盲信道估计方法 |
CN103209145A (zh) * | 2013-04-15 | 2013-07-17 | 深圳市力合微电子有限公司 | 一种基于ofdm调制的电力线载波信道检测方法 |
-
2017
- 2017-06-18 CN CN201710460772.7A patent/CN107332797B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101056293A (zh) * | 2006-04-13 | 2007-10-17 | 中兴通讯股份有限公司 | 一种正交频分复用***的半盲信道估计方法 |
US20100054314A1 (en) * | 2006-12-27 | 2010-03-04 | Abb Technology Ag | Initialization of and modem for an ofdm data transmission |
CN102821078A (zh) * | 2009-12-31 | 2012-12-12 | 宁波大学 | 基于ofdm信号循环平稳特性的盲信道估计方法 |
CN102364902A (zh) * | 2011-10-08 | 2012-02-29 | 中国科学院上海微***与信息技术研究所 | 一种用于限幅ofdm***的信道估计方法 |
CN103209145A (zh) * | 2013-04-15 | 2013-07-17 | 深圳市力合微电子有限公司 | 一种基于ofdm调制的电力线载波信道检测方法 |
Non-Patent Citations (1)
Title |
---|
***: "《基于导频的OFDM信道估计技术研究》", 《CNKI优秀硕士学位论文全文库》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108649997A (zh) * | 2018-04-19 | 2018-10-12 | 国网重庆市电力公司电力科学研究院 | 一种多入多出电力线通信窄带噪声的自适应消除方法、*** |
CN108737299A (zh) * | 2018-05-11 | 2018-11-02 | 重庆邮电大学 | 一种低复杂度的lmmse信道估计方法 |
CN108737301A (zh) * | 2018-05-23 | 2018-11-02 | 南通大学 | 一种基于b样条神经网络的宽带通信发射机指纹估计方法 |
CN108737301B (zh) * | 2018-05-23 | 2020-08-07 | 南通大学 | 一种基于b样条神经网络的宽带通信发射机指纹估计方法 |
CN109617847A (zh) * | 2018-11-26 | 2019-04-12 | 东南大学 | 一种基于模型驱动深度学习的无循环前缀ofdm接收方法 |
CN109617847B (zh) * | 2018-11-26 | 2021-04-06 | 东南大学 | 一种基于模型驱动深度学习的无循环前缀ofdm接收方法 |
CN111490951B (zh) * | 2020-03-05 | 2023-03-07 | 西安宇飞电子技术有限公司 | 一种基于多相fft的信号估计和信道内插方法及装置 |
CN111490951A (zh) * | 2020-03-05 | 2020-08-04 | 熊军 | 一种基于多相fft的信号估计和信道内插方法及装置 |
CN112565141A (zh) * | 2020-10-22 | 2021-03-26 | 浙江工商大学 | 一种采用时域滤波的ofdm***抗多径时变信道混合干扰方法 |
CN112565141B (zh) * | 2020-10-22 | 2022-07-19 | 浙江工商大学 | 一种采用时域滤波的ofdm***抗多径时变信道混合干扰方法 |
CN114143146A (zh) * | 2021-10-28 | 2022-03-04 | 西安交通大学 | 一种基于图信号方法的ofdm***信道估计***和方法 |
CN113973031A (zh) * | 2021-10-28 | 2022-01-25 | 重庆邮电大学 | 一种ofdm***的信道均衡方法 |
CN114143146B (zh) * | 2021-10-28 | 2023-08-15 | 西安交通大学 | 一种基于图信号方法的ofdm***信道估计***和方法 |
CN113973031B (zh) * | 2021-10-28 | 2023-11-03 | 重庆邮电大学 | 一种ofdm***的信道均衡方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107332797B (zh) | 2020-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107332797A (zh) | 一种电力线ofdm通信***中的信道估计方法 | |
Al-Naffouri et al. | A model reduction approach for OFDM channel estimation under high mobility conditions | |
CN103873422B (zh) | 水声正交频分复用***符号内多径干扰消除方法 | |
CN101222458B (zh) | Mimo-ofdm信道的低阶递推最小均方误差估计 | |
CN106506415B (zh) | 一种多用户mimo-ofdm***信道估计的方法 | |
CN113676289B (zh) | 一种基于变换域最大比合并的otfs调制信号检测方法 | |
WO2008113216A1 (fr) | Méthode d'évaluation d'un canal | |
CN101018219B (zh) | 一种空频信号处理方法 | |
CN110048972A (zh) | 一种水声正交频分复用信道估计方法及*** | |
CN105337906A (zh) | 信道估计方法及装置 | |
CN107026804A (zh) | Mimo‑ofdm***中基于指数平滑的信道估计方法 | |
US8873681B2 (en) | Process for suppressing intercarrier interference in a OFDM receiver | |
Zaier et al. | Blind channel estimation enhancement for MIMO-OFDM systems under high mobility conditions | |
CN109302240A (zh) | 基于双选择衰落信道的低复杂度osdm串行均衡方法 | |
Bhoyar et al. | Leaky least mean square (LLMS) algorithm for channel estimation in BPSK-QPSK-PSK MIMO-OFDM system | |
Jie et al. | An improved DFT-based channel estimation algorithm for MIMO-OFDM systems | |
Hajizadeh et al. | Channel Estimation in OFDM System Based on the Linear Interpolation, FFT and Decision Feedback | |
Kahlon et al. | Channel estimation techniques in MIMO-OFDM systems–review article | |
CN111212002A (zh) | 基于子空间算法的海洋水声ofdm信道的盲辨识方法 | |
Manzoor et al. | A novel noise power and SNR estimation in WiMAX systems | |
Thomas et al. | Null-subcarrier based channel estimation and mutual interference reduction in MIMO OFDM systems | |
CN109217954A (zh) | 基于双选择衰落信道的低复杂度osdm块均衡方法 | |
CN104301263B (zh) | 一种多带uwb***低复杂度信道估计方法及装置 | |
Kaur et al. | Improved channel estimation in transmit diversity MIMO-OFDM using particle swarm optimization | |
Yoshida et al. | Degenerated-inverse-matrix-based channel estimation for OFDM systems |
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 |