广义逆矩阵求法19页PPT
第4章 矩阵的广义逆
![第4章 矩阵的广义逆](https://img.taocdn.com/s3/m/9b83f0e19e3143323868930f.png)
定义 3 设 A 为一个 m n 复矩阵,若有一个 n m 复矩阵 G 存在, 使( 1 )成立,即 AGA A ,则称 G 为 A 的一个 {1}-广义逆,记为
G A{1} 或 G A{1} ,也称 G 为 A 的一个减号广义逆,记为 G A , 即有 AA A A . (5)
A为列满秩
7
推论 设 A C mn , 则
(1) A左可逆的充要条件是 N ( A) {0};
( 2) A右可逆的充要条件是 R( A) C m .
证 充分性:N ( A) {0}
rank ( A) n
必要性: A左可逆
Ax 0只有零解
A为列满秩
1 ALபைடு நூலகம்A En
x N ( A)
由于 M-P 的 4 个方程都各有一定的解释,并且应用起来各有方 便之处,所以出于不同的目的,常常考虑满足部分方程的 G ,总之, 按照定义 2 可推得,满足 1 个,2 个,3 个,4 个 M-P 方程的广义逆 矩阵共有 15 类,即
1 2 3 4 C4 C4 C4 C4 15 .
使得
AGb b ( b R( A))
m n
则称G为A的广义逆矩阵 , 记为G A .
定理1设 A C
, 则A 存在广义逆矩阵A 的
充要条件是存在 G C nm , 使其满足AGA A
14
定理1 设 A C
m n
, 则A 存在广义逆矩阵A 的
nm
充要条件是存在 G C
15
由AGA A可得: AGAx0 Ax0 b 即,AGb b, 说明x Gb是方程 Ax b 的解. G是A的减号逆 , G A . m n nm 设 A C , 且 A C 是A的一个广义 推论 1 逆矩阵A , 则
高等数学逆矩阵ppt课件
![高等数学逆矩阵ppt课件](https://img.taocdn.com/s3/m/36d958e0af1ffc4fff47ac68.png)
268.
例7: 设方阵A满足矩阵方程 A2–A–2E = O, 证明: A, A+2E 都可逆, 并求它们的逆矩阵.
证明: 由 A2–A–2E=O, 得 A(A–E)=2E,
则
A
1
(
A
E
)
A1 E,
故A可逆, 且A-1 = 1 ( A E ).
2
2
又由 A2–A–2E=O, 得 (A+2E)(A–3E)+4E=O,
1 3
2, A12
2 3
1 3
3, A13
2 3
2 4
2,
同理可得 A21 6, A22 6, A23 2,
A31 4, A32 5, A33 2. 所以,
A
2 3
2
6 6
2
4 5
,
2
故
A1
|
1 A A|
1 3
1
2
3 3
1
5
122.
7
例3: 下列矩阵A,B是否可逆? 若可逆, 求其逆矩阵.
由伴随矩阵的性质: AA*= A*A = | A | E, 知
当| A | 0时,
A 1 A 1 A A E, | A| | A|
按逆矩阵的定义得, A1
1
A .
| A|
当| A | = 0 时, 称A为奇异矩阵, 否则称A为非奇异
矩阵.
4
由此可得, A是可逆矩阵的充分必要条件是A为非 奇异矩阵.
§2.3 逆 矩 阵
一、逆矩阵的概念和性质
在数的运算中, 当数 a 0 时, 有 aa-1 = a-1a = 1.
其中 a1 1 为a 的倒数, 或称a的逆(元). a
矩阵广义逆求法
![矩阵广义逆求法](https://img.taocdn.com/s3/m/32a0f7e319e8b8f67c1cb92b.png)
推论:设A Cmn, 1 )当r n时(列满秩),A+ =(AH A)1 AH, r 2)当r m时(行满秩),A+ =AH (AAH )1.
证明:当r n时, A AI n BC为A的满秩分解,由定理可得 1) A + =(A H A) 1 A H; 2)当r m时, A I m A BC ,
54 75 33 1 32 43 1 BH 234 130 182 78 34 53 15
3 1 2)因B列满秩,B =(B B) B = (1,3) (1,3) (1,3). 10 1
注:当r ( A) 1时,非零特征值只有1个,则A =
+
1
1
AH .
此时,设A=(aij ) mn , r ( AH A) r ( A) 1, 可知 1 n 0 H tr ( A A) 1 =tr j 1 0
a
i 1
利用奇异值分解求A +的简化步骤: 1)求出A H A的r个非零特征值1, , r , i 0; 2)求出A H A对应于特征值1, , r的标准正交特征向量
1 , , r .令V1 =(1 , , r );
1-1 H H + 3)则A V1 V1 A r-1
1 2 1 解:A 0 1 -1 0 1 -1 1 2 1 B= 2 5 ,C= 0 4 9 5 1 0 3 -3 4 0 1 -1 4 ,令 4 0 0 0 0 0 3 -3 . 1 -1 4
上述定理可简化:
Sr 定理3:设A C ,在A奇异值分解A=U 0 令V=(x1 , , x r ,x r 1 , ,x n )=(V1 ,V2 ).则
第八章 矩阵的广义逆
![第八章 矩阵的广义逆](https://img.taocdn.com/s3/m/7c8af1a2b0717fd5360cdc7d.png)
第八章矩阵的广义逆前言初等变换和标准形初等变换和标准形举例
§8.1 广义逆矩阵减号逆的概念
减号逆存在定理及求法减号逆存在定理及求法续
关于减号逆公式的注一个减号逆确定所有减号逆1减号逆的主要性质续减号逆的主要性质续
减号逆的主要性质续左逆与右逆的概念矩阵左逆与右逆的求法自反广义逆的概念
自反广义逆的存在与唯一性自反广义逆的唯一性自反广义逆与左(右)逆的关系用满秩分解求自反广义逆
自反广义逆的求法自反广义逆的求法续§8.2 伪逆矩阵
伪逆的存在性求伪逆举例
伪逆的唯一性
伪逆的性质
⎞
⎛−101求伪逆举例
§8.3 广义逆与线性方程组
一般矩阵方程有解的条件一般矩阵方程的通解
用减号逆求解相容线性方程组举例相容线性方程组的最小模解0130
−
相容方程组最小模解的充要条件
相容方程组最小模解的充要条件续
求相容方程组最小模解举例
Ax,即‖Ax-b‖>0.
不相容方程组的最小二乘解
R(A)
Ax 0
不相容方程组的最小二乘解举例用广义逆求最小二乘解定义8.3.2:线性方程组Ax=b 的一个最佳最小二乘
矩阵方程的最小二乘解。
矩阵的广义逆及其应用.ppt
![矩阵的广义逆及其应用.ppt](https://img.taocdn.com/s3/m/90f1ef0002768e9950e7386a.png)
第五章 矩阵的广义逆
§1 广义逆矩阵
(6) 若F是列满秩矩阵,则 F (F H F )1 F H
(7) 若G是行满秩矩阵,则 G GH (GGH )1
(8) 若矩阵A的满秩分解为A FG,则有 A G F ;
高等工程数学 理学院 杨文强
第五章 矩阵的广义逆
第五章 矩阵的广义逆
§1 广义逆矩阵 一、矩阵的广义逆
设A Rnn,对于线性方程组 Ax b,当A可逆时, 方程组有唯一解:x A1b.
若矩阵 A不可逆时,如何求解方程组 Ax b?
更一般,当矩阵 A Rmn不是方阵时,如何讨论 方程组 Ax b的解, 其中x Rn,b Rm ? 为了分析和解决上述问题,引入广义逆的概念.
高等工程数学 理学院 杨文强
第五章 矩阵的广义逆
§1 广义逆矩阵
定理2:设A Rmn,b Rm,x Rn,若性方程组 Ax b 是相容的,即方程组Ax b 有解,则其
通解为: x Ab (In A A)t,t是任意n 1向量. 证明:首先证明t Rn,x Ab (In A A)t是 方程组的解,然后证明方程组的任一解x,均可 表示成x Ab (In A A)t的形式.
A
1
1
1
2
(3)(1)3
0
3 3 2 4
0
1 2 4
0
1
2
0 4 8
高等工程数学 理学院 杨文强
第五章 矩阵的广义逆
§1 广义逆矩阵
1
A
0
0
1 2 4 (1)(2)2 1 1 0 0
第六章广义逆矩阵
![第六章广义逆矩阵](https://img.taocdn.com/s3/m/a3dd2716227916888486d711.png)
第六章广义逆矩阵§6.1 投影矩阵一、投影算子与投影矩阵v设L和M都是C n的子空间,且LÅM=C n.于是任意xÎC n都可唯一分解为x=y+z,yÎL,zÎM,称y是x沿着M到L的投影.v定义将任意xÎC n变为沿着M到L的投影的变换称为沿着M到L的投影算子,记为PL,M ,即PL,Mx=y。
v显然,R(P L,M)=L,N(P L,M)=M.v投影算子P L,M是一个线性算子。
v定义投影算子P L,M在C n的基e1,…,e n下的矩阵称为投影矩阵.记为P。
L,Mv幂等矩阵:A2=Av引理设AÎC n×n是幂等矩阵,则N(A)=R(I-A)。
证明:A2=AÞA(I-A)=OÞ对任意xÎR(I-A),存在yÎC n,x=(I-A)y,必有Ax=0。
故R(I-A)ÌN(A)Þdim R(I-A)£dim N(A)=n-dim R(A)即rank(I-A)£n-rank A。
考虑到I=A+(I-A)Þn£rank A+rank(I-A)有rank(I-A)=n-rank A,使得dim R(I-A)=n-dim R(A)=dim N(A),即得N(A)=R(I-A)。
v定理:P为投影矩阵的充要条件是P为幂等矩阵为投影矩阵,则对任意xÎC n有证明:设P=PL,MP2L,M x = P L,M (P L,M x) = P L,M y = y = P L,M x故P为幂等矩阵。
反之,设P为幂等矩阵n则对任意xÎC有x=x-Px+Px=(I-P)x+Px,其中(I-P)xÎN(P),PxÎR(P),使得C n=N(P)+R(P)。
设zÎN(P)∩R(P),由于N(P)=R(I-P)故存在u,vÎC n使得z=Pu=P2u=P(I-P)v Þz=Pu=(I-P)v=0故N(P)∩R(P)={0}。
矩阵的广义逆 ppt课件
![矩阵的广义逆 ppt课件](https://img.taocdn.com/s3/m/cc4fc660b8f67c1cfbd6b828.png)
意。
A-b为AX=b的特解, (In-A-A)Z为AX=0的通解.
E.H. Moore and Roger Penrose
二、Moore-Penrose (M-P) 广义逆
由Moore 1920年提出,1955年由Penrose独 立研究和发展。
1、 定义4.3 (P.98) 设矩阵 A Cmn,如果
A
L
1
的存在性
直观分析
➢
A
1 L
存在 矩阵A列满秩
➢
A
1 L
=
(AHA)–1AH
定理4.1(P.93) 设 A Cmn ,下列条件等价
1. A左可逆;
BA = In
2. A的零空间 N(A) = {0}; Ax = 0 x = BAx = 0
3. m n,秩(A) = n,即A是列满秩的; n-r(A) = 0
1. 矩阵A右可逆;
AC = Im
2. A的列空间 R(A) = Cm ; x = ACx
x R(A)
3. n m, 秩(A) = m, 即A是行满秩的r;(A)=dimR(A)
4. 矩阵 AAH 可逆,且
A
R
1
= AH(AAH)–1
r(AAH) = r(A)
讨论:可逆矩阵Ann的左、右逆和逆的关系 ➢ 可逆矩阵A的左、右逆就是矩阵A的逆A
GCnm ,使得
4.
矩阵AHA可逆,且
A
L
1
= (AHA)–1AH 。r(AHA) = r(A)
1 0
如前例 矩阵 A =
0
1
左可逆,AT右可
逆。如何求左或右逆? 2 1 可用行或列初等变换!
m r(A) C r(A )m ,n
广义逆矩阵求法
![广义逆矩阵求法](https://img.taocdn.com/s3/m/cdca1c826529647d272852cf.png)
广义逆矩阵
定理:设 阵方程
A 是数域 K 上一个s n 矩阵,则矩
AXA A
(1)
总是有解。如果 rank( A) r ,并且
( I nn A A) A A
A ( A ) 0
所以 X ( I nn A A) Z 是方程组 的通解。
AX 0
利用上述定理,可以得到非齐次线性方程组的 另一种形式的通解。
推论:设数域 K 是 n 元非齐次线性方程组 AX 有解,则它的通解为
即
I r B 1 Q P C D 先分析 Q 与 P 1 之间的关系。由已知 A ,
因此我们有
I r 0 1 0 0 Q P 1 分别把 Q , P 分块,设 行 Y1 }r Q Y2 }n r行
伪逆矩阵
定义:设 A C mn,若 A C nm ,且同时有
AA A A ,
H
A AA A
( AA ) AA ,
( A A) A A
H
则称 A 是 A 的伪逆矩阵。上述条件称为 Moore- Penrose 方程。 例:
1 设 A 0
取 B 0, D 0, C (0,,0, k Y ,0,,0)
1 i 2
则
Ir C
于是
0 1 Ir P C 0
1
逆矩阵PPT课件
![逆矩阵PPT课件](https://img.taocdn.com/s3/m/d2436e3faaea998fcc220eb9.png)
2、掌握逆矩阵的求法
学习重点:会判别逆矩阵是否存在;如何求逆矩阵 学习难点:熟练运用公式求逆矩阵
一、概念的引入
当数 在数的运算中, 其中 为 时, 有
的倒数, (或称
的逆);
在矩阵的运算中, 单位阵 相当于数的乘法运算中 的1, 那么,对于矩阵 , 如果存在一个矩阵 , 使得
证明: 若设
可得 所以 的逆矩阵是唯一的,即
(2)、
(3)、
证明
例
设
目前只能利用定义,用待定系数法解决!
解 则 设 是 的逆矩阵,
又因为
所以
显然当阶数大时,很繁!深切渴望好 方法!!!
三、逆矩阵的求法
定理 矩阵 可逆的充要条件是 ,且
牢 记 这 个 定 理
注:
现在有两种方法:待定系数法(略)和公式法。 例1 求方阵 的逆矩阵.
逆矩阵的计算方法
思考题
思考题解答
答Leabharlann 则矩阵 称为 的可逆矩阵或逆阵.
二、逆矩阵的概念和性质
1、 定义 对于 阶矩阵 ,如果有一个 阶矩阵 则说矩阵 是可逆的,并 ,使得 把矩阵 称为 的逆矩阵. 例 设
注意:
要同时成立!
现在要解决的问题:1. 方阵 满足什么条件时可逆? 2. 可逆时,逆阵怎样求?
2、性质 若 (1)、 是可逆矩阵,则 和 是 的逆矩阵是唯一的. 的可逆矩阵,则有
解
同理可得
用伴随阵求三阶以上 矩阵的逆阵计算量大
故
显然当阶数大时,还是有点很繁哦!继续 渴望好方法!!!
例2
解
用伴随阵求三阶以上矩阵的逆阵计算量大
例3
设
解
于是
广义逆矩阵
![广义逆矩阵](https://img.taocdn.com/s3/m/6b5d8404580216fc710afd1e.png)
1 0 0
10 0
1 0
0
0 1
0 1 0
0 0
0 1
把 y1 0,0,1,0T , y2 0,0,0,1T 扩充为R4 的一组标准正交基得:
y3 1,0,0,0T , y4 0,1,0,0T 再令 U y1, y2 , y3, y4 ,
设r 0,由满秩分解定理知,存在B Crmr ,C Crrn , 使得A BC
令X C H (CC H )1(BH B)1 BH
可以验证X满足广义逆矩阵方程
对于矩阵方程
几类弱逆
AXA A
(P1)
XAX X
AX H AX
(P2 ) ( P3 )
XAH XA (P4 )
则有唯一解 x A1b; 但当A是奇异方阵或长方矩阵时,它的解
不唯一,我们可以利用减号逆给出方程组的通解。
2)如果方程组相容,且其解有无穷多个,可求出具有极小范数的
解,即 min x , 其中 为欧氏范数,可以证明满足此条件 Axb
的解是唯一的,称为极小范数解。
3)若方程组不相容,则不存在通常意义下的解,但在许多实际 问题中,需要求出这样的解:
Al --最小二乘广义逆
A{1,4},它的形式记为
Am 最小范数广义逆
广义逆A-
A{1}是指仅满足第一个Penrose方程的广义逆,即若
AA-1A=A, 则记 A A{1}
说明: 1)利用初等行变换,可以求得A-
2)A的减号逆A-不唯一。
例:设
A
1 1
0 0,
1 0
第七章 广义逆矩阵
广义逆矩阵是逆矩阵的推广,与线性方程组的求解有密切 联系。给定一个线性方程组 Ax=b,当矩阵A可逆时,线性 方程组的解可表示为x=A-1 b
广义逆矩阵教案ppt课件
![广义逆矩阵教案ppt课件](https://img.taocdn.com/s3/m/f1d7f47e1a37f111f0855b86.png)
可逆方矩阵 .
返回
定 理 4设 A Cm n是左可 ,A L 1 逆 是 A 的 矩左 阵逆 矩 阵 ,则 方 程 Ax组 b有 解 的 充 要 条 件 是
(E mAL A 1)b0(1)
若(1)式成,则 立方程 Ax组 b有唯一解 x(AHA)AHb.
证 充分性: N(A){0} Ax0只有零解
ran(Ak)n A为列满秩
必要性: A左可逆
AL1AEn
返回
xN(A) xEnxAL 1(A)xAL100
N(A){0}
初等变换求左(右)逆矩阵:
(1)P(AEm)E0n G *
(2)EAnQEGm *0
返回
例 1 设矩阵 A为 1 2
A 0 1 0 0
证 : 必要性: 设 x0是 方A程 xb的 组解
(AL 1 A )A ( 0)x (AL 1 A )bA(AL1A)x0 AEnx0
Ax0 b
(EmAL A1)b0
返回
充分性: (EmAL A1)b0 x0 AL1b Ax0AAL1b b
唯一性:
设 x0,x1Βιβλιοθήκη A xb的解A (x 1 x 0 ) A 1 x A 0 x 0 x1x00
1 0 1
2
0
0
0
0
1
1 0 0
0 1 2
1
2
1
0 1 0
0
0
1
1 0 0 0 1 0
1
2
3
0 1 2
0
0
1
1 0 0
0 1 0
1
逆矩阵的计算ppt课件
![逆矩阵的计算ppt课件](https://img.taocdn.com/s3/m/8283cd2cae1ffc4ffe4733687e21af45b207fe52.png)
上页 下页 返回
26
用(8)代入(10),得 X = B( AX ) = ( BA )X
即有 BA = E。 于是有AB = BA = E。 具有这种性质的矩阵A称为是可逆的,而矩阵B 称为 矩阵A 的逆矩阵。
上页 下页 返回
27
证 AT ( A1 )T ( A1 A)T ET E.
上页 下页 返回
8
当| A | 0时, 定义 A0 E, Ak ( A1 )k ,
其中 k 为正整数。
当| A | 0,, 为整数时,有
A A A ,( A ) A .
上页 下页 返回
9
例9
1 求方阵 A 2
2 2
13 的 逆 阵.
上页 下页 返回
11
解 于是
A1
1 3
3 3
2 1
1
X A1CB1
52 , 21
B 1
3 5
21,
1 3
2 1
3 3 1
5221
1 2 3
3 0 1
3 5
21
1 0
0
1 2 2
3 5
21
2 10 10
1 4. 4
上页 下页 返回
12
矩阵的运算小结
一、已定义过的运算:
★矩阵与矩阵的加、减法; ★矩阵与数的乘积; ★矩阵与矩阵的乘积; ★方阵的行列式; ★逆矩阵; ★矩阵的转置。
Ex.4
设A
2 0
0 3
00 , 求A的 逆 矩 阵.
0 0 4
解 因| A| 24 0,故A可逆. 又
A11 12, A22 8, A33 6, Aij 0(i, j 1,2,3,且i j),
矩阵的逆及其求法PPT课件
![矩阵的逆及其求法PPT课件](https://img.taocdn.com/s3/m/02701466ec3a87c24128c464.png)
(A
E )
4
0
2
.
8 1 3
6 1 1
所以
(A
E )1
1 2
(A
E )
1 2
4 8
0 1
2
,
3
6 1 1 2 1 1 8 1
B
1 2
4 8
0 1
2
2
3 2
6 1
4
3
1 2
4 8
2 1
故
6
1 2
1
2
A B 4 7 3 .
2
3
11
2 2
1
2
.
充分发挥其作用,有必要对它进一步探讨。
定理3 A可逆 A 行 E Pm P2P1A E
Pm P2P1E A1 E 行 A1
求 A1方法 :( A E) 行(E A1)
21
第21页/共36页
例7 求下列矩阵的逆
矩阵 1 0 1
1. A 2 1 0 3 2 5
解:
1 0 1 1 0 0
A A1 E 1 0 , 因此 A 0.
充分性.设 A 0, 由定理 2.1 知
AA A A A E.
故有 A( 1 A* ) ( 1 A* )A E .
A
A
9
第9页/共36页
由逆矩阵定义知,A 可逆,且其逆为
A1 1 A* . A
定理 2.2 不仅给出了判断矩阵可逆的方法, 还给出了求解逆矩阵的一种方法 .
2
第2页/共36页
一、逆矩阵的概念
定义1 设 A 为 n 阶方阵,如有 n 阶方阵 B ,使 AB = BA = E .
则称 A 为可逆阵,B 为 A 的逆阵,记作 B A1 .
矩阵理论课件-第三章 矩阵的广义逆
![矩阵理论课件-第三章 矩阵的广义逆](https://img.taocdn.com/s3/m/9d8fc5aa6aec0975f46527d3240c844768eaa05c.png)
注2:由定理2知
A In
I
m
初等变换
PAQ Q
P
Ir 0 Q
0 0
P
1 0 -1 1
例:设A=
0
2
2
2 ,求A{1}.
-1 4 5 3
解:由
A I4
I3 0
初等变换
I2 0 Q
0 0
P ,这里
0
1 0 1 1
1 0
P=
0
1/ 2
1 2
0 0 1
,
这里只是给出了A{1}的一个构造性描述,在使用上并不直接, 因为还要求出一个A(1).
推论2:方程组(1)相容的充要条件是AA(1)b b,且其通解为 x=A(1)b+(I-A(1)A) y, y Cn任意.
证明:定理1中,取D=b Cm,B=1即得.
注1:因为A+ A{1},故Ax=b相容时,通解为 x=A+b+(I-A+A) y, y Cn.
证明:由A的奇异值分解(r(A)=r),有A=V
Sr 0
0 0
U
H,其中
Sr diag{1, , r},i 0,U和V是酉阵.
令G=U
Sr1 0
AGA=V
Sr 0
0 0
VH
,
可以验证G满足方程1)-4).如第1)3)方程
0 0
U
H
U
Sr1 0
0 0
VH
V
Sr 0
X=A(1) DB(1) +Y-A(1)AYBB(1).(2) 其中Y Cnq为任意.
证明::若AA(1)DB(1)B D,令X=A(1)DB(1)则满足AXB=D. :若AXB=D有解,则D=AXB=AA(1) AXBB(1) B=AA(1) DB(1) B.