矩阵束分解算法及在通信信号处理中的应用

摘要

矩阵束的三角化与对角化分解，如广义Schur分解(GSD)、广义特征值分解(GEVD)与广义奇异值分解(GSVD)，是矩阵计算的重要内容，并在无线通信与信号处理等工程领域中具有广泛的应用。
　　 无线与移动通信的空时接收信号处理问题，常可归结为相关矩阵束的分解。这类矩阵束中常含有阵元移位与时滞互相关矩阵等非厄米特非正定复矩阵。强相关或相干信号常会导致相关矩阵接近奇异，即具有坏条件数。近奇异非对称正定的矩阵束分解的有效计算方法的研究具有一定难度，但对于无线空时通信信号处理问题具有重要意义。在矩阵束分解算法及应用于无线通信的研究中，论文取得了下列创新成果：
　　 ①提出非对称矩阵束(A,B)的广义Schur分解(GSD)的交叉迭代算法。将该算法应用于改进强相关阵列信号波达方向估计的ESPRIT子空间方法的计算实现，提高了方向角度的估计精度和分辨概率。
　　 提出的交叉迭代算法通过在矩阵A,B间交换各自QR和RQ分解的酉因子，并乘上原有三角因子，实现矩阵束(A,B)的循环迭代更新，并最终收敛至三角矩阵束。论文证明了该算法的收敛性同直接计算单个矩阵AB-1Schur分解的QR迭代算法收敛性等价。然而，由于包含矩阵求逆与乘积运算，若直接计算AB-1，其条件数可能很大(即使A,B都是良态的)，在有限字长条件下，QR迭代精度无法保证。交叉迭代算法绕开了矩阵直接求逆与乘积运算，得到的结果具有和单个矩阵A,B条件数大致相当的精度。该算法的速度不受B秩亏损的影响。
　　 与矩阵束GSD的QZ算法相比，交叉迭代算法中两个矩阵A,B的QR(或RQ)分解及乘积更新过程在计算上相互独立，仅需交换酉因子，具有环型并行计算结构；交叉迭代算法的一步迭代计算主要集中于QR(与RQ)分解及相反过程(酉旋转矩阵与上三角阵的乘积)，硬件和软件的实现都易于模块化与并行化。
　　 ②建立了自适应GSD和GEVD的递推交叉迭代算法，并用于解决快速时变环境中未知CDMA系统自适应盲辨识问题。
　　 在实际的自适应应用中，根据实时更新的输入数据递推生成输出。与批处理的QZ算法相比，基于QR(与RQ)分解的交叉迭代算法更适合自适应实现。自适应秩一更新矩阵束GSD的一步交叉迭代的计算代价为O(n2)。论文提出了自适应GSD和GEVD的递推交叉迭代算法，实现对广义特征值与广义特征向量的自适应更新，即论文提出的广义子空间追踪(GSST)。与子空间追踪不同，广义子空间追踪不要求相关矩阵是厄米特阵，无需对广义特征向量的正交性约束，更适用于辨识一般非列正交矩阵，如CDMA系统多用户扩频序列矩阵。
　　 借助天线阵列，可以在CDMA系统期望用户扩频波形未知(如截听、灾害等)情况下，实现对用户有效扩频波形的盲辨识和接收信号的检测。针对多径信道的快速时变性和计算的低复杂性要求，论文提出了时变环境中未知CDMA系统自适应盲辨识广义子空间追踪方法。分别基于递推交叉迭代和递推更新lanczos迭代两种算法进行相关矩阵束的自适应广义特征值分解，利用广义特征向量自适应辨识未知期望用户信号的有效扩频序列，并实现接收信号的盲多用户检测。仿真表明了GSST算法的有效性，并比较了递推交叉迭代和递推Lanczos迭代两种GSST实现形式的收敛性能。
　　 ③提出了典型相关解(CCD)精确计算的正切CCD算法。应用于相关噪声环境中的多径CDMA盲多用户检测，改善了接收信号的检测输出性能。
　　 典型相关分析(CCA)通过分析两个随机变量之间的典型相关性，寻找其中的公共模式。它综合考虑两组输入数据的互相关性及每组数据的自相关性。典型相关分解(CCD)实现两路线性变换数据之间典型相关性的最大化，从而得到该数据相关部分(即信号)的最优估计。作为三个矩阵乘积的奇异值分解，CCD的计算易受到矩阵坏条件数和元素误差的影响。基于对数据矩阵Gram-Schmidt正交化预处理的已有算法亦存在误差累积效应，在矩阵规模较大时尤为明显。
　　 在结合矩阵束广义奇异值分解与乘积奇异值分解精确计算的正切算法基础上，论文建立了精确计算典型相关分解的正切CCD算法。该算法无需正交化预处理，不存在误差累积效应，应用于改进具有多接收阵元的有色噪声环境中多径CDMA盲多用户检测的子空间方法，可准确辨识信号与噪声子空间。仿真表明多径信道与用户有效特征波形盲估计及盲多用户检测性能得到明显改善；相对于原有方法，基于正切CCD的子空间方法的信号检测结果对矩阵规模大小和条件数不敏感，在高负载系统条件下性能优势更为明显。

目录

文摘

英文文摘

声明

1绪论

1.1计算——第三种科学方法

1.2矩阵计算

1.3移动通信的发展与趋势

1.4无线信道

1.5 移动通信中的阵列信号处理

1.6CDMA与多用户检测

2矩阵的Schur分解

2.1 Schur分解

2.2 QR迭代法

2.2.1幂法

2.2.2正交迭代法

2.2.3 QR迭代法

2.3 Hessenberg阵规约

2.4实用QR算法

2.4.1降阶

2.4.2带位移QR迭代

2.4.3单步位移策略

2.4.4双位移策略

2.4.5双隐式位移策略

2.5附录

2.5.1子空间之间的距离

2.5.2定理2.2之证明

3矩阵束特征值分解的交叉迭代算法与阵列信号DOA估计

3.1矩阵束的分解与计算

3.1.1矩阵束的分解

3.1.2 Ax-λBx的QZ算法

3.1.3广义特征向量与广义不变子空间的计算

3.2交叉迭代算法

3.3 Hessenberg——三角型矩阵束交叉迭代

3.4交叉迭代的的位移加速

3.5广义特征值分解及其性质

3.6由GSD到GEVD

3.7应用举例--ESPRIT方法

4递推交叉迭代广义子空间追踪与阵列CDMA盲辨识

4.1秩一修正矩阵束的广义Schur分解

4.2秩一修正GEVD

4.3借助天线阵列的多径CDMA系统的盲辨识

4.3.1系统模型

4.3.2盲信道辨识与多用户检测的子空间方法的几个预备性引理

4.3.3不使用扩频码的盲多用户检测

4.4 CDMA盲辨识子空间算法的自适应形式

5典型相关分解的正切算法与相关噪声盲多用户检测

5.1奇异值分解的精确数值计算

5.2广义奇异值分解

5.2.1对称正定问题

5.2.2广义奇异值分解

5.2.3广义奇异值分解的实际算法

5.3乘积奇异值分解

5.3.1乘积奇异值分解问题

5.3.2乘积奇异值分解的精确计算

5.3.3三矩阵乘积奇异值分解的精确计算

5.4典型相关分解精确计算的正切算法

5.5 CDMA多用户检测

5.5.1 同步CDMA信号的线性接收机

5.5.3盲多用户检测的直接法

5.5.4盲多用户检测的子空间方法

5.5.5多径信道中的盲多用户检测

5.6相关噪声盲多用户检测的正切CCD算法与仿真

6结论与展望

致谢

参考文献

附录