《矩阵的广义逆优秀课件.ppt》由会员分享，可在线阅读，更多相关《矩阵的广义逆优秀课件.ppt（19页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、矩阵的广义逆第1页，本讲稿共19页矩阵的广义逆矩阵的广义逆概述概述：矩阵的逆：矩阵的逆：A A n n n n ，B B n n n n ，B A=A B=I,B A=A B=I,则则则则B=A B=A 1 1 广义逆的目标：广义逆的目标：逆的推广逆的推广对一般的矩阵对一般的矩阵对一般的矩阵对一般的矩阵 A A mm n n可建立部分逆的性质。可建立部分逆的性质。可建立部分逆的性质。可建立部分逆的性质。当矩阵当矩阵当矩阵当矩阵A A n n n n可逆时，广义逆与逆相一致。可逆时，广义逆与逆相一致。可逆时，广义逆与逆相一致。可逆时，广义逆与逆相一致。可以用广义逆作求解方程组可以用广义逆作求解方
2、程组可以用广义逆作求解方程组可以用广义逆作求解方程组AX=bAX=b的理论分析。的理论分析。的理论分析。的理论分析。第2页，本讲稿共19页 4.1 矩阵的左逆与右逆矩阵的左逆与右逆一、满秩矩阵和单侧逆一、满秩矩阵和单侧逆1、左逆和右逆的定义、左逆和右逆的定义定义定义定义定义4 4.1 1(P P.93)A A C C m m n n，B B C C n n mm，BA=IBA=In n，则称矩阵则称矩阵则称矩阵则称矩阵B B 为矩阵为矩阵为矩阵为矩阵A A 的的的的左左左左逆，记为逆，记为逆，记为逆，记为 B=B=。例题例题1 矩阵矩阵A的左逆的左逆A=A=。A A C C m m n n ，
3、C C C C n n mm ，AC=IAC=Imm，则称矩阵则称矩阵则称矩阵则称矩阵C C 为为为为 矩阵矩阵矩阵矩阵A A 的的的的右右右右逆，记为逆，记为逆，记为逆，记为 C=C=。第3页，本讲稿共19页 2、左逆和右逆存在的条件、左逆和右逆存在的条件 的存在性的存在性直观分析直观分析直观分析直观分析 存在存在矩阵矩阵A A列满秩列满秩 =（A AH HA A）1A AH H 定理定理定理定理4 4.1.1(P.93)设设设设A A C C mm n n ，下列条件等价，下列条件等价，下列条件等价，下列条件等价1.1.A A左可逆左可逆左可逆左可逆2.2.A A的零空间的零空间的零空间的
4、零空间N N（A A）=0=0。3.3.mm n n，秩（秩（秩（秩（A A）=n=n，即矩阵即矩阵即矩阵即矩阵A A是列满秩的。是列满秩的。是列满秩的。是列满秩的。4.4.矩阵矩阵矩阵矩阵A AHH A A可逆。可逆。可逆。可逆。例题例题例题例题2 2 求矩阵求矩阵求矩阵求矩阵A=A=的左逆。的左逆。的左逆。的左逆。第4页，本讲稿共19页矩阵右逆的存在性矩阵右逆的存在性矩阵右逆的存在性矩阵右逆的存在性定理定理定理定理4.24.2 (P P.9494)A A C C m m n n ，则下列条件等价：则下列条件等价：则下列条件等价：则下列条件等价：1.1.矩阵矩阵矩阵矩阵A A右可逆。右可逆。
5、右可逆。右可逆。2.2.A A的列空间的列空间的列空间的列空间R R（A A）=C=Cmm3.3.n n mm ，秩（秩（秩（秩（A A）=m=m，A A是行满秩的。是行满秩的。是行满秩的。是行满秩的。4.4.矩阵矩阵矩阵矩阵A AA AH H 可逆可逆可逆可逆 =A AH H（AAAAH H）1 1讨论：可逆矩阵讨论：可逆矩阵讨论：可逆矩阵讨论：可逆矩阵A An n n n的左、右逆和逆的关系的左、右逆和逆的关系的左、右逆和逆的关系的左、右逆和逆的关系 可逆矩阵可逆矩阵可逆矩阵可逆矩阵A A的左、右逆就是矩阵的左、右逆就是矩阵的左、右逆就是矩阵的左、右逆就是矩阵A A的逆的逆的逆的逆A A
6、A A 1 1=（A AHHA A）1 1A AHH=A AH H（AAAAH H）1 1第5页，本讲稿共19页二、单侧逆和求解线性方程组二、单侧逆和求解线性方程组AX=b讨论讨论AX=b AX=b 有解与左、右逆存在的关系。有解与左、右逆存在的关系。有解与左、右逆存在的关系。有解与左、右逆存在的关系。借助于左、右逆求借助于左、右逆求借助于左、右逆求借助于左、右逆求AX=bAX=b的形如的形如的形如的形如X=BbX=Bb的解。的解。的解。的解。1、右可逆矩阵、右可逆矩阵定理定理定理定理4 4 4 4 (P.95)1.A A C C m m n n右可逆，则右可逆，则右可逆，则右可逆，则 b b
7、 C Cmm，AX=bAX=b有解。有解。有解。有解。2.X=b X=b 是方程组是方程组是方程组是方程组AX=bAX=b的解。的解。的解。的解。第6页，本讲稿共19页二、单侧逆和求解线性方程组二、单侧逆和求解线性方程组AX=b2 2、左可逆矩阵、左可逆矩阵、左可逆矩阵、左可逆矩阵求解分析：求解分析：求解分析：求解分析：定理定理定理定理4 4 3 3 (P P.9494)设矩阵设矩阵设矩阵设矩阵A A C C m m n n左可逆，左可逆，左可逆，左可逆，B B是矩阵是矩阵是矩阵是矩阵A A的的的的任何一个左逆，则任何一个左逆，则任何一个左逆，则任何一个左逆，则1.1.AX=bAX=b有形如有
8、形如有形如有形如X=BbX=Bb的解的充要条件是的解的充要条件是的解的充要条件是的解的充要条件是 (I Imm ABAB )b=0 b=0 ()()2.2.当当当当()()式成立时，方程组的解是惟一的，而且惟一解是式成立时，方程组的解是惟一的，而且惟一解是式成立时，方程组的解是惟一的，而且惟一解是式成立时，方程组的解是惟一的，而且惟一解是X=X=（A AH HA A）1 1A AH H b b证明：证明：证明：证明：讨论讨论讨论讨论：对任何满足式对任何满足式对任何满足式对任何满足式()的左逆的左逆的左逆的左逆B B，X=BbX=Bb都是方程组的都是方程组的都是方程组的都是方程组的 解，如何解释
9、方程组的解是惟一的？解，如何解释方程组的解是惟一的？解，如何解释方程组的解是惟一的？解，如何解释方程组的解是惟一的？第7页，本讲稿共19页 4.2 广义逆矩阵广义逆矩阵思想思想：用公理来定义广义逆。用公理来定义广义逆。一、减号广义逆一、减号广义逆定义定义定义定义4 4.2 2 (P.95)A A C C m m n n ，如果，如果，如果，如果，GG C C n n mm使得，使得，使得，使得，AGA=AAGA=A，则矩阵则矩阵则矩阵则矩阵GG为的为的为的为的A A减号广义逆。或减号广义逆。或减号广义逆。或减号广义逆。或11逆。逆。逆。逆。A A的减号逆集合的减号逆集合的减号逆集合的减号逆集合
10、A1=AA1=A1 1 1 1，A A2 2 1 1，A Ak k 1 1 例题例题例题例题1 1 A A C C n n n n可逆，则可逆，则可逆，则可逆，则A A 1 1 A1A1；A A单侧可逆，则单侧可逆，则单侧可逆，则单侧可逆，则A A 1 1L L A1A1；A A 1 1R R A1A1。减号逆的求法：减号逆的求法：减号逆的求法：减号逆的求法：定理定理定理定理4 4.5.5（P P.9595）减号逆的性质：减号逆的性质：减号逆的性质：减号逆的性质：定理定理定理定理4 4.6 6 (P.96)第8页，本讲稿共19页二、二、二、二、Moore-Penrose（M-P）广义逆广义逆由
11、由由由Moore 1920Moore 1920年提出，年提出，年提出，年提出，19551955年由年由年由年由PenrosePenrose发展。发展。发展。发展。1 1、定义定义定义定义4 4.3 3(P.98)设矩阵设矩阵设矩阵设矩阵A A C C m m n n ，如果，如果，如果，如果 GG C C n n m m，使得，使得，使得，使得1.1.AGA=A AGA=A 2.2.GAG=G GAG=G3.3.（AGAG）H H=AG=AG4.4.（GAGA）H H =GA=GA 则称则称则称则称GG为为为为A A的的的的M-PM-P广义逆，记为广义逆，记为广义逆，记为广义逆，记为G=G=A
12、 A+。A A1 1=A=A+；A A11L L=（A AHHA A）11A AHH=A=A+；A A 11R R=A=AHH（AAAAHH）11=A=A+；若若若若 A A+，则则则则A A+是是是是 A1 A1。例题例题例题例题2 2 讨论原有的逆的概念和讨论原有的逆的概念和讨论原有的逆的概念和讨论原有的逆的概念和M-PM-P广义逆的关系。广义逆的关系。广义逆的关系。广义逆的关系。第9页，本讲稿共19页3、M-PM-P广义逆的存在性及其求法广义逆的存在性及其求法广义逆的存在性及其求法广义逆的存在性及其求法定理定理定理定理4 4.8 8（P P.9999）任何矩阵都有任何矩阵都有任何矩阵都有
13、任何矩阵都有M-PM-P广义逆。广义逆。广义逆。广义逆。求法求法求法求法：设设设设A A满秩分解满秩分解满秩分解满秩分解A=BCA=BC，则则则则A A+=C=CH H（CCCCH H ）1 1（B BH H B B）1 1B BH H。（定理定理定理定理4 4.9.9）设）设）设）设A A奇异值分解奇异值分解奇异值分解奇异值分解 ：，则，则，则，则2、M-P 广义逆的惟一性广义逆的惟一性定理定理定理定理4 4.9.9(P.98)如果如果如果如果A A有有有有M-PM-P广义逆，则广义逆，则广义逆，则广义逆，则A A的的的的 M-P M-P广义逆是惟一的。广义逆是惟一的。广义逆是惟一的。广义逆
14、是惟一的。第10页，本讲稿共19页例题例题1 求下列特殊矩阵的广义逆；求下列特殊矩阵的广义逆；零矩阵零矩阵零矩阵零矩阵0 0；1 1阶矩阵阶矩阵阶矩阵阶矩阵(数数数数)a a；对角矩阵对角矩阵对角矩阵对角矩阵 例题例题例题例题3 3 设设设设 ，求求求求A A+。0 0+mm n n =0=0 n n mm 例题例题2 设向量设向量 的的M-P广义逆。广义逆。.第11页，本讲稿共19页4、M-P广义逆的性质广义逆的性质定理定理定理定理4.12 4.12(P.100)：则：则：则：则A A满足下列性质：满足下列性质：满足下列性质：满足下列性质：1.1.（A A+）+=A=A2.2.（A A+）H
15、 H=（A A HH ）+3.3.（A A）=+A A+4.4.A A列满秩，则列满秩，则列满秩，则列满秩，则A A+=（A A HH A A）1A A HH ，A A行满秩，则行满秩，则行满秩，则行满秩，则A A+=A=AHH （AAAAHH）1。5.5.A A有满秩分解：有满秩分解：有满秩分解：有满秩分解：A=BCA=BC，则则则则A A+=C=C+B B+。A A+与与与与A A1 1 性质的差异比较：性质的差异比较：性质的差异比较：性质的差异比较：（ABAB）11=B=B 11 A A 11 ，一般不成立一般不成立一般不成立一般不成立（ABAB）+=B=B+A A+。（只有满秩分解成立
16、）只有满秩分解成立）只有满秩分解成立）只有满秩分解成立）（A A11）k k=（A Ak k）1 1，但不成立（但不成立（但不成立（但不成立（A A+）k k=（A Ak k）+第12页，本讲稿共19页 4.3 投影变换投影变换（为讨论（为讨论（为讨论（为讨论A A+的应用做准备）的应用做准备）的应用做准备）的应用做准备）问题：问题：问题：问题：逆在什么情形下是有用的？逆在什么情形下是有用的？逆在什么情形下是有用的？逆在什么情形下是有用的？一、投影变换和投影矩阵一、投影变换和投影矩阵一、投影变换和投影矩阵一、投影变换和投影矩阵定义定义定义定义4 4.4 4（P P.101101）设设设设C C
17、n n=L=L MM ，向量向量向量向量x x C Cn n,x=y+z x=y+z，y y L L，z z M M，如果线性变换如果线性变换如果线性变换如果线性变换 ：C C n nC Cn n ，（x x）=y=y，则称则称则称则称 为从为从为从为从 C Cn n 沿子空间沿子空间沿子空间沿子空间MM到子空间到子空间到子空间到子空间L L的投影变换。的投影变换。的投影变换。的投影变换。投影变换的矩阵投影变换的矩阵投影变换的矩阵投影变换的矩阵R R（）=L=L；N N（）=M=M，C Cn n=R=R（）N N（）L L和和和和MM是是是是 的不变子空间；的不变子空间；的不变子空间；的不变子
18、空间；L L=I=I；MM=0=0投影的矩阵和变换性质投影的矩阵和变换性质投影的矩阵和变换性质投影的矩阵和变换性质：1.1.定理定理定理定理4 4.13.13（P P.101101）是投影是投影是投影是投影 是幂等变换是幂等变换是幂等变换是幂等变换2.2.推论推论推论推论：为投影变换的充要条件是变换矩阵是为投影变换的充要条件是变换矩阵是为投影变换的充要条件是变换矩阵是为投影变换的充要条件是变换矩阵是 幂等矩阵幂等矩阵幂等矩阵幂等矩阵第13页，本讲稿共19页二、正交投影和正交投影矩阵二、正交投影和正交投影矩阵1.正交投影的定义：正交投影的定义：定义定义定义定义4 4.5 5（P P.103103
19、）设设设设 ：C C n nC Cn n 是投影变换，是投影变换，是投影变换，是投影变换，C C n n=R R（）N N（），），），），如果如果如果如果 R R （）=N=N（），），），），则称为正交投影。则称为正交投影。则称为正交投影。则称为正交投影。2 2 正交投影矩阵正交投影矩阵正交投影矩阵正交投影矩阵定理定理定理定理4 4.14.14（P P.103103）是正交投影是正交投影是正交投影是正交投影 投影矩阵投影矩阵投影矩阵投影矩阵A A满足：满足：满足：满足：A A 2 2=A=AA AHH=A=A例题例题例题例题1 1 设设设设WW是是是是C C n n 的子空间，证明的子空间
20、，证明的子空间，证明的子空间，证明 存在到存在到存在到存在到WW的投影变的投影变的投影变的投影变换，换，换，换，使使使使R R（）=W=W。第14页，本讲稿共19页3、正交投影的性质、正交投影的性质定理定理定理定理4 4.16.16（P P.104104）设设设设WW是是是是C C n n的子空间，的子空间，的子空间，的子空间，x x0 0 C C n n，x x 0 0 W W，如果如果如果如果 是空间是空间是空间是空间C C n n向空间向空间向空间向空间WW的正交投影，则的正交投影，则的正交投影，则的正交投影，则含义：含义：含义：含义：点点点点 （x x0 0）是空间是空间是空间是空间
21、W W 中与点中与点中与点中与点x x 0 0距离最近的点。距离最近的点。距离最近的点。距离最近的点。第15页，本讲稿共19页4、A+A与与AA+的性质的性质定理定理4.15（P.104）A A+A A的性质：的性质：的性质：的性质：（A A+A A）2 2=A=A+A A，（，（，（，（A A+A A）HH=A =A+A A C C n n=R=R（A A+）N N（A A）R R （A A+）=N=N（A A）A A A A+的性质：的性质：的性质：的性质：（A A A A+）2 2=A A=A A+，（，（，（，（A A A A+）HH=A A =A A+C C mm=R=R（A A）N
22、 N（A A+）R R （A A+）=N=N（A A）A A+A A 是正交投影，将向量是正交投影，将向量是正交投影，将向量是正交投影，将向量 x x 投影到空间投影到空间投影到空间投影到空间R R（A A+）中。中。中。中。A A A A+是正交投影，将向量是正交投影，将向量是正交投影，将向量是正交投影，将向量 x x 投影到空间投影到空间投影到空间投影到空间R R（A A）中。中。中。中。含义：含义：含义：含义：第16页，本讲稿共19页4.4 最佳最小二乘解最佳最小二乘解一、最佳最小二乘解一、最佳最小二乘解A mn X n 1=b m1有解有解有解有解b b R R（A A）无解无解无解无
23、解b b R R（A A）1 1、AX=bAX=b的最佳最小二乘解的最佳最小二乘解的最佳最小二乘解的最佳最小二乘解定义定义定义定义4 4.6.6（P P.105105）u u 是最小二乘解是最小二乘解是最小二乘解是最小二乘解 x x0 0是最佳最小二乘解是最佳最小二乘解是最佳最小二乘解是最佳最小二乘解2 2、AX=bAX=b的最佳最小二乘解的计算的最佳最小二乘解的计算的最佳最小二乘解的计算的最佳最小二乘解的计算定理定理定理定理4 4.17.17 设方程组设方程组设方程组设方程组AX=bAX=b，则则则则A A+b b 是是是是AX=b AX=b 的最佳最小二乘解的最佳最小二乘解的最佳最小二乘解
24、的最佳最小二乘解。例题例题例题例题1 1 （P P.106106，eg8eg8）第17页，本讲稿共19页例题例题2、设，、设，=1.证明证明 R（A）2.在列空间在列空间R（A）上找一点上找一点X0，X0距离距离 最近。最近。第18页，本讲稿共19页 二、最佳拟合曲线二、最佳拟合曲线问题问题：在实际问题中，已知变量：在实际问题中，已知变量X和变量和变量Y之间之间存在函数关系存在函数关系Y=F（X），），但不知道但不知道F（X）的具的具体形式，由观察和实验数据寻求经验公式体形式，由观察和实验数据寻求经验公式:Y=f（X），），使得误差最小。使得误差最小。例题例题1（P P.107107，eg9e
25、g9）一组实验数据一组实验数据一组实验数据一组实验数据（1 1，2 2），（），（），（），（2 2，3 3），），），），（3 3，5 5），（），（），（），（4 4，7 7）的分布呈直线趋势，求最佳拟合直线。）的分布呈直线趋势，求最佳拟合直线。）的分布呈直线趋势，求最佳拟合直线。）的分布呈直线趋势，求最佳拟合直线。方法：方法：将误差向量表示为将误差向量表示为将误差向量表示为将误差向量表示为 e e =A=A b b，求方程组求方程组求方程组求方程组 A A b=0b=0的最小二乘解的最小二乘解的最小二乘解的最小二乘解 ，由，由，由，由 给出拟合参数。给出拟合参数。给出拟合参数。给出拟合参数。第19页，本讲稿共19页