CN103873405B - 一种基于apsk调制信号幅值的恒模盲均衡方法 - Google Patents

一种基于apsk调制信号幅值的恒模盲均衡方法 Download PDF

Info

Publication number
CN103873405B
CN103873405B CN201410073586.4A CN201410073586A CN103873405B CN 103873405 B CN103873405 B CN 103873405B CN 201410073586 A CN201410073586 A CN 201410073586A CN 103873405 B CN103873405 B CN 103873405B
Authority
CN
China
Prior art keywords
mrow
apsk
msub
signal
msup
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.)
Active
Application number
CN201410073586.4A
Other languages
English (en)
Other versions
CN103873405A (zh
Inventor
闫朝星
周三文
卢满宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aerospace Long March Launch Vehicle Technology Co Ltd
Beijing Institute of Telemetry Technology
Original Assignee
Aerospace Long March Launch Vehicle Technology Co Ltd
Beijing Institute of Telemetry Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Aerospace Long March Launch Vehicle Technology Co Ltd, Beijing Institute of Telemetry Technology filed Critical Aerospace Long March Launch Vehicle Technology Co Ltd
Priority to CN201410073586.4A priority Critical patent/CN103873405B/zh
Publication of CN103873405A publication Critical patent/CN103873405A/zh
Application granted granted Critical
Publication of CN103873405B publication Critical patent/CN103873405B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明涉及一种基于APSK调制信号幅值的恒模盲均衡方法,首先初始化长度为L的滤波器系数W(n)为W0(n),对长度为L的输入信号X(n)进行滤波后得到y(n),然后将y(n)与每个星座环的目标函数求差得到|y(n)|2‑Rk,并求得||y(n)|2‑Rk|取最小值时对应的k值,据此得到误差函数e(n),最后,用误差e(n)对滤波器系数W(n)进行更新,并输入新的接收信号,然后进入到下一次迭代,本发明方法改善了幅度相移键控(APSK)调制系统接收机的均衡器收敛后信号的码间干扰性能,大幅提高了信号质量。

Description

一种基于APSK调制信号幅值的恒模盲均衡方法
技术领域
本发明属于数字无线通信传输技术领域,涉及一种卫星移动通信系统中幅度-相移键控调制信号的恒模盲均衡方法,尤其涉及一种基于APSK调制信号幅值的恒模盲均衡方法。
背景技术
在宽带卫星通信系统中,除了传输信道的噪声,信号经历的信道还有地面站设备、星上的功率放大器等。卫星通信系统中的带通滤波器、发射机以及其他通信设备中都不同程度地存在着非线性相位响应。这些非线性响应对于宽带卫星通信中常用的幅度相移键控(APSK)调制方式的影响不可忽略,码间干扰(ISI)被认为是在移动无线通信信道中实现高速率数据传输的主要障碍,而均衡技术正是对付码间干扰的一项技术。
均衡技术能够很大程度上减少码间干扰,改善非线性失真引起的频谱扩展、星座扭曲以及信道畸变等问题,提高系统的传输速率。均衡器的基本结构包括横向滤波器结构、判决反馈结构等。根据所采样的数据不同,均衡的算法还可以分为数据辅助均衡和盲均衡。盲均衡可以分为基于梯度算法的Bussgang盲均衡和基于二阶统计量和高阶统计量的盲均衡等。
由D.N.Godard提出的恒模算法(CMA)均衡方法对相移键控(PSK)调制最佳,设a(n),n=1,2,...为发送的调制符号数据,接收信号x(n)经过均衡后为y(n),CMA算法的代价函数J(n)表示为:J(n)=E[(|y(n)|p-Rp)2],其中,Rp=E[|a(n)|2p]/E[|a(n)|p],p为设计参数。
CMA算法均衡器的滤波器系数W(n)的长度为L,该系数组可表示为W(n)={w0(n),w1(n),...,wL-1(n)}。W(n)的更新方法可表示为:
其中,表示对代价函数求偏导,μ为系数更新的步长,X(n)为L个接收信号的集合,表示为X(n)={x(n-L+1),...,x(n-1),x(n)}。特别的,当p=2代价函数J(n)表示为:J(n)=E[(|y(n)|2-R2)2],其中,
所以,滤波器系数的更新表示为:
CMA算法对于具有恒模的PSK调制时最佳的,但是对正交幅度调制(QAM)信号,CMA算法收敛时的误差信号并不为0。由K.N.Oh提出的多模算法(MMA)均衡方法把QAM信号星座点划分到同相(I)正交(Q)两路等模信号,更适用于QAM调制信号。但是对于宽带卫星通信中常采用的APSK调制,其星座点幅值不具有QAM调制的I、Q路分立的特点,而MMA算法均衡后的码间干扰性能要劣于CMA算法。
发明内容
本发明的目的在于克服现有技术的上述不足,提供一种基于APSK调制信号幅值的恒模盲均衡方法,该方法改善了幅度相移键控(APSK)调制系统接收机的均衡器收敛后信号的码间干扰性能,大幅提高了信号质量。
本发明的上述目的主要是通过如下技术方案予以实现的:
一种基于APSK调制信号幅值的恒模盲均衡方法,包括如下步骤:
步骤(一)、设2M阶APSK调制星座点中,第k圈星座点对应圆环的半径幅值为Ak,k=1,2,…,M-2,且M为正整数,M≥4,第k圈星座点对应的误差函数ek(n)表示为:
ek(n)=y(n)(|y(n)|2-Rk)
则误差函数e(n)表示为:
其中:
y(n)为通过APSK调制系统接收端的均衡器滤波后的信号;
arg{min||y(n)|2-Rk|}表示为||y(n)|2-Rk|值最小时所对应的k值,即kmin
Rk为第k圈星座点对应的目标函数,即:
其中:ak(n)为按照幅值为Ak的圆环上星座点进行映射的信号值;
步骤(二)、对于2M阶APSK调制信号,设a(n),n=1,2,...为APSK调制系统发送的调制符号数据,经过APSK调制系统的传输信道后接收信号为x(n),通过APSK调制系统接收端的均衡器滤波后的信号为y(n),则所述均衡器中滤波器系数的迭代更新方程为:
W(n+1)=W(n)-μe(n)X*(n)
其中:
X(n)为APSK调制系统的接收矢量,即L个接收信号x(n)的集合,表示为X(n)={x(n-L+1),...,x(n-1),x(n)};
W(n)为APSK调制系统接收端的均衡器滤波系数,长度为L,表示为W(n)={w0(n),w1(n),...,wL-1(n)};
μ为系数更新的步长;
L为正整数,L≥1。
在上述基于APSK调制信号幅值的恒模盲均衡方法中,步骤(二)中滤波器系数迭代更新过程中所述APSK调制系统接收端的均衡器滤波系数W(n)的初始值W0(n)的长度为L,L取奇数,表示为W0(n)=[00...010...00],即中间第(L+1)/2个数为1,其他为0。
本发明与现有技术相比具有如下有益效果:
(1)、本发明提出一种全新的基于APSK调制信号幅值的恒模盲均衡方法,通过确定接收信号幅度的平方与目标函数之间的差值绝对值的最小值对应的k值,从而确定误差函数,与现有技术采用固定的误差函数值相比,能够更有效提取非线性信道误差,从而大大改善了幅度相移键控(APSK)调制系统接收机的均衡器收敛后信号的码间干扰性能,提高了信号质量;
(2)、本发明基于APSK调制信号幅值的恒模盲均衡方法比传统算法更充分地利用了幅度相移键控(APSK)调制星座点所包含的幅度信息,能够更有效提取非线性信道误差;
(3)、本发明所设计的均衡方法在接收机中具有易于实现的结构,可以有较低的实现复杂度;且本发明不需要改变接收系统结构,仅需要改进提取误差信息的函数,不改变原系统中的滤波器系数更新方程和滤波过程,不影响整个系统其他单元的正常工作。
附图说明
图1为本发明基于APSK调制信号幅值的恒模盲均衡方法原理图;
图2为本发明16APSK调制的星座点分析图;
图3为本发明均衡前经过非线性信道的16APSK信号星座图,SNR=20dB;
图4为本发明均衡后的16APSK信号星座图,SNR=20dB;
图5为本发明基于APSK调制信号幅值的恒模盲均衡算法的性能,16APSK,SNR=20dB。
具体实施方式
下面结合附图和具体实施实例对本发明作进一步详细的描述。
如图1所示为本发明基于APSK调制信号幅值的恒模盲均衡方法原理图,首先,初始化长度为L的滤波器系数W(n)为W0(n),对长度为L的输入信号X(n)进行滤波后得到y(n)。然后,将y(n)与每个星座环的目标函数求差得到|y(n)|2-Rk,并求得||y(n)|2-Rk|取最小值时对应的k值,e(n)=ek(n),据此得到误差函数e(n)。最后,用误差e(n)对滤波器系数W(n)进行更新,并输入新的接收信号,然后进入到下一次迭代,具体实现方法如下:
步骤(一)、设2M阶APSK调制星座点中,第k圈星座点对应圆环的半径幅值为Ak,k=1,2,…,M-2,且M为正整数,M≥4,第k圈星座点对应的误差函数ek(n)表示为:
ek(n)=y(n)(|y(n)|2-Rk)
则误差函数e(n)表示为:
其中:
y(n)为通过APSK调制系统接收端的均衡器滤波后的信号;
arg{min||y(n)|2-Rk|}表示为||y(n)|2-Rk|值最小时所对应的k值,即kmin
Rk为第k圈星座点对应的目标函数,即:
其中:ak(n)为按照幅值为Ak的圆环上星座点进行映射的信号值。
步骤(二)、对于2M阶APSK调制信号,设a(n),n=1,2,...为发送的调制符号数据,经过APSK调制系统的传输信道后接收信号为x(n),通过APSK调制系统接收端的均衡器滤波后的信号为y(n),则所述均衡器中滤波器系数的迭代更新方程为:
W(n+1)=W(n)-μe(n)X*(n)
其中:
X(n)为APSK调制系统的接收矢量,即L个接收信号x(n)的集合,表示为X(n)={x(n-L+1),...,x(n-1),x(n)};
W(n)为APSK调制系统接收端的均衡器滤波系数,长度为L,表示为W(n)={w0(n),w1(n),...,wL-1(n)};
μ为系数更新的步长;
L为正整数,L≥1。
其中APSK调制系统接收端的均衡器滤波系数W(n)迭代更新过程中初始值W0(n)的长度为L,L取奇数,表示为W0(n)=[00...010...00],即中间第(L+1)/2个数为1,其他为0。
对于2M阶APSK调制信号,取M=4的16APSK为例,如图2所示为本发明16APSK调制的星座点分析图,设a(n),n=1,2,...为发送的调制符号数据,经过信道后接收信号为x(n),通过均衡器滤波后的信号为y(n)。均衡器中滤波器系数的更新方程为:
W(n+1)=W(n)-μe(n)X*(n)
其中:
e(n)为误差函数,X(n)为接收信号,X*(n)为X(n)的共轭,μ为系数更新的步长,W(n)为APSK调制系统接收端的均衡器滤波系数,滤波系数迭代过程中的初始值设为W0(n)。X(n)为L个接收信号x(n)的集合,表示为X(n)={x(n-L+1),...,x(n-1),x(n)};W(n)的长度为L,该系数组可表示为W(n)={w0(n),w1(n),...,wL-1(n)}。
设16APSK调制星座点中,第k圈星座点对应圆环的半径幅值为Ak,k=1,2。则误差函数e(n)可表示为:
ek(n)=y(n)(|y(n)|2-Rk)
则误差函数e(n)可表示为:
其中:
arg{min||y(n)|2-Rk|}表示为||y(n)|2-Rk|值最小时所对应的k值;Rk为第k圈星座点对应的目标函数。
本发明实施例的效果可以通过以下仿真进一步得到说明。
仿真条件如下,随机产生数据,进行16APSK调制,如图2所示,对应于DVB-S2标准定义的信道编码为3/4码率时16APSK调制的内圈星座点:对应于外圈星座点:滤波器系数长度取L=11,初始值可设为为W0(n)=[00000100000],仿真采用的信道响应为:hB=[-0.005-0.004j,0.009+0.030j,-0.024-0.104j,0.854+0.520j,-0.218+0.273j,0.049-0.074j,-0.016+0.020j],SNR=20dB。经过此信道的16APSK调制信号如图3所示(图3为本发明均衡前经过非线性信道的16APSK信号星座图,SNR=20dB)。按照如上所述均衡过程对接收信号进行均衡,其中步长为μ=10-3,得到均衡收敛后的星座点如图4所示(图4为本发明均衡后的16APSK信号星座图,SNR=20dB)。用码间干扰指标定量地分析均衡的性能,其中残留码间干扰(ISI)为滤波器系数与信道响应的函数,
其中,h(i)为非线性信道的系数,w(i)为均衡中的滤波器系数,为两者卷积后模值的最大值。
设计的均衡算法性能如图5所示(图5为本发明基于APSK调制信号幅值的恒模盲均衡算法的性能,16APSK,SNR=20dB)。由于16APSK信号的幅值只有两种情况,故MMA算法的优势不再存在,CMA均衡算法收敛后的ISI性能比MMA算法高约4dB,而本发明所改进的CMA算法收敛后的ISI性能比CMA算法提高大于5dB。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (1)

1.一种基于APSK调制信号幅值的恒模盲均衡方法,其特征在于包括如下步骤:
步骤(一)、设2M阶APSK调制星座点中,第k圈星座点对应圆环的半径幅值为Ak,k=1,2,…,M-2,且M为正整数,M≥4,第k圈星座点对应的误差函数ek(n)表示为:
ek(n)=y(n)(|y(n)|2-Rk)
则误差函数e(n)表示为:
<mrow> <mi>e</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>e</mi> <msub> <mi>k</mi> <mi>min</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>k</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mi>arg</mi> <mo>{</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mo>|</mo> <mo>|</mo> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <msub> <mi>R</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>}</mo> </mrow>
其中:
y(n)为通过APSK调制系统接收端的均衡器滤波后的信号;
arg{min||y(n)|2-Rk|}表示为||y(n)|2-Rk|值最小时所对应的k值,即kmin
Rk为第k圈星座点对应的目标函数,即:
<mrow> <msub> <mi>R</mi> <mi>k</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>E</mi> <mo>&amp;lsqb;</mo> <mo>|</mo> <msub> <mi>a</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>4</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mi>E</mi> <mo>&amp;lsqb;</mo> <mo>|</mo> <msub> <mi>a</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>=</mo> <mo>|</mo> <msub> <mi>A</mi> <mi>k</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>;</mo> </mrow>
其中:ak(n)为按照幅值为Ak的圆环上星座点进行映射的信号值;
步骤(二)、对于2M阶APSK调制信号,设a(n),n=1,2,...为APSK调制系统发送的调制符号数据,经过APSK调制系统的传输信道后接收信号为x(n),通过APSK调制系统接收端的均衡器滤波后的信号为y(n),则所述均衡器中滤波器系数的迭代更新方程为:
W(n+1)=W(n)-μe(n)X*(n)
其中:
X(n)为APSK调制系统的接收矢量,即L个接收信号x(n)的集合,表示为X(n)={x(n-L+1),...,x(n-1),x(n)};X*(n)为X(n)的共轭;
W(n)为APSK调制系统接收端的均衡器滤波系数,长度为L,表示为W(n)={w0(n),w1(n),...,wL-1(n)};
μ为系数更新的步长;
L为正整数,L≥1;
所述步骤(二)中滤波器系数迭代更新过程中所述APSK调制系统接收端的均衡器滤波系数W(n)的初始值W0(n)的长度为L,L取奇数,表示为W0(n)=[00...010...00],即中间第(L+1)/2个数为1,其他为0。
CN201410073586.4A 2014-02-28 2014-02-28 一种基于apsk调制信号幅值的恒模盲均衡方法 Active CN103873405B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410073586.4A CN103873405B (zh) 2014-02-28 2014-02-28 一种基于apsk调制信号幅值的恒模盲均衡方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410073586.4A CN103873405B (zh) 2014-02-28 2014-02-28 一种基于apsk调制信号幅值的恒模盲均衡方法

Publications (2)

Publication Number Publication Date
CN103873405A CN103873405A (zh) 2014-06-18
CN103873405B true CN103873405B (zh) 2017-08-29

Family

ID=50911540

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410073586.4A Active CN103873405B (zh) 2014-02-28 2014-02-28 一种基于apsk调制信号幅值的恒模盲均衡方法

Country Status (1)

Country Link
CN (1) CN103873405B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112565123B (zh) * 2020-11-06 2023-01-10 北京遥测技术研究所 一种联合环路锁定检测的uqpsk信号盲均衡方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101030954A (zh) * 2006-03-01 2007-09-05 澜起科技(上海)有限公司 正交振幅调制(qam)接收器的盲均衡器
CN101854317A (zh) * 2010-05-13 2010-10-06 深圳大学 一种盲均衡器及盲均衡处理方法
WO2012041399A1 (en) * 2010-09-28 2012-04-05 Telefonaktiebolaget L M Ericsson (Publ) Adaptive equaliser with asynchronous detection and inhibit signal generator

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101030954A (zh) * 2006-03-01 2007-09-05 澜起科技(上海)有限公司 正交振幅调制(qam)接收器的盲均衡器
CN101854317A (zh) * 2010-05-13 2010-10-06 深圳大学 一种盲均衡器及盲均衡处理方法
WO2012041399A1 (en) * 2010-09-28 2012-04-05 Telefonaktiebolaget L M Ericsson (Publ) Adaptive equaliser with asynchronous detection and inhibit signal generator

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
适用于高阶QAM信号的水声信道修正盲均衡算法;宁小玲,刘忠,罗亚松,付学志,杨泗杰;《重庆邮电大学学报》;20111031;第23卷(第5期);正文第517页 *

Also Published As

Publication number Publication date
CN103873405A (zh) 2014-06-18

Similar Documents

Publication Publication Date Title
CN103873404B (zh) 高阶正交幅度调制系统中基于i/q路幅值的多模盲均衡方法
CN105245478B (zh) 一种基于qam调制方式的自适应均衡算法
CN103546402B (zh) 一种发送信号的方法、装置和系统
CN103312640A (zh) 一种联合信道估计与iq不平衡补偿的方法
CN109246039A (zh) 一种基于双向时域均衡的软信息迭代接收方法
EP3146690B1 (en) Method and apparatus for recovering time-domain hybrid modulated qam signals
CN107566307A (zh) 盲均衡装置及方法、数据调制系统及方法
CN102437978B (zh) 数字微波均衡方法和装置
CN106230757A (zh) 基于预编码的fbmc系统实数域均衡方法
CN106253922B (zh) 一种基于盲均衡算法的pcm/fm的抗多径干扰遥测接收方法
CN106161304A (zh) 一种联合信道估计的发射端iq不平衡补偿方法
CN103873405B (zh) 一种基于apsk调制信号幅值的恒模盲均衡方法
JP2017092739A (ja) 伝送装置、受信装置、伝送システム及び変調方法
JPH06500914A (ja) 信号点、等化器利得などを調整する装置および方法
CN106656879B (zh) 一种高速高阶变步长自适应均衡方法
Maneekut et al. Hybrid probabilistic and geometric shaping for 64-QAM optical fiber transmission with maximum aposterior probability detection
CN107171990A (zh) 基于分数间隔与迭代算法的单载波信号时域均衡方法
CN102571663B (zh) 微波通信数据传输方法及装置
Yin The CPDA detector for the MIMO OCDM system
Makki et al. Comparative study on LMS and CMA channel equalization
CN112565123B (zh) 一种联合环路锁定检测的uqpsk信号盲均衡方法
JP2002280941A (ja) 適応等化器、受信装置およびタップ係数算出方法
WO2016095261A1 (zh) 一种分布式iq不平衡估计与抑制方法
CN105991488A (zh) 应用于ftn调制中的减小状态数的维特比解调方法
KR100579889B1 (ko) 부호 순열 필터링을 이용한 적응 등화 장치 및 방법

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