中科院矩阵分析课件
矩阵分析课件(2-1)

定义2.1.6. 称(2.1.3)式中右边的
对角形 - 矩阵为A ()的smith标准形, 称d1 ( ),d 2 ( ), ..., d r ( )为A ()的不变 因子。
上述定理说明,任何 - 矩阵( A )都 与它的smith标准形等价。
引入符号:
ri:表示矩阵的第i行;Ci:表示矩阵的第i列;
如果A1 ( )中至少有一个元素不为0,不妨设 左上角元素不为0,对A1 ( )重复A( )的讨论过 程, ; 最后将对角形中主对角线上元素变为 首系为1,因此,有 d1 ( )
. d r ( ) A( ) , 0 . 0 其中d( 的且d ( | di ( j )是首1 i ) 1 ) ( j 1, 2, , r; i 1, 2, , r - 1).
r (或 Ci) rj (或C j ):表示互换矩阵的第i , j i 两行(或两列);
Cri (或CCi ):表示矩阵的第i行(或列)乘常数C;
() ri rj:表示将矩阵的第i行乘上 ()后
加到第j行上;
()C i C j:表示将矩阵的第i列乘上 ()后
加到第j列上。
验证可知: 1 P ( i , j ) P ( i , j ),P ( i (c )) P ( i ( )), c -1 P ( i , j( () )) P ( i , j(- () )).
-1 -1
与线性代数中的证明类似,可以证明:
定理2.1.2: 对一个m n的 - 矩阵A ()
定理2.1.4 任意一个非零的n阶 - 矩阵A( ) 都等价于一个对角矩阵,即 d1 ( ) ... d r ( ) A( ) (2.1.3) 0 ... 0 其中r 1, d i ( )是首系为1的多项式且 d i ( ) | d i 1 ( ),(i 1, 2...r - 1)。
矩阵分析第4章ppt课件.ppt

从而
A
P 1
Ir 0
D 0
Q
1
P 1
Ir 0
I
r
D Q 1
BC
其中
B
P 1
Ir
0
Crmr ,
C Ir
D
Q
1
C rn r
A BIr D Q1 B D Q1 AQ B D
所以B是A中r 个线性无关的列
例 :分别求下面三个矩阵的满秩分解
1 2 1 0 1 2
(1)
0 0
0 1 1 (2) A 2 0 0
解: (1)由于
1 2
AAH 0 0
0 0
1 2
0 0
0 0
5 AAH 0
0
0 0 0
0 0 0
显然 AAH 的特征值为5,0,0,所以 A 的
奇异值为 5
(2)由于
0 2
AAH
0 2
1 0
1 0
1 1
0 0
AAH
2 0
0 4
显然 AAH 的特征值为 2,4,所以 A 的
2
0
cnn
2
n Unnn ,
c11 c21
cn1
R
c22
cn
2
cnn
显然矩阵 R 是一个正线上三角矩阵。
A是列满秩也有
注:Ar 1 2
r
c11 c21
cr1
1 2
r
c22
cr
2
crr
UrR 矩阵 R 是一个正线上三角矩阵
下面考虑分解的唯一性。设有两种分解式
A UR UR
1
2
1 2
1
1 2
矩阵分析第5章课件

第五章 向量与矩阵范数 前言
• 向量与矩阵范数是向量与矩阵的一个重要数 字特征---用它可以建立向量集或矩阵集的 拓扑结构,从而便于研究向量或矩阵序列,向 量或矩阵级数的收敛性质.因此,这一章的理 论在数值分析及其它领域中十分有用. • 本章是本课程重点内容之一.所有5节都要认 真学好.最后一节(矩阵幂级数)是研究矩阵 函数的重要工具.
Holder不等式与Minkowski不等式
• 下面两个不等式对本章的理论推导十分有用 • Holder不等式:对任意给定p>1和q=p/(p-1) (>1,即(1/p)+(1/q)=1)及任意ak,bk0成立 k=1nakbk (k=1nakp)1/p(k=1nbkp)1/p. (C-S不等式为其(p=2时)特例) • Minkowski不等式:对任意给定p1成立 (k=1n|ak+bk|p)1/p (k=1n|ak|p)1/p+(k=1n|bk|p)1/p
ACmn 定义 ‖A‖= maxi,k|aik| 则‖A‖显然是向量范数(向量的无穷大范数),但它 不是矩阵范数,反例如下:
1 1 1 1 1 2 A 1 1 , B 0 1 , AB 1 2
矩阵分析课件-第六章

cos A B=cosA cos B sin A sin B
dt
dt
d cos At=A sin At=-sin At A
dt
6 det eA=etrA,其中trA是A 的迹
7 cos A= 1 eiA+e-iA ,sinA= 1 eiA-e-iA
2
2i
8 sin2 A+cos2 A=E,sin -A=-sin A,cos -A =cosA
9当AB=BA时,有sin A B=sin A cos B cos A sin B
D
i
其中
D
J
i
=
D/ i Di
1
di-1
!Ddi-1
i
D/ i
Di
dixdi
设D = E-A =-1 p1 -2 p2 -s ps
i
j, i
j
, pi是i的代数重复度;
pi
d
,
i
D i =D i = =Ddi-1 i =0, D Ji =0,
故:D A=0.
f (k) j =p(k) j ,j=1,2, ,s;k=0,1, ,dj-1
即f x与p 在A的影谱上有相同的值, 则矩阵函数f A定义为:
f A=pA 称p 为f A的定义多项式。
定理6.2.1:设A
Cnn,J为A的若当标准形,P
Cnn n
且A=PJP-1,函数 f x 在A的影谱上有定义,
ln E+A的幂级数展开式见p201
&6.4 矩阵指数函数与矩阵三角函数
由定理5.5.3知:对任意n阶方阵A
e
At=
k=0
Aktk, k!
sin At= k=0
中科院学习课件 矩阵分析与应用 9 Determinants

Since σ (1, 2) = +1 and σ (2, 1) = −1, we obtain the familiar formula a11 a12 a21 a22 = a11 a22 bin | UCAS
7 / 23
Determinants | Determinants
Li Bao bin | UCAS 2 / 23
Determinants | Introduction
These men had something else in common — their ideas concerning the solution of linear systems were never adopted by the mathematical community of their time, and their discoveries quickly faded into oblivion. Eventually the determinant was rediscovered, and much was written on the subject between 1750 and 1900. During this era, determinants became the major tool used to analyze and solve linear systems, while the theory of matrices remained relatively undeveloped. The study and use of determinants eventually gave way to Cayley.s matrix algebra, and today matrix and linear algebra are in the main stream of applied mathematics, while the role of determinants has been relegated to a minor backwater position. Nevertheless, it is still important to understand what a determinant is and to learn a few of its fundamental properties. Our goal is not to study determinants for their own sake, but rather to explore those properties that are useful in the further development of matrix theory and its applications.
矩阵分析第4章课件

矩阵满秩分解不唯一;但同一矩阵的两个满
秩分解的因式矩阵之间存在密切的关系( 见P153,定理4.1.2).
ACrmn r=rank A min{m,n} A的秩等于它的行秩、列秩或行列式秩。A的行( 列)秩是它的最大线性无关组的行(列)数;A 的行列式秩是它的非0子式的最大阶数。 A=BC rank A rank B & rank A rank C
1
初等变换与初等矩阵性质
①3类初等矩阵都是可逆的(行列式不为0). ②将A依次作初等矩阵P1,…,Pr对应的行(列)初等变
换等价于左(右)乘A以可逆矩阵Pr,…,P1(P1,…,Pr).
③可适当选第一类初等矩阵的乘积P使PA(AP)的 行(列)是A的行(列)的任意排列.可适当选第三类 初等矩阵P(i,j(k))中的k使P(i,j(k))A的(i,j) 元变为0.可适当选第二类初等矩阵P(i(k))中的k 使P(i(k))A的非零(i,i)元变为1.综合起来推出: Er 0 存在初等矩阵的乘积P和Q,使 PAQ= 0 0 m n 其中r=rank A.一般地,ACr 都 Er 0 存在m,n阶可逆阵P和Q使 PAQ=
a11 a1n AB ann
b11 b1n a11b11 * bnn annbnn
a11 a1n 1/ a11 * 1 1 A , aii 0 det A 0 A det A a 1/ a nn nn
1 C11 1 2 C21 1 C22 2 n Cn1 1 Cn 2 2 ... Cnn n
矩阵分析第一章课件.ppt

性判别矩阵
C
V
CA
CAn
1
是列满秩的。
例 5:设
A
0 1
1 0
,
B
1 1
1 1
由于矩阵
B
AB
1 1
1 1 1 1 1 1
是行满秩的,所以相应的系统是可控制的。
二 矩阵理论在生物数学中的应用
在化的花瓣中存在一种特殊的生物模式。几乎所有 花,其花瓣数都是一种有规律的级数。例如百合花 的花瓣有3瓣;毛茛属的植物有5瓣花;许多翠雀属 的植物有8瓣花;万寿菊的花瓣有13瓣;紫菀属的植 物有21瓣花;大多数的雏菊有34,55,89 瓣花。 另外,在向日葵的花盘内葵花籽的螺旋式排列中也 可以发现类似的排列模式,同时植物的叶序中也存 在此种现象。这就是著名的Fibonacci级数模式。我 们称下面的数列
x2
4 3 , x3
1 3 , x4
2 3
同样可解出在第二组基下的坐标为
y1 1, y2 1, y3 1, y4 4
由此可以看出:一个向量在不同基底下的坐标是不相 同的。
基变换与坐标变换
设 1,2 , ,n(旧的)与 1, 2, , n (新的) 是 n 维线性空间V 的两组基底,它们之间的关系为
注意: 通过上面的例子可以看出线性空间的基底并不 唯一,但是维数是唯一确定的。利用维数的定义线性 空间可以分为有限维线性空间和无限维线性空间。目 前,我们主要讨论有限维的线性空间。
例 4 在4维线性空间 R22 中,向量组
0 1
1 1
,
1 1
0 1
,
1 0
1 1
,
1 1
中科院矩阵分析_第一章

矩阵的代数性质1•矩阵是线性映射的表示:线性映射的相加表示为矩阵的相加线性映射的复合表示为矩阵的相乘2•矩阵是一种语言,它是表示复杂系统的有力工具。
学习矩阵理论的重要用途之一就是学会用矩阵表示复杂系统的关系,培养根据矩阵推演公式的能力是学习矩阵论的目的之一。
定义一个矩阵有几种方式:可以通过定义矩阵的每一个元素来定义一个矩阵,也可以通过矩阵具有的性质来定义一个矩阵。
如:对称矩阵可以定义为:a ij=a ji也可以定义为:(x, Ay)=(Ax,y),还可以定义为:Ax= f(x),其中f(x)=x T Ax/2,即它对向量x的作用相当于函数f(x)在x处的梯度。
3. 矩阵可以表示为图像矩阵的大小可以表示为图像。
反之,一幅灰度图像本身就是矩阵。
图像压缩就是矩阵的表示问题•这时矩阵相邻元素间有局部连续性,既相邻的元素的值大都差别不大。
4. 矩阵是二维的(几何性质)矩阵能够在二维的纸张和屏幕等平面媒体上表示,使得用矩阵表示的问题显得简单清楚,直观,易于理解和交流。
很多二元关系很直观的就表示为矩阵,如关系数据库中的属性和属性值,随机马尔科夫链的状态转移概率矩阵,图论中的有向图或无向图的矩阵表示等。
第一章:线性空间和线性变换1. 线性空间集合与映射集合是现代数学最重要的概念,但没有严格的定义。
集合与其说是一个数学概念,还不如说是一种思维方式,即用集合(整体)的观点思考问题。
整个数学发展的历史就是从特殊到一般,从个体到整体的发展历程。
集合的运算及规则,两个集合的并、交运算以及一个集合的补;集合中元素没有重合,子集,元素设S, S'为集合映射:为一个规则:S S',使得S中元素a和S'中元素对应,记为a'= (a),或:a a'.映射最本质的特征在于对于S中的任意一个元素在S'中仅有唯一的一个元素和它对应。
映射的原象,象;映射的复合。
满射,单射,—映射。
若S'和S相同,则称为变换。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
矩阵分析及其应用3.1 矩阵序列定义3.1 设矩阵序列{A (k)},其中A (k)=()(k ij a )∈C m ⨯n ,当k →∞,)(k ija →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证明收敛。
相反,由于∑∑∑+=+=+=-≤≤nm k n m k n m k k k k k k 11212)1(11)sin( < 1/m从而只要l 充分大,则当m, n > l 时就有ε≤∑+=nm k k k 12)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|<M 或等价的|| A(k)||<M’|)(kij定理:有界的矩阵序列{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 <13.2矩阵级数定义3.4设矩阵序列{A (k)},其中A (k)=()(k ij a )∈C n ⨯n ,由它们 形成的无穷和 A (0)+A (1)+…+A (k)+…称为矩阵级数, 记为 ()0k 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那么称矩阵级数()0k k A∞=∑收敛,且和为S, 记为S=()0k k A∞=∑不收敛的矩阵级数称为发散的。
显然()0k k A∞=∑=S 是指()0k ij ij k a s ∞==∑,∀ i ,j即矩阵级数收敛是指它的每个分量所构成的数项级数收敛。
性质:矩阵级数()0k k A∞=∑收敛的充要条件为对任意向量x, 向量级数()0k k Ax ∞=∑收敛。
定义3.6 设矩阵级数()0k k A∞=∑的每个分量)(k ij a 所构成的数项级数()0k ij k a∞=∑绝对收敛,则称矩阵级数()0k k A ∞=∑绝对收敛。
关于绝对收敛,我们有如下的定理:性质1. 绝对收敛的()0k k A∞=∑交换求和次序不改变其绝对收敛性和极限值。
性质2. 矩阵级数()0k k A∞=∑绝对收敛的充要条件为正项级数()0||||k k A∞=∑收敛。
性质3. 如果矩阵级数()0k k A ∞=∑(绝对)收敛,那么()0k k PA Q ∞=∑也是(绝对)收敛,且有 ()0k k PAQ ∞=∑=P (()0k k A ∞=∑)Q 性质4. 设C n ⨯n 的两个矩阵级数S 1: A (1)+A (2)+…+A (k)+…S 2: B (1)+B (2)+…+B (k)+…都绝对收敛,其和分别为A 和B.则矩阵级数S 3: A (1) B (1)+ [A (1) B (2)+ A (2) B (1)]+…+[ A (1) B (k)+ A (2) B (k -1) +…+A (k) B (1)]+…绝对收敛且和为AB.证明:由于S 1: A (1)+A (2)+…+A (k)+…绝对收敛的充要条件为 正项级数||A (1)||+||A (2)||+…+||A (k)||+…收敛且与排列无关。
我们证明的思路是证明正项级数:||A (1) B (1)||+ ||A (1) B (2)+ A (2) B (1)||+…+||A (1) B (k)+ A (2) B (k -1) +…+A (k) B (1)||+…收敛。
引用魏氏定理,我们仅需验证下列正项级数: ||A (1)||⋅||B (1)||+ { ||A (1) ||⋅||B (2)||+ ||A (2) ||⋅||B (1)||}+…+{||A (1) ||⋅||B (k)||+ ||A (2) ||⋅||B (k -1) ||+…+||A (k) ||⋅||B (1)||}+… 收敛。
这由题设正项级数||A (1)||+||A (2)||+…+||A (k)||+…和正项级数 ||B (1)||+||B (2)||+…+||B (k)||+…的收敛性可得。
定理3.4 幂级数 I+A+A 2+…+A k +…收敛的充要条件为 A 的谱半径ρ(A)<1, 收敛时其和为(I -A)-1。
若有矩阵范数||.||使得||A||<1,则||(I -A)-1- (I+A+A 2+…+A k )||≤||A||k+1/(1-||A||)证明: 必要性. 由于I+A+A 2+…+A k +…收敛,从而 S (k)= I+A+A 2+…+A k 收敛。
记T (k)= I+A+A 2+…+A k+1, A k+1=T (k)- S (k)收敛,且T (k)- S (k) →0,这样我们有A k →0,从而ρ(A)<1.充分性:设ρ(A)<1,(I -A)-1存在,由于I+A+A 2+…+A k =(I -A)-1 -(I -A)-1 A k+1因A k →0,所以I+A+A 2+…+A k +…→ (I -A)-1.又因为(I -A)-1 - (I+A+A 2+…+A k )= (I -A)-1 A k+1从而||(I -A)-1 - (I+A+A 2+…+A k )||=|| (I -A)-1 A k+1||设B=(I -A)-1A k+1,从而(I -A)B=A k+1即B=AB+ A k+1,从而||B||≤ ||A||⋅||B||+ ||A k+1||≤ ||A||⋅||B||+ ||A||k+1因为矩阵范数||.||使得||A||<1,所以||B||≤||A||k+1/(1-||A||)成立。
定理3.6 设幂级数 ∑∞==0)(k k k z cz f 的收敛半径为r ,如果方阵A 满足ρ(A)< r , 则矩阵幂级数∑∞==0)(k k k A c A f 是绝对收敛的;如果ρ(A) > r ,∑∞=0k k kA c 是发散的。