基于四元数的直接判决并行恒模均衡算法

李森; 徐明莹; 张璐; 邓明旭

doi:10.11999/JEIT221413

基于四元数的直接判决并行恒模均衡算法

doi: 10.11999/JEIT221413

大连海事大学信息科学技术学院大连 116026

基金项目: 国家重点研发计划(2019YFE0111600)，国家自然科学基金(51939001, 61971083)，辽宁省振兴人才计划(XLYC2002078)，鹏城实验室重点项目(PCL2021A03-1)

详细信息

作者简介:
李森：女，教授，研究方向为非高斯信号处理、通信信号处理、统计信号处理

徐明莹：女，硕士生，研究方向为4元数理论、信道均衡算法

张璐：女，硕士生，研究方向为信道均衡算法、阵列信号处理

邓明旭：女，硕士生，研究方向为信道均衡算法、频谱感知

通讯作者:
李森　listen@dlmu.edu.cn

中图分类号: TN911.72
计量
- 文章访问数: 343
- HTML全文浏览量: 217
- PDF下载量: 41
- 被引次数: 1
出版历程
- 收稿日期: 2022-11-09
- 修回日期: 2023-07-05
- 网络出版日期: 2023-07-13
- 刊出日期: 2023-10-31

Concurrent Decision Directed and Constant Modulus Equalization Algorithm Based on Quaternion

College of Information Science and Technology, Dalian Maritime University, Dalian 116026, China

Funds: The National Key Research and Development Program of China (2019YFE0111600), The National Natural Science Foundation of China (51939001,61971083), Liao Ning Revitalization Talents Program ( XLYC2002078) , The Major Key Project of PCL(PCL2021A03-1)

摘要

摘要: 近年来4元数理论成为各界学者研究的热点并被应用到许多领域。该文基于4元数自适应滤波算法对正交极化信道均衡问题进行研究。为了解决4元数恒模(QCMA)算法的相位模糊问题，该文将QCMA算法与4元数最小均方算法(QLMS)相结合提出一种4元数直接判决并行恒模(QCMA+DD-QLMS) 算法。基于广义哈密顿实演算(GHR)的梯度运算规则对新的算法做了理论推导并进行MATLAB实验仿真，仿真结果表明该文所提算法不但能够解决QCMA算法相位模糊问题，还具有更小的稳态均方误差(MSE)。
- 4元数 /
- 正交极化信道均衡 /
- 4元数最小均方算法 /
- 4元数恒模算法 /
- 4元数直接判决并行恒模算法
Abstract: In recent years, quaternion theory has become a research hotspot for scholars and has been applied to many fields. In this paper, the quadrature polarization channel equalization problem is studied based on quaternion adaptive filtering algorithm. In order to solve the phase ambiguity problem of Quaternion Constant Modulus Algorithm (QCMA), a concurrent quaternion Direct Decision constant modulus algorithm (QCMA+DD-QLMS) is proposed by combining QCMA algorithm with Quaternion Least Mean Square (QLMS) algorithm. Based on the gradient operation rules of Generalized Hamilton-Real (GHR), the new algorithm is theoretically deduced and simulated by MATLAB. The simulation results show that the algorithm proposed in this paper can not only solve the phase ambiguity problem of the QCMA, but also has smaller steady-state Mean Square Error (MSE).

HTML全文

1. 引言

提高无线通信系统的容量一直是无线通信领域研究的焦点。传统的极化分集方案是指在发射端和接收端同时使用一对具有正交极化方向的天线，以相同的频率不同的极化方向传输同样的信息，利用极化方向相互正交的天线发出的信号呈现出的不相关衰落特性来减轻多径效应，从而在发射机和接收机之间保持可靠的通信链路和可接受的服务质量。但是这种传统的极化分集方案并不能实现对频谱资源的有效利用，为了进一步提高频谱利用率， Isaeva等人^[1]提出了在相互正交的两个极化天线上分别传输不同信息的4维调制(Quaternion Orthogonal Design, QOD) 技术。4维信号在复杂的多径、反射、折射等传播过程中其极化方向的变化会导致两个不同极化方向的接收天线之间产生干扰，因此在应用QOD技术提高通信系统频谱利用率的同时，也需要设计高效的信号处理算法来实现4维信号的均衡和干扰抑制。

4元数是爱尔兰数学家Hamilton在1843年首先提出的，因其本身固有的结构特点成为处理3维或4维信号相关问题的一个重要工具。例如最近的一些研究成果成功地将4元数应用在彩色图像处理^[2]、风速预测^[3]和矢量传感器阵列信号处理^[4]等方面。近年来，人们对4元数神经网络也越来越感兴趣^[5]，提出了4元数卷积神经网络和4元数递归神经网络等^[6]。最早将4元数应用到自适应滤波算法中的是Took和Mandic^[7]，他们提出了4元数最小均方误差(Quaternion Least Mean Square, QLMS)算法。针对4元数不同的梯度求导规则，研究人员又陆续提出了QLMS算法的其他表达形式以及改进算法。例如：Xu等人^[8]提出了基于广义哈密顿实演算(Generalized Hamilton-Real, GHR)梯度求解规则的QLMS算法；并在此基础上进一步提出了能够统一处理4元数圆和非圆信号的自适应滤波宽线性(Widely Linear, WL)QLMS算法^[9]。为了解决QLMS算法收敛速度慢的问题， Mengüç等人^[10]用4元数误差信号的最小峰度函数代替4元数最小均方误差函数，提出了4元数最小平均峰度(Quaternion Least-Mean Kurtosis, QLMK)算法并进一步提出了宽线性4元数最小平均峰度(Widely Linear Quaternion Least-Mean Kurtosis WL-QLMK)算法^[11]。文献[12]使用误差熵作为代价函数，基于信息理论学习技术提出了4元数核最小误差熵(Quaternion Kernel Minimum Error Entropy, QKMEE)自适应算法，解决了QLMS算法在非高斯脉冲噪声环境下性能下降的问题。

以上的算法都是针对单输入单输出(Single Input Single Output, SISO)系统而言的，文献[13]探讨了利用多对双极化天线同时发射和接收垂直与水平极化电磁波的双极化多输入多输出(Multiple Input Multiple Output, MIMO)系统，并使用4元数描述了双极化MIMO系统的信道模型，相较于双极化SISO传输系统，双极化MIMO系统在信道容量及信号传输效率方面有很大的提高。

4元数理论和4元数自适应滤波算法的研究进展为解决4维调制信号的均衡和干扰抑制问题奠定了理论基础。文献[14]用4元数有限冲击响应滤波器来模拟正交极化信道对传输信号产生的影响。文献[15]将2维调制的4QAM(Quadrature Amplitude Modulation, QAM)方案扩展到正交极化天线系统中的QOD的16Q²AM(Quaternion QAM, Q²AM)方案，并利用4元数4个部分之间彼此相互正交的特性将16Q²AM调制信号表示为4元数形式，在此基础上利用QLMS算法提出了一种基于训练信号的均衡方案解决了4元数有限冲击响应信道引起的码间干扰问题。文献[16]又进一步提出了一种基于4元数的恒模盲均衡(Quaternion Constant Modulus Algorithm, QCMA)算法。仿真结果表明QLMS和QCMA两种算法都较好地实现了正交极化信道的均衡。由于QLMS算法是基于训练信号的，因此会降低通信系统的传输效率。QCMA盲均衡算法虽然可以提高通信系统的传输效率，但是由于QCMA算法的代价函数与相位无关导致均衡后的星座图会产生随机角度的相位偏移，使均衡器输出的稳态均方误差(Mean Square Error, MSE)过大。为了解决这个问题，参考在2维调制方案中的并行盲均衡算法^[17-20]的思想，本文将直接判决(Direct Decision, DD)^[21]的QLMS算法和QCMA算法相结合，提出了一种改进的4元数直接判决并行恒模(QCMA+DD-QLMS)算法，利用QCMA算法能使眼图尽快张开和直接判决的QLMS算法收敛速度快、MSE小的优点，在实现均衡输出无相位模糊的同时还具有更小的MSE。

2. 4元数基本理论和4维调制(QOD)技术

2.1 4元数基本理论

4元数 $q \in \mathbb{Q}$ ( $\mathbb{Q}$ 表示4元数域)的通用定义为

$q = {q_{\rm{a}}} + {q_{\rm{b}}}{\boldsymbol{i}} + {q_{\rm{c}}}{\boldsymbol{j}} + {q_{\rm{d}}}{\boldsymbol{k}},{\text{ }}{q^ * } = {q_{\rm{a}}} - {q_{\rm{b}}}{\boldsymbol{i}} - {q_{\rm{c}}}{\boldsymbol{j}} - {q_{\rm{d}}}{\boldsymbol{k}}$

(1)

其中， ${q^ * }$ 代表 $q$ 的共轭， ${q_{\rm{a}}},{q_{\rm{b}}},{q_{\rm{c}}},{q_{\rm{d}}} \in \mathbb{R}$ ( $\mathbb{R}$ 表示实数域)。 ${\boldsymbol{ i}}$ , ${\boldsymbol{ j}}$ , ${\boldsymbol{k}}$ 是虚部单位向量，它们之间满足式(2)的关系

${\boldsymbol{ij}} = {\boldsymbol{k}},{\text{ }}{\boldsymbol{jk}} = {\boldsymbol{i}},{\text{ }}{\boldsymbol{ki}} = {\boldsymbol{j}},{\text{ }}{\boldsymbol{ji}} = - {\boldsymbol{k}},{\text{ }}{\boldsymbol{kj}} = - {\boldsymbol{i}},{\text{ }}{\boldsymbol{ik}} = - {\boldsymbol{j}}$

(2)

4元数在理论上可以理解为对复数的扩展，所以它与复数有着相似的运算法则，但4元数乘法不满足交换律，即： ${q_1}{q_2} \ne {q_2}{q_1}$ 。不过4元数的共轭满足 $q{q^ * } = {q^ * }q$ ，因此，4元数的模值表达式为： ${\left| q \right|^2} = q{q^ * }$ 。

对于4元数的梯度求导，本文采用基于GHR的运算规则，GHR演算简化了4元数的乘积规则和链式规则，能够高效计算4元数的导数，GHR导数分为左导数和右导数，本文主要关注左导数，它的求导规则为^[8]

$\left. \begin{gathered} \frac{{\partial f}}{{\partial q}} = \frac{1}{4}\left( {\frac{{\partial f}}{{\partial {q_{\rm{a}}}}} - \frac{{\partial f}}{{\partial {q_{\rm{b}}}}}{\boldsymbol{i}} - \frac{{\partial f}}{{\partial {q_{\rm{c}}}}}{\boldsymbol{j}} - \frac{{\partial f}}{{\partial {q_{\rm{d}}}}}{\boldsymbol{k}}} \right) \\ \frac{{\partial f}}{{\partial {q^ * }}} = \frac{1}{4}\left( {\frac{{\partial f}}{{\partial {q_{\rm{a}}}}} + \frac{{\partial f}}{{\partial {q_{\rm{b}}}}}{\boldsymbol{i}} + \frac{{\partial f}}{{\partial {q_{\rm{c}}}}}{\boldsymbol{j}} + \frac{{\partial f}}{{\partial {q_{\rm{d}}}}}{\boldsymbol{k}}} \right) \end{gathered} \right\}$

(3)

其中， $f$ 是关于4元数 $q$ 的4元数函数， ${{\partial f}}/{{\partial {q_{\rm{a}}}}}, {{\partial f}}/{{\partial {q_{\rm{b}}}}}, {{\partial f}}/{{\partial {q_{\rm{c}}}}},{{\partial f}}/{{\partial {q_{\rm{d}}}}}$ 是 $f$ 对 $q$ 的每个部的偏导。

特殊地，对于共轭的导数为

$\frac{{\partial q}}{{\partial q}} = \frac{{\partial {q^ * }}}{{\partial {q^ * }}} = 1{\kern 1pt} ,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \frac{{\partial q}}{{\partial {q^ * }}} = \frac{{\partial {q^ * }}}{{\partial q}} = - \frac{1}{2}$

(4)

设 $g$ 也是关于 $q$ 的4元数函数，对于4元数函数相乘时的求导GHR演算中引入了一种新的求导规则^[9]

$\frac{{\partial \left( {fg} \right)}}{{\partial q}} = f\frac{{\partial g}}{{\partial q}} + \frac{{\partial f}}{{\partial {q^g}}}g{\text{ }}\frac{{\partial \left( {fg} \right)}}{{\partial {q^ * }}} = f\frac{{\partial g}}{{\partial {q^ * }}} + \frac{{\partial f}}{{\partial {q^{g * }}}}g$

(5)

其中 ${{\partial f}}/{{\partial {q^g}}}$ , ${{\partial f}}/{{\partial {q^{g * }}}}$ 也是一种GHR左导数，但只对某些特殊的函数形式才有定义，具体内容可以参考文献[9]中给出的表1。

2.2 4维调制(QOD)技术

为了进一步提高频谱利用率，在双极化天线的基础上提出了QOD技术，信号通过QOD技术调制后经由正交极化天线发射将产生4维信号。在2维调制技术中基带信号可以采用复数表示，例如4QAM调制时基带信号的实部和虚部可以分别取1或–1，构成4种组合；16QAM调制时基带信号的实部和虚部可以分别取1,–1,3,–3中的一个值，构成16种组合。在QOD技术中，可以用4元数来表示4维基带信号^[15]。例如：4元数的4个部分别取1或–1，则构成16种不同的组合，称这种调制方式为16Q²AM；同理，4元数的4个部分分别取1,–1,3,–3中的一个值，则一共有256种组合，称这种调制方式为256Q²AM。

3. 4元数自适应均衡算法

4元数值无线通信系统模型如所示，4维信号 $s\left( n \right)$ 经4元数有限冲击响应信道 $h\left( n \right)$ 传输后在均衡器输入端和4元数加性噪声 $n\left( n \right)$ 一起进入到4元数均衡器中。4元数均衡器的输入信号矢量为 ${\boldsymbol{x}}\left( n \right) = {\left[ {x\left( n \right),x\left( {n - 1} \right), \cdots ,x\left( {n - L + 1} \right)} \right]^{\rm T}}$ ，其中 $x\left( n \right)$ 可以表示为

图 1 4元数值无线通信系统模型

下载: 全尺寸图片幻灯片

$x\left( n \right) = h\left( n \right) * s\left( n \right) + n\left( n \right)$

(6)

其中 $*$ 表示卷积运算。均衡器的输出为

$y\left( n \right) = {{\boldsymbol{w}}^{\rm T}}\left( n \right){\boldsymbol{x}}\left( n \right)$

(7)

其中， ${\boldsymbol{w}}\left( n \right) = {\left[ {{w_0}\left( n \right),{w_1}\left( n \right), \cdots ,{w_{L - 1}}\left( n \right)} \right]^{\rm T}}$ 为4元数值的均衡器权矢量，L为均衡器的长度。

3.1 4元数最小均方误差(QLMS)算法

QLMS算法的代价函数为

${J_{{\rm{QLMS}}}}\left( n \right) = \min {\rm{E}}\left[ {{{\left| {d\left( n \right) - y\left( n \right)} \right|}^2}} \right]$

(8)

其中， $d\left( n \right)$ 为均衡器期望输出的4维信号，也就是已知的训练序列。

QLMS算法权矢量 ${\boldsymbol{w}}\left( n \right)$ 的更新公式为^[14]

${\boldsymbol{w}}\left( {n + 1} \right) = {\boldsymbol{w}}\left( n \right) - \mu \left( {d\left( n \right) - y\left( n \right)} \right){{\boldsymbol{x}}^ * }\left( n \right)$

(9)

其中， $\mu > 0$ 为QLMS算法迭代步长。

3.2 4元数恒模盲均衡(QCMA)算法

QCMA算法的代价函数为

${J_{{\rm{QCMA}}}}\left( n \right) = \min {\rm{E}}\left[ {{{\left( {{{\left| {y\left( n \right)} \right|}^2} - {R_{\rm{c}}}} \right)}^2}} \right]$

(10)

其中， ${R_{\rm{c}}} = {{{\rm{E}}\left[ {{{\left| {s\left( n \right)} \right|}^4}} \right]} /{{\rm{E}}\left[ {{{\left| {s\left( n \right)} \right|}^2}} \right]}}$ 是一个正实数。QCMA算法的权矢量更新公式为^[15]

${{\boldsymbol{w}}_{\rm{c}}}\left( {n + 1} \right) = {{\boldsymbol{w}}_{\rm{c}}}\left( n \right) - {\mu _{\rm{c}}}\left[ {{{\left| {y\left( n \right)} \right|}^2} - {R_{\rm{c}}}} \right]y\left( n \right){{\boldsymbol{x}}^ * }\left( n \right)$

(11)

其中， ${{\boldsymbol{w}}_{\rm{c}}}\left( n \right) = {\left[ {{w_{{\rm{c}},0}}\left( n \right),{w_{{\rm{c}},1}}\left( n \right), \cdots ,{w_{{\rm{c}},L - 1}}\left( n \right)} \right]^{\rm T}}$ 为QCMA算法的权矢量， ${\mu _{\rm{c}}} > 0$ 为QCMA算法的迭代步长。

迭代步长大小决定了每次更新后均衡器抽头系数的变化程度，进而影响算法的收敛速度和稳定性。具体而言，在一定范围内，算法的迭代步长越大，抽头系数更新的幅度就越大，算法收敛速度也会更快。然而，如果步长过大则可能导致算法不稳定，甚至出现震荡现象。相反，如果步长过小则算法收敛速度会变慢。因此需要在保证收敛的前提下选择合适的迭代步长。

3.3 4元数直接判决并行恒模均衡算法

QLMS算法在实际应用中，必须有训练序列才能完成均衡，会降低通信系统的传输效率。QCMA盲均衡算法虽然不需要训练序列，但是它的均衡结果会产生随机角度的相位偏移，导致MSE较大。为了克服QCMA算法稳态误差大和QLMS算法需要训练序列的问题，本文将基于判决引导的最小均方误差(Decision-Directed QLMS, DD-QLMS)算法与QCMA算法相结合提出了一种4元数直接判决并行恒模均衡(QCMA+DD-QLMS)算法，首先利用QCMA算法张开眼图，然后利用DD-QLMS算法进一步减小MSE。4元数直接判决并行恒模均衡器的原理图如所示，该方案由一个QCMA均衡器和一个同时工作的DD-QLMS均衡器组成，设 ${{\boldsymbol{w}}_{\rm{d}}}\left( n \right) = {\left[ {{w_{{\rm{d}},0}}\left( n \right),{w_{{\rm{d}},1}}\left( n \right), \cdots ,{w_{{\rm{d}},L - 1}}\left( n \right)} \right]^{\rm T}}$ 为 DD-QLMS算法的权矢量，则直接判决并行恒模均衡器的权矢量可以表示为 ${\boldsymbol{w}}\left( n \right) = {{\boldsymbol{w}}_{\rm{c}}}\left( n \right) + {{\boldsymbol{w}}_{\rm{d}}}\left( n \right)$ ，均衡器的输出为

图 2 4 元数直接判决并行恒模均衡算法原理图

下载: 全尺寸图片幻灯片

${y_1}\left( n \right) = {{\boldsymbol{w}}_{\rm{c}}^{\rm T}}{\left( n \right)}{\boldsymbol{x}}\left( n \right) + {{\boldsymbol{w}}_{\rm{d}}^{\rm T}}{\left( n \right)}{\boldsymbol{x}}\left( n \right)$

(12)

其中，QCMA算法权矢量 ${{\boldsymbol{w}}_{\rm{c}}}\left( n \right)$ 的更新方程按式(11)计算。

在QLMS算法中均衡器期望输出信号 $d\left( n \right)$ 为已知的训练序列，而在DD-QLMS算法中均衡器期望输出信号 $d\left( n \right)$ 为对并行均衡器输出信号直接判决后的信号 $g\left[ {{y_1}\left( n \right)} \right] = \mathop {\arg \min }\limits_{{S_k} \in S} \left[ {{{\left| {{y_1}\left( n \right) - {s_k}} \right|}^2}} \right]$ ，其中 $S$ 为4维调制信号的星座图集合， $g\left[ {{y_1}\left( n \right)} \right]$ 表示在输入信号星座图集合中找到与均衡器输出 ${y_1}\left( n \right)$ 最接近的信号。因此DD-QLMS算法的代价函数定义为

${J_{{\rm{DD}} {\text{-}} {\rm{QLMS}}}}\left( n \right) = \min {\rm{E}}\left[ {{{\left| {g\left[ {{y_1}\left( n \right)} \right] - {y_1}\left( n \right)} \right|}^2}} \right]$

(13)

则其权矢量 ${{\boldsymbol{w}}_{\rm{d}}}\left( n \right)$ 的更新方程

${{\boldsymbol{w}}_{\rm{d}}}\left( {n + 1} \right) = {{\boldsymbol{w}}_{\rm{d}}}\left( n \right) - {\mu _{\rm{d}}}\left\{ {g\left[ {{y_1}\left( n \right)} \right] - {y_1}\left( n \right)} \right\}{{\boldsymbol{x}}^ * }\left( n \right)$

(14)

${\mu _{\rm{d}}} > 0$ 为DD-QLMS算法迭代步长。由于DD-QLMS算法中的期望信号为并行均衡器输出的判决信号，为了保证算法的有效性在QCMA+DD-QLMS并行算法中上述权矢量的更新过程只有在利用QCMA算法将眼图打开后才工作。定义利用QCMA算法进行权矢量更新后整个均衡器的输出为

${\hat y_1}\left( n \right) = {\boldsymbol{w}}_{\rm{c}}^{\rm T}\left( {n + 1} \right){\boldsymbol{x}}\left( n \right) + {\boldsymbol{w}}_{\rm{d}}^{\rm T}\left( n \right){\boldsymbol{x}}\left( n \right)$

(15)

通过判断QCMA算法权矢量更新前后均衡器的判决输出是否相同来判断眼图是否打开，即

$\begin{split} \varepsilon \left(m\right)&=\varepsilon \left\{g\left[{\widehat{y}}_{1}\left(n\right)\right]-g\left[{y}_{1}\left(n\right)\right]\right\}\\ & =\left\{\begin{array}{l}1,m=0+0_i+0_j+0_k\text{，}眼图打开\\ 0,m\ne 0+0_i+0_j+0_k\text{，}眼图未打开\end{array}\right. \end{split}$

(16)

如果眼图没有打开则继续利用QCMA算法来打开眼图，DD-QLMS算法不发挥作用；如果眼图已经打开则进一步利用DD-QLMS算法减小均衡器的输出误差。因此最终在QCMA+DD-QLMS并行算法中DD-QLMS算法的权矢量更新公式为

$\begin{split} {{\boldsymbol{w}}_{\rm{d}}}\left( {n + 1} \right) = & {{\boldsymbol{w}}_{\rm{d}}}\left( n \right) - {\mu _{\rm{d}}}\varepsilon \left\{ {g\left[ {{{\hat y}_1}\left( n \right)} \right] - g\left[ {{y_1}\left( n \right)} \right]} \right\}\\ & \cdot\left\{ {g\left[ {{y_1}\left( n \right)} \right] - {y_1}\left( n \right)} \right\}{x^ * }\left( n \right)\\[-10pt] \end{split}$

(17)

4. 仿真实验

设信道脉冲响应 ${\boldsymbol{h}}\left( n \right)$ 为4抽头的4元数FIR滤波器，其系数的具体取值如表1所示^[16]。

表 1 4-抽头4元数信道脉冲响应

抽头数	实部	i部	j部	k部
0	–1.1696	–0.1135	0.0958	–0.9583
1	0.4989	–0.3639	–0.2515	–0.6047
2	1.0571	–0.3791	1.4957	0.5183
3	0.4911	0.8653	0.1552	–0.8966

下载: 导出CSV

| 显示表格

MSE是衡量盲均衡算法的一个重要性能指标，MSE的定义为 ${\rm{MSE}} = {{\displaystyle\sum\nolimits_{l = 1}^m {{{\left( {g\left[ {{y_l}\left( n \right)} \right] - {y_l}\left( n \right)} \right)}^2}} }/m}$ ，其中 $m$ 代表Mentecarlo运行次数，在本文中取 $m = 300$ 。由于4维调制信号有4维坐标，无法直接画出4维调制信号的星座图，所以本文中以2维星座图或3维星座图的形式来表示4维调制信号。当以2维星座图的形式来表示4维调制信号时共有两个2维星座图，1个2维星座图表示 ${\boldsymbol{i}}$ 和实部，1个2维星座图表示 ${\boldsymbol{j}}$ 和 ${\boldsymbol{k}}$ 部。当以3维星座图的形式来表示4维调制信号时共有4张图，每张图中包括4元数4个部中的任意3个部。

实验1　输入信号为16Q²AM信号，噪声采用4元数值加性高斯白噪声，信噪比设为25 dB，均衡器长度 $L = 15$ ，采用中心抽头系数初始化的方法。 QCMA算法的迭代步长设为 ${\mu _{\rm{c}}} = 0.001\;3$ ^[16]， QCMA+DD-QLMS算法的迭代步长首先设为情形1(C1)中的参数： ${\mu _{\rm{c}}} = 0.001\;3$ , ${\mu _{\rm{d}}} = 0.001\;2$ 。对QCMA和QCMA+DD-QLMS算法进行实验仿真。均衡前的星座图如图3所示，各个极化方向上的信号已经完全无法分辨。图4和图5分别给出了利用QCMA和QCMA+DD-QLMS算法均衡后的2维星座图，图6和图7分别给出了利用QCMA和QCMA+DD-QLMS算法均衡后3维星座图。从均衡后的星座图中可以看出，由于QCMA算法的代价函数中只考虑了信号的模值，而没有考虑到信号的相位信息，因此其均衡后的星座图出现了明显的相位偏移现象，而QCMA+DD-QLMS并行算法中先用QCMA算法来张开眼图，然后在DD-QLMS算法中利用并行均衡器输出的直接判决信号作为参考信号找出输出信号中与之最接近的点，既考虑了幅值又考虑了相位，因此均衡后的星座图没有出现相位模糊。

图 3 16Q²AM调制时均衡前3维星座图

下载: 全尺寸图片幻灯片

图 4 16Q²AM 调制时 QCMA 算法均衡后 2 维星座图

下载: 全尺寸图片幻灯片

图 5 16Q²AM 调制时 QCMA+DD-QLMS 算法均衡后 2 维星座图

下载: 全尺寸图片幻灯片

图 6 16Q²AM 调制时 QCMA 算法均衡后 3 维星座图

下载: 全尺寸图片幻灯片

图 7 16Q²AM 调制时 QCMA+DD-QLMS 算法均衡后 3 维星座图

下载: 全尺寸图片幻灯片

为了验证步长对算法收敛性及稳态性能的影响，在中除了给出上述两种情况下的MSE收敛曲线还增加了QCMA+DD-QLMS算法在情形2(C2) ${\mu _{\rm{c}}} = 0.002$ , ${\mu _{\rm{d}}} = 0.001\;2$ 时的收敛曲线。从图8的对比结果图可看出，对比QCMA算法，QCMA+DD-QLMS算法具有更小的稳态MSE。而对于具有不同迭代步长的QCMA+DD-QLMS算法，迭代步长与算法的收敛速度是呈正向关系的，即迭代步长越大算法的收敛速度越快；而迭代步长与算法的稳态性能是呈反向关系的，即迭代步长越大算法的稳态性能越差、稳态MSE越大。

图 8 16Q²AM 调制时QCMA和 QCMA+DD-QLMS算法MSE收敛曲线对比图

下载: 全尺寸图片幻灯片

实验2　输入信号为256Q²AM信号，噪声同样采用4元加性高斯白噪声，信噪比设为45 dB，两个算法的迭代步长分别设为 ${\mu _{\rm{c}}} = 0.000\;005$ , ${\mu _{\rm{d}}} = 0.000\;08$ ，均衡器长度 $L = 15$ ，采用中心抽头系数初始化的方法，对QCMA算法和QCMA+DD-QLMS进行实验仿真。均衡前的星座图如图9所示，各个极化方向上的信号同样也是完全无法分辨。图10和图11分别给出了利用QCMA和QCMA+DD-QLMS算法均衡后的2维星座图，图12和图13分别给出了利用QCMA和QCMA+DD-QLMS算法均衡后的3维星座图，图14给出了两种算法的MSE收敛曲线对比图。

图 9 256Q²AM调制时均衡前3维星座图

下载: 全尺寸图片幻灯片

图 10 256Q²AM调制时QCMA算法均衡后2维星座图

下载: 全尺寸图片幻灯片

图 11 256Q²AM调制时QCMA+DD-QLMS算法均衡后2维星座图

下载: 全尺寸图片幻灯片

图 12 256Q²AM调制时QCMA算法均衡后3维星座图

下载: 全尺寸图片幻灯片

图 13 256Q²AM调制时QCMA+DD-QLMS算法均衡后3维星座图

下载: 全尺寸图片幻灯片

图 14 256Q²AM调制时QCMA和QCMA+DD-QLMS算法MSE收敛曲线对比图

下载: 全尺寸图片幻灯片

从仿真结果的均衡后星座图可看出，两种算法均完成了均衡任务，但是QCMA算法均衡后的星座图出现明显的相位偏移现象，收敛后的MSE也较大，而QCMA+DD-QLMS算法均衡后的信号都聚集在各自星座点的周围，并没有出现相位模糊的现象且得到了更小的稳态MSE。但QCMA+DD-QLMS算法是将QCMA算法与DD-QLMS算法并行操作，很显然它的计算复杂度要高于QCMA算法。

5. 结论

QCMA算法虽然实现简单，计算复杂度低，但是均衡后会出现相位模糊现象且MSE较大。为了解决这个问题，本文提出一种基于QCMA算法和DD-QLMS算法相结合的QCMA+DD-QLMS并行算法。所提出的算法充分利用了QCMA算法能较好地打开眼图的特性和DD-QLMS算法MSE小的优点。对QCMA和QCMA+DD-QLMS算法进行了仿真验证，从均衡后的星座图及算法的MSE收敛曲线来看，虽然QCMA+DD-QLMS算法较QCMA算法的计算复杂度高，但是QCMA+DD-QLMS算法很好地解决了QCMA算法的相位模糊问题，使得均衡后的输出具有更小的MSE。

图 1 4元数值无线通信系统模型

下载: 全尺寸图片幻灯片

图 2 4 元数直接判决并行恒模均衡算法原理图

下载: 全尺寸图片幻灯片

图 3 16Q²AM调制时均衡前3维星座图

下载: 全尺寸图片幻灯片

图 4 16Q²AM 调制时 QCMA 算法均衡后 2 维星座图

下载: 全尺寸图片幻灯片

图 5 16Q²AM 调制时 QCMA+DD-QLMS 算法均衡后 2 维星座图

下载: 全尺寸图片幻灯片

图 6 16Q²AM 调制时 QCMA 算法均衡后 3 维星座图

下载: 全尺寸图片幻灯片

图 7 16Q²AM 调制时 QCMA+DD-QLMS 算法均衡后 3 维星座图

下载: 全尺寸图片幻灯片

图 8 16Q²AM 调制时QCMA和 QCMA+DD-QLMS算法MSE收敛曲线对比图

下载: 全尺寸图片幻灯片

图 9 256Q²AM调制时均衡前3维星座图

下载: 全尺寸图片幻灯片

图 10 256Q²AM调制时QCMA算法均衡后2维星座图

下载: 全尺寸图片幻灯片

图 11 256Q²AM调制时QCMA+DD-QLMS算法均衡后2维星座图

下载: 全尺寸图片幻灯片

图 12 256Q²AM调制时QCMA算法均衡后3维星座图

下载: 全尺寸图片幻灯片

图 13 256Q²AM调制时QCMA+DD-QLMS算法均衡后3维星座图

下载: 全尺寸图片幻灯片

图 14 256Q²AM调制时QCMA和QCMA+DD-QLMS算法MSE收敛曲线对比图

下载: 全尺寸图片幻灯片

表 1 4-抽头4元数信道脉冲响应

抽头数实部 i部 j部 k部

0 –1.1696 –0.1135 0.0958 –0.9583
1 0.4989 –0.3639 –0.2515 –0.6047
2 1.0571 –0.3791 1.4957 0.5183
3 0.4911 0.8653 0.1552 –0.8966

下载: 导出CSV

参考文献(21)

[1]	ISAEVA O M and SARYTCHEV V A. Quaternion presentations polarization state[C]. The 2nd Topical Symposium on Combined Optical-Microwave Earth and Atmosphere Sensing, Atlanta, US, 1995: 195–196.
[2]	周珂, 吴成茂, 李昌兴. 基于四元数傅里叶变换的盲彩色图像质量评价[J]. 激光与光电子学进展, 2020, 57(18): 181021. doi: 10.3788/LOP57.181021 ZHOU Ke, WU Chengmao, and LI Changxing. Quality assessment of blind color images using quaternion Fourier transform[J]. Laser &Optoelectronics Progress, 2020, 57(18): 181021. doi: 10.3788/LOP57.181021
[3]	SAOUD L S, AL-MARZOUQI H, and DERICHE M. Wind speed forecasting using the stationary wavelet transform and quaternion adaptive-gradient methods[J]. IEEE Access, 2021, 9: 127356–127367. doi: 10.1109/ACCESS.2021.3111667
[4]	JIANG Mengdi, LIU Wei, and LI Yi. Adaptive beamforming for vector-sensor arrays based on a reweighted zero-attracting quaternion-valued LMS algorithm[J]. IEEE Transactions on Circuits and Systems II:Express Briefs, 2016, 63(3): 274–278. doi: 10.1109/TCSII.2015.2482464
[5]	GAUDET C J and MAIDA A S. Deep quaternion networks[C]. 2018 International Joint Conference on Neural Networks (IJCNN), Rio de Janeiro, Brazil, 2018: 1–8.
[6]	田新宇. 基于四元数的神经网络算法研究[D]. [硕士论文], 电子科技大学, 2022. TIAN Xinyu. Research on neural network algorithm based on quaternion[D]. [Master dissertation], University of Electronic Science and Technology of China, 2022.
[7]	TOOK C C and MANDIC D P. The quaternion LMS algorithm for adaptive filtering of hypercomplex processes[J]. IEEE Transactions on Signal Processing, 2009, 57(4): 1316–1327. doi: 10.1109/TSP.2008.2010600
[8]	XU Dongpo, JAHANCHAHI C, TOOK C C, et al. Quaternion derivatives: The GHR calculus[EB/OL]. https://arxiv.org/abs/1409.8168, 2014.
[9]	MENGÜÇ E C. Novel quaternion-valued least-mean kurtosis adaptive filtering algorithm based on the GHR calculus[J]. IET Signal Processing, 2018, 12(4): 487–495. doi: 10.1049/iet-spr.2017.0340
[10]	MENGÜÇ E C, ACIR N, and MANDIC D P. Widely linear quaternion-valued least-mean kurtosis algorithm[J]. IEEE Transactions on Signal Processing, 2020, 68: 5914–5922. doi: 10.1109/TSP.2020.3029959
[11]	OGUNFUNMI T and SAFARIAN C. The quaternion stochastic information gradient algorithm for nonlinear adaptive systems[J]. IEEE Transactions on Signal Processing, 2019, 67(23): 5909–5921. doi: 10.1109/TSP.2019.2944757
[12]	STERN S and FISCHER R F H. Quaternion-valued multi-user MIMO transmission via dual-polarized antennas and QLLL reduction[C]. 2018 25th International Conference on Telecommunications (ICT), Saint Malo, France, 2018: 63–69.
[13]	WYSOCKI B J, WYSOCKI T A, and SEBERRY J. Modeling dual polarization wireless fading channels using quaternions[C]. Joint IST Workshop on Mobile Future, 2006 and the Symposium on Trends in Communications. SympoTIC '06, Bratislava, Slovakia, 2006: 68–71.
[14]	LIU Wei. Antenna array signal processing for quaternion-valued wireless communication systems[C]. 2014 IEEE Benjamin Franklin Symposium on Microwave and Antenna Sub-systems for Radar, Telecommunications, and Biomedical Applications (BenMAS), Philadelphia, USA, 2014: 1–3.
[15]	LIU Wei. Channel equalization and beamforming for quaternion-valued wireless communication systems[J]. Journal of the Franklin Institute, 2017, 354(18): 8721–8733. doi: 10.1016/j.jfranklin.2016.10.043
[16]	曾乐雅, 许华, 王天睿. 自适应切换双模盲均衡算法[J]. 电子与信息学报, 2016, 38(11): 2780–2786. doi: 10.11999/JEIT160099 ZENG Leya, XU Hua, and WANG Tianrui. Dual mode blind equalization algorithm based on adaptive switching[J]. Journal of Electronics &Information Technology, 2016, 38(11): 2780–2786. doi: 10.11999/JEIT160099
[17]	XU Tongtong, XIANG Zheng, YANG Hua, et al. Concurrent modified constant modulus algorithm and decision directed scheme with Barzilai-Borwein method[J]. Frontiers in Neurorobotics, 2021, 15: 699221. doi: 10.3389/fnbot.2021.699221
[18]	LI Sen, WANG Yan, and LIN Bin. Concurrent blind channel equalization in impulsive noise environments[J]. Chinese Journal of Electronics, 2013, 22(4): 741–746.
[19]	胡婉如, 梅如如, 崔健, 等. 一种基于CMA和DDLMS算法的双模式盲均衡算法[J]. 电讯技术, 2021, 61(1): 83–88. doi: 10.3969/j.issn.1001-893x.2021.01.013 HU Wanru, MEI Ruru, CUI Jian, et al. A dual-mode blind equalization algorithm based on CMA and DDLMS algorithm[J]. Telecommunication Engineering, 2021, 61(1): 83–88. doi: 10.3969/j.issn.1001-893x.2021.01.013
[20]	陈康. 室内可见光通信自适应均衡技术的研究[D]. [硕士论文], 长春理工大学, 2019. CHEN Kang. Research on adaptive equalization technology of indoor visible light communication[D]. [Master dissertation], Changchun University of Science and Technology, 2019.
[21]	CHEN S, COOK T B, and ANDERSON L C. A comparative study of two blind FIR equalizers[J]. Digital Signal Processing, 2004, 14(1): 18–36. doi: 10.1016/j.dsp.2003.04.001

施引文献

期刊类型引用(0)
其他类型引用(1)

资源附件(0)

访问统计

图(14) / 表(1)

计量

文章访问数: 343
HTML全文浏览量: 217
PDF下载量: 41
被引次数: 1

1. 引言
2. 4元数基本理论和4维调制(QOD)技术
2.1 4元数基本理论
2.2 4维调制(QOD)技术
3. 4元数自适应均衡算法
3.1 4元数最小均方误差(QLMS)算法
3.2 4元数恒模盲均衡(QCMA)算法
3.3 4元数直接判决并行恒模均衡算法
4. 仿真实验
5. 结论

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

基于四元数的直接判决并行恒模均衡算法

doi: 10.11999/JEIT221413

通讯作者:
李森　listen@dlmu.edu.cn

计量

Concurrent Decision Directed and Constant Modulus Equalization Algorithm Based on Quaternion

1. 引言

2. 4元数基本理论和4维调制(QOD)技术

2.1 4元数基本理论

2.2 4维调制(QOD)技术

3. 4元数自适应均衡算法

3.1 4元数最小均方误差(QLMS)算法

3.2 4元数恒模盲均衡(QCMA)算法

3.3 4元数直接判决并行恒模均衡算法

4. 仿真实验

5. 结论

期刊类型引用(0)

其他类型引用(1)

计量

目录

1. 引言

2. 4元数基本理论和4维调制(QOD)技术

2.1 4元数基本理论

2.2 4维调制(QOD)技术

3. 4元数自适应均衡算法

3.1 4元数最小均方误差(QLMS)算法

3.2 4元数恒模盲均衡(QCMA)算法

3.3 4元数直接判决并行恒模均衡算法

4. 仿真实验

5. 结论

留言板

基于四元数的直接判决并行恒模均衡算法

doi: 10.11999/JEIT221413

通讯作者: 李森 listen@dlmu.edu.cn

计量

出版历程

Concurrent Decision Directed and Constant Modulus Equalization Algorithm Based on Quaternion

1. 引言

2. 4元数基本理论和4维调制(QOD)技术

2.1 4元数基本理论

2.2 4维调制(QOD)技术

3. 4元数自适应均衡算法

3.1 4元数最小均方误差(QLMS)算法

3.2 4元数恒模盲均衡(QCMA)算法

3.3 4元数直接判决并行恒模均衡算法

4. 仿真实验

5. 结论

期刊类型引用(0)

其他类型引用(1)

计量

出版历程

目录

1. 引言

2. 4元数基本理论和4维调制(QOD)技术

2.1 4元数基本理论

2.2 4维调制(QOD)技术

3. 4元数自适应均衡算法

3.1 4元数最小均方误差(QLMS)算法

3.2 4元数恒模盲均衡(QCMA)算法

3.3 4元数直接判决并行恒模均衡算法

4. 仿真实验

5. 结论

通讯作者:
李森　listen@dlmu.edu.cn