一、命题
1.1 命题描述
A
A
A是一个实对称矩阵,则
A
A
A可以相似对角化,则
A
A
A可以正交相似对角化,即
∃
正
交
矩
阵
Q
,
对
角
矩
阵
D
,
s
.
t
.
A
=
Q
T
D
Q
\exist 正交矩阵Q,对角矩阵D ,s.t. A = Q^TDQ
∃正交矩阵Q,对角矩阵D,s.t.A=QTDQ
我们要证明的是 任给正整数
p
p
p,
∃
对
称
矩
阵
B
,
s
.
t
.
A
=
B
p
\exist 对称矩阵B,s.t. \ A = B^p
∃对称矩阵B,s.t. A=Bp,此时记作
B
=
A
1
/
p
B = A^{1/p}
B=A1/p
1.2 联系
实对称矩阵 A A A可以分解为 A = B T B A = B^TB A=BTB是易证的,一般的线性代数书里面讨论的也是这种分解
但是能否保证 B B B对称呢?答案是肯定的,后面会给出一个证明(很简单)。
笔者是在学习数理统计的过程中发现书上直接将协方差矩阵分解为两个矩阵乘积才开始纠结这个问题。
1.3 例子
在MATLAB中可以直接使用^
运算符进行求解,例如:
>> A
A =
2.0000 0.3000
0.3000 2.0000
>> A^(1/2)
ans =
1.4102 0.1064
0.1064 1.4102
>> ans^2
ans =
2.0000 0.3000
0.3000 2.0000
或者可以在Mathematica中可以用 MatrixPower[A, 1/3]
求矩阵的幂,我们会发现,求得的矩阵都是对称阵
二、证明
设矩阵A可分解为 A = Q T D Q A = Q^TDQ A=QTDQ
D D D为对角矩阵,定义 D b = d i a g ( d 11 b , . . . , d n n b ) D^b = diag(d_{11}^b,...,d_{nn}^b) Db=diag(d11b,...,dnnb),也就是给对角线上每个元素求b次幂
则有
因此 A 1 / p = Q T D 1 / p Q A^{1/p} = Q^TD^{1/p}Q A1/p=QTD1/pQ
显然 Q T D 1 / p Q Q^TD^{1/p}Q QTD1/pQ是对称矩阵,证明结束
因此在MATLAB中,如果对于一个对称矩阵求幂次,得到的一定是对称矩阵