《矩阵分析》
- 格式:ppt
- 大小:582.50 KB
- 文档页数:22
矩阵分析方法及应用论文矩阵分析方法是一种应用矩阵论和线性代数的数学工具,用于研究和解决与矩阵相关的问题。
矩阵可以用于描述线性变换、矢量空间和方程组等数学对象。
矩阵分析方法可以应用于多个领域,包括数学、物理、工程、计算机科学等。
在以下回答中,我将简要介绍矩阵分析方法的基本原理和一些应用,并提供一些相关论文的例子。
首先,让我们来了解一下矩阵分析的基本原理。
矩阵是一个由数值排列成的矩形数组,可以表示为一个m×n的矩阵,其中m表示行数,n表示列数。
矩阵的元素可以是实数或复数。
通过矩阵分析,我们可以研究矩阵的性质、运算规则和应用。
矩阵乘法是矩阵分析中最基本的操作之一。
当两个矩阵相乘时,第一个矩阵的列数必须等于第二个矩阵的行数。
矩阵乘法的结果是一个新的矩阵,其行数等于第一个矩阵的行数,列数等于第二个矩阵的列数。
矩阵乘法可以表示线性变换和矢量的线性组合等概念。
另一个重要的矩阵分析方法是特征值和特征向量的计算。
矩阵的特征值是矩阵与一个非零向量之间的一个简单乘法关系。
特征向量是与特征值对应的非零向量。
特征值和特征向量在物理、工程和计算机科学等领域中有广泛的应用,例如图像处理、机器学习和数据压缩等。
矩阵分析方法在多个领域有着广泛的应用。
下面是一些矩阵分析方法的应用领域及相应的论文例子:1. 图像处理:矩阵分析方法在图像处理中被广泛应用,例如图像压缩和恢复。
论文例子:《基于矩阵分解的图像压缩算法研究》、《基于矩阵分析方法的图像恢复技术研究》。
2. 数据处理:矩阵分析方法在数据挖掘和机器学习中起着重要作用,例如矩阵分解和矩阵推荐系统。
论文例子:《基于矩阵分解的矩阵推荐系统研究》、《基于矩阵分析的数据挖掘技术研究》。
3. 信号处理:矩阵分析方法在信号处理中具有广泛的应用,例如语音信号处理和音频编码。
论文例子:《基于矩阵分析方法的语音信号处理技术研究》、《基于矩阵分解的音频编码算法研究》。
4. 控制系统:矩阵分析方法在控制系统设计和分析中具有重要作用,例如状态空间表示和线性二次型控制器设计。
矩阵分析的应用
1、商品细分:商品细分矩阵分析是一种从市场上容易得到的数据,根据客户的不同需求,确定不同的属性,并将属性进行技术分析,从而得出市场消费者对产品的需求以及品牌的相对优势,从而帮助商家分析出满足客户需求的产品细分结构。
2、客户关系管理:矩阵分析可以帮助企业分析其客户的需求特点和关系,根据客户的不同行业、地理位置、企业规模等特点来确定客户群体,从而制定科学的客户关系管理策略,提高企业的客户关系管理水平。
3、绩效考核:矩阵分析的强大分析功能可以帮助企业分析销售团队的绩效,研究其团队绩效评估指标,比如业绩贡献、潜在客户开发情况、拜访状况等,从而实现企业员工绩效考核的客观、准确、合理的目标管理。
;。
P25⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=∴+-=-=+-=⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-==-=-=∈∈∀+=-+-=--+=+-+=+∈∀11-11-11-00011-11-11-000),,(),,()()(0110101111011011100110011101)()3(0100,0010,10012)()()()()(,)()()()()(,)1.(1421121121121121122121121221121121121111211211212211212121212121该基下的对应的矩阵为同理:变换的像分别求上一组基的线性以取这样的一组基这是一个三维空间,可可以写为)对于空间(的线性变换是根据定义可知,设设E E E E E E T E E E T E E E T E E E T E E E a a a a W W W T X T B X X B BX X B X T FW X X T X T B X X B B X X B B X B X X B X B BX X X X B X X T WX X T T T T TT TT TTT T T T λλλλλλ()()()()()()()()()()()()()()123123123123-1123123123123123123123123-1123-1123115.,,,,,,,,101110-121,,=,,,,,=,,,,,,,,,,,,,,=,T A T B A P P T T P T P AP P AP B P APηηηηηηεεεεεεεεεηηηηηηεεεεεεηηηηηηηηηεεεεεεηη==⎛⎫⎪= ⎪ ⎪⎝⎭=⎡⎤⎣⎦=⎡⎤⎣⎦===解:由题意知:其中,设则则由()()()23-1123123-11-1,=,,,,-110100010100010=100010=110010=1101-1100110100110101010101001110110110101-12111P P B P AP ηηηηεεε-⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪=== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭得到-111132⎛⎫ ⎪ ⎪ ⎪⎝⎭1.16(1)证明:()()()()()()()221223131212122T f t T f t x x x x t t x x t t +=+++++++⎡⎤⎡⎤⎣⎦⎣⎦ Q ()()()22123231312T x x t x t x x x x t x x t ⎡⎤++=+++++⎣⎦()()2123011,,1011,,110Tx x x t t ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦∴()()22121213112232T f t f t T x x t x t x x t x t ⎡⎤+=+++++⎡⎤⎣⎦⎣⎦()()2212123122T x x t t x t t ⎡⎤=++++⎣⎦()()221231212,,2,,TT x x x t t t t ⎡⎤=++⎢⎥⎣⎦()()221231212011,,1012,,110Tx x x t t t t ⎡⎤⎢⎥=++⎢⎥⎢⎥⎣⎦()()()()()2223131212122x x x x t t x x t t =+++++++∴()()12T f t f t +=⎡⎤⎣⎦()()12T f t T f t +⎡⎤⎡⎤⎣⎦⎣⎦ ()()2123T f t T x x t x t λλλλ=++⎡⎤⎣⎦()()2123,,,,T T x x x t t λλλ⎡⎤=⎢⎥⎣⎦()()2123011,,101,,110Tx x x t t λλλ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦()()()2231312x x x x t x x t λλλ=+++++()T f t λ=⎡⎤⎣⎦ ∴T 是[]3F x 的线性变换 (2)解: ()()2123T f t T xx tx t=++⎡⎤⎣⎦ ()()()21231x T x T t x T t =++()()()()2212311T f t x t t x t x t =+++++⎡⎤⎣⎦∴()21T t t =+;()21T t t =+;()21T t t =+∴()()220111,,1011,,110T t t t t ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦∴T 在基21,,t t 下的矩阵A 为011101110⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(3)解:()()211112111E A λλλλλλ---=--=-+--1232;1λλλ===-()112=1,1,1Tλξ=时,可以求得特征向量()()2323==1,1,0=1,0,1TTλλξξ=---1时,可以求得特征向量,故111=110101P ⎡⎤⎢⎥-⎢⎥⎢⎥-⎣⎦()()21231,,t t P ∂∂∂=令,,()()2221111,,1101011,1,1t t t t t t ⎡⎤⎢⎥=-⎢⎥⎢⎥-⎣⎦=++--则T 在基1∂=21t t ++,2∂=1t -,3∂=21t -下的矩阵为对角矩阵.P45第二章 内积空间练习题1.解:(1)Q ()11221x y x y αβ,=++,∴()11221x y x y λαβλλ,=++。
第1章 线性空间和线性变换(详解)1-1 证:用iiE 表示n 阶矩阵中除第i 行,第i 列的元素为1外,其余元素全为0的矩阵.用ij E (,1,2,,1)i j i n <=-表示n 阶矩阵中除第i 行,第j 列元素与第j 行第i 列元素为1外,其余元素全为0的矩阵.显然,iiE ,ijE 都是对称矩阵,iiE 有(1)2n n -个.不难证明iiE ,ijE 是线性无关的,且任何一个对称矩阵都可用这n+(1)2n n -=(1)2n n +个矩阵线性表示,此即对称矩阵组成(1)2n n +维线性空间.同样可证所有n 阶反对称矩阵组成的线性空间的维数为(1)2n n -.评注:欲证一个集合在加法与数乘两种运算下是一个(1)2n n +维线性空间,只需找出(1)2n n +个向量线性无关,并且集合中任何一个向量都可以用这(1)2n n +个向量线性表示即可.1-2解: 11223344x x x x ααααα=+++令解出1234,,,x x x x 即可.1-3 解:方法一 设11223344x x x x =+++A EE E E即123412111111100311100000x x x x ⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤=+++⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦故12341231211203x x x x x x x x x x +++++⎡⎤⎡⎤=⎢⎥⎢⎥+⎣⎦⎣⎦于是12341231,2x x x x x x x +++=++=1210,3x x x +==解之得12343,3,2,1x x x x ==-==-即A 在1234,,,E E E E 下的坐标为(3,3,2,1)T--.方法二 应用同构的概念,22R ⨯是一个四维空间,并且可将矩阵A 看做(1,2,0,3)T,1234,,,E E E E 可看做(1,1,1,1),(1,1,1,0),(1,1,0,0),(1,0,0,0)T T T T .于是有1111110003111020100311000001021000300011⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥→⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦因此A 在1234,,,E E E E 下的坐标为(3,3,2,1)T--.1-4 解:证:设112233440k k k k αααα+++=即1234123412313412411111110110110110k k k k k k k k k k k k k k k k k ⎡⎤⎡⎤⎡⎤⎡⎤+++⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦+++++⎡⎤==⎢⎥++++⎣⎦于是12341230,0k k k k k k k +++=++= 1341240,0k k k k k k ++=++=解之得12340k k k k ====故1234,,,αααα线性无关. 设123412341231341241111111011011011a b x x x x c d x x x x x x x x x x x x x ⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤=+++⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦+++++⎡⎤=⎢⎥++++⎣⎦于是12341230,0x x x x x x x +++=++= 1341240,0x x x x x x ++=++=解之得122,x b c d a x a c =++-=-34,x a d x a b =-=-1234,,,x x x x 即为所求坐标.1-5 解:方法一 (用线性空间理论计算)32312233410()121,,,021,1,(1),(1)p x x x x x y y x x x y y ⎡⎤⎢⎥⎢⎥⎡⎤=+=⎣⎦⎢⎥⎢⎥⎣⎦⎡⎤⎢⎥⎢⎥⎡⎤=---⎣⎦⎢⎥⎢⎥⎣⎦又由于23231,1,(1),(1)111101231,,,00130001x x x x x x ⎡⎤---⎣⎦⎡⎤⎢⎥-⎢⎥⎡⎤=⎣⎦⎢⎥-⎢⎥⎣⎦于是()p x 在基231,1,(1),(1)x x x ---下的坐标为11234111113012306001306000122y y y y -⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦方法二 将3()12p x x =+根据幂级数公式按1x -展开可得32323()12(1)(1)(1)(1)(1)(1)(1)2!3!36(1)6(1)2(1)p x x p p p p x x x x x x =+''''''=+-+-+-=+-+-+- 因此()p x 在基231,1,(1),(1)x x x ---下的坐标为[]3,6,6,2T.评注:按照向量坐标定义计算,第二种方法比第一种方法更简单一些.1-6 解:①设[][]12341234,,,,,,=ββββααααP将1234,,,αααα与1234,,,ββββ代入上式得20561001133611001121011010130011⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥=⎢⎥⎢⎥--⎢⎥⎢⎥-⎣⎦⎣⎦P 故过渡矩阵1100120561100133601101121001110131122223514221915223112822-⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥=⎢⎥⎢⎥--⎢⎥⎢⎥-⎣⎦⎣⎦⎡⎤---⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦P②设1212343410(,,,)10y y y y ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦ξββββ将1234,,,ββββ坐标代入上式后整理得11234792056181336027112111310130227y y y y -⎡⎤-⎢⎥⎢⎥⎡⎤⎡⎤⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦评注:只需将,i i αβ代入过渡矩阵的定义[][]12341234,,,,,,=ββββααααP 计算出P .1-7 解:因为12121212{,}{,}{,,,}span span span +=ααββααββ由于秩1212{,,,}3span =ααββ,且121,,ααβ是向量1212,,,ααββ的一个极大线性无关组,所以和空间的维数是3,基为121,,ααβ. 方法一 设1212{,}{,}span span ∈ξααββ,于是由交空间定义可知123411212111011030117k k k k -⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥+++=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦解之得1222122,4,3(k l k l l l l =-==-为任意数)于是11222[5,2,3,4]T k k l =+=-ξαα(很显然1122l l ββ=+ξ)所以交空间的维数为1,基为[5,2,3,4]T-.方法二 不难知12121212{,}{,},{,}{,}span span span span ''==ααααββββ其中2213[2,2,0,1],[,2,1,0]3TT ''=--=-αβ.又12{,}span 'αα也是线性方程组 13423422x x x x x x =-⎧⎨=-⎩ 的解空间.12{,}span 'ββ是线性方程组13423413232x x x x x x ⎧=-+⎪⎨⎪=-⎩ 的解空间,所以所求的交空间就是线性方程组1342341342342213232x x x x x x x x x x x x =-⎧⎪=-⎪⎪⎨=-+⎪⎪=-⎪⎩ 的解空间,容易求出其基础解系为[5,2,3,4]T-,所以交空间的维数为1,基为[5,2,3,4]T -.评注:本题有几个知识点是很重要的.12(1){,,,}n span ααα的基底就是12,,,nααα的极大线性无关组.维数等于秩12{,,,}n ααα.1212(2){,}{,}span span +ααββ1212{,,,}span =ααββ.(3)方法一的思路,求交1212{,}{,}span span ααββ就是求向量ξ,既可由12,αα线性表示,又可由12,ββ线性表示的那部分向量.(4)方法二是借用“两个齐次线性方程组解空间的交空间就是联立方程组的解空间”,将本题已知条件改造为齐次线性方程组来求解.1-8解:(1):解出方程组1234123420510640x x x x x x x x ---=⎧⎨---=⎩(Ⅰ)的基础解系,即是1V 的基, 解出方程组123420x x x x -++=(Ⅱ)的基础解系,即是2V 的基;(2): 解出方程组1234123412342051064020x x x x x x x x x x x x ---=⎧⎪---=⎨⎪-++=⎩的基础解系,即为12V V ⋂的基;(3):设{}{}1121,,,,,k l V span V span ααββ==,则11,,,,,k l ααββ的极大无关组即是12V V +的基. 1-9解:仿上题解.1-10解: 仿上题解.1-11 证:设210121()()()0k k l l l l --++++=ξξξξA AA①用1k -A从左侧成①式两端,由()0k=ξA 可得10()0k l -=ξA因为1()0k -≠ξA,所以00l =,代入①可得21121()()()0k k l l l --+++=ξξξA AA②用2k -A从左侧乘②式两端,由()0k=ξA可得00l =,继续下去,可得210k l l -===,于是21,(),(),,()k -ξξξξA AA线性无关.1-12 解:由1-11可知,n 个向量210,(),(),,()n -≠ξξξξAAA线性无关,它是V 的一个基.又由21212121[,(),(),,()][(),(),,()][(),(),,(),0]000010000100[,(),(),,()]0000010n n n n n n----⨯==⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ξξξξξξξξξξξξξξA A A AA A A A AAA AA 所以A 在21,(),(),,()n -ξξξξA AA下矩阵表示为n 阶矩阵000100001000000010⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦评注:n 维线性空间V 中任何一组n 个线性无关的向量组都可以构成V 的一个基,因此21,(),(),,()n -ξξξξA AA是V 的一个基.1-13证: 设()()()111,,,,,,,,,,,r s m r s A A ξξξββααα==设11,,,,,,r r s ξξξξξ是的极大无关组,则可以证明11,,,,,,r r s ααααα是的极大无关组.1-14 解:(1)由题意知123123[,,][,,]=ααααααA A123123111[,,][,,]011001⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦βββααα设A在基123,,βββ下的矩阵表示是B ,则11111123111011103011001215001244346238--⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎡⎤⎢⎥=---⎢⎥⎢⎥⎣⎦B P AP (2)由于0A ≠,故0=AX 只有零解,所以A的核是零空间.由维数定理可知A 的值域是线性空间3R .1-15解:已知()()2323,,,,A αααααα=11A(1) 求得式()()2323,,,,P εεεααα=11中的过渡矩阵P ,则1B P AP -=即为所求;(2)仿教材例1.5.1.(见<矩阵分析>史荣昌编著.北京理工大学出版社.) 1-16解:设()23,,A ααα=1,则{}23(),,;()R A span N A ααα=1就是齐次方程组0Ax = 的解空间. 1-17证:由矩阵的乘法定义知AB BA 与的主对角线上元素相等,故知AB BA 与的迹相等;再由1-18 题可证. 1-18证:对k 用数学归纳法证。
矩阵分析实验报告学院:电气学院专业:控制工程姓名:XXXXXXXX学号:211208010001矩阵分析实验报告实验题目利用幂法求矩阵的谱半径实验目的与要求1、 熟悉matlab 矩阵实验室的功能和作用;2、 利用幂法求矩阵的谱半径;3、 会用matlab 对矩阵分析运算。
实验原理理念谱半径定义:设n nA C⨯∈,1λ,2λ,3λ, ,j λ, n λ是A 的n 个特征值,称()max ||j jA ρλ=为关于A 的谱半径。
关于矩阵的谱半径有如下结论:设n nA C⨯∈,则(1)[]()()kkA A ρρ=;(2)22()()()H H A A AA A ρρ==。
由于谱半径就是矩阵的主特征值,所以实验换为求矩阵的主特征值。
算法介绍定义:如果1λ是矩阵A 的特征值,并且其绝对值比A 的任何其他特征值的绝对值大,则称它为主特征值。
相应于主特征值的特征向量1V 称为主特征向量。
定义:如果特征向量中最大值的绝对值等于单位值(例如最大绝对值为1),则称其为是归一化的。
通过形成新的向量'12=c n V (1/)[v v v ],其中c=v 且1max {},j i n i ≤≤=v v 可将特征向量 '12n [v v v ]进行归一化。
设矩阵A 有一主特征值λ,而且对应于λ有唯一的归一化特征向量V 。
通过下面这个称为幂法(power method )的迭代过程可求出特征对λ,V ,从下列向量开始:[]'0=111X (1)用下面递归公式递归地生成序列{}k X :k k Y AX =k+111k k X Y c +=(2)其中1k c +是k Y 绝对值最大的分量。
序列{}k X 和{}k c 将分别收敛到V 和λ:1lim k X V =和lim k c λ= (3)注:如果0X 是一个特征向量且0X V ≠,则必须选择其他的初始向量。
幂法定理:设n ×n 矩阵A 有n 个不同的特征值λ1,λ2,···,,λn ,而且它们按绝对值大小排列,即:123n λλλλ≥≥≥⋅⋅⋅≥ (4)如果选择适当的X 0,则通过下列递推公式可生成序列{[()()()]}12k kk k n X x x x '=⋅⋅⋅和{}k c : k k Y AX = (5)和:111k k k X Y c ++=(6)其中: ()1k k j c x +=且{}()()1max k k j i i nx x ≤≤=(7)这两个序列分别收敛到特征向量V 1和特征值λ1。
矩阵分析在控制系统中的应用摘要:详细综述了LMI 在控制系统中的发展现状和应用,主要涉及了不确定系统的鲁棒性能和鲁棒稳定性、不确定系统的鲁棒控制器设计、LMI 在时滞系统中的应用及存在的问题、不确定系统的鲁棒滤波应用状况、不确定系统的模型验证应用等,并分析了基于LMI 方法的变结构控制、极点配置、模糊控制等其它相关内容。
给出了上述控制问题的LMI 描述及相关求解方法,最后并指出了LMI 进一步的应用研究方向。
主题词: 线性; 矩阵; 控制系统; 控制器1 引言在过去的10 余年内,由于LMI 的优良性质和数学的规范以及解法的突破,使其在控制系统分析和设计方面得到了广泛的重视和应用。
研究者发现许多控制问题均可描述为LMI 问题[1~4 ] ,并呈现继续增长的趋势。
本文对LMI 在控制系统中的发展和现状进行综述,着重讨论LMI 在不确定控制系统中的应用研究成果、现状以及发展。
2 线性矩阵不等式LMI 一般形式为F ( x) ≡F0 + Σmi =1xi F i > 0 (1)其中x ∈Rm ———变量; F i = F Ti ∈Rn×n 是给定的。
显然式(1) 表明矩阵F( x) 是正定的。
式(1) 的另一个含义是集合{ x/ F( x) > 0} 是凸的。
LMI 问题可描述为:给定F( x) > 0 ,找到x,使得f ( x) > 0 ,或证明LMI F( x) 是不可解的。
动态系统分析的LMI 方法可以追溯到100 多年以前。
1890 年Lyapunov 在出版他的被称为Lyapunov 理论的著作中,提出微分方程Ûx( t) = Ax ( t) (2)稳定,当且仅当存在对称正定矩阵P = P T > 0 ,使得下面的不等式成立A T P + PA < 0 (3)同时Lyapunov 也指出这样的LMI 可以精确求解。
20 世纪40 年代,前苏联科学家Lur’e、Postnikov 及其它学者将Lyapunov 方法应用于控制工程中的一些典型的问题,尤其是当执行机构具有非线性时的系统稳定性,虽然他们没有形成精确的矩阵不等式,但是所提出的稳定性准则具有LMI的雏形。
《矩阵分析》(第3版)史荣昌,魏丰.第一章课后习题答案第1章 线性空间和线性变换(详解)1-1 证:用ii E 表示n 阶矩阵中除第i 行,第i 列的元素为1外,其余元素全为0的矩阵.用ij E (,1,2,,1)i j i n <=-L 表示n 阶矩阵中除第i 行,第j 列元素与第j 行第i 列元素为1外,其余元素全为0的矩阵.显然,ii E ,ij E 都是对称矩阵,ii E 有(1)2n n -个.不难证明ii E ,ij E 是线性无关的,且任何一个对称矩阵都可用这n+(1)2n n -=(1)2n n +个矩阵线性表示,此即对称矩阵组成(1)2n n +维线性空间.同样可证所有n 阶反对称矩阵组成的线性空间的维数为(1)2n n -.评注:欲证一个集合在加法与数乘两种运算下是一个(1)2n n +维线性空间,只需找出(1)2n n +个向量线性无关,并且集合中任何一个向量都可以用这(1)2n n +个向量线性表示即可.1-2解: 11223344x x x x ααααα=+++令 解出1234,,,x x x x 即可.1-3 解:方法一 设11223344x x x x =+++A E E E E即123412111111100311100000x x x x ⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤=+++⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦ 故12341231211203x x x x x x x x x x +++++⎡⎤⎡⎤=⎢⎥⎢⎥+⎣⎦⎣⎦于是12341231,2x x x x x x x +++=++=1210,3x x x +==解之得12343,3,2,1x x x x ==-==-即A 在1234,,,E E E E 下的坐标为(3,3,2,1)T --.方法二 应用同构的概念,22R ⨯是一个四维空间,并且可将矩阵A 看做(1,2,0,3)T ,1234,,,E E E E 可看做(1,1,1,1),(1,1,1,0),(1,1,0,0),(1,0,0,0)T T T T .于是有1111110003111020100311000001021000300011⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥→⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦因此A 在1234,,,E E E E 下的坐标为(3,3,2,1)T --.1-4 解:证:设112233440k k k k αααα+++=即1234123412313412411111110110110110k k k k k k k k k k k k k k k k k ⎡⎤⎡⎤⎡⎤⎡⎤+++⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦+++++⎡⎤==⎢⎥++++⎣⎦于是12341230,0k k k k k k k +++=++=1341240,0k k k k k k ++=++=解之得12340k k k k ====故1234,,,αααα线性无关. 设123412341231341241111111011011011a b x x x x c d x x x x x x x x x x x x x ⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤=+++⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦+++++⎡⎤=⎢⎥++++⎣⎦于是12341230,0x x x x x x x +++=++= 1341240,0x x x x x x ++=++=解之得122,x b c d a x a c =++-=-34,x a d x a b =-=-1234,,,x x x x 即为所求坐标.1-5 解:方法一 (用线性空间理论计算)32312233410()121,,,021,1,(1),(1)p x x x x x y y x x x y y ⎡⎤⎢⎥⎢⎥⎡⎤=+=⎣⎦⎢⎥⎢⎥⎣⎦⎡⎤⎢⎥⎢⎥⎡⎤=---⎣⎦⎢⎥⎢⎥⎣⎦又由于23231,1,(1),(1)111101231,,,00130001x x x x x x ⎡⎤---⎣⎦⎡⎤⎢⎥-⎢⎥⎡⎤=⎣⎦⎢⎥-⎢⎥⎣⎦于是()p x 在基231,1,(1),(1)x x x ---下的坐标为11234111113012306001306000122y y y y -⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦方法二 将3()12p x x =+根据幂级数公式按1x -展开可得32323()12(1)(1)(1)(1)(1)(1)(1)2!3!36(1)6(1)2(1)p x x p p p p x x x x x x =+''''''=+-+-+-=+-+-+- 因此()p x 在基231,1,(1),(1)x x x ---下的坐标为[]3,6,6,2T.评注:按照向量坐标定义计算,第二种方法比第一种方法更简单一些. 1-6 解:①设[][]12341234,,,,,,=ββββααααP将1234,,,αααα与1234,,,ββββ代入上式得20561001133611001121011010130011⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥=⎢⎥⎢⎥--⎢⎥⎢⎥-⎣⎦⎣⎦P 故过渡矩阵1100120561100133601101121001110131122223514221915223112822-⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥=⎢⎥⎢⎥--⎢⎥⎢⎥-⎣⎦⎣⎦⎡⎤---⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦P②设1212343410(,,,)10y y y y ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦ξββββ将1234,,,ββββ坐标代入上式后整理得11234792056181336027112111310130227y y y y -⎡⎤-⎢⎥⎢⎥⎡⎤⎡⎤⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦评注:只需将,i i αβ代入过渡矩阵的定义[][]12341234,,,,,,=ββββααααP计算出P .1-7 解:因为12121212{,}{,}{,,,}span span span +=ααββααββ由于秩1212{,,,}3span =ααββ,且121,,ααβ是向量1212,,,ααββ的一个极大线性无关组,所以和空间的维数是3,基为121,,ααβ.方法一 设1212{,}{,}span span ∈ξααββI ,于是由交空间定义可知123411212111011030117k k k k -⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥+++=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦解之得1222122,4,3(k l k l l l l =-==-为任意数)于是11222[5,2,3,4]T k k l =+=-ξαα(很显然1122l l ββ=+ξ)所以交空间的维数为1,基为[5,2,3,4]T -. 方法二 不难知12121212{,}{,},{,}{,}span span span span ''==ααααββββ其中2213[2,2,0,1],[,2,1,0]3TT ''=--=-αβ.又12{,}span 'αα也是线性方程组13423422x x x x x x =-⎧⎨=-⎩ 的解空间.12{,}span 'ββ是线性方程组13423413232x x x x x x ⎧=-+⎪⎨⎪=-⎩ 的解空间,所以所求的交空间就是线性方程组1342341342342213232x x x x x x x x x x x x =-⎧⎪=-⎪⎪⎨=-+⎪⎪=-⎪⎩ 的解空间,容易求出其基础解系为[5,2,3,4]T -,所以交空间的维数为1,基为[5,2,3,4]T -.评注:本题有几个知识点是很重要的.12(1){,,,}n span αααL 的基底就是12,,,n αααL 的极大线性无关组.维数等于秩12{,,,}n αααL .1212(2){,}{,}span span +ααββ1212{,,,}span =ααββ.(3)方法一的思路,求交1212{,}{,}span span ααββI 就是求向量ξ,既可由12,αα线性表示,又可由12,ββ线性表示的那部分向量.(4)方法二是借用“两个齐次线性方程组解空间的交空间就是联立方程组的解空间”,将本题已知条件改造为齐次线性方程组来求解.1-8解:(1):解出方程组1234123420510640x x x x x x x x ---=⎧⎨---=⎩(Ⅰ)的基础解系,即是1V 的基, 解出方程组123420x x x x -++=(Ⅱ)的基础解系,即是2V 的基; (2): 解出方程组1234123412342051064020x x x x x x x x x x x x ---=⎧⎪---=⎨⎪-++=⎩的基础解系,即为12V V ⋂的基;(3):设{}{}1121,,,,,k l V span V span ααββ==L L ,则11,,,,,k l ααββL L 的极大无关组即是12V V +的基. 1-9解:仿上题解.1-10解: 仿上题解.1-11 证:设210121()()()0k k l l l l --++++=ξξξξL A AA①用1k -A从左侧成①式两端,由()0k=ξA可得10()0k l -=ξA因为1()0k -≠ξA,所以00l =,代入①可得21121()()()0k k l l l --+++=ξξξL A A A②用2k -A从左侧乘②式两端,由()0k=ξA可得00l =,继续下去,可得210k l l -===L ,于是21,(),(),,()k -ξξξξL A AA 线性无关.1-12 解:由1-11可知,n 个向量210,(),(),,()n -≠ξξξξL A AA线性无关,它是V 的一个基.又由21212121[,(),(),,()][(),(),,()][(),(),,(),0]000010000100[,(),(),,()]00000010n n n n n n----⨯==⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ξξξξξξξξξξξξξξL L L L L L L M M M M L LA A A AA A A A AAA A A 所以A在21,(),(),,()n -ξξξξL A AA下矩阵表示为n 阶矩阵00001000010000000010⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦L L L M M M M L L评注:n 维线性空间V 中任何一组n 个线性无关的向量组都可以构成V 的一个基,因此21,(),(),,()n -ξξξξL A A A是V 的一个基.1-13证: 设()()()111,,,,,,,,,,,r s m r s A A ξξξββααα==L L L L L 设11,,,,,,r r s ξξξξξL L L 是的极大无关组,则可以证明11,,,,,,r r s αααααL L L 是的极大无关组. 1-14 解:(1)由题意知123123[,,][,,]=ααααααA A123123111[,,][,,]011001⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦βββααα 设A在基123,,βββ下的矩阵表示是B ,则11111123111011103011001215001244346238--⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎡⎤⎢⎥=---⎢⎥⎢⎥⎣⎦B P AP (2)由于0A ≠,故0=AX 只有零解,所以A的核是零空间.由维数定理可知A的值域是线性空间3R .1-15解:已知()()2323,,,,A αααααα=11A(1) 求得式()()2323,,,,P εεεααα=11中的过渡矩阵P ,则1B P AP -=即为所求; (2)仿教材例1.5.1.(见<矩阵分析>史荣昌编著.北京理工大学出版社.) 1-16解:设()23,,A ααα=1,则{}23(),,;()R A span N A ααα=1就是齐次方程组0Ax = 的解空间. 1-17证:由矩阵的乘法定义知AB BA 与的主对角线上元素相等,故知AB BA 与的迹相等;再由1-18 题可证. 1-18证:对k 用数学归纳法证。