关于复向量空间、內积、范数,零矩阵的判定条件等请参考:
线性代数基础知识(上)
线性代数基础知识(下)
广义逆矩阵
- 左逆矩阵与右逆矩阵
- {1}逆
- Moore-Penrose广义逆
- {1,4}逆与极小范数解
- {1,3}逆与最小二乘解
- PM逆与极小范数最小二乘解
左逆矩阵与右逆矩阵
我们知道,只有方阵才有逆矩阵,且可逆方阵对方阵是有限制条件的,只有行列式不为零的方阵才可逆。可逆方阵给我们解线性方程组带来了很大的方便:设
Ax=y(A称为
Ax=y的系数矩阵)是关于x的方程,若A是方阵且A可逆,则
x=A−1y。然而对于一般的
m×n系数矩阵A,有没有简洁的办法来求解该线性方程组呢?我们注意到,实际上只需有一个矩阵L,如果我们用L左乘
Ax=y的两端,能够“抵消”掉A,得到
x=Ly,那么就得到了解。显然这样的L应满足
LAx=x,考虑其中一种最简单的情形:
LA=I,这就引出了左逆矩阵的概念:
- 定义:设A是
m×n矩阵,若存在
n×m矩阵L,使得
LA=I,则称L是A的一个左逆矩阵
先来看一下左逆矩阵存在的条件:
- 定理:设A是m×n矩阵,若A的左逆矩阵L存在,则
m⩾n
证明:(反证法)
假设
m<n,则可设
L=[Um×mV(n−m)×m],
A=[Rm×mSm×(n−m)],由
LA=[URVRUSVS]=[IOOI],得
UR=I,US=O,故U可逆,故
S=O,这与
VS=I是矛盾的。故假设不成立,得证。
- 定理:设A是m×n矩阵,则A的左逆矩阵L存在的充要条件为A是列满秩矩阵
证明:
必要性:若
LA=I,则
n=r(I)=r(LA)⩽r(A),故
r(A)=n,即A是列满秩的。
充分性:若A是列满秩的,则由
r(AHA)=r(A)=n知,
AHA是满秩方阵,即
AHA可逆。设
L=(AHA)−1AH,则
LA=(AHA)−1AHA=I,故A的左逆矩阵存在。
这说明并非所有
m×n矩阵A都有左逆矩阵。既然有左逆矩阵,自然也有右逆矩阵的概念。
- 定义:设A是
m×n矩阵,若存在
n×m矩阵R,使得
AR=I,则称R是A的一个右逆矩阵
- 定理:设A是m×n矩阵,若A的右逆矩阵R存在,则
m⩽n
- 定理:设A是m×n矩阵,则A的右逆矩阵R存在的充要条件为A是行满秩矩阵
{1}逆
我们讨论了左逆、右逆的概念及其存在的条件,显然它们不是求解线性方程组的最有力工具,左逆矩阵只能用来求解系数矩阵列满秩的线性方程组。现在我们从一般的线性方程组出发,探究什么样的矩阵可以用来求解任意的线性方程组。
对于一般的线性方程组
Am×nx=y,如果它有解,我们不妨作最简单的假定,即假定它的解具有
x=By的形式。现在我们来探究一下这样的矩阵B究竟是什么:
- 定理:命题(1)“
∀y∈Cm,若关于x的线性方程组
Am×nx=y有解,则
x=By是它的一个解,其中B只与A有关”成立的充要条件为命题(2)“矩阵B满足
ABA=A”成立
证明:
必要性:
∀z∈Cm,当
y=Az时,关于x的线性方程组
Ax=y=Az必有解(显然z就是它的一个解)。若命题(1)成立,则命题“
∀z∈Cm,关于x的线性方程组
Ax=Az的一个解是
x=BAz”也成立。
(注意,这里只是对命题1中的条件“若线性方程组
Am×nx=y有解”作了一个转化,当y任取时,
Am×nx=y是不一定有解的,而当y具有形式
y=Az时,
Am×nx=y就有解了,此时命题(1)后半部分的断言依然成立)
故“
∀z∈Cm,ABAz=Az”成立。则
∀z∈Cm,(ABA−A)z=0,由零矩阵的判定条件知
ABA−A=O,故
ABA=A。
充分性:若ABA=A,则
∀y∈Cm,若关于x的线性方程组
Am×nx=y有解,那么可设
x=x0是它的一个解,则有
y=Ax0,因为
ABy=ABAx0=Ax0=y,所以
x=By也是它的一个解。故命题(1)成立。
上述定理说明,我们期望找到的矩阵B其实就是满足
ABA=A的矩阵B,我们把满足该条件的矩阵B称为A的一个广义逆矩阵,更确切地,B称为A的一个{1}逆(因为还有其他类型的广义逆矩阵):
- 定义:设A是一个
m×n矩阵,若存在一个
n×m矩阵B使得
ABA=A,则称B是A的一个{1}逆,记作
B=A(1)。通常,将A的全体{1}逆的集合写作
A{1}。
任意一个
m×n矩阵A的{1}逆一定存在吗?下面的定理给出了肯定的回答:
- 定理:设
A∈Crm×n(即A是一个秩为r的
m×n复矩阵)。若
r=0,则
A{1}=Cn×m;若
r>0,存在可逆矩阵P、Q使得
PAQ=[IrOOO],则
A{1}={Q[IrL21L12L22]P∣∣∣∣L12∈Cr×(m−r),L21∈C(n−r)×r,L22∈C(n−r)×(m−r)}
证明:只证
r>0的情况。由于
A=P−1[IrOOO]Q−1,故任意矩阵
X=Q[IrL21L12L22]P,易验证
AXA=A,故
X∈A{1}。任意矩阵
A(1)∈A{1},设
A(1)=Q[L11L21L12L22]P,由
AA(1)A=A得
[IrOOO][L11L21L12L22][IrOOO]=[IrOOO],故
L11=Ir,即
A(1)=Q[IrL21L12L22]P。得证。
上述定理不仅证明了任意矩阵A都有{1}逆,还给出了{1}逆的求法,即对A进行初等变换,求出变换所用的可逆矩阵P、Q,就可以求出A的任意一个{1}逆。
- 推论:设
A∈Crm×n,则A的{1}逆唯一的充要条件为
r=m=n>0
【注】此时A的{1}逆就是A的逆矩阵。
回到求解线性方程组的问题上来。我们已经知道任取A的一个{1}逆
A(1), 若
Ax=y有解,则
x=A(1)y一定是它的一个解。那什么条件下
Ax=y才有解?它的通解又是什么?如果我们将
x=A(1)y代入原方程,就得到
AA(1)y=y,这是在原方程有解的条件下得到的结论。然而,如果
AA(1)y=y,这不就说明
x=A(1)y是原方程的一个解吗?这就得到了线性方程组有解的充要条件。
- 定理:线性方程组
Am×nx=y有解的充要条件是存在A的一个{1}逆
A(1)使得
AA(1)y=y
证明:
必要性:若
Am×nx=y有解,则
y=Ax=AA(1)Ax=AA(1)y。
充分性:若
AA(1)y=y,则
x=A(1)y是原方程的一个解,故原方程有解。
- 定理:线性方程组
Am×nx=y有解的充要条件是任意A的一个{1}逆
A(1)都有
AA(1)y=y
证明:同上。
关于
Ax=y的通解(这里是指
Ax=y的所有解都要能写成通解的形式,不仅仅是找到
Ax=y的很多个解,还要找到
Ax=y的全部解),有以下定理:
- 定理:若
Am×nx=y有解,则任取A的一个{1}逆
A(1),
Ax=y的通解是
x=A(1)y+(I−A(1)A)z,其中
z∈Cn任取
证明:
将
x=A(1)y+(I−A(1)A)z,z∈Cn代入原方程,可得
Ax=AA(1)y+A(I−A(1)A)z=y+(A−AA(1)A)z=y,可见
x=A(1)y+(I−A(1)A)z,z∈Cn都是原方程的解。
任取原方程的一个解
x0,则
Ax0=y,令
z=x0,则
x=A(1)y+(I−A(1)A)z=A(1)y+(I−A(1)A)x0=x0+A(1)y−A(1)y=x0可见
x=A(1)y+(I−A(1)A)z,z∈Cn还包括了原方程的所有解。得证。
上述定理说明,A的任意一个{1}逆都能给出
Ax=y的解的完整表述。这说明{1}逆就是解线性方程组的一个完备的工具。细心的朋友可能已经发现,既然A的任意一个{1}逆
A(1)都能写出原方程的全部解,而我们又知道
x=A(1)y必然是原方程的解,那么我们可以推断:若给定A的一个{1}逆G,那么对于A的任意一个{1}逆M,必存在
z∈Cm,使得
My=Gy+(I−GA)z。现在的问题是,对于任意的
z∈Cm,是否存在A的一个{1}逆M,使得
My=Gy+(I−GA)z成立?就是说,集合
{My∣M∈A{1}}是否包含了原方程的所有解?当y=0时,只有A列满秩时,
{My∣M∈A{1}}才包含了原方程的所有解(因为此时原方程只有零解)。那么当
y̸=0时呢?在解决这个问题之前,我们先将{1}逆这个工具运用到更一般的矩阵方程上。
- 定理:关于
Xm×n的矩阵方程
AXB=D有解的充要条件为存在(或任意)A的一个{1}逆
A(1)和B的一个{1}逆
B(1)使得
D=AA(1)DB(1)B成立;任取A的一个{1}逆
A(1)和B的一个{1}逆
B(1),若
AXB=D有解,则其通解为
X=A(1)DB(1)+Y−A(1)AYBB(1),其中
Y∈Cm×n任取
证明:
若
AXB=D有解,则
D=AXB=AA(1)AXBB(1)B=AA(1)DB(1)B若
D=AA(1)DB(1)B,则
X=A(1)DB(1)是原方程的一个解。
将
X=A(1)DB(1)+Y−A(1)AYBB(1)代入原方程得
AXB=AA(1)DB(1)B+AYB−AA(1)AYBB(1)B=D+AYB−AYB=D故
X=A(1)DB(1)+Y−A(1)AYBB(1)都是原方程的解。
任取原方程的一个解
X0,则
AX0B=D,令
Y=X0,则
X=A(1)DB(1)+Y−A(1)AYBB(1)=A(1)DB(1)+X0−A(1)AX0BB(1)=X0+A(1)DB{1}−A(1)DB(1)=X0故
X=A(1)DB(1)+Y−A(1)AYBB(1)还包含了原方程的所有解。得证。
现在,我们考虑关于M的矩阵方程
AMA=A,解该方程,得到如下结论:
- 定理:给定
Am×n的一个{1}逆G,则
M=G+Y−GAYAG,Y∈Cm×n给出了A的全部{1}逆
证明:解方程
AMA=A,可得方程的通解为
M=GAG+Z−GAZAG,Z∈Cm×n令
Y=Z−G,得
M=GAG+Y+G−GAYAG−GAGAG=G+Y+GAG−GAG−GAYAG=G+Y−GAYAG得证。
这说明A的所有{1}逆都可以用A的某个给定的{1}逆表达出来。现在回到问题:集合
{My∣M∈A{1}}是否包含了
Ax=y的所有解?有了上面的结论的铺垫,我们现在可以解决这个问题:
- 定理:若
Am×nx=y,y̸=0有解,则其通解是
x=My,M∈A{1}
证明:
设G是A的一个{1}逆,则由
Ax=y有解知
AGy=y,且原方程的通解是
x=Gy+(I−GA)z。设M是A的任意一个{1}逆,则存在矩阵Y,使得
M=G+Y−GAYAG。问题转化为,对任意
z∈Cn,能否找到Y,使得
Gy+(I−GA)z=My=Gy+Yy−GAYAGy=Gy+(I−GA)Yy成立。即是否存在Y使得
(I−GA)(Yy−z)=0。显然,只要找到Y满足
Yy=z即可。因为
y̸=0,故容易验证
y(1)=(yHy)−1yH是y的一个{1}逆,且满足
y(1)y=1。由于
z=zy(1)y,故关于Y的矩阵方程
Yy=z有解,且
Y=zy(1)就是它的一个解。综上,找到了Y,即找到了M使得
My=Gy+(I−GA)z,故得证。
实际上,设x是
Ax=y的一个解,则有
x=Gy+(I−GA)x=Gy+(I−GA)xy(1)y=(G+(I−GA)xy(1))y容易验证
G+(I−GA)xy(1)是A的一个{1}逆,故上述定理成立。
在机器学习中,线性回归模型是最基础也最简单的模型之一,在对实际数据进行拟合时,往往是不可能做到完全拟合的。故对于一般的线性回归问题,我们往往考虑其最小二乘解(或者完全等价地,最小化线性回归的代价函数,即均方误差函数)。而且我们通常不希望解的范数太大,故还需考虑其极小范数解。有了上面的结论,我们在寻找这些特殊解时,就可以把目标定在寻找特殊的{1}逆上。那么都有哪些特殊的{1}逆呢?我们先来看一下其中性质优良的Penrose-Moore广义逆。
Penrose-Moore广义逆
Penrose于1955年提出了Penrose-Moore条件,满足这些条件中的任何一个的矩阵G都可以称为A的一个广义逆矩阵,它们分别是:
- AGA=A
- GAG=G
- AG是共轭对称的
- GA是共轭对称的
{1}逆是满足第一个条件的广义逆矩阵,这也是{1}逆的记法的来源。类似的还有{1,2}逆、{1,3,4}逆等等。共有
24−1=15类广义逆矩阵,其中得到重要应用的有{1,2}逆(自反广义逆矩阵)、{1,2,3}逆(正规化广义逆矩阵)、{1,2,4}逆(弱广义逆矩阵)、{1,2,3,4}逆(Penrose-Moore广义逆)等。
定义:设
A∈Cm×n,
G∈Cn×m,若G满足如下四个条件,则称G是A的Penrose-Moore广义逆矩阵,简称PM逆,记为
G=A+:
-
AGA=A
-
GAG=G
-
(AG)H=AG
-
(GA)H=GA
PM逆不仅在数学规划中有着重要的应用,还在概率统计、数值分析、系统控制、博弈论、信号处理和网络理论等领域有着广泛的应用。这是因为PM逆具有着非常优良的数学性质,使得其在各个领域的理论分析中占有着重要的地位。现在,我们就来看看PM逆都有哪些优良的性质。
- 定理:任意矩阵
A∈Crm×n,A的PM逆是存在且唯一的
证明:
存在性:当
r=0时,易验证
On×m是A的一个PM逆。
当
r>0时,存在A的满秩分解
A=KL,其中
K∈Cm×r是列满秩矩阵,
L∈Cr×n是行满秩矩阵。由于
r(KHK)=r(K)=r以及
r(LLH)=r(L)=r,故
KHK和
LLH是满秩方阵。故
KHKLLH是可逆方阵。设
G=LH(KHKLLH)−1KH,现在证明G是A的一个PM逆:
AGA=KLLH(KHKLLH)−1KHKL=K(LLH)(LLH)−1(KHK)−1(KHK)L=KL=A
GAG=LH(KHKLLH)−1KHKLLH(KHKLLH)−1KH=LH(LLH)−1(KHK)−1(KHK)(LLH)(LLH)−1(KHK)−1KH=LH(KHKLLH)−1KH=G
GA=LH(KHKLLH)−1KHKL=LH(LLH)−1(KHK)−1(KHK)L=LH(LLH)−1L
(GA)H=LH((LLH)−1)HL=LH(LLH)−1L=GA
AG=KLLH(KHKLLH)−1KH=K(LLH)(LLH)−1(KHK)−1KH=K(KHK)−1KH
(AG)H=K((KHK)−1)HKH=K(KHK)−1KH=AG这就证明了A的PM逆的存在性。
唯一性:设
X,Y分别是A的一个PM逆,则
X=XAX=(XA)HX=AHXHX=(AYA)HXHX=AHYHAHXHX=(YA)H(XA)HX=YAXAX=YAX=Y(AX)H=YXHAH=YXH(AYA)H=YXHAHYHAH=Y(AX)H(AY)H=YAXAY=YAY=Y这就证明了PM逆的唯一性。得证。
该定理不仅证明了PM逆的存在性和唯一性,还给出了求PM逆的一种算法:满秩分解法。
- 定理:
r(A)=r(A+)=r(AA+)=r(A+A)=r(AA+A)=r(A+AA+)
证明:
因为
r(A)=r(AA+A)⩽r(AA+)⩽r(A+)
r(A+)=r(A+AA+)⩽r(A+A)⩽r(A)
r(AA+)⩽r(A)
r(A+A)⩽r(A+)所以
r(A)=r(A+)=r(AA+)=r(A+A)=r(AA+A)=r(A+AA+)得证。
PM逆有以下列出的一些计算性质(通过定义容易验证):
-
(A+)+=A
-
(AT)+=(A+)T
-
(AH)+=(A+)H
-
(kA)+=k1A+,k∈C,k̸=0
- 若A是n阶
(n⩾2)方阵,则
(A∗)+=(A+)∗
- 一般
(AB)+̸=B+A+,但
(AHA)+=A+(AH)+且
(AAH)+=(AH)+A+
-
(A+A)+=A+A且
(AA+)+=AA+(通过验证
A+A和
AA+是幂等矩阵易得)
-
(I−A+A)+=I−A+A且
(I−AA+)+=I−AA+(通过验证
I−A+A和
I−AA+是幂等矩阵易得)
- 若U、V为酋矩阵,则
(UAV)+=VHA+UH
从PM逆的特性上看,PM逆可能是最接近逆矩阵的广义逆了(唯一性、秩的关系、计算性质等)。当方阵A可逆时,
A+就是A的逆矩阵。此外,如果对矩阵
Am×n作一些限定,则能得到类似于逆矩阵的如下性质:
- 定理:
A+A=In的充要条件为
Am×n是列满秩矩阵
证明:
必要性:若
A+A=I,则由
r(A)=r(A+A)=n知A是列满秩矩阵
充分性:若A是列满秩矩阵,则由
r(A+A)=r(A)=n知
A+A是满秩方阵,用
(A+A)−1左乘
A+AA+A=A+A,即得
A+A=I
【注】实际上,当A列满秩时,A的左逆矩阵存在,
L=(AHA)−1AH就是A的一个左逆矩阵。容易验证该左逆矩阵就是A的PM逆。
- 定理:
AA+=Im的充要条件为
Am×n是行满秩矩阵
证明:与上同理。
【注】实际上,当A行满秩时,A的右逆矩阵存在,
R=AH(AAH)−1就是A的一个右逆矩阵。容易验证该右逆矩阵就是A的PM逆。
PM逆作为一种特殊的{1}逆,当然可以像{1}逆那样表达线性方程组的解的结构:
- 定理:线性方程
Am×nx=y有解的充要条件为
y=AA+y,若它有解,则通解为
x=A+y+(I−A+A)z,z∈Cn
- 定理:关于
Xm×n的矩阵方程
AXB=D有解的充要条件为
D=AA+DB+B,若它有解,则通解为
X=A+DB++Y−A+AYBB+,Y∈Cm×n
借助PM逆,我们还能解决方程解的唯一性问题:
- 定理:若线性方程
Am×nx=y有解,则解唯一的充要条件为A是列满秩矩阵
证明:
必要性:考虑方程的通解
x=A+y+(I−A+A)z,z∈Cn,显然若方程的解唯一,则必有
∀z∈Cn,(I−A+A)z=0。故由零矩阵的判定条件知
A+A=I,故A列满秩。
充分性:若A列满秩,则
A+A=I,故方程的通解
x=A+y+(I−A+A)z=A+y,可见方程的解是唯一的。
如果A列满秩,则A的左逆矩阵存在,
L=(AHA)−1AH就是A的一个左逆矩阵,且有
LAx=x=Ly。所以上面的定理的结论可以再强一些:
-
定理:线性方程
Am×nx=y的解存在且唯一的充要条件为A是列满秩矩阵
-
定理:若关于
Xm×n的矩阵方程
AXB=D有解,则解唯一的充要条件为A列满秩且B行满秩
证明:
必要性:考虑通解
X=A+DB++Y−A+AYBB+,Y∈Cm×n,若解唯一,则
∀Y∈Cm×n,Y=A+AYBB+。下面分情况讨论:
若
m⩾n,则可取Y列满秩,由
r(Y)=r(A+AYBB+)⩽r(B)得,B行满秩。故
BB+=I,故
∀Y∈Cm×n,Y=A+AY即
(I−A+A)Y=O。则
∀z∈Cm,(I−A+A)z=0。由零矩阵的判定条件得,
A+A=I,故A列满秩。
若
m⩽n,则可取Y行满秩,由
r(Y)=r(A+AYBB+)⩽r(A)得,A列满秩。故
A+A=I,故
∀Y∈Cm×n,Y=YBB+即
Y(I−BB+)=O。则
∀zT∈Cn,z(I−BB+)=0。由零矩阵的判定条件得,
BB+=I,故B行满秩。
综上,无论何种情况,若方程的解唯一,则A列满秩且B行满秩。
充分性:若A列满秩且B行满秩,则
A+A=I且
BB+=I。故方程的通解
X=A+DB++Y−A+AYBB+=A+DB++Y−Y=A+DB+,可见方程的解是唯一的。得证。
当A列满秩、B行满秩时,
L=(AHA)−1AH和
R=BH(BBH)−1分别是A和B的左逆矩阵和右逆矩阵,则
LAXBR=X=LDR。故上述定理的结论可以再强一些:
- 定理:矩阵方程
AXB=D的解存在且唯一的充要条件为A列满秩且B行满秩
至此,关于PM逆如何计算的问题,我们只提到了一种方法:满秩分解法。实际上,求解PM逆还有很多行之有效的方法。这里介绍两例:
-
奇异值分解法:设矩阵
A∈Cm×n,A的奇异值分解为
A=UΣVH,其中U、V均为酋矩阵,
Σ是对角矩阵。则
A+=VΣ+UH
-
Greville递推法
【注】若初始列向量
a1=0,图中公式不适用,此时这样计算:
a1+=a1T
接下来,我们探讨极小范数解和最小二乘解,在探讨这个问题之前,先看下相容方程组的定义:
- 定义:关于x的线性方程组
Ax=y称为相容方程组或一致方程,当且仅当该方程组有解
说白了就是给有解方程组起个名称而已。显然齐次线性方程组都是相容的,即都是一致方程。
{1,4}逆与极小范数解
什么是极小范数解?顾名思义,就是范数最小的解。既然极小范数解是一个解,我们当然要在一致方程的范畴中讨论这个问题(后面会讨论非一致方程的极小范数最小二乘解,此时极小范数是另一个含义,不过与此大体相同)。
- 定义:一致方程
Am×nx=y的极小范数解定义为
x^满足
∣∣x^∣∣=Ax=ymin∣∣x∣∣
我们已经知道,一致方程
Am×nx=y当
y̸=0时的通解为
x=A(1)y,A(1)∈A{1}。显然,当
y=0时,
Ax=y唯一的极小范数解是
x=0;当
y̸=0时,我们可以寻求一类特殊的{1}逆M使得
x=My是极小范数解。现在,我们来看A的{1}逆M应该满足什么样的条件:
- 定理:若
M∈A{1},则命题(1)“
∀y∈Cm,x=My是一致方程
Am×nx=y的极小范数解”成立的充要条件为命题(2)“
MA=A+A”成立
证明:
命题(1)等价于“
∀b∈Cm,x=MAb是方程
Am×nx=Ab的极小范数解”。
∀b∈Cm,方程
Am×nx=Ab的通解为
x=A+Ab+(I−A+A)z,z∈Cn。因为
M∈A{1},故必存在
z^∈Cn,使得
MAb=A+Ab+(I−A+A)z^。设有实值函数
f(z)=∣∣A+Ab+(I−A+A)z∣∣2,z∈Cn。
必要性:由命题(1)知,f在
z=z^处取得最小值。则对
∀v∈Cn,实值函数
gv(t)=f(z^+tv),t∈R和
hv(t)=f(z^+itv),t∈R都在
t=0处取得最小值(i是虚数单位)。由极值的必要条件,
dtdgv∣t=0=0且
dtdhv∣t=0=0。以
dtdgv为例:
dtdgv=dtd∣∣A+Ab+(I−A+A)(z^+tv)∣∣2=dtd∣∣MAb+t(I−A+A)v∣∣2=dtd(∣∣MAb∣∣2+t2∣∣(I−A+A)v∣∣2+2tRe{<MAb,(I−A+A)v>})=2t∣∣(I−A+A)v∣∣2+2Re{<MAb,(I−A+A)v>}【注】
Re{}是指实部,
Im{}是指虚部。
当
t=0时,
dtdgv=2Re{<MAb,(I−A+A)v>}=0。同理,当
t=0时,
dtdhv=2iIm{<MAb,(I−A+A)v>}=0。故
<MAb,(I−A+A)v>=0,即
bH(MA)H(I−A+A)v=0。由零矩阵的判定条件得
bH(MA)H(I−A+A)=O,再由零矩阵的判定条件得
(MA)H(I−A+A)=O。两端取共轭转置得
(I−A+A)MA=MA−A+A=O,即
MA=A+A。
充分性:若
MA=A+A,则
(MA)H(I−A+A)=O。因为
f(z)=∣∣A+Ab+(I−A+A)z∣∣2=∣∣MAb+(I−A+A)(z−z^)∣∣2=∣∣MAb∣∣2+∣∣(I−A+A)(z−z^)∣∣2+2Re{<MAb,(I−A+A)(z−z^)>}=∣∣MAb∣∣2+∣∣(I−A+A)(z−z^)∣∣2+2Re{bH(MA)H(I−A+A)(z−z^)}=∣∣MAb∣∣2+∣∣(I−A+A)(z−z^)∣∣2⩾∣∣MAb∣∣2对任意
z∈Cn成立,故
x=MAb是原方程的极小范数解。
我们找到了M应满足的条件:
MA=A+A,这个条件看上去并不那么“显然”,接下来对它作个等价转化:
- 定理:若
M∈A{1},则命题(1)“
∀y̸=0,y∈Cm,x=My是一致方程
Am×nx=y的极小范数解”成立的充要条件为命题(2)“
(MA)H=MA”成立
证明:只需证明若
M∈A{1},则
MA=A+A⟺(MA)H=MA。
⇒:
由于
A+A是共轭对称的,故
(MA)H=MA。
⇐:
先用PM逆的定义证明
MAA+=A+:
A(MAA+)A=AA+A=A
(MAA+)A(MAA+)=M(AMA)A+=MAA+
A(MAA+)=AA+,因为
AA+是共轭对称的,故
A(MAA+)是共轭对称的
(MAA+)A=MA,因为
MA是共轭对称的,故
(MAA+)A是共轭对称的
这就证明了
MAA+=A+。用A右乘该式,就得到
MA=A+A。
上面的定理告诉我们M就是A的{1,4}逆,而且在定理的证明过程中容易发现{1,4}逆的一个充要条件:
- 定理:矩阵M是A的一个{1,4}逆的充要条件为
MA=A+A
证明:根据前面的定理已有的结论,只需证明若
MA=A+A,则
AMA=A即可。用A左乘
MA=A+A,就得到
AMA=A。得证。
至此,我们基本上已经解决了极小范数解的问题。现在还剩下一个问题:极小范数解是唯一的还是有多个?实际上,有了前面的结论的支撑,我们很容易证明极小范数解是唯一的:
- 定理:一致方程的极小范数解是唯一的
证明:我们已经知道一致方程
Ax=y的全部极小范数解是
x=My,M∈A{1,4}。现在任取
M1,M2∈A{1,4},由于
Ax=y是一致方程,故
AM1y=AM2y=y。
M1y=(M1A)M2y=(A+A)M2y=(M2A)M2y=M2y这就证明了极小范数解是唯一的。
{1,3}逆与最小二乘解
在线性回归问题中,在给定的数据集下一般是没有精确解的,那么这时我们就希望找到的解“尽量精确”。如何做到尽量精确呢?注意到线性方程组
Ax=y实际上是说,
Ax这个向量与
y这个向量之间无偏差,故对于非一致方程(不相容方程组、无解线性方程组)
Ax=y,我们可以用
Ax和
y之间的偏差有多大来反映“解”
x有多精确。常采用的“偏差”是
Ax和
y之间的欧式距离,即
∣∣Ax−y∣∣。这就引出了最小二乘解的定义:
- 定义:非一致方程
Am×nx=y的最小二乘解
x^定义为
x^满足
∣∣Ax^−y∣∣=x∈Cnmin∣∣Ax−y∣∣
在求解极小范数解的问题时,我们一开始就已经知道极小范数解具有
x=A(1)y的形式。在最小二乘问题中,我们并不好知道最小二乘解会具有怎样的形式。不过,注意到实际上最小二乘解的定义对任何线性方程组都是适用的(对于一致方程,容易看出它的精确解与最小二乘解是等价的),我们可以将一致方程的情形下的最小二乘解看作是一个特例(这里只是猜测),这样我们可以推测非一致方程的最小二乘解会具有类似于
x=A(1)y的形式(至少某一些是)。不过,在搞清一致和非一致情形下的最小二乘解究竟是否有这样的关系前,我们不能妄加限制非一致情形下的最小二乘解的形式。所以就像最初研究{1}逆一样,我们假定
x=My是
Ax=y的最小二乘解,来探究一下这样的矩阵M具有怎样的性质。
- 定理:命题(1)“
∀y∈Cm,
x=My是方程组
Am×nx=y的最小二乘解”成立的充要条件为命题(2)“
AHAM=AH”成立
证明:
对
∀y∈Cm,设有实值函数
f(x)=∣∣Ax−y∣∣2,x∈Cn。
必要性:记
x^=My,若命题(1)成立,则f在
x=x^处取得最小值。故对
∀v∈Cn,实值函数
gv(t)=f(x^+tv),t∈R和
hv(t)=f(x^+itv),t∈R(i是虚数单位)都在
t=0处取得最小值。由极值的必要条件,
dtdgv∣t=0=0且
dtdhv∣t=0=0。以
dtdgv为例:
dtdgv=dtd∣∣A(x^+tv)−y∣∣2=dtd∣∣(AM−I)y+tAv∣∣2=dtd(∣∣(AM−I)y∣∣2+t2∣∣Av∣∣2+2tRe{<(AM−I)y,Av>})=2t∣∣Av∣∣2+2Re{<(AM−I)y,Av>}当
t=0时,
dtdgv=2Re{<(AM−I)y,Av>}=0。同理,当
t=0时,
dtdhv=2iIm{<(AM−I)y,Av>}=0。故
<(AM−I)y,Av>=0,即
yH(AM−I)HAv=0。由零矩阵的判定条件得,
yH(AM−I)HA=O,再由零矩阵的判定条件得
(AM−I)HA=O。两端取共轭转置得
AH(AM−I)=O,即
AHAM=AH。
充分性:若
AHAM=AH,则
(AM−I)HA=O。因为
f(x)=∣∣Ax−y∣∣2=∣∣(AM−I)y+A(x−My)∣∣=∣∣(AM−I)y∣∣2+∣∣A(x−My)∣∣2+2Re{<(AM−I)y,A(x−My)>}=∣∣(AM−I)y∣∣2+∣∣A(x−My)∣∣2+2Re{yH(AM−I)HAv}=∣∣(AM−I)y∣∣2+∣∣A(x−My)∣∣2⩾∣∣AMy−y∣∣2对
∀x∈Cn成立。故
x=My是原方程的最小二乘解。
同样地,条件
AHAM=AH看起来不是那么“显然”,所以我们对它作个等价转化:
- 定理:命题(1)“
∀y∈Cm,
x=My是方程组
Am×nx=y的最小二乘解”成立的充要条件为命题(2)“
AMA=A且
(AM)H=AM”成立
证明:
只需证明
AHAM=AH⟺AMA=A∧(AM)H=AM。
⇒:
用A右乘
AHAM=AH得
AHAMA=AHA,即
AHA(MA−I)=O。因为
r(AHA)=r(A),故由零矩阵的判定条件得
A(MA−I)=O,即
AMA=A。用
MH左乘
AHAM=AH得
(AM)HAM=(AM)H,两端取共轭转置得
(AM)HAM=AM,故
(AM)H=AM。
⇐:
AHAM=AH(AM)H=(AMA)H=AH。
上面的定理告诉我们,M就是A的{1,3}逆。现在,像研究一致方程的通解一样,我们还能找到最小二乘解的“通解”:
- 定理:线性方程组
Ax=y的全部最小二乘解为
x=A(1,3)y+(I−A(1,3)A)z,z∈Cn
证明:
因为
∣∣A(A(1,3)y+(I−A(1,3)A)z)−y∣∣=∣∣AA(1,3)y+(A−AA(1,3)A)z−y∣∣=∣∣AA(1,3)y−y∣∣,故
x=A(1,3)y+(I−A(1,3)A)z,z∈Cn都是原方程的最小二乘解。
任取原方程的任一最小二乘解
x,由于
∣∣Ax−y∣∣2=∣∣AA(1,3)y−y+A(x−A(1,3)y)∣∣2=∣∣AA(1,3)y−y∣∣2+∣∣A(x−A(1,3)y∣∣2+2Re{<AA(1,3)y−y,A(x−A(1,3)y)>}=∣∣AA(1,3)y−y∣∣2+∣∣A(x−A(1,3)y∣∣2=∣∣AA(1,3)y−y∣∣2故
∣∣A(x−A(1,3)y∣∣2=0即
Ax=AA(1,3)y。该方程的通解为
x=A(1,3)y+(I−A(1,3)A)z,z∈Cn,故
x=A(1,3)y+(I−A(1,3)A)z,z∈Cn包含了原方程的所有最小二乘解。综上,得证。
在上述定理的证明过程中,容易看出
Ax=y的全部最小二乘解就是
Ax=AA(1,3)y的通解。此外,最小二乘解与正规方程有着紧密的联系:
- 定理:
Ax=y的全部最小二乘解是正规方程
AHAx=AHy的通解
证明:
只需证明方程
Ax=AA(1,3)y与
AHAx=AHy等价即可。
⇒:
用
AH左乘
Ax=AA(1,3)y得
AHAx=AHAA(1,3)y=((AA(1,3))HA)Hy=(AA(1,3)A)Hy=AHy
⇐:
与上同理可得
AHAx=AHAA(1,3)y,即
AHA(x−A(1,3)y)=O。因为
r(AHA)=r(A),故由零矩阵的判定条件得
A(x−A(1,3)y)=O,即
Ax=AA(1,3)y。
最小二乘解是唯一的吗?由于
Ax=y的全部最小二乘解就是一致方程
Ax=AA(1,3)y的通解,故我们只需看
Ax=AA(1,3)y的解是否唯一即可。由PM逆一节的相关结论知道,
Ax=AA(1,3)y的解唯一的充要条件是A是列满秩矩阵。所以,当A不是列满秩矩阵时,
Ax=y的最小二乘解就不唯一。
至此,我们对最小二乘解有了大致的了解。现在回到线性回归问题上,我们为了使非一致方程
Ax=y的解“尽量精确”,可以取它的最小二乘解
x=A+y+(I−A+A)z,z∈Cn(注意,PM逆是{1,3}逆的特例)。显然,如果要用这个通式计算方程的最小二乘解,取
z=0时的解
x=A+y是最好计算的。那么问题来了,这个解有没有什么特殊之处?注意到它的形式和一致方程的极小范数解的形式一样,我们可以揣度它应该是非一致方程
Ax=y的所有最小二乘解中范数最小的。下面我们就证明这个结论。
PM逆与极小范数最小二乘解
-
定义:线性方程组
Am×nx=y的极小范数最小二乘解
x^定义为
x^首先是原方程的一个最小二乘解,
x^还满足
∣∣x^∣∣=z∈Cnmin{A+y+(I−A+A)z}
-
定理:线性方程组
Am×nx=y唯一的极小范数最小二乘解是
x=A+y
证明:
∀z∈Cn,==⩾∣∣A+y+(I−A+A)z∣∣2∣∣A+y∣∣2+∣∣(I−A+A)z∣∣2+2Re{<A+y,(I−A+A)z>}∣∣A+y∣∣2+∣∣(I−A+A)z∣∣2∣∣A+y∣∣2当且仅当
(I−A+A)z=0时取等号,故A^+y是原方程唯一的极小范数最小二乘解。
总结
这篇博客用尽可能少的额外概念对广义逆矩阵作了一个大体的介绍,广义逆矩阵是描述线性方程组解的结构最有力的工具之一,故博客还利用广义逆矩阵提供了解决机器学习中的线性回归问题的一个途径。解决线性回归问题的另一个途径是最小化线性回归的代价函数:利用矩阵微积分可以证明使得代价函数最小的参数向量
x满足正规方程
ATAx=ATy,其中A是数据集的设计矩阵(每一行是一个样本的特征向量,每个样本的第一个特征均为1,即A的第一列都是1),y是标记向量。若A列满秩,则
x=(ATA)−1ATy。通常,如果不存在线性相关的特征列且样本数量大于特征数量的话,A就是列满秩的,故对于大多数情况
x=(ATA)−1ATy是可行的。