广义逆矩阵A+

前言

非奇异的"方阵"才能叫"可逆矩阵",奇异的"方阵"叫"不可逆矩阵"!这是基础线性代数中的概念。现代应用中为了解决各种线性方程组(系数矩阵是非方阵方阵为奇异),将逆矩阵的概念推广到"不可逆方阵"和"长方形矩阵"上,从而产生了"广义逆矩阵"的概念!归根结底逆矩阵的初衷不变:为了解线性方程组。有了广义逆矩阵之后,可以给出线性方程组(包括相容的和矛盾的方程组)各种"解"的统一描述

所以,广义逆矩阵把最基础、最原始的方阵逆矩阵情况也包含其中。

广义逆矩阵的概念

定义1:设矩阵A\in C^{m\times n},如果有矩阵X\in C^{\color{red}{n}\times \color{red}{m}}满足下列4个Penrose方程

  • AXA = A
  • XAX = X
  • (AX)^{H} = AX
  • (XA)^{H} = XA

中的几个或者全部,则称矩阵XA的广义逆矩阵!如果4个方程全都满足,则称矩阵X为矩阵AMoore-Penrose逆。显然,如果A可逆矩阵(方阵),那么X = A^{-1}一定4个方程都满足。

按照上面的定义,可以分为满足1个、2个、3个、4个方程的广义逆矩阵,共有15种:

C^{1}_{4} + C^{2}_{4} + C^{3}_{4} + C^{4}_{4} = 15

定理1:设A\in C^{m\times n},则AMoore-Penrose逆存在且唯一

定义2:我们一般把AMoore-Penrose逆记为A^{+},也称为A的加号逆。

按满足条件的不同,15种广义逆中最常用的是下面5个:

A\{1\}、A\{1,2\}、A\{1,3\}、A\{1,4\}、A^{+}

不同的广义逆在不同的特定问题中可能有各自特殊的用途,但是作为一定存在且满足条件方程数最多加号逆,是最重要的!其在解线性方程组中的应用是本文关注的重点!

加号逆求法

求加号逆的方法有两种:奇异值分解(SVD)满秩分解(FRD)。其中FRD的实现要简单的多!下面给出用满秩分解求加号逆的定理:

定理2:设A\in C^{m\times n}_{r},它的满秩分解为:A = FG,则加号逆的求法为:

A^{+} = \color{red}{G^{H}(GG^ \color{blue}{H})^{-1}(F^{ \color{blue}{H}}F)^{-1}F^H}

2点说明:

2条推理:

  • 若A为行满秩,则A^{+} = A^{H}(AA^H)^{-1}
  • 若A为列满秩,则A^{+} = (A^HA)^{-1}A^{H}

3条重要性质:

  • (A^{+})^{+} = A
  • (A^{+})^{H} = (A^{H})^{+}\quad (A^{+})^{T} = (A^{T})^{+}
  • rank(A)=rank(A^{+})=rank(AA^{+})=rank(A^{+}A)

加号逆求解线性方程组

首先回顾一下任意线性方程组的解的判断:

不论正方形还是长方形方程组Ax=b,都可以用""来判断方程解的情况:

Ax=b → \begin{cases} \color{red}{无解} & R(A) < R(A,b) \\ 唯一解 & R(A) = R(A,b) = n \\ \color{red}{无穷解} & R(A) = R(A,b) < n \end{cases}

说明:

  • 其中有解的方程组又可称为"相容方程组",无解的方程组又可称为"矛盾方程组"!
  • 相容方程组包括唯一解和无穷解的情况;
  • 矛盾方程组的最小二乘解一般是不唯一的!所有最小二乘解中的极小(2)范数解是唯一的!所以我们一般求矛盾方程组的唯一极小范数最小二乘解

下面给出4条关于加号逆在解线性方程组中的作用:

  • Ax=b有解(或相容)的充分必要条件为:AA^{+} = b
  • x = \color{red}{A^{+}b} + (I-A^{+}A)y\quad (y \in C^{n}任意)是相容方程组Ax=b通解;或是矛盾方程组Ax=b全部最小二乘解
  • x_{0} = \color{red}{A^{+}b}是相容方程组Ax=b通解中唯一极小范数解;或是矛盾方程组Ax=b唯一极小范数最小二乘解
  • 若相容线性方法组Ax=b唯一解(方阵),那么这个解就是:x_{0} = \color{red}{A^{+}b}

下面对上面的4条说明进行举例验证:


例1(方程无解):用广义逆矩阵方法判断下面的线性方程组是否有解。如果有唯一解,则求这个唯一解;如果有通解,则求通解和通解中的唯一极小范数解;如果无解,求全部的最小二乘解和其中唯一的极小范数最小二乘解:

\begin{cases} 2x_1 + 4x_2 + x_3 + x_4 = 10 \\ x_1 + 2x_2 - x_3 + 2x_4 = 6 \\ -x_1 - 2x_2 - 2x_3 + x_4 = -7 \end{cases}

首先用""来判断方式解的情况:

rank(A) = 2 < rank([A,b]) = 3

所以是无解的!下面我们再用加号逆来判断一下解的情况:

A^{+} = \left( \begin{matrix} 0.0606 & 0.0303 & -0.0303 \\ 0.1212 & 0.0606 & -0.0606 \\ 0.0303 & -0.1515 & -0.1818 \\ 0.0303 & 0.1818 & 0.1515 \end{matrix} \right)

更具加号广义逆的方程解的判断:

AA^{+}b = (11,5,-6)^{T} \neq b

所以原方程是无解的!也说明用判断和用广义逆判断结果是一致的。

既然方程无解,那么我们可以求它的全部最小二乘解和其中的唯一的极小范数最小二乘解,根据上面4条的说明,先求全部最小二乘解

x = \color{red}{A^{ +}b} + (I-A^{+}A)y = \left( \begin{matrix} 1 \\ 2 \\ \frac{2}{3} \\ \frac{1}{3} \end{matrix} \right)+ \frac{1}{11}\left( \begin{matrix} 9 & -4 & -1 & -1 \\ -4 & 3 & -2 & -2 \\ -1 & -2 & 5 & 5 \\ -1 & -2 & 5 & 5 \end{matrix} \right) \left( \begin{matrix} y_1 \\ y_2 \\ y_3 \\ y_4 \end{matrix} \right)

其中y_1,y_2,y_3,y_4 \in C 任意 ,其中唯一的极小范数最小二乘解为:

x_0 = \color{red}{A^{+}b} = (1,2,\frac{2}{3},\frac{1}{3})^{T}


例2(方程无穷解):求下面这个有无穷解的线性方程组的通解和唯一极小范数解

\begin{cases} x_1 - x_2 - x_3 + 4x_4 = 0 \\ x_1 - x_2 + x_3 - 3x_4 = 1 \\ x_1 - x_2 - 2x_3 + 3x_4 = -0.5 \end{cases}

还是先用判断一下解的情况:

rank(A) = 3 = rank([A,b]) < n = 4

所以方程组有无穷解。下面用广义加号逆来求通解

x = \color{red}{A^{+}b} + (I-A^{+}A)y = \left( \begin{matrix} 0.25 \\ -0.25 \\ 0.5 \\ 0 \end{matrix} \right)+ \frac{1}{11}\left( \begin{matrix} 0.5 & 0.5 & 0 & 0 \\ 0.5 & 0.5 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{matrix} \right) \left( \begin{matrix} y_1 \\ y_2 \\ y_3 \\ y_4 \end{matrix} \right)

其中y_1,y_2,y_3,y_4 \in C 任意 ,其中唯一的极小范数解为:

x_0 = \color{red}{A^{+}b} = (0.25,-0.25,0.5,0)^{T}


例3(方程有唯一解)

\begin{cases} 5x_1 + 2x_2 + 0x_3 = 3 \\ -x_1 + 4x_2 + 2x_3 = -3 \\ 2x_1 - 3x_2 + 10x_3 = 15 \end{cases}

还是先用""判断一下方程解的情况:

rank(A) = 3 = rank([A,b]) = n

所以方程有唯一解,下面我们用广义加号逆方法来求这个唯一解:

x_0 = \color{red}{A^{+}b} = (1,-1,1)^{T}


最后补充说明4点:根据很好理解

  • 方程个数 > 未知数个数,一般是无解(有无穷解的可能);
  • 方程个数 < 未知数个数,一般是无穷解(有无解的可能);
  • 方程个数 = 未知数个数,3种可能都有!但只有方阵才可能有唯一解出现;
  • 不论长长方形、宽长方形、正方形;不论原方程组是有无穷解、唯一解、无解;广义加号逆都有适用的公式!万能

综上,可以把广义加号逆当成是基础的方阵逆完全拓展!即广义加号逆适用于各种形状的矩阵求逆!并且:当矩阵为方阵且可逆时,广义加号逆 = 方阵直接逆inv(A)

本文相关程序下载地址

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,547评论 6 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,399评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,428评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,599评论 1 274
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,612评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,577评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,941评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,603评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,852评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,605评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,693评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,375评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,955评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,936评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,172评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 43,970评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,414评论 2 342

推荐阅读更多精彩内容