![阵列信号处理及MATLAB实现(第2版)](https://wfqqreader-1252317822.image.myqcloud.com/cover/854/41202854/b_41202854.jpg)
2.1 矩阵代数的相关知识
2.1.1 特征值与特征向量
令,
,若标量λ和非零向量e满足方程
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-37-3.jpg?sign=1739592613-e613bJGU9ccpHBp3LWWrqWEJ8tSTVHNJ-0-78ab9a78cd5d0db0db8c24386ccd0c00)
则称λ是矩阵A的特征值,e是与λ对应的特征向量。特征值与特征向量总是成对出现,称(λ,e)为矩阵A的特征对,特征值可能为零,但是特征向量一定非零。
2.1.2 广义特征值与广义特征向量
令,
,若标量λ和非零向量e满足方程
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-37-6.jpg?sign=1739592613-XF588KjLT5onBik4dZPWna884GXYOwni-0-a0d89f09e90a5b2af8b85fdcf3ee6f71)
则称λ是矩阵A相对于矩阵B的广义特征值,e是与λ对应的广义特征向量。如果矩阵B非满秩,那么λ可以是任意值(包括零)。当矩阵B为单位阵时,式(2-2)就称为普通的特征值问题,因此式(2-2)可以看成对普通特征值问题的推广。
2.1.3 矩阵的奇异值分解
对于复矩阵,称
的n个特征根
的算术根
(i=1,2,…,n)为A的奇异值。若记
,其中
是A的全部非零奇异值,则称m×n矩阵为A的奇异值矩阵。
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-38-7.jpg?sign=1739592613-5BkB2hCMNFMqyLb3QSl5KtPbyW3ibI8q-0-38f6c63755beb765e6add0c98a2e740e)
奇异值分解定理:对于m×m维矩阵A,则分别存在一个m×n维酉矩阵U和一个n×n维酉矩阵V,使得
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-38-8.jpg?sign=1739592613-DhrlgqKgaIlBsWQF54U7orGcRN1XaDu9-0-12056e911a937892b02606daa3374a19)
其中,上标H表示矩阵的共轭转置。
2.1.4 Toeplitz矩阵
定义:具有2n-1个元素的n阶矩阵
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-38-9.jpg?sign=1739592613-wdUAASRKNHG48F6Bh85wY5nSGoOM9NES-0-951c884ae6f7d1bf4f93e962645c7876)
称为Toeplitz矩阵,简称T矩阵。
T矩阵也可简记为
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-38-10.jpg?sign=1739592613-G40ygpGMylkMxC8ZZ8ZafRCkkl9OmMGY-0-45356718cc25eb0cd88763f5bca796fe)
式中,记号中的“1”和“n”表示矩阵A元素的下标,i,j=1,2,…,n。T矩阵完全由第1行和第1列的2n-1个元素确定。可见,Toeplitz矩阵中位于任意一条平行于主对角线的元素全都是相等的,且关于副对角线对称。
2.1.5 Hankel矩阵
定义:具有以下形式的n+1阶矩阵
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-39-1.jpg?sign=1739592613-emkMJa0LFpVOVBKmuKbGlS83OwUTmTUH-0-97aecb9cdde335aa3d6896bec7e36583)
称为Hankel矩阵或正交对称矩阵(Orthosymmetric Matrix)。
可见,Hankel矩阵完全由其第1行和第n列的2n+1个元素确定。其中,所有垂直于主对角线的直线上有相同的元素。
2.1.6 Vandermonde矩阵
定义:具有以下形式的m×n阶矩阵
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-39-2.jpg?sign=1739592613-2rlaJX4ggsfu8LpsM49vnoKCKYAeNX4h-0-cc36d49f46aa85b2d8219f5b4db7843f)
称为Vandermonde矩阵。如果ai≠aj,那么V(a1,a2,…,an)是非奇异的。
2.1.7 Hermitian矩阵
如果矩阵An×n满足
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-39-3.jpg?sign=1739592613-E9wtcWWKchrJSv4KIz8cNLn1gUWmxGr0-0-86e782e95f67b7b41c5c0fa62e35e3c4)
则A称为Hermitian矩阵。Hermitian矩阵有以下主要性质:
(1)所有特征值都是实数。
(2)对应于不同特征值的特征向量相互正交。
(3)Hermitian矩阵可分解为的形式,这一分解称为谱定理,也就是矩阵A的特征值分解定理,其中
,
是由特征向量构成的酉矩阵[1]。
2.1.8 Kronecker乘积
定义:p×q矩阵A和m×n矩阵B的Kronecker乘积记作A⊗B,它是一个pm×qn矩阵,定义为
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-40-3.jpg?sign=1739592613-faSllTwx01Iskr8j3xndEEASl3PEWKvV-0-bee3c30f2d479f25024cd76de170cc9f)
Kronecker乘积有一个重要的性质,即:,
,
,以下等式成立:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-40-7.jpg?sign=1739592613-YMuFY0IHnz1eoqguUWn6H5SAYq8uAF4V-0-d47cd4574f90f99c701e1cdd36b291a8)
其中,vec(·)为向量化算子,,且vec(A)具有如下形式:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-40-9.jpg?sign=1739592613-DWAFrTMUNKZjsN23qBE4GZPfneZNkIWf-0-ebd82a9fee8ed6ec8a4099db9445f6f1)
Kronecker乘积具有如下性质:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-40-10.jpg?sign=1739592613-W9bQlnBC8wIZUqjjQQC9XvYHdImMerFR-0-67335202885045388a886bb28f1618ad)
2.1.9 Khatri-Rao乘积
考虑两个矩阵A(I×F)和B(J×F),它们的Khatri-Rao乘积A⊙B为一个IJ×F维矩阵,其定义为
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-41-1.jpg?sign=1739592613-QeezhDudtefcvZQNBoD5FQ895tLXHmgm-0-f6c4e56c192ea44d64eaa8abfb82b633)
其中,aF为A的第f列,bF为B的第f列,即Khatri-Rao乘积是列向量的Kronecker乘积。
Khatri-Rao乘积具有如下性质:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-41-2.jpg?sign=1739592613-8D1BvxRtU3x6GinjOAuufFYSflS4sOBq-0-2e0423f1b34a9ff9b33fb8fcb15b28c6)
令,Khatri-Rao乘积具有如下性质:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-41-4.jpg?sign=1739592613-VnXNpjLwgPRh6rafLx6Kmt3HhJ0YFuYW-0-8771fd573d48bb4cc878aecf77f00931)
其中,unvec(·)是矩阵化算子,它是vec(·)的逆运算,具有以下形式:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-41-5.jpg?sign=1739592613-eJVLDPv0sGSDjHcNnKyYpdrCPsS8IKdD-0-2384fdc60b97d9fda2c5c78e04d2dfdc)
而diag(x)表示一个对角矩阵,其元素为向量x中的元素。
2.1.10 Hadamard乘积
矩阵和
的Hadamard乘积定义为
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-42-1.jpg?sign=1739592613-0oHrZKfR7pQGDrbv2Wq4ejWSTxF0guR5-0-37bea13d29da4d4b7ef3bd0bc697ddd4)
2.1.11 向量化
通常,张量和矩阵用向量来表示比较方便,定义矩阵的向量化为[1,2]
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-42-3.jpg?sign=1739592613-U2J3GMyT6L5iyerJ3VsCgYZ8htArLA7q-0-2d0e4edf3729df24b1eee478f7ed30b1)
式中,vec算子用于将矩阵Y的所有列堆积成一个向量;重塑(reshape)是向量化的逆函数,它将一个向量转化成一个矩阵。例如,可定义为(使用MATLAB表示法并类似于MATLAB中的reshape函数):
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-42-5.jpg?sign=1739592613-wrJK2DXIRhe2iGFUlllxNxAlHChap3op-0-db5c1478e98b4031b824839c22a14e9e)
类似地,定义张量的向量化为相应的模展开矩阵
。例如,三阶张量
的向量化可写成如下形式:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-42-9.jpg?sign=1739592613-eIN97YwFtrHwMDArm1JEK6tACpwgKZHk-0-8e48bfafa7ef004fad26dcb23253904e)
vec算子的基本性质包括:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-42-10.jpg?sign=1739592613-0Miiu54aQy2CEeJ7lT7FDmg8sgcbHEU6-0-5342ae90c7c4539a52f0592198773d25)