Gram 矩阵及其主要性质
文章目錄
- Gram 矩陣
- 6 大性質
Gram 矩陣
假設 AAA 是一個 m×nm\times nm×n 階矩陣,
AAA 由列向量 αi\mathbf{\alpha}_iαi? 表示, 即
A=[α1α2?αn]A=\begin{bmatrix}\mathbf{\alpha}_1 & \mathbf{\alpha}_2 &\cdots & \mathbf{\alpha}_n \end{bmatrix}A=[α1??α2????αn??]
則
G=ATA=[α1Tα2T?αnT][α1α2?αn]=[α1Tα1α1Tα2?α1Tαnα2Tα1α2Tα2?α2Tαn???αnTα1αnTα2?αnTαn]\begin{aligned} G &= \, A^{\mathsf T}A \\[3pt] &= \begin{bmatrix} \mathbf{\alpha}_1^{\mathsf T} \\ \mathbf{\alpha}_2^{\mathsf T} \\ \vdots \\ \mathbf{\alpha}_n^{\mathsf T} \end{bmatrix} \begin{bmatrix} \mathbf{\alpha}_1 & \mathbf{\alpha}_2 & \cdots & \mathbf{\alpha}_n \end{bmatrix} \\[3pt] & = \begin{bmatrix} \mathbf{\alpha}_1^{\mathsf T}\mathbf{\alpha}_1 & \mathbf{\alpha}_1^{\mathsf T}\mathbf{\alpha}_2 & \cdots & \mathbf{\alpha}_1^{\mathsf T}\mathbf{\alpha}_n \\ \mathbf{\alpha}_2^{\mathsf T}\mathbf{\alpha}_1 & \mathbf{\alpha}_2^{\mathsf T}\mathbf{\alpha}_2 & \cdots &\mathbf{\alpha}_2^{\mathsf T}\mathbf{\alpha}_n \\ \vdots & \vdots & & \vdots \\ \mathbf{\alpha}_n^{\mathsf T}\mathbf{\alpha}_1 & \mathbf{\alpha}_n^{\mathsf T}\mathbf{\alpha}_2 & \cdots & \mathbf{\alpha}_n^{\mathsf T}\mathbf{\alpha}_n \end{bmatrix} \end{aligned} G?=ATA=??????α1T?α2T??αnT????????[α1??α2????αn??]=??????α1T?α1?α2T?α1??αnT?α1??α1T?α2?α2T?α2??αnT?α2??????α1T?αn?α2T?αn??αnT?αn?????????
AAA 由行向量 βiT\mathbf{\beta}_i^{\mathsf T}βiT? 表示, 即
A=[β1Tβ2T?βmT]A=\begin{bmatrix}\mathbf{\beta}_1^{\mathsf T} \\ \mathbf{\beta}_2^{\mathsf T} \\ \vdots \\ \mathbf{\beta}_m^{\mathsf T} \end{bmatrix}A=??????β1T?β2T??βmT????????
則
G=AAT=[β1Tβ2T?βmT][β1β2?βm]=[β1Tβ1β1Tβ2?β1Tβmβ2Tβ1β2Tβ2?β2Tβm???βmTβ1βmTβ2?βmTβm]\begin{aligned} G &= \, AA^{\mathsf T} \\[3pt] &= \begin{bmatrix}\mathbf{\beta}_1^{\mathsf T} \\ \mathbf{\beta}_2^{\mathsf T} \\ \vdots \\ \mathbf{\beta}_m^{\mathsf T} \end{bmatrix} \begin{bmatrix} \mathbf{\beta}_1 & \mathbf{\beta}_2 & \cdots & \mathbf{\beta}_m \end{bmatrix} \\[3pt] & = \begin{bmatrix} \mathbf{\beta}_1^{\mathsf T}\mathbf{\beta}_1 & \mathbf{\beta}_1^{\mathsf T}\mathbf{\beta}_2 & \cdots & \mathbf{\beta}_1^{\mathsf T}\mathbf{\beta}_m \\ \mathbf{\beta}_2^{\mathsf T}\mathbf{\beta}_1 & \mathbf{\beta}_2^{\mathsf T}\mathbf{\beta}_2 & \cdots &\mathbf{\beta}_2^{\mathsf T}\mathbf{\beta}_m \\ \vdots & \vdots & & \vdots \\ \mathbf{\beta}_m^{\mathsf T}\mathbf{\beta}_1 & \mathbf{\beta}_m^{\mathsf T}\mathbf{\beta}_2 & \cdots & \mathbf{\beta}_m^{\mathsf T}\mathbf{\beta}_m \end{bmatrix} \end{aligned} G?=AAT=??????β1T?β2T??βmT????????[β1??β2????βm??]=??????β1T?β1?β2T?β1??βmT?β1??β1T?β2?β2T?β2??βmT?β2??????β1T?βm?β2T?βm??βmT?βm?????????
6 大性質
下面只考慮列向量 Gram 矩陣
(1) G=ATAG = \, A^{\mathsf T}AG=ATA 是對稱矩陣
GT=(ATA)T=ATA=GG^{\mathsf T } = \, (A^{\mathsf T}A)^{\mathsf T} = \, A^{\mathsf T}A = G GT=(ATA)T=ATA=G
(2) 對于實矩陣 AAA rank(ATA)=rank(A)\mathrm{rank} (A^{\mathsf T}A) = \mathrm{rank} (A)rank(ATA)=rank(A)
證明 {Ax=0ATAx=0\begin{cases} A\mathsf{x} = 0 \\ A^{\mathsf T}A\mathbf{x} = 0 \end{cases}{Ax=0ATAx=0? 同解即可.
證明過程詳見經典例題(第3小問)
(3) 若 ATA=0A^{\mathsf T}A=0ATA=0, 則 A=0A = 0A=0
由上面性質
rank(ATA)=rank(A)=rank(0)=0\begin{aligned} \mathrm{rank} (A^{\mathsf T}A) &= \mathrm{rank} (A) \\ &= \mathrm{rank} \ (0) = 0 \end{aligned}rank(ATA)?=rank(A)=rank?(0)=0?
(4) 對于實矩陣 AAA, 則 ATAA^{\mathsf T}AATA 是半正定矩陣
xTATAx=(Ax)TAx≥0\mathbf{x}^{\mathsf T}A^{\mathsf T}A\mathbf{x} = (A\mathbf{x})^{\mathsf T}A\mathbf{x} \geq 0 xTATAx=(Ax)TAx≥0
(5) 對于任意 nnn 階實對稱半正定矩陣 MMM, 存在矩陣 AAA 使得 M=ATAM=A^{\mathsf T}AM=ATA 成立.
因為矩陣 MMM 實對稱, 所以 MMM 可以正交對角化, 即M=QΛQTM = Q\Lambda Q^{\mathsf T}M=QΛQT 又因為矩陣 MMM 半正定, 所以其特征值 $\lambda_i \geq 0 $, 所以可記 Λ1/2=diag(λ1,…,λn)\Lambda^{1/2} = \mathrm{diag} ( \sqrt{\lambda_1}, \dots, \sqrt{\lambda_n})Λ1/2=diag(λ1??,…,λn??) 且 KaTeX parse error: Expected 'EOF', got '}' at position 29: …2}Q^\{\mathsf T}? 則可得
M=QΛQT=(Λ1/2QT)TΛ1/2QT=ATA\begin{aligned} M &= Q\Lambda Q^{\mathsf T} \\ &= (\Lambda^{1/2}Q^{\mathsf T})^{\mathsf T}\Lambda^{1/2}Q^{\mathsf T} \\ &= A^{\mathsf T}A \end{aligned}M?=QΛQT=(Λ1/2QT)TΛ1/2QT=ATA?
(6) 若 A=[α1α2?αn]A=\begin{bmatrix}\mathbf{\alpha}_1 & \mathbf{\alpha}_2 &\cdots & \mathbf{\alpha}_n \end{bmatrix}A=[α1??α2????αn??] 列滿秩, 則 ATAA^{\mathsf T}AATA 正定
- 由性質 (2), 知 rank(ATA)=rank(A)=n\mathrm{rank} (A^{\mathsf T}A) = \mathrm{rank} (A) = nrank(ATA)=rank(A)=n
- 因為 Ax=0A\mathbf{x}=0Ax=0 只有零解, 結合性質 (4), 對于非零 x∈Rn\mathbf{x}\in \mathbb{R}^nx∈Rn
xTATAx=(Ax)TAx>0\mathbf{x}^{\mathsf T}A^{\mathsf T}A\mathbf{x} = (A\mathbf{x})^{\mathsf T}A\mathbf{x} > 0 xTATAx=(Ax)TAx>0
原文鏈接
[1] matnoble.me/posts/gram
[2] 關注我吧
總結
以上是生活随笔為你收集整理的Gram 矩阵及其主要性质的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 拳王虚拟项目公社:虚拟资源自动化挣钱,空
- 下一篇: MySQL 数值拼接字符串