CN103716265B - 一种相位噪声补偿抑制改进方法 - Google Patents
一种相位噪声补偿抑制改进方法 Download PDFInfo
- Publication number
- CN103716265B CN103716265B CN201410007166.6A CN201410007166A CN103716265B CN 103716265 B CN103716265 B CN 103716265B CN 201410007166 A CN201410007166 A CN 201410007166A CN 103716265 B CN103716265 B CN 103716265B
- Authority
- CN
- China
- Prior art keywords
- matrix
- phase noise
- estimated
- channel
- noise
- 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
- Noise Elimination (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明属于无线通信技术领域，具体涉及在无线通信系统中通过运用迭代的方法来实现相位噪声的补偿抑制。本发明是先通过信道估计序列，估计出等效离散时域信道冲激响应。然后通过插值估计出相位噪声的公共相位误差（CPE），最后通过迭代的方法实现相位噪声补偿，提高系统的可靠性，减小误码率。
Description
技术领域
本发明属于无线通信技术领域，具体涉及在无线通信系统中通过运用迭代的方法来实现相位噪声的补偿抑制。
背景技术
随着人们对无线通信需求的爆炸式增长，频谱资源越发的显得宝贵。相应地，这就要求我们更好的处理诸如本地振荡器的相位噪声问题。
在实际的通信系统中，由模拟前端产生的信号失真将会引起系统性能下降。由于射频电路的面积、功耗和成本的权衡，由模拟前端产生的信号失真不可避免。无线通信的模拟前端一项重要的线性失真是由压控振荡器或者锁相环的随机相位噪声引起的。
相位噪声，实际上是对频率源频率稳定度的一种表征。通常情况下，频率稳定度分为长期频率稳定度和短期频率稳定度。所谓短期频率稳定度，是指由随机噪声引起的相位起伏或频率起伏。至于因为温度、老化等引起的频率慢漂移，则称之为长期频率稳定度。通常主要考虑的是短期稳定度问题，可以认为相位噪声就是短期频率稳定度，只不过是一个物理现象的两种不同表示方式。对于振荡器，频率稳定度是它在整个规定的时间范围内产生相同频率的一种量度。如果信号频率存在瞬时的变化，不能保持不变，那么信号源就存在着不稳定性，起因就是相位噪声。
如果没有相位噪声，那么振荡器的整个功率都应集中在中心频率处。但相位噪声的出现将振荡器的一部分功率扩展到相邻的频率中去，产生了边带。相位噪声通常定义为在某一给定偏移频率处的dBc/Hz值，其中，dBc是以dB为单位的该频率处功率与总功率的比值。一个振荡器在某一偏移频率处的相位噪声定义为在该频率处1Hz带宽内的信号功率与信号的总功率比值。
相位噪声的补偿或者抑制算法研究大部分集中在OFDM系统中。总体来说，有基于反馈和无反馈的相位噪声补偿抑制算法。在单载波通信系统中，相位噪声不仅会产生公共相位误差(CPE)，在频率选择信道下还会产生码间串扰(ISI)。
发明内容
本发明的目的在于提供一种相位噪声补偿抑制改进方法,提高信号传输的可靠性，降低误码率。
为了方便地描述本发明的内容，首先对本发明中所使用的属于定义进行说明：
特殊字(UW，Unique Word)。为了在接收端进行同步或参数估计等，在发送端发送的具有某些特定特性的、对接收端已知的特殊序列。
本发明的目的通过如下步骤实现：
S1、利用信道估计序列实现信道估计，得到等效时域信道的冲激响应估计值
S2、接收端通过迭代对相位噪声进行补偿，包括：
S21、设置UW的长度大于等效时域，将所述UW加到要传输的数据序列中去，其中，UW为特殊字，是为了在接收端进行同步或参数估计等，在发送端发送的具有某些特定特性的、对接收端已知的特殊序列；
S22、第i个UW通过信道时的接收信号可以近似表示为：
则过传输数据序列中的UW，采取插值的方法估计出第i个UW与第i+1个UW之间传输的第i个数据块的相位噪声常数与a的比：其中，a是相位噪声常数；
S23、对于接收信号而言，去除CP后接收信号频域可以用矩阵形式表示为：YN×1＝AN×NHN×NXN×1+WN×1，其中AN×N为相位噪声频域构成的toeplitz矩阵，HN×N为估计信道的频域构成的对角矩阵，XN×1为传输数据频域构成的矩阵，YN×1为接收信号的频域构成的矩阵，WN×1为高斯白噪声频域构造的矩阵；
S24、接收信号时域表示为：yN×1＝diag(pN×1)hN×NxN×1+wN×1，其中，yN×1为接收信号的时域构成的矩阵，diag(pN×N)为相位噪声时域构成的对角矩阵，hN×N为信道的时域构成的toeplitz矩阵，xN×1为传输数据时域构成的矩阵，wN×1为噪声时域构造成的矩阵。为了降低复杂度，通过构造插值矩阵PN×N，使pN×1＝PN×Ncs×1。我们把作为初始条件构造A1,N×N，利用YN×1＝AN×NHN×NXN×1+WN×1及降噪处理估计出然后利用yN×1＝diag(pN×1)hN×NxN×1+wN×1以及pN×1元素的模值为1估计出
S25、利用构造A2,N×N，利用YN×1＝AN×NHN×NXN×1+WN×1及降噪处理估计出再利用yN×1＝diag(pN×1)hN×NxN×1+wN×1以及pN×1元素的模值为1估计出不断的迭代下去，直到：
||yN×1-diag(pi+1,N×1)hN×Nxi+1,N×1||2≤||yN×1-diag(pi+1,N×1)hN×Nxi,N×1||2，收敛为止，从而实现传输数据块估计
本发明的有益效果是：本发明是先通过信道估计序列，估计出等效离散时域信道冲激响应，然后通过插值估计出相位噪声的公共相位误差(CPE)，最后通过迭代的方法实现相位噪声补偿，提高系统的可靠性，减小误码率。
附图说明
图1是本发明使用的单载波频域均衡系统模型图
图2是本发明发送端物理层发送数据结构图；
图3是本发明使用的相位噪声统计模型图；
图4是本发明实现相位噪声估计补偿抑制的流程图；
图5是本发明算法误比特率(BER)性能曲线图。
具体实施方式
下面结合附图来说明本发明的具体实施方式：
S1、利用信道估计序列实现信道估计，包括：
S11、信道估计序列是由一些已知符号构成的序列，例如802.11.ad标准中单载波信道估计序列是[-Gb128,-Ga128,Gb128,-Ga128,-Gb128,Ga128,-Gb128,-Ga128,-Gb128]，其中Ga128和Gb128是格雷序列构成。
S12、在接收端我们可以得到时域的信号表示为矩阵的形式：yN×1＝AN×NhN×NxN×1+wN×1，其中，yN×1是一个N×1列向量的形式，其受到相位噪声和高斯白噪声的影响。AN×N是一个N×N的对角矩阵由相位噪声构成，hN×N是由等效时域信道冲激响应构成的托普利兹(toeplitz)矩阵，xN×1是由传输数据构成的N×1列向量，wN×1是N×1的噪声向量。
S13、我们既可以利用一些传统信道估计法，例如基于序列相关的信道估计、最小二乘法(LS)信道估计等等，也可以利用一些比较好的新颖的方法来实现信道估计，如正交匹配追踪算法(OMP，Orthogonal Matching Pursuit)。
S14、我们以频域最小二乘法LS信道估计为例(等效时域信道冲激响应长度小于N)，FN×N是归一化的N×N傅立叶矩阵，0≤k≤N-1，0≤i≤N-1。频域接收信号矩阵表示为YN×1＝aXN×NHN×1+WN×1,其中，YN×1＝FN×NyN×1，XN×N是一个对角矩阵，其主对角元素是由时域传输数据xN×1经过N点fft变换得到的XN×1，HN×1是一个N×1矩阵，HN×1＝FN×NhN×1，hN×1是由等效时域信道冲激响应h后面添0构成的N×1矩阵。WN×1是时域噪声wN×1经过N点归一化fft变化后得到。在这里我们认为相位噪声是一个常数a，所以再经过N点归一化ifft得到然后通过设定门限值得到估计的等效时域信道冲激响应这里a是相位噪声常数。
S2、接收端通过迭代对相位噪声进行补偿，包括：
S21、通过信道估计我们得到了估计的等效时域信道的冲激响应这里我们以块状数据传输为例，我们之前已经定义UW为已经知道的序列，我们把UW加到要传输的数据序列中去，实现去除码间干扰。
S22、我们假定UW长度大于等效时域信道长度，这里我们知道传输的数据序列中有UW，UW通过信道实质是等效时域信道冲激响应h与UW卷积，但是其同时受到相位噪声和高斯白噪声的影响，这里我们认为UW序列中受到的相位噪声为一常数(不同UW受到的相位噪声常数不同，这里可以认为相位噪声常数实际是相位噪声的公共相位误差CPE)，第i个UW通过信道时的接收信号可以近似表示为：这样我们可以通过传输数据序列中的UW，采取插值的方法估计出第i个UW与第i+1个UW之间传输的第i个数据块的相位噪声常数与a的比：其中，a是相位噪声常数。
S23、对于接收信号而言，去除CP后接收信号频域可以用矩阵形式表示为：YN×1＝AN×NHN×NXN×1+WN×1，其中AN×N为相位噪声频域构成的toeplitz矩阵，HN×N为估计信道的频域构成的对角矩阵，XN×1为传输数据频域构成的矩阵，YN×1为接收信号的频域构成的矩阵，WN×1为高斯白噪声频域构造的矩阵。相应地，接收信号时域表示为：yN×1＝diag(pN×1)hN×NxN×1+wN×1，其中，yN×1接收信号的时域构成的矩阵，diag(pN×N)为相位噪声时域构成的对角矩阵，hN×N为信道的时域构成的toeplitz矩阵，xN×1为传输数据时域构成的矩阵，wN×1为噪声时域构造成的矩阵。为了降低复杂度，通过构造插值矩阵PN×N，使pN×1＝PN×Ncs×1。我们把作为初始条件构造A1,N×N，利用YN×1＝AN×NHN×NXN×1+WN×1及降噪处理估计出然后利用yN×1＝diag(pN×1)hN×NxN×1+wN×1以及pN×1元素的模值为1估计出接着再利用构造A2,N×N，利用YN×1＝AN×NHN×NXN×1+WN×1及降噪处理估计出再利用yN×1＝diag(pN×1)hN×NxN×1+wN×1以及pN×1元素的模值为1估计出不断的迭代下去，直到：
||yN×1-diag(pi+1,N×1)hN×Nxi+1,N×1||2≤||yN×1-diag(pi+1,N×1)hN×Nxi,N×1||2
收敛为止，从而实现传输数据块估计补偿了相位噪声，提高了误码率。
图1是本发明使用的单载波频域均衡系统模型图。
图2是本发明发送端物理层发送数据结构图。
我们把信道估计序列(CES)放在要发送的实际数据之前，在实际要传输的数据序列之间插入已知序列UW，传输是以数据块的形式进行传输。这样UW可以作为等效的循环前缀使用(即CP)。此例所用的UW是格雷序列，UW扩展后的数据块长度Nb＝512，数据块CP长度Ncp＝64。设发送用户数据为：xi＝[xi(0),xi(1),...xi(NS-1)]T，Ns是用户数据长度为448。用户数据间插入的长度Ncp的UW为：u＝[u(0),u(1),...,u(Ncp-1)]T。UW扩展后的第i数据块表示为xi＝[xi T,uT]T，其长度Nb＝Ns+Ncp。等效离散信道脉冲响应(CIR，channel impulseresponse)为：h＝[h(0),h(1),...h(L-1)]T，L是信道长度(L≤Ncp)。假设其它条件理想，第i个数据块受到相位噪声矩阵Ai的影响，及AWGN向量的加性影响，其中它的第m个元素wi(m)是均值为0，方差为的AWGN。接收信号为其中h0是Nb×Nb的toeplitz矩阵，它的第一列是h＝[h(0),h(1),...h(L-1),0,0,...,0]T，yi＝[yi(0),yi(1),...,yi(Nb-1)]T。
图3是本发明使用的相位噪声统计模型图。
相位噪声一般用它的功率谱密度(PSD)来表征，通信标准IEEE802.15.3c和IEEE802.11ad给出了关于相位噪声PSD的“一个极点/零点”模型：
式中，f表示偏移载波中心的频率。PSD(0)是一个常数，在IEEE802.15.3c和IEEE802.11ad标准分别取值-87dbc/Hz和-90dbc/Hz。fp＝1MHz是极点频率，fZ＝100MHz是零点频率。
我们把平稳高斯色相位噪声看作是零均值的高斯白噪声通过低通滤波器得到的带限高斯色噪声。该低通滤波器传输函数的模平方可以分解为：
从上式中选取一个稳定系统的传输函数,得到模拟低通滤波器的传输函数H(s)：
其中s＝jΩ＝j2πf，a＝1/2πfZ，b＝1/2πfp，fp和fZ分别是的极点频率和零点频率。我们利用双线性变换法实现模拟滤波器转化为数字滤波器：
其中常数C＝π(fp+fz)/tan(π(fp+fz)/2fs)使在(fp+fz)/2这个频率点处的模拟频率转换到数字频率成线性关系。由此可以得到数字滤波器的传输函数：
图4是本发明实现相位噪声补偿抑制的流程图；
图5是使用图1的系统模型、图2的数据传输结构、图3的相位噪声模型结构和图4的算法流程，应用到具体的通信系统中，这里选取单载波频域均衡系统作为一个实例，仿真得到的是本发明算法在单载波频域系统中的误比特率(BER)性能曲线图。其中，图5表示在IEEE 802.11.ad信道标准定义的视距(LOS)信道模型中不同比特信噪比Eb/N0(dB)的性能曲线图。本例的仿真系统是属于高频高速超宽带通信系统，它主要仿真参数是：载波频率为60GHz，符号速率为1.76Gbps，QPSK调制,相位噪声为-80dbc/Hz@1MHz。
Claims (1)
1.一种相位噪声补偿抑制改进方法，其特征在于，包括以下步骤：S1、利用信道估计序列实现信道估计，得到等效时域信道的冲激响应估计值
S2、接收端通过迭代对相位噪声进行补偿，包括：
S21、设置UW的长度大于等效时域信道冲激相应长度，将所述UW加到要传输的数据序列中去，其中，UW为特殊字，是为了在接收端进行同步或参数估计等，在发送端发送的具有某些特定特性的、对接收端已知的特殊序列；
S22、第i个UW通过信道时的接收信号可以近似表示为：
则过传输数据序列中的UW，采取插值的方法估计出第i个UW与第i+1个UW之间传输的第i个数据块的相位噪声常数与a的比：其中，a是相位噪声常数；
S23、对于接收信号而言，去除CP后接收信号频域可以用矩阵形式表示为：YN×1＝AN× NHN×NXN×1+WN×1，其中AN×N为相位噪声频域构成的toeplitz矩阵，HN×N为估计信道的频域构成的对角矩阵，XN×1为传输数据频域构成的矩阵，YN×1为接收信号的频域构成的矩阵，WN×1为高斯白噪声频域构造的矩阵；
S24、接收信号时域表示为：yN×1＝diag(pN×1)hN×NxN×1+wN×1，其中，yN×1为接收信号的时域构成的矩阵，diag(pN×N)为相位噪声时域构成的对角矩阵，hN×N为信道的时域构成的toeplitz矩阵，xN×1为传输数据时域构成的矩阵，wN×1为噪声时域构造成的矩阵，为了降低复杂度，通过构造插值矩阵PN×N，使pN×1＝PN×Ncs×1，我们把作为初始条件构造A1,N×N，利用YN×1＝AN×NHN×NXN×1+WN×1及降噪处理估计出然后利用yN×1＝diag(pN×1)hN×NxN×1+wN×1以及pN×1元素的模值为1估计出
S25、利用构造A2,N×N，利用YN×1＝AN×NHN×NXN×1+WN×1及降噪处理估计出再利用yN×1＝diag(pN×1)hN×NxN×1+wN×1以及pN×1元素的模值为1估计出不断的迭代下去，直到：||yN×1-diag(pi+1,N×1)hN×Nxi+1,N×1||2≤||yN×1-diag(pi+1,N×1)hN×Nxi,N×1||2，收敛为止，从而实现传输数据块估计
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410007166.6A CN103716265B (zh) | 2014-01-07 | 2014-01-07 | 一种相位噪声补偿抑制改进方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410007166.6A CN103716265B (zh) | 2014-01-07 | 2014-01-07 | 一种相位噪声补偿抑制改进方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103716265A CN103716265A (zh) | 2014-04-09 |
| CN103716265B true CN103716265B (zh) | 2017-05-03 |
Family
ID=50408859
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410007166.6A Expired - Fee Related CN103716265B (zh) | 2014-01-07 | 2014-01-07 | 一种相位噪声补偿抑制改进方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103716265B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104022983B (zh) * | 2014-05-09 | 2018-02-23 | 国家电网公司 | 一种ofdm系统中的cpe抑制方法 |
| CN104954305A (zh) * | 2015-05-31 | 2015-09-30 | 电子科技大学 | 无线通信系统下的一种相位噪声改进估计方法 |
| CN105915476A (zh) * | 2016-04-15 | 2016-08-31 | 电子科技大学 | 一种基于贝叶斯的相噪补偿方法 |
| CN107819544B (zh) * | 2016-09-13 | 2020-12-18 | 张贻和 | 一种降低信道误码率的方法 |
| CN108259126B (zh) * | 2016-12-29 | 2020-12-15 | 华为技术有限公司 | 相位噪声参考信号传输方法和装置 |
| CN106850469B (zh) * | 2016-12-30 | 2020-08-25 | 上海华为技术有限公司 | 一种相位噪声估计的方法及相关装置 |
| CN107231216B (zh) * | 2017-07-04 | 2019-09-27 | 电子科技大学 | 基于gamp算法的相位噪声补偿抑制方法 |
| CN108521311B (zh) * | 2018-04-24 | 2020-11-27 | 成都吉纬科技有限公司 | 一种基于格雷序列的信噪比估计方法 |
| CN109194594B (zh) * | 2018-08-28 | 2020-12-18 | 电子科技大学 | 一种基于连续载波聚合的相位噪声抑制方法 |
| CN110191070A (zh) * | 2019-05-31 | 2019-08-30 | 北京邮电大学 | 一种相位噪声抑制方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101394381A (zh) * | 2007-09-19 | 2009-03-25 | 中国科学院微电子研究所 | 一种用于相位噪声补偿的自适应均衡器 |
| CN102130879A (zh) * | 2011-03-10 | 2011-07-20 | 上海交通大学 | 正交频分复用系统中相位噪声消除方法 |
| CN103338171A (zh) * | 2013-06-27 | 2013-10-02 | 北京大学 | 一种基于频域信道估计的接收端均衡方法和系统 |
-
2014
- 2014-01-07 CN CN201410007166.6A patent/CN103716265B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101394381A (zh) * | 2007-09-19 | 2009-03-25 | 中国科学院微电子研究所 | 一种用于相位噪声补偿的自适应均衡器 |
| CN102130879A (zh) * | 2011-03-10 | 2011-07-20 | 上海交通大学 | 正交频分复用系统中相位噪声消除方法 |
| CN103338171A (zh) * | 2013-06-27 | 2013-10-02 | 北京大学 | 一种基于频域信道估计的接收端均衡方法和系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103716265A (zh) | 2014-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103716265B (zh) | 一种相位噪声补偿抑制改进方法 | |
| CN1802831B (zh) | 用于ofdm信号的自适应相位补偿的方法和装置 | |
| Vihriala et al. | On the waveforms for 5G mobile broadband communications | |
| CN105830380B (zh) | 高效短训练域序列生成方法、信号发送、接收方法及装置 | |
| CN100385824C (zh) | 一种mimo-ofdm系统的自适应信道估计方法 | |
| CN103685096B (zh) | 一种基于最优导频的mimo-ofdm系统信道估计方法 | |
| CN103179062B (zh) | Sc-fde系统低复杂度信道估计下的相位噪声抑制方法 | |
| CN113765581A (zh) | 基于压缩感知与波束对齐的ris快时变信道估计方法 | |
| CN103746947A (zh) | 一种相位噪声的估计方法 | |
| CN108156101B (zh) | 一种mimo-scfde系统联合迭代信道估计与迭代均衡方法 | |
| CN104639490A (zh) | 一种频率相关iq不平衡与信道的联合估计与补偿方法 | |
| Zhang et al. | An adaptive matching pursuit algorithm for sparse channel estimation | |
| CN105915476A (zh) | 一种基于贝叶斯的相噪补偿方法 | |
| CN104486266A (zh) | 一种基于mimo-ofdm系统的信道估计方法及装置 | |
| CN104917711A (zh) | 无线通信系统下的一种相噪补偿改进方法 | |
| CN106332095A (zh) | 基于级联频域均衡的超奈奎斯特传输方法 | |
| CN108768905B (zh) | 一种基于混沌成型和混沌匹配滤波器的无线通信方法 | |
| CN107332606B (zh) | 基于双采样的leo系统差分空时正交频分复用编码方法 | |
| CN104954305A (zh) | 无线通信系统下的一种相位噪声改进估计方法 | |
| CN105991489A (zh) | 利用频域过采样来实现信道均衡的方法 | |
| Dinis et al. | A hybrid ARQ scheme for faster than Nyquist signaling with iterative frequency-domain detection | |
| Darsena et al. | LTV equalization of CPM signals over doubly-selective aeronautical channels | |
| Guangrong et al. | Fractionally spaced equalization algorithms in 60GHz communication system | |
| Yang et al. | Semi-blind channel estimation for OFDM systems | |
| CN103414667B (zh) | 一种基于二维离散导频的ofdm自适应信道估计方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| 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: 20170503 Termination date: 20200107 |