中科院矩阵分析课件
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
矩阵分析及其应用
3.1 矩阵序列
定义3.1 设矩阵序列{A (k)},其中A (k)=()(k ij a )∈C m ⨯n ,当k →∞,
)(k ij
a →a ij 时,称矩阵序列{A (k)}收敛,并称矩阵A=(a ij )为矩 阵序列{A (k)}的极限,或称{A (k)}收敛于A, 记为
A A k k =∞
→)(lim 或 A (k)→ A 不收敛的矩阵序列称为发散的。
由定义,矩阵序列A (k) 发散的充要条件为存在ij 使
得数列)
(k ij a 发散。
类似地,我们可以定义矩阵收敛的Cauchy 定义
定义3.1' 矩阵序列{A (k)}收敛的充要条件为
对任给ε>0 存在N(ε), 当 k , l ≥ N(ε) 时有
||A (k)-A (l )|| < ε
其中||.||为任意的广义矩阵范数。 例1 ⎪⎪⎪⎪⎭
⎫ ⎝⎛-=∑=-n k n n k k e n n 12)()sin()1sin(11A 如果直接按定义我们因为求不出A (n )的极限从而
很难应用定义3.1证明收敛。 相反,由于∑∑∑+=+=+=-≤≤n
m k n m k n m k k k k k k 112
12
)1(11)sin( < 1/m
从而只要l 充分大,则当m, n > l 时就有
ε≤∑+=n
m k k k 1
2)sin( 这样A (l ) 收敛。
定理3.1 A (k)→ A 的充要条件为
||A (k) -A||→0
证明:利用广义矩阵范数的等价性定理,仅对∞范数可以证明。
即 c 1 ||A (k) -A||∞ ≤ ||A (k) -A||≤ c 2 ||A (k) -A||∞
性质0 若A (k)→ A , 则 ||A (k)|| → ||A|| 成立。
性质1. 设A(k)→ A m⨯n,B(k)→ B m⨯n, 则
α⋅ A(k)+β⋅ B(k)→α⋅ A+β⋅ B,∀α,β∈C
性质2. 设A(k)→ A m⨯n,B(k)→ B n⨯l, 则
A(k)⋅B(k)→ A⋅B
证明:由于矩阵范数地等价性,我们可以只讨论相容的
矩阵范数。
||A(k)⋅B(k)-A⋅B|| ≤ || A(k)⋅B(k)-A⋅B(k)||+||AB(k)- A⋅B||
≤ || A(k)-A||⋅||B(k)||+||A||⋅||B(k)-B||
注意||B(k)||→||B||,则结论可得。
特别地有
性质2’. A(k)→ A的充要条件为
A(k) x→Ax, 对任意x成立
或者y H A(k) x→ y H Ax, 对任意x,y成立.
(在无穷维空间中称为弱收敛,但在有限维空间中
和一般收敛性定义是等价的)
对于Hermite(对称)矩阵我们有如下的定理:
设A(k),k=1,2,…,和A都为Hermite矩阵,那么
A(k)→ A的充要条件为
x H A(k) x→x H Ax, 对任意x成立
推论:设A(k),k=1,2,…, 为半正定的Hermite矩阵,且单调减少,即A(k)和A(k)-A(k+1)为半正定Hermite矩阵,那么A(k)有极限.
性质3设A(k)和A都为可逆矩阵,且A(k)→A,则
(A(k))-1→A-1
证明:因为A-1⋅ (A(k)) →I. 所以存在K,当k >K时有
||I- A-1⋅ (A(k))||<1/2
我们有(A(k))-1= A-1+( I- A-1⋅ (A(k))) (A(k))-1
从而||(A(k))-1||≤||A-1||+||( I- A-1⋅ (A(k)))||⋅|| (A(k))-1||
当k>K时,有
||(A(k))-1||≤||A-1||+1/2⋅|| (A(k))-1||
即||(A(k))-1||≤2⋅||A-1||
因为A-1- (A(k))-1= A-1 (A(k)- A) (A(k))-1
从而|| A-1- (A(k))-1||≤||A-1||⋅||A(k)- A||⋅||(A(k))-1||
(当k>K时) ≤||A-1||⋅||A(k)- A||⋅2||A-1||
(当k→∞时) →0
由定理3.1有(A(k))-1→ A-1
定义3.2矩阵序列{A(k)}称为有界的,如果存在常数
M>0,使得对一切k都有
a| |)(k ij 定理:有界的矩阵序列{A (k)}一定有收敛的子列。 定义3.3 设A 为方阵,且当k →∞时有A k →0,则称A 为 收敛矩阵。 定理3.2(迭代法基本定理) A k →0的充要条件为谱半径 ρ(A)<1. 证明:必要性:设A k →0,证明ρ(A)<1. 对A 的任意特征值λ和相应的特征向量x 有 λx =Ax . 这样我们有A k x =λk x 从而有|λ|k ⋅||x||=||A k x||≤||A k ||⋅||x|| 从而有|λ|k ≤||A k ||→0 这样有|λ|<1, 由于λ为A 的任意特征值, 所以ρ(A)<1, 即必要性得证。 充分性。已知ρ(A)<1,证明A k →0. 取ε=(1-ρ(A))/2 >0,由定理2.10有,存在某种相容的 矩阵范数||.||M 使得 ||A||M < ρ(A)+ ε <1 从而||A k || M ≤(||A||M )k <(ρ(A)+ ε)k 所以当k →∞有||A k || M →0, 从而A k →0. 定理3.3 A k →0的充分条件为存在矩阵范数||.||M 使得 ||A||M <1 3.2矩阵级数 定义3.4设矩阵序列{A (k)},其中A (k)=()(k ij a )∈C n ⨯n ,由它们 形成的无穷和 A (0)+A (1)+…+A (k)+…称为矩阵级数, 记为 ()0 k k A ∞=∑,即有 ()0k k A ∞=∑= A (0)+A (1)+…+A (k)+… 定义3.5 记S (N )=()0N k k A =∑,称其为矩阵级数()0k k A ∞=∑的部分和. 如果矩阵序列{S (N)}收敛,且有极限S, 即有 S (N)→S