矩阵论——讲稿

合集下载

矩阵讲义全

矩阵讲义全

本课程的说明:矩阵分析理论是在线性代数的基础上推广的(数学是在已有的基础理论上模仿,推广而发展的。

要大胆猜想,小心证明!) 矩阵分析理论的组成:四部分:一、基础知识(包括书上的前三章内容)重点、难点:约当标准形与多项式矩阵,矩阵的分解等; 二、矩阵分析(第四章:矩阵函数及其应用)重点、难点:范数,矩阵幂级数,微分方程组; 三、矩阵特征值的估计(第五章)重点、难点:Gerschgorin 圆盘定理;广义逆矩阵; 四、非负矩阵(第六章)(注:不讲)重点、难点:基本不等式,素矩阵,随机矩阵等。

§1 线性空间与度量空间一、线性空间: 1.数域:Df 1:若复数的一个非空集合P 含有非零的数,且其中任意两数的和、差、积、商(除数不为0)仍在这个集合中,则称数集P 为一个数域 eg 1:Q (有理数),R (实数),C (复数),Z (整数),N (自然数)中哪些是数域?哪些不是数域? 2.线性空间— 设P 是一个数域,V 是一个非空集合,若满足:<1> 可加性—指在V 上定义了一个二元运算(加法)即:V ∈∀βα, 经过该运算总存在唯一的元素V ∈γ与之对应,称γ为α与β的和,记βαγ+= 并满足:① αββα+=+② )()(γβαγβα++=++ ③ 零元素—=有θαθααθ+∈∀∈∃Vt s V .(线性空间必含θ)。

④ αβαβθβααβ-+∈∀∈∃=记的负元素为=有对V V<2> 数积:(数乘运算)—在P 与V 之间定义了另一种运算。

即V P k ∈∈∀α,经该运算后所得结果,仍为V 中一个唯一确定的元素(存在唯一确定的元素V ∈δ与之对应),称δ为k 与α的乘积。

记为αδk =并满足:① αα=⋅1② P l k ∈∀, αα)()(kl l k = ③ P l k ∈∀, αααl k l k +=+)( ④ γβα∈∀, βαβαk k k +=+)(则称V 为数域P 上的线性空间(向量空间)记为)...(∙+P V 习惯上V 中的元素—向量, θ—零向量, 负元素—负向量结论:可以证明,线性空间中的零向量是唯一的,负元素也是唯一的,且有:θα=⋅0 θθ=⋅k αα-=⋅-)1( )(βαβα-+=-eg2:}{阶矩阵是n m A A V ⨯= P —实数域R按照矩阵的加法和数与矩阵的乘法,就构成实数域R 上的线性空间,记为:n m R ⨯同样,若V 为n 维向量,则可构成R 上的n 维向量空间n R —线性空间。

矩阵论课件

矩阵论课件
矩阵论
第二章
第一节
矩阵与约当标准形
矩阵
第二节 不变因子及初等因子
第三节 约当标准形 第四节 凯莱—哈米尔顿定理 最小多项式
4 December 2014 河北科技大学
机动 目录 上页 下页 返回 结束
矩阵论
第一节
定义 设 P
矩阵
为数域, 为数字,P[ ] [ ]为关于 中的元素(数)为元素的矩
4 December 2014
河北科技大学
机动 目录 上页 下页 返回 结束
定理 设 矩阵 A( ) aij
阵,且 rank( A( )) r ,则


矩阵论
m n
为非零的多项式矩
A( )
d1 ( ) d 2 ( ) r ( ) 0 J ( ) 0 0 d ( ) r 0 0 0 diag d1 ( ), d 2 ( ), , d r ( ), 0, , 0 --称为 A( )的 Smith (史密斯)标准形.
矩阵论
Dn ( ) a ;
n
Dn1 ( )
4 December 2014
D1 ( ) 1.
河北科技大学
机动 目录 上页 下页 返回 结束
矩阵论
定义 把 矩阵 A( ) 的每个次数大于零的不变因子
在复数域 [ ]中分解成标准分解式,即分解成首项 系数为1的互不相同的一次因式方幂的乘积,所有 这些一次因式 的方幂 ( 相同的必须按出现次数 计 算) ,称为 A( )的初等因子.
[ ]中分解成标准分解式,所有出现的一次因式的
标准形)
方幂就是 A( )的全部初等因子.

研究生矩阵论第1讲 线性空间

研究生矩阵论第1讲 线性空间

矩阵论1、意义随着科学技术的发展,古典的线性代数知识己不能满足现代科技的需要,矩阵的理论和方法业巳成为现代科技领域必不可少的工具.有人认为:“科学计算实质就是矩阵的计算”.这句话概括了矩阵理论和方法的重要性及其使用的广泛性.因此,学习和掌握矩阵的基本理论和方法,对于理、工科研究生来说是必不可少的数学工具.2、内容《矩阵论》和工科《线性代数》课程在研究矩阵的内容上有较大的差异:线性代数:研究行列式、矩阵的四则运算(加、减、乘、求逆 ) 以及第一类初等变换 (非正交的)、对角标准形 (含二次型) 以及n阶线性方程组的解等基本内容.矩阵论:研究矩阵的几何理论(线性空间、线性算子、内积空间等)、第二和第三类初等变换(正交的)、分析运算(矩阵微积分和级数)、矩阵的范数和条件数、广义逆和分解、若尔当标准形以及几类特殊矩阵和特殊运算等,内容十分丰富.3、方法在研究的方法上,矩阵论和线性代数也有很大的不同:线性代数:引入概念直观,着重计算.矩阵论:着重从几何理论的角度引入矩阵的许多概念和运算,把矩阵看成是线性空间上线性算子的一种数量表示.深刻理解它们对将来正确处理实际问题有很大的作用.第1讲 线性空间内容: 1.线性空间的概念;2.基变换和坐标变换;3.子空间和维数定理;4.线性空间的同构线性空间和线性变换是矩阵分析中经常用到的两个极其重要的概念,也是通常几何空间概念的推广和抽象,线性空间是某类客观事物从量的方面的一个抽象.§1 线性空间的概念1. 群,环,域代数学是用符号代替数(或其它)来研究数(或其它)的运算性质和规律的学科,简称代数.代数运算:假定对于集A 中的任意元素a 和集B 中的任意元素b ,按某一法则和集C 中唯一确定的元素c 对应,则称这个对应为A 、B 的一个(二元)代数运算.代数系统:指一个集A 满足某些代数运算的系统.1.1群定义1.1 设V 是一个非空集合,在集合V 的元素之间定义了一种代数运算,叫做加法,记为“+”.即,对V 中给定的一个法则,对于V 中任意元素βα,,在V 中都有惟一的一个元ν和他们对应,称ν为βα,的和,记为βαν+=.若在“+”下,满足下列四个条件,则称V 为一个群.1)V 在“+”下是封闭的.即,若,,V ∈βα有 V ∈+βα;2) V 在“+”下是可结合的.即,)()(γβαγβα++=++ ,V ∈γ;3)在V 中有一个元e ,若,V ∈β有 βββ=+=+e e ;e 称为单位元;4)对于,V ∈β有 e =+=+αββα.称α为β的逆元.注:对V 任意元素βα,,都有αββα+=+,则称V 为交换群或阿贝尔群.1.2 环定义1.2 设V 是一个非空集合,在集合V 的元素之间定义了两种代数运算,分别叫做加法、乘法,记为“+”和“*”.即,对V 中给定的一个法则,对于V 中任意元素α,β,在V 中都有惟一的一个元ν和他们对应,称ν为α,β的和和积,记为βαν+=(βαν*=).满足下列三个条件,则称V 为一个环. 1)V 在“+”下是阿贝尔群;2) V 在“*”下是可结合的.即,)()(νβανβα**=**;3)乘法对加法满足左、右分配律,即对于V 中任意元素α,β,ν,有 βνανβαν**)(*+=+,νβνανβα*+*=*+)(.注:对V 任意元素βα,,都有αββα*=*,则称V 为交换环.1.3 域定义 1.3 设V 满足环的条件,且在对“加法”群中去除单位元的集合对于“乘法”满足交换群的条件,则称V 为域.例:有理数集对于通常的数的加法和乘法运算构成域,称之为有理数域.最常见的数域有有理数域Q 、实数域R 、复数域C .实数域和复数域是工程上较常用的两个数域.此外,还有其它很多数域.如{}.,2)2(Q b a b a Q ∈+=,不难验证,)2(Q 对实数四则运算封闭的,所以)2(Q 也是一个数域.而整数集合Z 就不是数域. 数域有一个简单性质,即所有的数域都包含有理数域作为它的一部分.特别,每个数域都包含整数0和1. 2. 线性空间定义 1.4 设V 是一个非空集合,P 是一个数域.在集合V 的元素之间定义了一种代数运算,叫做加法,记为“+”:即,给出了一个法则对于V 中任意元素βα,,在V 中都有惟一的一个元ν和他们对应,称ν为βα,的和,记为βαν+=.在数域P 和集合V 的元素之间还定义了一种代数运算,称为数量乘法(数乘),记为“∙”:即,对于数域P 中任一数k 和V 中任一元α,在V 中都有惟一的一个元δ和它们对应,称δ为k 和α的数乘,记为αδ∙=k .如果加法和数乘这两种运算在V 中是封闭的,且满足如下八条规则:⑴ 交换律αββα+=+;⑵ 结合律)()(γβαγβα++=++ ,V ∈γ;⑶ V V ∈∃∈∀0,α,有αα=+0,(0称为零元素);⑷ V V ∈∃∈∀βα,,有 0=+βα,(β称为的α负元素,记为α-); ⑸ P V ∈∈∀1,α,有 αα=∙1;⑹ αα∙=∙∙)()(kl l k ,P l k ∈,;⑺ ααα∙+∙=∙+l k l k )(;⑻ βαβα∙+∙=+∙k k k )(,则称集合V 为数域P 上的线性空间.当数域P 为实数域时,V 就称为实线性空间;P 为复数域,V 就称为复线性空间.例 1.按通常向量的加法和数乘运算,由全体实n 维向量组成的集合,在实数域R 上构成一个实线性空间,记为n R ;由全体复n 维向量组成的集合,在复数域C 上构成—个复线性空间,记为n C .例 2.按照矩阵的加法及数和矩阵的乘法,由数域P 上的元素构成的全体n m ⨯矩阵所成的集合,在数域P 上构成一个线性空间,记为n m P ⨯.而其中秩为)0(>r r 的全体矩阵所成的集合rR 则不构成线性空间,为什么?(事实上,零矩阵r R O ∉).例3.按通常意义的函数加法和数乘函数,闭区间[]b a ,上的连续函数的全体所成的集合,构成线性空间[]b a C ,.例4. 设+R ={全体正实数},其“加法”及“数乘”运算定义为xy y x =+, k x x k = 。

矩阵论矩阵分析

矩阵论矩阵分析

第三章 矩阵分析在此之前我们只研究了矩阵的代数运算,但在数学的许多分支和工程实际中,特别是涉及到多元分析时,还要用到矩阵的分析运算.本章首先讨论矩阵序列的极限和矩阵级数,然后介绍矩阵函数和它的计算,最后介绍矩阵的微积分,以及矩阵分析在解微分方程组和线性矩阵方程中的应用.§3.1 矩阵序列 定义 3.1 设有Cm n⨯中的矩阵序列{}()k A ,其中()()()k k ij m nAa ⨯=.若()lim (1,2,,;1,2,,)k ij ij k a a i m j n →+∞=== ,则称矩阵序列{}()k A 收敛于()ij m n A a ⨯=,或称A 为矩阵序列{}()k A 的极限,记为()lim k k A A →+∞=或()()k A A k →→+∞不收敛的矩阵序列称为发散. 由定义可见,Cm n⨯中一个矩阵序列的收敛相当于mn 个数列同时收敛.因此,可以用初等分析的方法来研究它.但同时研究mn 个数列的极限未免繁琐.与向量序列一样,可以利用矩阵范数来研究矩阵序列的极限. 定理 3.1 设()k A,C (012)m n A k ,,,⨯∈= .则()lim k k AA →+∞=的充分必要条件是()lim 0k k A A →+∞-=,其中 是C m n ⨯上的任一矩阵范数.证 先取Cm n⨯上矩阵的G-范数.由于()()()()1=1k k k ij ij ij ij Gi,jm nk ijiji j a a a a A Aaa =-≤-=-≤-所以()lim k k A A →+∞=的充分必要条件是()lim 0k Gk A A→+∞-=.又由范数的等价性知,对C m n⨯上任一矩阵范数 ,存在正常数α,β,使得()()()k k k GGAAAA AA αβ-≤-≤-故()lim 0k Gk AA→+∞-=的充分必要条件是()lim 0k k A A →+∞-=.证毕推论 设()k A,C(012)m nA k ,,,⨯∈= ,()lim k k A A →+∞=.则()lim k k A A →+∞=其中 是Cm n⨯上任一矩阵范数.证 由()()k k AA A A -≤-即知结论成立.证毕需要指出的是,上述推论的相反结果不成立.如矩阵序列()1(1)112k k A k ⎛⎫- ⎪=+ ⎪ ⎪⎝⎭不收敛.但()Flim lim k k x A →+∞== 收敛的矩阵序列的性质,有许多与收敛数列的性质相类似. 定理3.2 设()lim k k AA →+∞=,()lim k k B B →+∞=,其中()k A ,()k B ,A ,B 为适当阶的矩阵,α,β∈C .则 (1)()()lim ()k k k AB A B αβαβ→+∞+=+;(2) ()()lim k k k A BAB →+∞=;(3)当()k A与A 均可逆时,()11lim ()k k AA --→+∞=.证 取矩阵范数 ,有()()()()()()()()()()()()()()()k k k k k k k k k k k k k A B A B A A B B A B ABA B A B A B AB A B B A A Bαβαβαβ+-+≤-+--=-+-≤-+-由定理3.1和推论知(1)和(2)成立.因为()1()k A -,1A -存在,所以()lim det det 0k k AA →+∞=≠,又有()lim adj adj k k A A →+∞=.于是()()11()adj adj lim ()lim det det k k k k k A AA A A A--→+∞→+∞=== 证毕 定理3.2(3)中条件()k A与A 都可逆是不可少的,因为即使所有的()k A可逆也不能保证A一定可逆.例如()11111k Ak ⎛⎫+ ⎪= ⎪ ⎪⎝⎭对每一个()k A都有逆矩阵()1()1k kk A k k --⎛⎫=⎪-+⎝⎭,但()11lim 11k k A A →+∞⎛⎫== ⎪⎝⎭而A 是不可逆的. 在矩阵序列中,最常见的是由一个方阵的幂构成的序列.关于这样的矩阵序列有以下的概念和收敛定理. 定义3.2 设n nA C ⨯∈,若()lim 0k k A→+∞=,则称A 为收敛矩阵.定理3.3 设n nA C⨯∈,则A 为收敛矩阵的充分必要条件是ρ(A )<1.证 必要性.已知A 为收敛矩阵,则由谱半径的性质,有(())()k k k A A A ρρ=≤其中 是Cn n⨯上任一矩阵范数,即有lim (())0kk A ρ→+∞=,故ρ(A )<1.充分性.由于ρ(A )<1,则存在正数ε,使得ρ(A )+ε<1.根据定理2.14,存在C n n⨯上的矩阵范数m ,使得()1m A A ρε≤+<从而由kk m mAA ≤得lim 0kmk A →+∞=.故lim 0k k A →+∞=. 证毕推论 设n nA C ⨯∈.若对Cn n⨯上的某一矩阵范数 有1A <,则A 为收敛矩阵.例3.1 判断下列矩阵是否为收敛矩阵:(1)181216A -⎛⎫= ⎪-⎝⎭; (2)0.20.10.20.50.50.40.10.30.2A ⎛⎫⎪= ⎪ ⎪⎝⎭. 解 (1)可求得A 的特征值为156λ=,212λ=-,于是5()16A ρ=<,故A 是收敛矩阵; (2)因为10.91A =<,所以A 是收敛矩阵.§3.2 矩阵级数定义3.3 由Cm n⨯中的矩阵序列{}()k A 构成的无穷和(0)(1)()k A A A ++++ 称为矩阵级数,记为()k k A+∞=∑.对任一正整数N ,称()()NN k k SA ==∑为矩阵级数的部分和.如果由部分和构成的矩阵序列{}()N S收敛,且有极限S ,即()lim N N SS →+∞=,则称矩阵级数()0k k A +∞=∑收敛,而且有和S ,记为()k k S A+∞==∑不收敛的矩阵级数称为发散的.如果记()()()k k ij m n Aa ⨯=,()ij m n S s ⨯=,显然()0k k S A +∞==∑相当于()(1,2,,;1,2,,)k ij ij k a s i m j n +∞====∑即mn 个数项级数都收敛. 例3.2 已知()1π24(0,1,)10(1)(2)k kk A k k k ⎛⎫⎪ ⎪== ⎪ ⎪++⎝⎭研究矩阵级数()k A+∞∑的敛散性.解 因为k 00()()001π2410(1)(2)1π1242341012N Nk Nk k N k N k k N N S A k k N ====⎛⎫⎪ ⎪== ⎪ ⎪++⎝⎭⎛⎫⎛⎫-- ⎪ ⎪⎝⎭ ⎪= ⎪- ⎪+⎝⎭∑∑∑∑所以()4π2lim 301N N S S →+∞⎛⎫⎪= ⎪ ⎪⎝⎭故所给矩阵级数收敛,且其和为S . 定义3.4 设()()()C (0,1,)k k m n ij m n Aa k ⨯⨯=∈= .如果mn 个数项级数()0(1,2,,;1,2,,)k ijk ai m j n +∞===∑ 都绝对收敛,即()k ijk a +∞=∑都收敛,则称矩阵级数()k k A+∞=∑绝对收敛.利用矩阵范数,可以将判定矩阵级数是否绝对收敛转化为判定一个正项级数是否收敛的问题.定理3.4 设()()()C(0,1,)k k m nij m nAa k ⨯⨯=∈= .则矩阵级数()0k k A +∞=∑绝对收敛的充分必要条件是正项级数()0k k A +∞=∑收敛,其中 是C m n ⨯上任一矩阵范数.证 先取矩阵的1m -范数.若1()k k m A +∞=∑收敛,由于1()()()11(1,2,,;1,2,,)mnk k k ijij i j m aa A i m j n ==≤===∑∑从而由正项级数的比较判别法知()k ijk a+∞=∑都收敛,故()k k A+∞=∑绝对收敛.反之,若()k k A+∞=∑绝对收敛,则()0k ijk a+∞=∑都收敛,从而其部分和有界,即()0(1,2,,;1,2,,)Nk ijijk aM i m j n =≤==∑ 记,max ij i jM M =,则有1()()()0011110()()NNmnm n Nk k k ijijk k i j i j k m AaamnM =========≤∑∑∑∑∑∑∑故1()k k m A +∞=∑收敛.这表明()k k A+∞=∑绝对收敛的充分必要条件是1()k k m A +∞=∑收敛.由矩阵范数的等价性和正项级数的比较判别法知,1()k k m A+∞=∑收敛的充分必要条件是()0k k A +∞=∑收敛,其中 是C m n ⨯上任一矩阵范数. 证毕利用矩阵级数收敛和绝对收敛的定义,以及数学分析中的相应结果,可以得到以下一些结论.定理3.5 设()k k AA +∞==∑,()0k k B B +∞==∑,其中()k A ,()k B ,A ,B 是适当阶的矩阵,则(1)()()0()k k k AB A B +∞=+=+∑;(2)对任意λ∈C ,有()k k AA λλ+∞==∑;(3)绝对收敛的矩阵级数必收敛,并且任意调换其项的顺序所得的矩阵级数仍收敛,且其和不变; (4)若矩阵级数()k k A+∞=∑收敛(或绝对收敛),则矩阵级数()k k PAQ +∞=∑也收敛(或绝对收敛),并且有()()0()(3.1)k k k k PAQ P A Q+∞+∞===∑∑(5)若()k k A+∞=∑与()k k B+∞=∑均绝对收敛,则它们按项相乘所得的矩阵级数(0)(0)(0)(1)(1)(0)(0)()(1)(1)(()()(3.2)k k k A B AB A B A B A B A B -++++++++ 也绝对收敛,且其和为AB . 证 只证(4)和(5).若()0k k A+∞=∑收敛,记()()0NN k k SA ==∑,则()lim N N S A →+∞=.从而()()00lim(lim)NNk k N N k k PAQ P AQ PAQ →+∞→+∞====∑∑可见()k k PAQ +∞=∑收敛,且式(3.1)成立.若()k k A+∞=∑绝对收敛,则由定理3.4知()k k A +∞=∑收敛,但()()()k k k PA Q P AQ Aα≤≤其中α是与k 是无关的正数,从而()k k PAQ +∞=∑收敛,即()k k PAQ +∞=∑绝对收敛.当()k k A+∞=∑和()k k B+∞=∑绝对收敛时,由定理3.4知()k k A+∞=∑和()0k k B +∞=∑收敛,设其和分别为1σ与2σ,从而它们按项相乘所得的正项级数(0)(0)(0)(1)(1)(0)(0)()(1)(1)()(0)()()k k k A B A B A B ABABAB-++++++++也收敛,其和为12σσ.因为(0)()(1)(1)()(0)(0)()(1)(1)()(0)k k k k k k A B A B A B ABABAB--+++≤+++所以矩阵级数(3.2)绝对收敛.记()()1NN k k SA==∑,()()2NN k k SB ==∑,()(0)()(1)(1)()(0)3()NN k k k k SA B A B A B -==+++∑则()()()(1)()(2)(1)(2)()()(1)()()123N N N N N N N N N S S S A B A B A B A B A B --=++++++又记()()1NN k k Aσ==∑,()()2NN k k B σ==∑,()(0)()(1)(1)()(0)3()NN k k k k A B A B A B σ-==+++∑显然()()()()()()123123N N N N N N S S S σσσ-≤-故由()()12lim N N N S S AB →+∞=和()()()123lim ()0N N N N σσσ→+∞-=,得()3lim N N S AB →+∞=证毕下面讨论一类特殊的矩阵级数——矩阵幂级数. 定义3.5 设n nA C⨯∈,C(0,1,)k a k ∈= .称矩阵级数kk k a A+∞=∑为矩阵A 的幂级数.利用定义来判定矩阵幂级数的敛散性,需要判别2n 个数项级数的敛散性,当矩阵阶数n 较大时,这是很不方便的,且在许多情况下也无此必要.显然,矩阵幂级数是复变量z 的幂级数0kk k a z+∞=∑的推广.如果幂级数kk k a z+∞=∑的收敛半径为r ,则对收敛圆z r <内的所有z ,kk k a z+∞=∑都是绝对收敛的.因此,讨论kk k a A+∞=∑的收敛性问题自然联系到kk k a z+∞=∑的收敛半径.定理3.6 设幂级数kk k a z+∞=∑的收敛半径为r ,Cn nA ⨯∈.则(1)当ρ(A )<r 时,矩阵幂级数0kk k a A+∞=∑绝对收敛;(2)当ρ(A )>r 时,矩阵幂级数kk k a A+∞=∑发散.证 (1)因为ρ(A )<r ,所以存在正数ε,使得ρ(A )+ε<r .根据定理2.14,存在Cn n⨯上的矩阵范数m ,使得m ()A A r ρε≤+<从m m(())kk k k k k a A a A a A ρε≤≤+而由于幂级数(())kkk aA ρε+∞=+∑收敛,故矩阵幂级数0k k k a A +∞=∑绝对收敛.(2)当ρ(A )>r 时,设A 的,n 个特征值为12,,,n λλλ ,则有某个l λ满足l r λ>.由Jordan 定理,存在n 阶可逆矩阵P ,使得11112(10)i n n P AP J λδδδλλ--⎛⎫⎪⎪== ⎪⎪⎝⎭代表或而kk k a J+∞=∑的对角线元素为0(1,2,,)k k jk a j n λ+∞==∑ .由于0k k lk a λ+∞=∑发散,从而0k k k a J +∞=∑发散.故由定理 3.5(4)知,kkk a A+∞=∑也发散. 证毕推论 设幂级数kkk a z +∞=∑的收敛半径为r ,C n n A ⨯∈.若存在C n n ⨯上的某一矩阵范数 使得A r <,则矩阵幂级数0kk k a A+∞=∑绝对收敛.例3.3 判断矩阵幂级数018216kkk k+∞=-⎛⎫ ⎪-⎝⎭∑的敛散性. 解 令181216A -⎛⎫= ⎪-⎝⎭.例3.1中已求得5()6A ρ=.由于幂级数0kk kz +∞=∑的收敛半径为r =1,故由ρ(A )<1知矩阵幂级数kk kA+∞=∑绝对收敛.最后,考虑一个特殊的矩阵幂级数. 定理3.7 设Cn nA ⨯∈.矩阵幂级数kk A+∞=∑(称为Neumann 级数)收敛的充分必要条件是ρ(A )<1,并且在收敛时,其和为1()I A --. 证 当ρ(A )<1时,由于幂级数kk z+∞=∑的收敛半径r =1,故由定理 3.6知矩阵幂级数0kk A+∞=∑收敛.反之,若kk A+∞=∑收敛,记0kk S A+∞==∑,()()0NN k k SA ==∑则()lim N N S S →+∞=.由于()(1)()(1)lim lim ()=lim lim N N N N N N N N N A S S S S O --→+∞→+∞→+∞→+∞==--故由定理3.3知ρ(A )<1.当kk A+∞=∑收敛时,ρ(A )<1,因此I -A 可逆,又因为()1()N N S I A I A +-=-所以()111()()N N S I A A I A -+-=---故()1lim ()N N S S I A -→+∞==- 证毕 例3.4 已知0.20.10.20.50.50.40.10.30.2A ⎛⎫ ⎪= ⎪ ⎪⎝⎭,判断矩阵幂级数0k k A +∞=∑的敛散性.若收敛,试求其和.解 因为10.91A =<,所以kk A+∞=∑收敛,且102814141()44624214202535k k A I A +∞-=⎛⎫ ⎪=-= ⎪ ⎪⎝⎭∑ §3.3 矩阵函数矩阵函数是以矩阵为变量且取值为矩阵的一类函数.本节介绍矩阵函数的定义和计算方法,并讨论常用矩阵函数的性质. 一、矩阵函数的定义 定义3.5 设幂级数0k k k a z +∞=∑的收敛半径为r ,且当z r <时,幂级数收敛于函数f (z ),即0()()kk k f z a zz r +∞==<∑如果Cn nA ⨯∈满足ρ(A )<r ,则称收敛的矩阵幂级数kk k a A+∞=∑的和为矩阵函数,记为f (A ),即0()(3.3)kk k f A a A+∞==∑根据这个定义,可以得到在形式上和数学分析中的一些函数类似的矩阵函数.例如,对于如下函数的幂级数展开式02120101e ()!(1)sin ()(21)!(1)cos ()(2)!(1)(1)(1)ln(1)(1)1kzk k k k k kk kk k k k z r k z zr k z zr k z z r z zr k +∞=+∞+=+∞=+∞-=+∞+===+∞-==+∞+-==+∞-==-+==+∑∑∑∑∑ 相应地有矩阵函数01e !kk A A k +∞==∑(C n n A ⨯∈) 210(1)sin (21)!kk k A A k +∞+=-=+∑ (C n n A ⨯∈)20(1)cos (2)!k kk A A k +∞=-=∑ (C n n A ⨯∈)1()k k I A A +∞-=-=∑ (ρ(A )<1)1(1)ln()1k k k I A A k +∞+=-+=+∑ (ρ(A )<1)称e A为矩阵指数函数,sin A 为矩阵正弦函数,cos A 为矩阵余弦函数.如果把矩阵函数f (A )的变元A 换成At ,其中t 为参数,则相应得到()()(3.4)kk k f At a At +∞==∑在实际应用中,经常需要求含参数的矩阵函数.二、矩阵函数值的计算以上利用收敛的矩阵幂级数的和定义了矩阵函数f (A ),在具体应用中,要求将f (A )所代表的具体的矩阵求出来,即求出矩阵函数的值.这里介绍几种求矩阵函数值的方法.以下均假设式(3.3)或式(3.4)中的矩阵幂级数收敛. 方法一 利用Hamilton-Cayley 定理利用Hamilton-Cayley 定理找出矩阵方幂之间的关系,然后化简矩阵幂级数求出矩阵函数的值.举例说明如下. 例3.5 已知0110A ⎛⎫=⎪-⎝⎭,求e At.解 可求得2det()1I A λλ-=+.由Hamilton-Cayley 定理知2A I O +=,从而2A I =-,3A A =-,4A I =,5A A =,…即2(1)k k A I =-,21(1)(1,2,)k k A Ak +=-=故243501e 1!2!4!3!5!cos sin (cos )(sin )sin cos Atk k k t t t t A t I t Ak t t t I t A t t +∞=⎛⎫⎛⎫==-+-+-+- ⎪ ⎪⎝⎭⎝⎭⎛⎫=+= ⎪-⎝⎭∑例3.6 已知4阶方阵A 的特征值为π,-π,0,0,求sin A ,cos A .解 因为2422det()(π)(π)πI A λλλλλλ-=-+=-,所以422πA A O -=.于是422πA A =,523πA A =,642πA A =,743πA A =,…即2222πkk A A -=,21223π(2,3,)k k A A k +-==故213223023321323332(1)1(1)sin π(21)!3!(21)!11(1)π3!π(21)!sin ππ1ππk k k k k k k k k A A A A Ak k A A A k A A A A +∞+∞+-==+∞+=--==-+++⎛⎫-=-+ ⎪+⎝⎭=+=-∑∑∑-22222022222(1)1(1)cos π(2)!2!(2)!cos π12ππk k k k k k A A I A Ak k I A I A +∞+∞-==--==-+=+=-∑∑-方法二 利用相似对角化 设C n nA ⨯∈是可对角化的,即存在C n n n P ⨯∈,使得112diag(,,,)n P AP A λλλ-== 则有11112112()()()diag(,,,)diag((),(),,())kkk k k k k k k kk k k k k n k k k n f A a A a P P P a P P a a a P P f f f P λλλλλλ+∞+∞+∞--===+∞+∞+∞-===-==Λ=Λ==∑∑∑∑∑∑同理可得112()diag((),(),,())n f At P f t f t f t P λλλ-=例3.7 已知460350361A ⎛⎫ ⎪=-- ⎪ ⎪--⎝⎭,求e At,cos A .解 可求得2det()(2)(1)I A λλλ-=+-,即A 的特征值为12λ=-,231λλ==.对应12λ=-的特征向量为T 1(1,1,1)p =-,对应231λλ==的两个线性无关的特征向量为T 2(2,1,0)p =-,T 3(0,0,1)p =.于是120110101P --⎛⎫ ⎪= ⎪ ⎪⎝⎭ 使得1211P AP --⎛⎫⎪= ⎪ ⎪⎝⎭故22212222e 2e e 2e 2e 0e e e e2e e 0e e e 2e 2e e tt t t t At tt t t t t t t t tt P P --------⎛⎫⎛⎫-- ⎪ ⎪==--⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭1cos(2)cos cos1cos12cos1cos 22cos12cos 20cos 2cos12cos 2cos10cos 2cos12cos 22cos1cos1A P--⎛⎫⎪= ⎪ ⎪⎝⎭--⎛⎫ ⎪=-- ⎪ ⎪--⎝⎭方法三 利用Jordan 标准形 设Cn nA ⨯∈,且C n nn P ⨯∈,使得121s J J P AP J J -⎛⎫ ⎪⎪== ⎪ ⎪⎝⎭ 其中×1(1,2,,)1i iii i i r rJ i s λλλ⎛⎫ ⎪⎪== ⎪⎪⎝⎭由定理1.12得111111001(1)01(1C C ()C ()()1!(1)!()1!()()1!(1)!i i i i i i i r k r k k i k i k ik k k k i i k i k k k k k ik i r r k k k i k k k k k tr r i f J t a J t a t t t r a t t t f f f r λλλλλλλλλλλλλλλλ--+-+∞+∞-==--+∞==--⎛⎫ ⎪⎪== ⎪ ⎪ ⎪⎝⎭⎛⎫' ⎪- ⎪⎪= ⎪⎪'⎪ ⎪⎝⎭'-=∑∑∑)()()()1!()i tt f t f f λλλλλ=⎛⎫ ⎪ ⎪⎪ ⎪⎪' ⎪ ⎪⎝⎭从而1010110011()()()()()k kk kk k k k k k k k k k k k k k k s k s f At a A t a PJP t a J tP a J t P P P a J t f J t P Pf J t +∞+∞-==+∞=+∞--=+∞=-==⎛⎫⎪ ⎪== ⎪⎪ ⎪ ⎪⎝⎭⎛⎫ ⎪= ⎪ ⎪⎝⎭∑∑∑∑∑例3.8 已知101120403A -⎛⎫⎪= ⎪ ⎪-⎝⎭,求e A,sin At .解 例1.9已求得100111210P ⎛⎫ ⎪=-- ⎪ ⎪⎝⎭,11112P AP J -⎛⎫⎪== ⎪ ⎪⎝⎭于是12222e e e 0ee e 3e e e 2e+e e 4e 03e A P P -⎛⎫⎛⎫ ⎪ ⎪==- ⎪⎪ ⎪ ⎪-⎝⎭⎝⎭-- 1sin cos sin sin sin 2sin 2cos 0cos sin 2cos sin 2sin 2cos sin sin 24cos 02cos sin t t t At P t Pt t t t t t t t t t t t t t t t t t t t -⎛⎫⎪= ⎪ ⎪⎝⎭-⎛⎫ ⎪=+---+ ⎪ ⎪-+⎝⎭根据Jordan 标准形理论可得 定理3.8 设Cn nA ⨯∈,1λ,2λ,…,n λ是A 的n 个特征值,则矩阵函数f (A )的特征值为1()f λ,2()f λ,…,()n f λ. 方法四 待定系数法 设Cn nA ⨯∈,且A 的特征多项式为1212()det()()()()(3.5)srr r s I A ψλλλλλλλλ=-=---其中1λ,2λ,…,s λ是A 的全部互异特征值,12s r r r n +++= .为计算矩阵函数()k kk k f At a A t +∞==∑,记0()k k k k f t a t λλ+∞==∑.将f (λt )改写为()(,)()(,)(3.6)f t q t r t λλψλλ=+其中q (λ,t )是含参数t 的λ的幂级数,r (λ,t )是含参数t 且次数不超过n -1的λ的多项式,即1110(,)()()()n n r t b t b t b t λλλ--=+++由Hamilton-Cayley 定理知ψ(A )=O ,于是由式(3.6)得1110()(,)()(,)()()()n n f At q A t A r A t b t Ab t A b t Iψ--=+=+++可见,只要求出()(0,1,,1)k b t k n =- 即可得到f (At ).注意到()()0(0,1,,1;1,2,,)l i i l r i s ψλ==-=将式(3.6)两边对λ求导,并利用上式,得d d ()(,)d d iil ll l f t r t λλλλλλλλ=== 即d d ()(,)(0,1,,1;1,2,,)(3.7)d d iil l li l l t t f r t l r i s μλλλμλμλ====-=由式(3.7)即得到以0()b t ,1()b t ,…,1()n b t -为未知量的线性方程组. 综上分析,用待定系数法求矩阵函数f (At )或f (A )的步骤如下: 第一步:求矩阵A 的特征多项式(3.5);第二步:设1110()n n r b b b λλλ--=+++ .根据()()()()(0,1,,1;1,2,,)i l l l i i tr t f l r i s λλλλ===-=或()()()()(0,1,,1;1,2,,)l l i i i r f l r i s λλ==-=列方程组求解0b ,1b ,…,1n b -;第三步:计算1110()(())()n n f At f A r A b A b A b I --==+++ 或.例3.9 已知101120403A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭,求e At,cos A .解 可求得2det()(1)(2)I A λλλ-=--.设2210()r b b b λλλ=++则由210212210(1)e (1)2e (2)42e tt t r b b b r b b t r b b b ⎧=++=⎪'=+=⎨⎪=++=⎩解得222120e e e 2e 2e 3e e 2e t t t t t t t t b t b t b t ⎧=--⎪=-++⎨⎪=-⎩于是2222210e 2e 0e e e e 2ee e e e 4e 02e e t tAt t t tt t t t t ttt t b A b A b I t t t t t ⎛⎫-⎪=++=-++-- ⎪ ⎪ ⎪-+⎝⎭而由21021210(1)cos1(1)2sin1(2)42cos 2r b b b r b b r b b b =++=⎧⎪'=+=-⎨⎪=++=⎩解得210sin1cos1cos 23sin12cos12cos 22sin1cos 2b b b =-+⎧⎪=-+-⎨⎪=+⎩从而22102sin1cos 20sin1cos 2sin1cos1cos 2cos 2sin1cos1cos 24sin102sin1cos1A b A b A b I +-⎛⎫ ⎪=++=-+--+ ⎪ ⎪-+⎝⎭如果求得矩阵A 的最小多项式,且其次数低于A 的特征多项式的次数,则计算矩阵函数就要容易一些.例3.10 已知311202113A -⎛⎫ ⎪=- ⎪ ⎪--⎝⎭,求e At ,sin A . 解 例1.9已求得A 的Jordan 标准形为2212J ⎛⎫⎪= ⎪ ⎪⎝⎭于是A 的最小多项式为2()(2)A m λλ=-.设10()r b b λλ=+由21021(2)2e (2)e t tr b b r b t ⎧=+=⎪⎨'==⎪⎩ 解得2120e (12)et t b t b t ⎧=⎪⎨=-⎪⎩ 于是2101e e 21221At t tt t b A b I t tt t t t +-⎛⎫ ⎪=+=-- ⎪ ⎪--+⎝⎭又由101(2)2sin 2(2)cos 2r b b r b =+=⎧⎨'==⎩ 解得10cos 2sin 22cos 2b b =⎧⎨=-⎩从而10sin 2cos 2cos 2cos 2sin 2cos 2sin 22cos 22cos 2cos 2cos 2sin 2cos 2A b A b I +-⎛⎫ ⎪=+=-- ⎪ ⎪--+⎝⎭三、常用矩阵函数的性质常用的矩阵函数有e A,sin A ,cos A ,它们有些性质与普通的指数函数和三角函数相同,但由于矩阵乘法不满足交换律,从而有些性质与一般指数函数和三角函数不相同. 定理3.9 对任意Cn n A ⨯∈,总有(1)sin(-A )=-sin A ,cos(-A )=cos A ; (2)i e cos isin AA A =+,i -i 1cos (e e )2A A A =+,i -i 1sin (e e )2iA A A =-. 证 (1)由sin A 与cos A 的矩阵幂级数形式直接得到;(2)i 221000i (1)(1)e i !(2)!(21)!cos isin k k k Ak k k k k k A A A k k k A A+∞+∞+∞+===--==++=+∑∑∑又有-i e cos()isin()cos isin A A A A A =-+-=- 从而i -i 1cos (e e )2A A A =+,i -i 1sin (e e )2iA A A =- 定理3.10 设A ,C n nB ⨯∈,且AB =BA ,则(1)ee e e e A BA B B A +==;(2)sin(A +B )=sin A cos B +cos A sin B ;(3)cos(A +B )=cos A cos B -sin A sin B .证 (1)0022011e e !!1()(2)2!1()e !A Bk k k k k A B k A B k k I A B A AB B A B k +∞+∞==+∞+=⎛⎫⎛⎫= ⎪⎪⎝⎭⎝⎭=++++++=+=∑∑∑(2)i()-i()i i -i -i i -i i -i i -i i -i 11sin()(e e )(e e e e )2i 2i 1111(e e )(e e )(e e )(e e )2i 222i sin cos cos sin A B A B A B A B A A B B A A B B A B A B A B+++=-=-=-+++-=+ 同理可证(3). 证毕在定理3.10中,取A =B ,即得 推论 对任意Cn nA ⨯∈,有22cos 2cos sin A A A =-,sin2A =2sin A cos A 值得注意的是,当AB ≠BA 时,ee e A BA B +=或e e e A B B A +=不成立.如取0010A ⎛⎫= ⎪⎝⎭,0100B ⎛⎫= ⎪⎝⎭,则0110A B ⎛⎫+= ⎪⎝⎭,00100100AB BA ⎛⎫⎛⎫=≠= ⎪ ⎪⎝⎭⎝⎭且10e 11A⎛⎫= ⎪⎝⎭,11e 01B ⎛⎫= ⎪⎝⎭,-1-1-1-1e+e e e 1e 2e e e+e A B+⎛⎫= ⎪⎝⎭-- 可见1121e e e e 1211A BB A ⎛⎫⎛⎫=≠= ⎪ ⎪⎝⎭⎝⎭e e e A B A B +≠,e e e A B B A +≠定理3.11 设Cn nA ⨯∈,则有(1)tr dete eA A=;(2)1(e )e A A --=.证 (1)设A 的特征值为1λ,2λ,…,n λ.则由定理3.8知,e A的特征值为1e λ,2e λ,…,e n λ,从而1212tr dete =e e e e e n n A A λλλλλλ++==…+…(2)由于tr dete =e0AA≠,所以e A 总是可逆的.又由定理3.10,得e e e e A A A A OI--===故1(e )e A A --=. 证毕需要指出的是,对任何n 阶方阵A ,e A总是可逆的,但sin A 与cos A 却不一定可逆.如取π00π/2A ⎛⎫=⎪⎝⎭,则00sin 01A ⎛⎫= ⎪⎝⎭,10cos 00A -⎛⎫= ⎪⎝⎭.可见sin A 与cos A 都不可逆.§3.4 矩阵的微分和积分在研究微分方程组时,为了简化对问题的表述及求解过程,需要考虑以函数为元素的矩阵的微分和积分.在研究优化等问题时,则要碰到数量函数对向量变量或矩阵变量的导数,以及向量值或矩阵值函数对向量变量或矩阵变量的导数.本节简单地介绍这些内容. 一、函数矩阵的微分和积分定义 3.6 以变量t 的函数为元素的矩阵()(())i j m n A t a t ⨯=称为函数矩阵,其中()(1,2,,;1,2,,)ij a t i m j n == 都是变量t 的函数.若t ∈[a ,b ],则称A (t )是定义在[a ,b )上的;又若每个()ij a t 在[a ,b ]上连续、可微、可积,则称A (t )在[a ,b ]上是连续、可微、可积的.当A (t )可微时,规定其导数为()(())ijm n A t a t ⨯''=或d d ()()d d ij m nA t a t t t ⨯⎛⎫= ⎪⎝⎭而当A (t )在[a ,b ]上可积时,规定A (t )在[a ,b ]上的积分为()()d ()d bb ijaam nA t t a t t ⨯=⎰⎰例3.11 求函数矩阵23sin cos ()2e 01t t t t t A t t t ⎛⎫⎪= ⎪ ⎪⎝⎭的导数. 解2cos sin 1d ()2ln 2e 2d 003t t t t A t t t t -⎛⎫⎪= ⎪ ⎪⎝⎭关于函数矩阵,有下面的求导法则.定理3.12 设A (t )与B (t )是适当阶的可微矩阵,则(1)d d d(()())()()d d d A t B t A t B t t t t+=+ (2)当λ(t )为可微函数时,有d d d (()())()()()()d d d t A t t A t t A t t t t λλλ⎛⎫=+ ⎪⎝⎭(3)d d d (()())()()()()d d d A t B t A t B t A t B t t t t ⎛⎫=+ ⎪⎝⎭; (4)当u =f (t )关于t 可微时,有d d()()()d d A u f t A u t u'= (5)当1()A t -是可微矩阵时,有111d d (())()()()d d A t A t A t A t t t ---⎛⎫=- ⎪⎝⎭证 只证(2)和(5).设()(())ij m n A t a t ⨯=,()(())ij n p B t b t ⨯=,则111d d (()())(()())d d d d ()()()()d d d d ()()()()d d nik kj m n k n nik kj ik kj k k m nA tB t a t b t t t a t b t a t b t t t A t B t A t B t t t ⨯===⨯=⎛⎫⎛⎫=+ ⎪ ⎪⎝⎭⎝⎭⎛⎫=+ ⎪⎝⎭∑∑∑由于1()()A t A t I-=,两边对t 求导,得11d d ()()()()d d A t A t A t A t O t t --⎛⎫+= ⎪⎝⎭从而111d d ()()()()d d A t A t A t A t t t ---⎛⎫=- ⎪⎝⎭证毕 定理3.13 设C n nA ⨯∈,则有(1)d e e e d AtAt At A A t ==; (2)dsin cos (cos )d At A At At A t ==;(3)dcos sin (sin )d At A At At A t=-=-.证 这里只证(1).(2)和(3)的证明与(1)类似.由0e !k Atkk t A k +∞==∑,并利用绝对收敛级数可以逐项求导,得101111d d e d d !(1)!e (1)!k k At k k k k k k Atk t t A A t t k k tA A A k -+∞+∞==-+∞-===-==-∑∑∑同样11111d e ==e d (1)!(1)!k k At k k At k k t t A A A A t k k --+∞+∞-==⎛⎫= ⎪--⎝⎭∑∑ 证毕根据定义和积分的有关性质,可得定理3.14 设A (t ),B (t )是区间[a ,b ]上适当阶的可积矩阵,A ,B 是适当阶的常数矩阵,λ∈C ,则 (1)(()())d ()d ()d bb baaaA tB t t A t t B t t +=+⎰⎰⎰;(2)()d ()d bba aA t t A t t λλ=⎰⎰;(3)()()d ()d bbaaA tB t A t t B =⎰⎰,()d ()d b baaAB t t A B t t =⎰⎰;(4)当A (t )在[a ,b ]上连续时,对任意t ∈(a ,b ),有()d ()d ()d t aA A t tττ=⎰(5)当A (t )在[a ,b]上连续可微时,有()d ()()baA t t A b A a '=-⎰以上介绍了函数矩阵的微积分概念及一些运算法则.由于d()d A t t仍是函数矩阵,如果它仍是可导矩阵,即可定义其二阶导数.不难给出函数矩阵的高阶导数11d d d ()()d d d k k k k A t A t t t t --⎛⎫= ⎪⎝⎭二、数量函数对矩阵变量的导数定义 3.7 设f (X )是以矩阵()ij m n X x ⨯=为自变量的mn 元函数,且(1,2,,;1,2,,)ijfi m j n x ∂==∂ 都存在,规定f 对矩阵变量X 的导数d d f X 为 1111d d ij m nm mn ff x x n f fX x ff x x ⨯∂∂⎛⎫ ⎪∂∂ ⎪⎛⎫∂ ⎪== ⎪ ⎪ ⎪∂⎝⎭∂∂ ⎪ ⎪∂∂ ⎪⎝⎭特别地,以T12(,,,)n x ξξξ= 为自变量的函数f (x )的导数T12d (,,,)d nf f f f x ξξξ∂∂∂=∂∂∂ 称为数量函数对向量变量的导数,即为在数学分析中学过的函数f 的梯度向量,记为grad f .例 3.12 设T 12(,,,)n a a a a = 是给定的向量,T 12(,,,)n x ξξξ= 是向量变量,且T T ()f x a x x a ==求d d f x. 解 因为1()nk kk f x a ξ==∑而(1,2,,)j jfa j n ξ∂==∂ 所以 TT 1212d (,,,)(,,,)d n nf f f f a a a a x ξξξ∂∂∂===∂∂∂ 例3.13 设()ij m n A a ⨯=是给定的矩阵,()ij n m X x ⨯=是矩阵变量,且()tr()f x Ax =求d d fX. 解 因为1()nikkj m m k AX ax ⨯==∑.所以11()tr()m nsk ks s k f X AX a x ====∑∑而(1,2,,;1,2,,)ijfi n j m x ∂==∂ 故T d ()d ji n m ij n mf f a A X x ⨯⨯⎛⎫∂=== ⎪ ⎪∂⎝⎭ 例 3.14 设()ij n n A a ⨯=是给定的矩阵,T 12(,,,)n x ξξξ= 是向量变量,且T ()f x x Ax =求d d f x. 解 因为T1111()()n nn ns sk ks sk k s k s k f x x Ax aa ξξξξ=======∑∑∑∑而1111,11,111()nj j j j jk k j jj j j j n nj k j n nsj s jk ks k fa a a a a a a a ξξξξξξξξξ--++===∂=+++++++∂=+∑∑∑所以1111111T T d d ()n ns s k k s k n nsn s nk k s k n f a a f x f a a A x Ax A A xξξξξξξ====∂⎛⎫⎛⎫+ ⎪ ⎪∂⎪ ⎪ ⎪== ⎪ ⎪ ⎪∂ ⎪ ⎪+ ⎪⎪∂⎝⎭⎝⎭=+=+∑∑∑∑ 特别地,当A 是对称矩阵时,有d 2d fAx x=例3.15 设()ij n n X x ⨯=是矩阵变量,且det X ≠0.证明1T ddet (det )()d X X X X-= 证 设ij x 的代数余子式为ij X .把det X 按等i 行展开,得1det nikik k X xX ==∑于是det ij ijX X x ∂=∂故 T1T 1Tddet det ()(adj )d ((det ))(det )()ij n n ij n nX X X X X x X X X X ⨯⨯--⎛⎫∂=== ⎪ ⎪∂⎝⎭== 三、矩阵值函数对矩阵变量的导数定义3.8 设矩阵()(())ij s t F X f X ⨯=的元素()(1,2,,;1,2,,)ij f X i s j t == 都是矩阵变量()ij m n X x ⨯=的函数,则称F (X )为矩阵值函数,规定F (X )对矩阵变量X 的导数d d FX为111d d 1FF x x n F X FF x x m mn ∂∂⎛⎫ ⎪∂∂ ⎪⎪= ⎪∂∂ ⎪ ⎪∂∂ ⎪⎝⎭ ,其中1111tij s stf f x x ij ij F x f f x x ij ij ∂∂⎛⎫ ⎪∂∂ ⎪∂⎪=⎪∂ ⎪∂∂ ⎪∂∂ ⎪⎝⎭即其结果为(ms )×(nt )矩阵. 作为特殊情形,这一定义包括了向量值函数对于向量变量的导数,向量值函数对于矩阵变量的导数,矩阵值函数对于向量变量的导数等.例3.16 设T12(,,,)n x ξξξ= 是向量变量,求T T d d d d x xx x=. 解 由定义,得T 1TT 2T 100010d d 001n nx x x I x x ξξξ⎛⎫∂ ⎪∂ ⎪⎛⎫ ⎪∂ ⎪⎪ ⎪===∂ ⎪ ⎪⎪ ⎪⎪⎝⎭ ⎪∂ ⎪∂⎝⎭同理可得T 12d ,,,d n n x x x x I x ξξξ⎛⎫∂∂∂== ⎪∂∂∂⎝⎭例3.17 设T1234(,,,)a a a a a =是给定向量,24()ij X x ⨯=是矩阵变量,求Td()d Xa X,d()d Xa X. 解 因为41121k k k n k k k x a Xa x a ==⎛⎫⎪ ⎪= ⎪ ⎪⎝⎭∑∑,44T 1211()(,)k k k k k k Xa x a x a ===∑∑ 所以T T TT T 13141112T T T T 2122232431243124()()()()d()d ()()()()00000000Xa Xa Xa Xa x x x x Xa XXa Xa Xa Xa xx x x a aa a a a a a ⎛⎫∂∂∂∂⎪∂∂∂∂ ⎪=⎪∂∂∂∂ ⎪ ⎪∂∂∂∂⎝⎭⎛⎫= ⎪⎝⎭而131411122122232412341234()()()()d()()()()()d 00000000Xa Xa Xa Xa x x x x Xa Xa Xa Xa Xa Xxx x x a a a a a a a a ∂∂∂∂⎛⎫⎪∂∂∂∂ ⎪=⎪∂∂∂∂ ⎪∂∂∂∂⎝⎭⎛⎫⎪ ⎪=⎪⎪⎝⎭§3.5 矩阵分析应用举例本节介绍矩阵函数及矩阵微积分的一些应用. 一、求解一阶线性常系数微分方程组在数学或工程技术中,经常要研究一阶常系数微分方程组1111122112211222221122d ()()()()()d d ()()()()()d d ()()()()()d n n n n n n n nn n n x t a x t a x t a x t f t t x t a x t a x t a x t f t t x t a x t a x t a x t f t t ⎧=++++⎪⎪⎪=++++⎨⎪⎪=++++⎪⎩满足初始条件0()(1,2,,)i ix t c i n ==的解.如果记T12(),(,,,)ij n n n A a c c c c ⨯==T 12()((),(),,())n x t x t x t x t = ,T 12()((),(),,())n f t f t f t f t =则上述微分方程组可写为0d ()()()(3.8)d ()x t Ax t f t tx t c⎧=+⎪⎨⎪=⎩因为d d ()(e ())e ()()e d d d ()e ()e ()d At At At At At x t x t A x t t t x t Ax t f t t -----=-+⎛⎫=-= ⎪⎝⎭将上式两边在[0t ,t ]上积分,得00d (e ())d e ()d d tt A A t t x f τττττττ--=⎰⎰ 即00e()e()e ()d tAt A A t x t x t f ττττ----=⎰于是微分方程组的解为00()()e e e ()d tA t t At A t x t c f τττ--+⎰=例3.18 求解微分方程组初值问题113212313123d ()()()1d d ()()2()1d d ()4()3()2d (0)1,(0)0,(0)1x t x t x t t x t x t x t tx t x t x t t x x x ⎧=-++⎪⎪⎪=+-⎪⎨⎪=-++⎪⎪⎪===⎩ 解 记123()10111120,0,()(),()140312()x t A c x t x t f t x t -⎛⎫⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪ ⎪====- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭⎝⎭则微分方程组可以写成式(3.8)的矩阵形式.例3.9已求得222e 2e 0e e e e 2e e e e e 4e 02e e t t tAt t t tt t t t t t t t t t t t t ⎛⎫-⎪=-++-- ⎪ ⎪-+⎝⎭依次计算下列各量e e e e e 2e t t At t t t t c t t ⎛⎫- ⎪= ⎪ ⎪-⎝⎭,00e 1e e ()d e 1e 2e 22e t t t A t tf d τττττττ-------⎛⎫⎛⎫- ⎪ ⎪=-=-+ ⎪⎪ ⎪ ⎪-⎝⎭⎝⎭⎰⎰, 0e 1e e ()d e 12e 2t t At A t tf τττ-⎛⎫- ⎪=-+ ⎪ ⎪-⎝⎭⎰故微分方程组的解为123e e e 1(2)e 1()()()e e 1(1)e 1()e 2e 2e 2(32)e 2t t t tt t t t t t t t t x t x t x t t t x t t ⎛⎫⎛⎫⎛⎫----⎛⎫⎪ ⎪ ⎪ ⎪==+-+=-+ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪----⎝⎭⎝⎭⎝⎭⎝⎭二、求解矩阵方程在控制论与系统理论中,要遇到形如AX +XB =F 的矩阵方程求解问题,这个矩阵方程也称为Lyapunov 方程.关于这个矩阵方程的解有如下结果. 定理3.15 给定矩阵方程 AX +XB =F (3.9) 其中Cm mA ⨯∈,Cn nB ⨯∈,Cm nF ⨯∈.如果A 和B 的所有特征值具有负实部(这种矩阵称为稳定矩阵),则该矩阵方程有惟一解e e d At Bt X F t +∞=-⎰证 记()e e At Bt Y t F =.则有Y (0)=F ,且d ()e e e e ()()(3.10)d At Bt At Bt Y t A F F B AY t Y t B t=+=+设12,,,m λλλ 是A 的m 个特征值,12,,,n μμμ 是B 的n 个特征值.根据利用Jordan 标准形求矩阵函数的方法(见§3.3)知,e At的元素是形如e (0)j tr t r λ≥的项的线性组合.因为A 的所有特征值j λ的实部是负的,所以lim eAtt O →+∞=.同理lim e Bt t O →+∞=.于是lim ()lim e e At Bt t t Y t F O →+∞→+∞==又由于e e At BtF 的元素是形如()e (0)i j tr t r λμ+≥的项的线性组合,且积分()0ed i j tr t t λμ+∞+⎰都存在,故积分e e d At Bt F t +∞⎰存在.对式(3.10)两边从0到+∞积分,得()()0()(0)()d ()d Y Y AY t t Y t t B +∞+∞+∞-=+⎰⎰即()()0()d ()d A Y t t Y t t B F+∞+∞-+-=⎰⎰这说明0e e d At Bt X F t +∞=-⎰是矩阵方程(3.9)的解.惟一性的证明见第七章. 证毕 推论1 设Cm mA ⨯∈,Cn nB ⨯∈,Cm nF ⨯∈,则矩阵微分方程d ()()()d (0)X t AX t X t B tX F⎧=+⎪⎨⎪=⎩的解为()e e At BtX t F =推论2 设A ,C n nF ⨯∈,且A 的所有特征值具有负实部,则矩阵方程HA X XA F+=-的惟一解为H 0ee d (3.11)A tAt X F t+∞=⎰如果F 是Hermite 正定矩阵,则解矩阵X 也是Hermite 正定矩阵.证 只需证明后一结论.当F 是Hermite 正定矩阵时,由式(3.11)可知X 是Hermite 矩阵.又对0Cnx ≠∈,由于eAt总是可逆的,所以e 0Atx ≠,于是HH H e e (e )(e )0A t At At At x F x x F x =>.从而HH 0(e )(e )d 0At At x Xx x F x t +∞=>⎰故X 是Hermite 正定矩阵. 证毕三、最小二乘问题 设Cm nA ⨯∈,C n b ∈.当线性方程组Ax =b 无解时,则对任意C nx ∈都有Ax -b ≠0.此时希望找出这样的向量0C n x ∈,它使2Ax b -达到最小,即022Clim (3.12)nx Ax b Ax b ∈-=-称这个问题为最小二乘问题,称0x 为矛盾方程组Ax =b 的最小二乘解.以下结论给出了当A ,b 分别是实矩阵和实向量时,Ax =b 的最小二乘解所满足的代数方程.定理3.16 设R m nA ⨯∈,R mb ∈,0R n x ∈.若0R n x ∈是Ax =b 的最小二乘解,则0x 是方程组TT(3.13)A Ax A b=的解.称式(3.13)为Ax =b 的法方程组.证 由于2T 2TTTTTT()()()f x Ax b Ax b Ax b x A Ax x A b b Ax b b=-=--=--+若0x 为Ax =b 的最小二乘解,则它应是f (x )的极小值点,从而d 0(3.14)d x f x=根据例3.12和例3.14,得T T d 22d fA Ax A b x=- 由式(3.14)即知T T00A Ax A b -=,故0x 是式(3.13)的解. 证毕 对于含约束条件的最小二乘问题,有如下的结果. 例3.19 设Rm nA ⨯∈,R m b ∈,Rk nB ⨯∈,R kd ∈,且Bx =d 有解.试求约束极小。

华中科技大学研究生矩阵论课件

华中科技大学研究生矩阵论课件

子空间的“和”为“直和”的充要–条件 :
定理1·8 设W=W1+W2,则下列各条等价:
(1)
W=W1W2
(2)
X W,X=X 1+X2的表
是惟一的
(3) W中零向量的表示是惟一的
(4)
dim W =dimW1+dimW2
.
26
例1
P12 eg18
例2 设在Rn×n中,子空间
W 1={A AT =A } , W2={B BT= –B }, 证明Rn×n=W1W2。
线性空间的一般性的观点:
线性空间的一般形式:
V(F),元素被统称为向量:, ,,
线性空间的简单性质(共性):
定理1 . 1:V(F)具有性质:
(1) V(F)中的零元素是惟一的。
(2) V(F)中任何元素的负元素是惟一的。
(3)数零和零元素的性质: 数0 0=0,k0=0,k =0 =0 或k=0
例3 子空间W的“直和补子空间”
.
27
1·2 内积空间
主题:定义内积的概念,借助于内积建立线性 空间的度量关系。
一、 欧氏空间和酉空间 1 几何空间中度量关系的定义基础 2 内积的定义 定义1·7 (P13) :要点 • 内积(,)是二元运算:Vn(F) F • (,)的公理性质 • (,)是任何满足定义的运算。 • 讨论(,1+2), (,k)
(II);{ 2 1 0 1 0 0 0 0 }
0
0
1
0
3
1
0
3
1. 求从基(I)到基(II)的过渡矩阵C。
2. 求向量 7 3 在基(II)的坐标Y。
1
2
§1.1 五、 子空间
概述:线性空间Vn(F)中,向量集合V可 以有集合的运算和关系: Wi V, W1W2, W1W2, 问题: 这些关系或运算的结果是否仍然为 线性空间 ?

北京化工大学《矩阵论》课件3-4,5

北京化工大学《矩阵论》课件3-4,5

df 。 dX
f ( X ) t r AX ( ) = a st x ts = a ji xij + a st xts
s 1 t 1
m
n
m
n
s 1 t 1 s j 或t i
因为
f a ji ,所以 xij
f df ) nm (a ji ) nm = AT =( x dX ij
b b a k 1 k 1 a
n
n
证毕 性质 10 性质 11
d dx

x a
A(t ) d t = A( x) ,其中 A(t ) 在 [a, b] 上连续, x [a, b] 。

b a
(t ) 在 [a, b] 上连续。 A(t ) d t A(b) A(a) ,要求 aij
a b b a a
b
, C

b a
A(t ) B d t ( A(t ) d t ) B , AB (t ) d t A B(t ) d t ,
a
b
其中第一式中 B 是常数矩阵,第二式中 A 是常数矩阵。
89

设 A(t ) = (aij (t )) mn ,B = (bij ) nl , 则 ( a ik (t )bkj ) d t = ( a ik (t ) d t )bkj 。
定义
若函数矩阵 A(t ) = (aij (t )) mn 的每个元素 aij (t ) 对变量 t 均可导(可微) ,
则称 A(t ) 为可导(或可微)矩阵,其导数矩阵为
d dt
A(t ) (或 A(t ) )= ( ddt aij (t )) mn
如果 A(t ) 的每个元素 aij (t ) 都在区间 [a, b] 上可积, 则称 A(t ) 是区间 [a, b] 上的可积 矩阵,其积分矩阵为

01_矩阵论_第一章

01_矩阵论_第一章

注记 3 线性空间的本质是线性运算。同一 个集合,若定义两种不同的线性运算,就构成 不同的线性空间;若定义的运算不是线性运算, 就不能构成线性空间。如前述的数学例子。
注记 4 抽象的线性空间的作用在于,由它 得出的一切结论对诸如上述线性空间的研究。
例 5 向量组 {e1 = (1, 0, 0, …, 0)T, e2 = (0, 1, 0, …, 0)T, …, en = (0, 0, …, 0, 1)T} 是 F n 的一组基, 则 dimF n = n。
例 6 求矩阵空间 R22 的维数与一组基。
解 任取矩阵 A,其中
a11 A a 21 a12 a22
a0 a 2 3 1 f x 1, x, x , x , a2 a 3


因为 所以 f(x) 在基 {1, x, x2, x3} 下的坐标为 (a0, a1, a2, a3)T。
在线性空间和线性变换的讨论中,矩阵是 一个重要的工具。特别地,矩阵还是线性变换 的便利表达方法。
§ 1.1 线性空间
一、线性空间的概念
在线性代数课程中,我们把有序数组称为 向量,把 n 维向量的全体所构成的集合 Rn 称为 n 维向量空间。一般地,如果 V 为非空的 n 维向 量的集合,且集合 V 对于向量加法及数乘两种 运算封闭,那么就称集合 V 为向量空间。
注记 1 有些教材中,向量空间与线性空间 表示的是同一个概念,但我们通常用向量空间 来表示某一数域上的以该数域中的 n 元有序数组 为元素构成的线性空间。
此外,从上述线性空间的例子中可以看到, 许多常见的研究对象都可以在线性空间中作为 向量来研究,只不过它们此时未必是有序数组 了。
注记 2 定义 1.1 中的加法和数乘运算分别 是V 中的一个二元运算以及数域 F 和 V 中元素 间的运算,它们已不再局限于数的加法、乘法 或者数值向量的加法、数乘概念。如上述的几 个例题。

矩阵论课件

矩阵论课件

第一讲线性空间一、线性空间的定义及性质[知识预备]★集合:笼统的说是指一些事物(或者对象)组成的整体集合的表示:枚举、表达式集合的运算:并(),交()另外,集合的“和”(+):并不是严格意义上集合的运算,因为它限定了集合中元素须有可加性。

★数域:一种数集,对四则运算封闭(除数不为零)。

比如有理数域、实数域(R)和复数域(C)。

实数域和复数域是工程上较常用的两个数域。

线性空间是线性代数最基本的概念之一,也是学习现代矩阵论的重要基础。

线性空间的概念是某类事物从量的方面的一个抽象。

1.线性空间的定义:设V是一个非空集合,其元素用x,y,z等表示;K是一个数域,其元素用k,l,m等表示。

如果V满足[如下8条性质,分两类](I)在V中定义一个“加法”运∈时,有唯一的和算,即当x,y V+∈(封闭性),且加法运算x y V满足下列性质(1)结合律()()++=++;x y z x y z(2)交换律 x y y x +=+;(3)零元律 存在零元素o ,使x +o x =;(4)负元律 对于任一元素x V ∈,存在一元素y V ∈,使x y +=o ,且称y 为x 的负元素,记为(x -)。

则有()x x +-= o 。

(II )在V 中定义一个“数乘”运算,即当x V ∈,k K ∈时,有唯一的kx V ∈(封闭性),且数乘运算满足下列性质(5)数因子分配律()+=+;k x y kx ky(6)分配律()+=+;k l x kx lx(7)结合律()()=;k lx kl x=;(8)恒等律1x x [数域中一定有1]则称V为数域K上的线性空间。

注意:1)线性空间不能离开某一数域来定义,因为同一个集合,如果数域不同,该集合构成的线性空间也不同。

(2)两种运算、八条性质数域K中的运算是具体的四则运算,而V中所定义的加法运算和数乘运算则可以十分抽象。

(3)除了两种运算和八条性质外,还应注意唯一性、封闭性。

唯一性一般较显然,封闭性还需要证明,出现不封闭的情况:集合小、运算本身就不满足。

矩阵论浙大研究生第3讲

矩阵论浙大研究生第3讲
向量的内积与范数
函数向量的内积与范数
矩阵论-矩阵代数基础
3
矩阵论-矩阵代数基础
4
1
内积与范数
向量的内积与范数 函数向量的内积与范数 两个函数向量的夹角定义
内积与范数
矩阵的内积与范数 向量内积与范数的推广 mn×1向量 拉长
随机向量的内积与范数 等价于 正交性
矩阵论-矩阵代数基础 5 矩阵论-矩阵代数基础 6
则秩1更新
( 1R 1x )( 1R 1x ) H 1 1x H R 1x ˆ (n) R ˆ (n 1) x(n)x H (n ) 为 R
R 1 ( n ) 1R 1 (n 1) g( n ) g H ( n ) ——更新公式 ˆ 1 (n )x(n ) ˆ 1 (n 1)x(n ) g(n) R g (n) 1R
内积与范数
矩阵的内积与范数 诱导范数
内积与范数
矩阵的内积与范数 诱导范数 例题
常用的诱导范数-p范数 绝对列和范数 绝对列和范数 p取不同的数 矩阵的最大奇异值 绝对行和范数
矩阵论-矩阵代数基础 7 矩阵论-矩阵代数基础 8
绝对行和范数
2
内积与范数
矩阵的内积与范数 “元素形式”范数
随机向量
概率密度函数 随机向量
Kronecker积
向量化与矩阵化
按列堆栈
按行堆栈
矩阵论-矩阵代数基础
27
矩阵论-矩阵代数基础
28
7
稀疏表示和压缩感知
一个含有大多数零元素的向量或矩阵称为稀疏向量 (sparse vector)或者稀疏矩阵(sparse matrix)
第一章习题
非稀疏
1.10,1.11(二选一) 1.14,1.17,1.19,1.37,1.82 下周三交

矩阵论—H-L定理精品PPT课件

矩阵论—H-L定理精品PPT课件

B1 B0 A a1E
B2 B1 A a2E

Bn1 Bn Bn1
2 A an1 A anE
E
以An , An1, , A, E 依次右乘③的第一式、第二式、
…、第n式、第n+1式,得
B0 An An B1 An1 B0 An a1 An1 B2 An2 B1 An1 a2 An
写在最后
成功的基础在于好的学习习惯
The foundation of success lies in good habits
24
谢谢聆听
·学习就是为了达到一定目的而努力去干, 是为一个目标去 战胜各种困难的过程,这个过程会充满压力、痛苦和挫折
Learning Is To Achieve A Certain Goal And Work Hard, Is A Process To Overcome Various Difficulties For A Goal
同理可得 gA( x) gB ( x). 又 gA( x), gB ( x)都是首1多项式, gA( x) gB ( x).
注:反之不然,即最小多项式相同的矩阵未必相似.
如:
1 1 0 0
1 1 0 0
A
0 0 0
1 0 0
0 1 0
0 0 2
,
B
0 0 0
1 0 0
0 2 0
0
0 2
的最小多项式皆为( x 1)2( x 2), 但A与B不相似.
再设
f ( ) n a n1 1
a a
n1
n
则, f ( )E nE a n1E a E a E ①
1
n1
n

B( )( E A) nB0 n1(B1 B0 A) n2(B2 B1A)

第三章 矩阵论

第三章 矩阵论

设n阶矩阵A的互异特征值 1 (n1重), 2 (n2 重)
, ,s (ns 重), ni n,A的特征多项式是
i 1 s
f ( ) I A ( 1 )n1 ( 2 )n2 ( s )ns
则A的最小多项式必有如下形式,
m( ) ( 1 )m1 ( 2 )m2 ( s )ms
定义 次数最低且首项系数为1的矩阵A的化 零多项式, 称为A的最小多项式, 记为 m( ) .
定理3.8 多项式 ( ) 是矩阵A的化零多项式 当且仅当 m( ) ( ) . 特别地, 有 m( ) f ( ) , 其 中 f ( )是A的特征多项式.
推论1 矩阵A的最小多项式是唯一的. 定理3.9 矩阵A的特征多项式、最小多项 式有相同的根.(重数可能不同)
1
的特征向量;
X 21 , X 22 , X 2l2是属于的特征值 2 的特征向量;
X k 1 , X k 2 , X klk 是属于的特征值 的特征向量; k

X11 , X12 , X1l1 , X 21 , X 22 , X 2 l2 ,, X k1 , X k 2 , X klk
定义 若A 经有限次初等变换后变为 B , 则称 A B 相抵.记为 A B . 与 相抵关系是 方阵的一种等价关系,具有 1.自反性 2.对称性 3.传递性
定理: B 的充要条件是存在两个可 A 逆矩阵 P 与 Q ,使得
AX X
例 已知三维线性空间V的基 1 , 2 , 3 , 线性变换T满足, T 1 1 2 2 2 3 T 2 2 1 2 2 3 T 3 2 1 2 2 3 求T的特征值与特征向量.

第二章矩阵论课件

第二章矩阵论课件
m m ki i , ki i , i 1 i 1
内积的四条规定可推出如下性质: 1° 2° 3° 4°


n n , l j j l j ( , j ) j1 j 1
n m m k i i , l j j j 1 i 1 i 1
x1 1 x2 2 xn n
y1 1 y2 2 yn n
由性质,知:
, xi y j i , j
n n
令:
aij i , j
a ji
i 1
j 1
i, j 1,2,, n
等价。 则 1 , 2 ,..., m为与
单位化: 令
i i , i 1,2,...,m, i
1 , 2 ,..., m 等价的标准正交向量组。
( n 定理2: 维内积空间必有标准正交基。n 0)
第二节
欧氏空间
定义8: 实数域上的内积空间称为欧几里得 [Euclid]空间, 简称欧氏空间。 由于欧氏空间是 实数域上的内积空间, 因而内积的共轭对称性 就成了对称性。 设 V 是 n 维欧氏空间,1 , 2 ,, n是 V 的一组基, 对于V 中两个向量:
i j ij
1 2 n
并设 设1 , 2 ,, n是 V 的标准正交基,
1 , 2 , n 1 , 2 , n A
即 , ,, , , A 由2)知 , ,, 也是 V 的标准正交基,
0
(2)
则 0 是 经单位化得到的单位向量。 定理1: [cauchy—schwarz不等式]对于内积 空间中任意向量 , ,有 ( , )

南航矩阵论第一章PPT课件

南航矩阵论第一章PPT课件

例如 , N Z, Z Q , Q R
显然有下列关系 :
(1) A A; A A; A
(2) A B且 B C
AC
.
9
机动 目录 上页 下页 返回 结束
定义 3 . 给定两个集合 A, B, 定义下列运算:
并集 A B x x A 或 x B 交集 A B x x A 且 x B
.
1
教材: 《矩阵论》,戴华编,科学出版社。 主要参考书:
1. 方保镕,周继东编,矩阵论,清华大学出版社,2004.
2. 刘慧等,矩阵论及应用,化学工业出版,2003. 3.程云鹏,矩阵论,西安工业大学出版,2000. 4. 罗家洪,矩阵分析引论,华南理工大学出版,2002.
.
2
第1章 线性空间与内积空间
X
f
Y
则称 f 为单射;
f (X)
若 f 既是满射又是单射, 则称 f 为双射 或一一映射.
.
29
机动 目录 上页 下页 返回 结束
设 f1 : A1 B1 f2 : A2 B2
则称映射 f1与f2相等, 如果 A1 A2 , B1 B2,并且对任意a A1有 f1(a) f2 (a)
为A B上的一个二元关系。
例 2:A={矩阵论五班学生}。
显然,均来自于南京的同学关系R是A上的 一个二元关系。
想一想: 在该例中还存在什么关系?
.
15
例 3:
A={ 张华,王兵,陈平,李兰
a1
a2
a3
a4
B={ 软件,硬件,自动化,遥感
b1
b2
b3
b4
则:R1={(a1,b1),(a1,b3),
.
38

第一章 矩阵论

第一章 矩阵论
A (a1 a2 an ) (a1 , a2 ,, an )
b1 b2 (3) 列矩阵 (又称为列向量) B bm
6
矩阵概念
(4)方阵中从左上角元素到右下角元素的元 素族称为主对角线.主对角线以外的元素都 是零的方阵称为对角矩阵,简称对角阵.记 为 a11 0 0 0 a22 0 diag a11 , a22 ,, ann 0 ann 0
第一章 矩阵论
1.1 矩阵的概念及其运算
1.2 矩阵的行列式
1.3 矩阵的初等变换与矩阵的秩
1.4 矩阵的逆
1.5 矩阵的分块、克莱默法则
1
1.1 矩阵的概念及其运算
1.1.1 矩阵的概念 1.1.2 矩阵的线性运算与乘法 1.1.3 矩阵的转置
2
矩阵概念
引例 产品分配问题:某厂向三个商店发送四个产品.
7
矩阵概念
(5)单位矩阵 主对角线上的元素全是1的对角阵
(6)上三角阵 主对角线下方所有元素均为零的 方阵; 下三角阵 主对角线上方所有元素均为零的 方阵.
1 0 E 0 0 0 1 0 0 1
a11 a12 a1n 0 a22 a2 n A 0 0 ann
a11 0 0 a21 a22 0 B an1 an 2 ann
8
1.1.2 矩阵的线性运算与乘法
线性运算
加法 :两同型矩阵
a11 a12 a21 a22 A am1 am 2 之和为 A B ( aij
b11 b12 b1n a1n b21 b22 b22 a 22 , B bm1 bm 2 bmn amn bij ) m n .

矩阵论简明教程整理全PPT课件

矩阵论简明教程整理全PPT课件

k
ei
e
H j
E ei , ej , k
第45页/共188页
Remark
det E u,v, det In uvH det 1 vHu
1 vHu (由n Im AB m In BA 得到)
第46页/共188页
四、其他特殊矩阵
1幂零矩阵:Ak 0, k : 某正整数; 2幂等矩阵:A2 A; 3 实对称正定矩阵:
a a jn 1 j1 2 j2
anjn
j1 j2 jn
第13页/共188页
二、块矩阵的行列式
1、设A Cmm , B Cmn , C Cnm , D Cnn , 则
1 A
0A
BA
0 AD
0D 0D CD
2 A B 1mn C D 1mn B A
CD
AB
DC
3 A B m A B
minrank A, rank B
第30页/共188页
推论1
设ACmn , B Cnk ,且AB 0,则
rank A rank B n
第31页/共188页
§1.4 特殊矩阵
一、 几类基本的特殊矩阵
1、零矩阵,单位矩阵 2、对角矩阵
a11
D
a22
diag
a11
,
a22
,
ann
第50页/共188页
§2.1 矩阵的特征值与特征向量
一、特征值与特征向量 1、定义 定义1
设ACnn ,若存在数 C和x Cn , x 0使得 Ax x
则称是A的特征值,x称为A属于的特征向量。
第51页/共188页
2、特征多项式 定义2
设ACnn , 称In A为A的特征矩阵,称detIn A 为A的特征多项式,称detIn A 0为A的特征方程。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

(Ⅱ) 定义的数乘运算封闭, 即
∀ x ∈V , ∀ k ∈ K , 对应唯一 元素(kx)∈V , 且满足 (5) 数对元素分配律: k( x + y) = kx + ky (∀y ∈V ) (6) 元素对数分配律: (k + l )x = kx + lx (∀l ∈ K ) (7) 数因子结合律: k(lx) = (kl )x (∀l ∈ K ) (8) 有单位数:单位数1∈ K , 使得 1x = x . 则称V 为 K 上的线性空间.
例 3 K = R 时, R n —向量空间;
R m×n —矩阵空间
第一章 线性空间与线性变换(第 1 节)
3
Pn[t]—多项式空间; C[a,b] —函数空间 K = C 时, Cn —复向量空间; Cm×n —复矩阵空间 例 4 集合 R + = {m m是正实数 } ,数域 R = {k k是实数 } .
0
a 12
a
22
ai
j1
I
S 2
=
{A
=
a11
0
0
a
22
a 11
, a22

R}
S 1
U
S 2
=
{A
=
a11 a21
a 12
a
22
aa 12 21
=
0,
ai
j

R}
S 1
+
S 2
=
{A
=
a11 a21
a 12
a 22
ai j ∈ R}
2.数域:关于四则运算封闭的数的集合.
2.减法运算:线性空间V 中, x − y = x + (− y) .
3.线性组合: x, xi ∈V , 若存在 ci ∈ K ,
使
x
=
c 1
x 1
+
L
+
cm
xm
,
则称 x

x 1
,L,
xm
的线性组合,或者
x
可由
x 1
,L,
xm
线性表示.
数乘封闭,(5)~(8)成立.故 R+ 是 R 上的线性空间. 例 5 集合 R 2 = {α = (ξ1 , ξ 2 ) ξ i ∈ R} ,数域 R .设 β = (η1 , η2 ), k ∈ R .
运算方式 1 加法: α + β = (ξ1 + η1 , ξ 2 + η2 ) 数乘: kα = (kξ1 , kξ2 )
(Ⅰ) 定义的加法运算封闭, 即
∀ x, y ∈V , 对应唯一 元素( x + y) ∈V , 且满足 (1) 结合律: x + ( y + z) = ( x + y) + z (∀z ∈V ) (2) 交换律: x + y = y + x
(3) 有零元: ∃θ ∈V , 使得 x + θ = x (∀x ∈V ) (4) 有负元: ∀x ∈V , ∃ (− x)∈V , 使得 x + (− x) = θ .
加法: m, n ∈ R + , m ⊕ n = mn
数乘: m ∈ R + , k ∈ R, k ⊗ m = m k
验证 R + 是 R 上的线性空间.
证 加法封闭,且(1)~(2)成立.
(3) m ⊕ θ = m ⇒ mθ = m ⇒ θ = 1 (4) m ⊕ (−m) = θ ⇒ m(−m) = 1 ⇒ (−m) = 1 m
=
S 2
时,称映射σ

S 上的变换. 1
例 2 S = { A = (ai j )n×n ai j ∈ R} (n ≥ 2) .
映射σ 1 :σ 1 ( A) = detA
(S → R)
变换σ 2 :σ 2 ( A) = (detA) In (S → S )
二、线性空间及其性质
1.线性空间:集合V 非空,给定数域 K ,若在V 中
矩阵论讲稿
讲稿编者: 张 凯 院
使用教材:《矩阵论》(第 2 版) 西北工业大学出版社 程云鹏 等编
辅助教材:《矩阵论导教导学导考》 《矩阵论典型题解析及自测试题》 西北工业大学出版社 张凯院 等编
课时分配:第一章 17 学时 第二章 5 学时 第三章 8 学时
第四章 8 学时 第五章 8 学时 第六章 8 学时
运算方式 2 加法: α ⊕ β = (ξ1 + η1 , ξ2 + η2 + ξ1η1 )
数乘:
k

=
(kξ1 ,
kξ 2
+
1 k(k 2

1)ξ
2 1
)
可以验证 R 2 (+ ⋅) 与 R2 (⊕ o) 都是 R 上的线性空间.
[注]
在 R2 (⊕ o) 中,
θ = (0,0) ,
−α
=
(−ξ1 , − ξ2
+
x 2

+
x 2
=
x 2

=
x 2
第一章 线性空间与线性变换(第 1 节)
4
例 6 在线性空间V 中,下列结论成立.
0x = θ :1x + 0x = (1 + 0)x = 1x ⇒ 0x = θ
kθ = θ : kx + kθ = k( x + θ ) = kx ⇒ kθ = θ
(−1)x = (− x) : (−1)x = (−1)x + [ x + (− x)] = [(−1)x + 1x] + (− x) = (− x)

S 2
∀b ∈
S 2

b∈
S 1
,
即S 2

S 1
交:
S 1
I
S 2
=
{a
a

S 1

a∈
S2 }
并:
S 1
U
S 2
=
{a
a

S 1

a

S 2
}
和: S 1
+
S 2
=
{a
=
a 1
+
a 2
a 1

S 1
,
a 2

S 2
}
例1
S 1
=
{A
=
a 11
a21
0
a
22
ai j ∈ R}
S 2
=
{A
=
a11
+
ξ
2 1
)

Th1 线性空间V 中的零元素唯一,负元素也唯一.

设θ 与θ 都是V
1
2
的零元素,
则θ1
= θ1
+θ2
= θ2
+ θ1
= θ2

x 1

x
2
都是
x
的负元素,
则由 x +
x 1
=θ及x+
x 2

可得
x 1
=
x 1

=
x 1
+
(x
+
x2 )
=
(
x 1
+
x) +
x 2
=
(x +
x 1
)
第一章 线性空间与线性变换(第 1 节)
1
第一章 线性空间与线性变换
§1.1 线性空间
一、集合与映射
1.集合:能够作为整体看待的一堆东西.
列举法:
S
=
{
a 1
,
a
2
,
a 3
,L}
性质法: S = { a a 所具有的性质 }
相等
(
S 1
=
S2 ) :指下面二式同时成立
∀a

S 1

a

S2 ,
即S 1
例如:实数域 R ,复数域 C ,有理数域 Q ,等等.
3.映射:设集合
S与 1
S 2
,若对任意的 a

S1 ,按照法则σ
,对应唯一的
第一章 线性空间与线性变换(第 1 节)
2
b ∈ S2 , 记作σ (a) = b.
称σ
为由
S到 1
S 2
的映射;称 b 为 a
的象,
a 为 b 的象源.
变换:当
S 1
相关文档
最新文档