矩阵方程xa=b例题解法

合集下载

【精选】数学实验一矩阵运算与Matlab命令24

【精选】数学实验一矩阵运算与Matlab命令24
B1=[b11 b12 b13 b14 ;b21 b22 b23 b24; b31 b32 b33 b34]
运行
17
矩阵的运算(矩阵的加减、数乘、乘积)
C=A1+B1 D=A1-B1 syms c, cA=c*A1 A2=A1(:,1:3), B1 G=A2*B1
18
矩阵的运算(矩阵的加减、数乘、乘积)
求解方程组Ax=b x=A\b 若A为可逆方阵, 输出原方程的解x; 若A为nxm(n>m)阵, 且A’A可逆,输出
原方程的最小二乘解x.
21
矩阵的运算(求解线性方程组)
求矩阵方程:
设A、B满足关系式:AB=2B+A,求B。 其中A=[3 0 1; 1 1 0; 0 1 4]。
取出A的1、3行和1、3列的交叉处元素 构成新矩阵A1
程序
A=[1 0 1 1 2;0 1 -1 2 3;

3 0 5 1 0;2 3 1 2 1],
vr=[1, 3];vc=[1, 3];
A1=A(vr, vc)
观察结果
26
分块矩阵(矩阵的标识)
将A分为四块,并把它们赋值到矩阵B 中,观察运行后的结果。
3
2
2

35 20 60 45
, B 10
15
50
40

20 12 45 20
将 表 格 写 成 矩 阵 形 式
6
计算
输入下面Matlab指令 A=[4 2 3;1 3 2;1 3 3;3 2 2], B=[35 20 60 45;10 15 50 40;20

3 0 5 1 0;2 3 1 2 1]

初等变换求解矩阵方程

初等变换求解矩阵方程

初等变换求解矩阵方程作者:张馨元来源:《读写算》2013年第45期【摘要】本文给出了求解矩阵方程AX=B,XA=B以及AXB=C的初等变换法.【关键词】初等变换初等矩阵单位矩阵矩阵方程1.引言矩阵方程是指含有未知矩阵的矩阵等式.本文主要研究了三种典型矩阵方程,即AX=B、XA=B和AXB=C的求解.当矩阵A,B可逆时,一般上述三种矩阵的计算结果是:X=A-1B、X=BA-1和X=A-1CB-1.也就是要计算矩阵方程需要先求相应的逆矩阵,然后再做矩阵的乘法运算.显然这样比较麻烦,本文运用初等变换的相关理论,给出这三种矩阵的解法.引理1 对矩阵A施行一次初等行变换就相当于在A的左边乘上相应的初等矩阵;对A施行一次初等列变换就相当于在A的右边乘上相应的初等矩阵.引理 2 设A为可逆矩阵,则A可单用初等行变换化为单位矩阵E,也可以单用初等列变换化为单位矩阵E.2.AX=B定理1 设矩阵方程AX=B中的矩阵A可逆,则此矩阵方程的解可以通过下列矩阵的初等变换得到:证明由可逆,根据引理1和2,存在初等矩阵,使,将这些初等矩阵去左乘矩阵方程,得,即,由上可见,如果用一系列的初等行变换把A化为E,同时把这些相同的初等行变换施加在B上, B就化成了X,即矩阵方程的解.3.XA=B定理2 设矩阵方程XA=B中的矩阵A可逆,则此矩阵方程的解可以通过下列矩阵的初等变换得到:证明由A可逆,根据引理1和2,存在初等矩阵,使,将这些初等矩阵去右乘矩阵方程XA=B,得列,即,由上可见,如果用一系列的初等列变换把A化为E,同时把这些相同的初等列变换施加在B上, B就化成了X,即矩阵方程的解.4.AXB=C利用矩阵方程AX=B和XA=B的解法,综合可得矩阵方程AXB=C的解法.定理3 假设矩阵方程AXB=C中的矩阵A和B均可逆.则此矩阵方程的解可以通过下列矩阵的初等变换得到:1)2)3)4)5.定理的应用例 1 解矩阵方程.解:设, .故.例2 设且满足XA=B,求X解:故.例3 设,,,求解矩阵方程AXB=C.解:法一:故.法二:故.上例的两个方法是运用定理3的(1)和(2)而得的,我们也可以运用定理3的(3)和(4)来求解上述矩阵方程.参考文献[1]北京大学数学系几何与代数教研室代数小组.高等代数.高等教育出版社.1988年[2]同济大学数学教研室.线性代数(第五版).高等教育出版社,2007年。

矩阵的初等变换与线性方程组

矩阵的初等变换与线性方程组

第三章 矩阵的初等变换与线性方程组说明与要求:上一章已经介绍了求解线性方程组的克莱姆法则.虽然克莱姆法则在理论上具有重要的意义,但是利用它求解线性方程组,要受到一定的限制.首先,它要求线性方程组中方程的个数与未知量的个数相等,其次还要求方程组的系数行列式不等于零.即使方程组具备上述条件,在求解时,也需计算n +1个n 阶行列式.由此可见,应用克莱姆法则只能求解一些较为特殊的线性方程组且计算量较大.本章讨论一般的n 元线性方程组的求解问题.一般的线性方程组的形式为⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++mn mn m m n n n n b x a x a x a b x a x a x a b x a x a x a 22112222212111212111 (I)方程的个数m 与未知量的个数n 不一定相等,当m =n 时,系数行列式也有可能等于零.因此不能用克莱姆法则求解.对于线性方程组(I ),需要研究以下三个问题:(1)怎样判断线性方程组是否有解?即它有解的充分必要条件是什么? (2)方程组有解时,它究竟有多少个解及如何去求解? (3)当方程组的解不唯一时,解与解之间的关系如何? 目的与要求:掌握矩阵的初等变换,能用初等变换化矩阵为行阶梯形、行最简形和标准型。

理解矩阵的秩概念、掌握用初等变换求矩阵的秩。

了解初等矩阵的概念,掌握用初等变换求逆矩阵的方法。

掌握用初等变换求解线性方程组。

本章重点:矩阵的初等变换;解线性方程组;秩;线性方程组解的判定. 。

本章难点:秩;线性方程组解的判定.§3.1 矩阵的初等变换在本章的§2.3节中给出了矩阵可逆的充分必要条件,并同时给出了求逆矩阵的一种方法——伴随矩阵法.但是利用伴随矩阵法求逆矩阵,当矩阵的阶数较高时计算量是很大的.这一节将介绍求逆矩阵的另一种方法——初等变换法.为此我们先介绍初等矩阵的概念,并建立矩阵的初等变换与矩阵乘法的联系.一. 初等变换定义下面三种变换称为矩阵的初等行变换:1.互换两行(记);2.以数乘以某一行(记);3.把某一行的倍加到另一行上(记)。

三种典型矩阵方程的简单解法

三种典型矩阵方程的简单解法

即对矩阵 … 施行初等列变换 ,当把 A 变成 E 时 ,B
B 就变成 X 。(f ) 式提供了一个具体解矩阵方程 XA = B 的
简单方法 。 例2 解下列矩阵方程 。
2005 年 6 月第 3 期 三种典型矩阵方程的简单解法 3 (i) X - 1 3 1 2 4 0 1 2 3 - 1 3 3 A - 1 = 3 1 2 4 2 1 0 0 1 2 = (0 2 3) ; 1 2 4 0 1 2 3 0 0 5 0 1 2 1 2 0 1 2 3 - 1 3 ,B = ( 0 2 3) 。 1 4 7 0 1 2 0 0 - 8 1 0 - 1 0 1 0 0 0 1
… ,
X
可得解矩阵方程 AXB = C 的简单解法 例3 解下列矩阵方程 。
1 (i) 2 3 2 2 4 3 1 X 3 2 5 1 3 1 1 = 2 3 0 ;

- 4

- 3 1 (ii) X 4 7

3 2 5 8 3
于是有 X = ( - 4 - 3 3) 。
1 4 2 5 8 2 5 8 3 1 3 6 1 2 1 5 3 2 。 3 1 4 1 4 7 2 1 5 3 2 。 3 0 - 3 - 6 0 - 6 - 20 6 = 0 1 1 7 1 4 A 7
- 1 - 1 X = PL PL - 1 … P1- 1B 。 证毕
,再左乘 B 即得 X。 ,再右乘 B 即得 X。
- 1 - 1
若 XA = B ,则有 XAA - 1 = BA - 1 ,即 X = BA - 1 。于是
- 1
(1) (2)
又若 AXB = C ,则有 A AXBB
CB
A B

矩阵方程的解法

矩阵方程的解法

矩阵方程的解法本文首先介绍了行对称矩阵的定义及性质,利用矩阵的广义逆,奇异值分解,给出了矩阵方程AX=B有行对称解的充分必要条件及有解时通解的表达式;并给出了矩阵方程解集合中与给定矩阵的最佳逼近解的表达式。

最后利用奇异值分解给出了矩阵方程有行对称解的充分必要条件及有解时通解的表达式。

矩阵方程问题是指在满足一定条件的矩阵集合中求矩阵方程的解的问题。

不同的约束条件,不同的矩阵方程,就导致了不同的约束矩阵方程问题。

约束矩阵方程问题在结构设计,参数识别,主成分分析,勘测,遥感,生物学,电学,固体力学,结构动力学,分子光谱学,自动控制理论,振动理论,循环理论等领域都有重要应用。

约束矩阵方程问题的内容非常广泛、约束矩阵方程问题又分为线性约束矩阵方程问题和非线性约束矩阵方程问题、有关线性约束矩阵方程问题的研究成果相当丰富、其中最简单的矩阵方程AX = B是研究最透彻的一类问题、求解线性矩阵方程一般会遇到两种情况:一是当矩阵方程有解时,如何求它的解及最佳逼近;二是当矩阵方程无解时,如何求它的最小二乘解。

对于本文所研究的AX=B、这两类简单矩阵方程,国内外学者已经作了大量研究。

都在相应的文献中对其进行了大量的研究,解决了求此方程的一些约束解和最小二乘解的问题。

自从针对工程应用领域提出了行对称矩阵概念之后,这方面研究已经取得了一些成果,如对行对称矩阵的一些性质,行对称矩阵的QR分解。

本文先对行对称矩阵进行介绍,再将行对称矩阵与约束矩阵方程结合起来,先研究了矩阵方程AX=B有行对称实矩阵解的充要条件,有解时,用奇异值分解及广义逆求出解及最佳逼近。

再对矩阵方程有行对称实矩阵解的充要条件进行了研究,利用奇异值分解得出了有解时的充要条件及解的表达式。

设表示全体n*m阶实矩阵集合,rank(A)表示矩阵A的秩,表示次对角线上元素全为1,其余元素全为0的方阵,即=,显然有成立。

表示n阶正交矩阵全体。

本文要讨论以下问题:问题1 给定矩阵A,B,求实行对称方阵X,使得AX=B。

线性代数笔记13——Ax=b的通解

线性代数笔记13——Ax=b的通解

线性代数笔记13——Ax=b的通解 关于最简⾏阶梯矩阵和矩阵秩,可参考《》 召唤⼀个⽅程Ax = b: 3个⽅程4个变量,⽅程组有⽆数解,现在要关注的是b1b2b3之间满⾜什么条件时⽅程组有解,它的解是什么? 在这个例⼦中可以马上看出,b1+b2 = b3,⼀般的⽅法是消元法化简: 化简到这⼀步就可以确定主元是x1和x3。

通过最后⼀⾏可知,b3 – b2 - b1 = 0。

b1b2b3可以是任意数,所以只要满⾜b3 – b2 - b1 = 0,⽅程组就有解。

这样的组合很多,可以很容易找到⼀个特解: 现在我们知道了b中三个分量的关系,并且还知道只有当 b属于A的列空间时有解。

通过上⼀章的⽅法可知,列空间的基就是主元所在的列: 到此为⽌回答了第⼀个问题,什么样的b才能使Ax = b有解。

现在需要回答另⼀个问题,Ax = b的所有解是什么? 可以先找出⼀个特解,⽅法是令所有⾃由元为0,然后解出主元: 已经找到了⼀个特解,那么⽅程组的其它解,也就是通解是什么呢? 假设Ax= 0的零空间的任意向量是x n,Ax = b有⼀个特解x p,那么有: ⼆者相加: 所以⽅程组的通解是x n + x p。

对于⽅程组的某解x p来说,x p与零空间内任意向量之和仍为解。

现在看看零空间: 综合特解,得到Ax = b的通解: 矩阵的秩和主元个数相同。

如果A是⼀个m⾏n列的矩阵,其主元的个数⼀定⼩于m,并且也⼩于n。

如果A的每⼀列都有主元,那么A是满秩矩阵,没有⾃由元,如果此时有解,则解是唯⼀的,就是特解,即x = x p,此时不需要求解零空间,零空间只包含零向量。

作者:我是8位的本⽂以学习、研究和分享为主,如需转载,请联系本⼈,标明作者和出处,⾮商业⽤途!扫描⼆维码关注公众号“我是8位的”。

第三章%20%20矩阵的初等变换与线性方程组[1]

第三章%20%20矩阵的初等变换与线性方程组[1]
14 September 2009 河北科大理学院
第三章 矩阵的初等变换与线性方程组
15
第11讲 线性方程组的解
对于线性方程组Ax=b, R(A)=r. 不妨设 B = ( A, b ) 的行最简形为
⎛1 ⎜0 ⎜M ⎜0 ⎜0 ⎜0 ⎜M ⎜0 ⎝ 0 1 M 0 0 0 M 0 L 0 b11 L 0 b21 M M L 1 br 1 L0 0 L0 0 M M L0 0 L b1,n− r L b2,n− r 2 M M ⎟ L br ,n− r d r ⎟ L 0 d r +1 ⎟ L 0 0 ⎟ M M ⎟ 0 ⎟ L 0 ⎠
r (1) A ~ B 的充分必要条件是存在m阶可逆阵P, 使得 PA = B. c (2) A ~ B 的充分必要条件是存在n阶可逆阵Q, 使得 AQ = B.
(3) A ~ B 的充分必要条件是存在m阶可逆阵P 及n阶可逆阵Q, 使得 PAQ = B.
r 推论 方阵A可逆的充分必要条件是 A ~ E
14 September 2009 河北科大理学院
第三章 矩阵的初等变换与线性方程组
7
四 利用初等变换求逆矩阵及相关问题
A −1 ?
r ( A , E ) ~ ( E , A −1 )
⎛ 0 −2 1 ⎞ 例1 设 A = ⎜ 3 0 −2 ⎟ , 求 A−1 . ⎜ −2 3 0 ⎟ ⎝ ⎠
14 September 2009
相容 不相容 d1 ⎞ 特解 通解 d ⎟
(cii ≠ 0, i = 1, 2, L , n) (1) 当 d r +1 ≠ 0 时,方程组无解; (2) 当 d r +1 = 0 时,方程组有解; 且 r = n 时,有唯一解, r < n 时,有无限多个解.

矩阵方程求解方法

矩阵方程求解方法

矩阵方程求解方法本文所述的矩阵方程是指形如Ax=b的方程,其中A是一个mxn的矩阵,称为方程的系数矩阵。

x和b是mx1的矩阵。

特别的,当b=0时,这种方程又称为其次方程。

本文将讨论这种矩阵的有解条件和求解方法。

矩阵方程的有解条件为了解释矩阵方程的有解条件,我们首先要熟悉一些概念。

一个矩阵方程的增广矩阵是系数矩阵A和b并在一起构成的矩阵,记作(A,b)。

假定, ,则矩阵方程的增广矩阵就是矩阵的秩定义为其行向量中极大线性无关组中包含向量的个数,等价的说法是,矩阵的秩是r,则矩阵通过行列初等变换,变换成左上角是一个r阶单位矩阵,其他都是0的矩阵。

矩阵A的秩记作r(A),其中r是英文单词rank的缩写。

有了这两个基本概念,我们就可以准确描述矩阵方程的有解条件了:矩阵方程Ax=b的有解条件是矩阵A的秩等于增广矩阵(A,b)的秩,也就是r(A)=r(A,b)。

证明很简单,既然矩阵A的秩是r,那么肯定可以找到两个可逆的矩阵P,Q,满足--1)其中Ir表示r阶单位矩阵。

应用到原来的方程,可以得到:--2)我们把Q-1x当作一个未知的变量,PAQ当作系数,这就构成一个新的矩阵方程。

而这个矩阵方程的左侧系数除了前r行是有1的之外,其余行是0。

为了它有解,Pb的后m-r行必须也是0。

这样(A,b)的秩必然是r。

必须注意到Q-1是可逆的,因此以Q-1x为未知变量的方程有解意味着以x 为未知变量的原方程也是有解的。

矩阵方程的解对于矩阵方程Ax=b,如果满足r(A)=r(A,b),则矩阵方程是有解的。

为了求它的解,我们首先把矩阵方程通过行列初等变换变化成前文2)式的形式,代入1)式后得到:--3)其中Q-1x和Pb是一个列向量,我们可以把它们分割成rx1和(n-r)x1的两个矩阵,分别记作x’1和x’2,及b’1和b’2。

则很显然我们可以得到:--4)很显然,b’2必须为0,因为展开后b’2等于0 x’1 +0 x’2 =0而由4式可以看出,x’1= b’1,x’2可以为任意向量。

体上的矩阵方程XA=B及其反问题

体上的矩阵方程XA=B及其反问题

体上的矩阵方程XA=B及其反问题
王卿文;朱淑花
【期刊名称】《山东师范大学学报:自然科学版》
【年(卷),期】1997(012)004
【摘要】改进了体上矩阵方程XA=B的一般解的求法,给出了此方程及其反问题有自共轭解,正定自共轭解的充要条件有其解集表达式。

推广了有关数域上线性方程及矩阵反问题及其一系列结果。

【总页数】3页(P447-449)
【作者】王卿文;朱淑花
【作者单位】昌潍师范专科学校数学系;昌潍师范专科学校数学系
【正文语种】中文
【中图分类】O151.21
【相关文献】
1.子矩阵约束下矩阵方程XA=B的对称非负定解 [J], 彭振赟
2.矩阵方程AX-XA=0的求解问题 [J], 陈邦考;胡志龙
3.实四元数体上矩阵方程(A*XA,B*XB)=(C,D)的自共轭半正定解及亚半正定解 [J], 李小平
4.矩阵方程XA=B的亚(半)正定自共轭四元数矩阵解 [J], 王卿文;杨昌兰
5.体上矩阵方程AXB=C的(反)次自共轭解 [J], 林春艳
因版权原因,仅展示原文概要,查看原文内容请购买。

矩阵方程的求解问题

矩阵方程的求解问题

矩阵方程的求解问题白秀琴(平顶山工业职业技术学院,基础部,河南 平顶山 467001)摘要:主要考察了矩阵方程的求解问题,给出了一般矩阵方程当系数矩阵满足不同条件时的两种求解方法。

关键词:矩阵;矩阵的逆;矩阵方程矩阵是线性代数中的最重要的部分,它贯穿于线性代数的始终,可以说线性代数就是矩阵的代数,矩阵是处理高等数学很多问题的有力工具。

矩阵方程是矩阵运算的一部分,这里我们主要讨论如何求解矩阵方程的问题。

掌握简单的矩阵方程的求法,对于求解复杂的矩阵方程有很大帮助。

简单的矩阵方程有三种形式:.,,C AXB C XA C AX ===如果这里的A 、B 都是可逆矩阵,则求解时需要找出矩阵的逆,注意左乘和右乘的区别。

它们的解分别为:.,,1111----===BAX CAX C A X例如,求解方程C AC =先考察A 是否可逆,如果A 可逆时,方程两边同时左乘1-A ,得,11C A AX A --=即,1C A X -=这里要注意只能左乘不能右乘,因为矩阵的乘法不满足交换律。

同样,对于方程,C XA =只能右乘1-A ,得,11--=CAXAA即.1-=CAX 而对于方程,C A X B =只能是左乘1-A 而右乘1-B ,得,1111----=CB A ACBBA 即.11--=CBA X看下面解矩阵方程例题: 例1:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡315432343122321X 解:先求出1-A ,则,111253232313431223211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-则 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=332123315432111253232313154321343122321X 例2:⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡212101343122321X解:先求出1-A ,则,111253232313431223211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-则 ⎥⎥⎦⎤⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡-=-27525120111253232312121013431223212121011X 例3:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡3154321325343122321X 解:先求出1-A ,则,111253232313431223211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡- ⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡-532113251则 ⎥⎦⎤⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=--532131543211125323231132531543234312232111X ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎦⎤⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=131148735331332123例4:解矩阵方程,2X A E AX +=+其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=161020101A ,E 是三阶单位方阵。

11.matlab求解矩阵方程

11.matlab求解矩阵方程

11.如何求解矩阵方程?
在一个矩阵方程中,未知元素是一个矩阵。

线性代数中的矩阵方程总可以写成如下两种形式中的一种
AX = B 或 XA = B
这里, X 是未知矩阵,而 A 、 B 则是已知的矩阵。

这两种类型的矩阵方程的解分别为
X = A -1B 或 X = BA -1
在 MATLAB 环境中,可以分别用下面两条命令求解这两类方程
X = A \B 或 X = B /A
注意, A 左除 B (用反斜杠)相当于用 A 的逆矩阵左乘矩阵 B ;
A 右除
B (用斜杠)相当于用 A 的逆矩阵右乘矩阵 B 。

例如求矩阵方程
[11−10221−10]X =[1−11121
]
⎤可用如下命令
A = [1 1 -1; 0 2 2; 1 -1 0];
B = [1 -1; 1 1; 2 1];
X = A \B
计算机运行后,显示数据结果
X =
1.8333 0.5000
-0.1667 -0.5000
0.6667 1.0000
这就是矩阵方程的解。

显然,X 是一个 3× 2 阶的矩阵。

再例如,求矩阵方程
X [11−10221−10]=[1−11110
]
可用如下命令
A = [1 1 -1;0 2 2; 1 -1 0];
B = [1 -1 1; 1 1 0];
X = B/A
计算机运行后,显示数据结果
X =
-0.3333 0.3333 1.3333
0.6667 0.3333 0.3333
显然,这一矩阵方程的解是一个2× 3 的矩阵。

矩阵方程xa=b例题解法

矩阵方程xa=b例题解法

矩阵方程xa=b例题解法
两种方法:
1、转换成AX=B 的形式。

XA=B 两边取转置得A^TX^T = B^T 对(A^T,B^T)用初等行变换化为(E,(A^T)^-1B^T) = (E,X^T)
2、构造分块矩阵A B 用初等列变换化为E BA^-1 = E X
注:不要先求A^-1,那样会多计算一次矩阵的乘法!
扩展资料:
对于矩阵方程,当系数矩阵是方阵时,先判断是否可逆。

如果可逆,则可以利用左乘或右乘逆矩阵的方法求未知矩阵,如果方阵不可逆或是系数矩阵不是方阵,则需要用矩阵的广义逆来确定矩阵方程有解的条件,进而在有解的情形求出通解。

举个例子:
1 3
2 ……
3
4 -1
2 6 5 * X = 8 8 3
-1 -3 1 ……-4 1 6
上列就是个矩阵方程。

元素是实数的矩阵称为实矩阵,元素是复数的矩阵称为复矩阵。

而行数与列数都等于n的矩阵称为n阶矩阵或n阶方阵。

n×n的方块矩阵A的一个特征值和对应特征向量是满足的标量以及非零向量。

其中v为特征向量,为特征值。

A的所有特征值的全体,叫做A的谱,记为。

矩阵的特征值和特征向量可以揭示线性变换的深层特性。

线性矩阵方程AX=B和XA=B的初等变换解法

线性矩阵方程AX=B和XA=B的初等变换解法
,
的 初 等变 换 解法
r
的 形 式 根据前 面 的 讨 论 对 其 增 广 矩 阵进 行 初 等行 变 换 即 得

这是 A
0 1
X 一 B

,
,
述 不 失 一般性 不 妨 设 A 的 前 列 是 线 性 无 关 的 那 么 由 矩 阵 初 等 变 换 的 有关 性 质 可
,
(A
,
B
) 一

:
( 1 )
.
B ) ~
.
.
r
(
A )
,
,
由此 可 知
r
有 r (
r
A

.
B
;
,
B

B
)
=
r
(A )

l
(A
,
B )
.

(
A
)
熟知的结论 如下
:
我 们 将其 作 为 本 文 的 引 理 列 出 线 性 方 程组 A
:
仿 照定 理
定理
2
r
的 证法 据 引 理 2 又 可 证 明 对 于矩 阵方 程 A X 一 B 设 ( A
,
:
A
,
,
,
B
:
,
)
A )




,

A X 一B
的初 等变 换 解法
,

,
A (X
,
,
X
:
.

:
X ) 一 (B
r
.
B

(
B )
Z

王晓峰著《线性代数》习题解答

王晓峰著《线性代数》习题解答

王晓峰著《线性代数》习题解答第一章1. 解下列方程组, 并在直角坐标系中作出图示.1)⎩⎨⎧=-=+21y x y x ;2)⎩⎨⎧=+=+5331y x y x ; 3)⎩⎨⎧=-=-2221y x y x .解: 1) 将第一个方程减去第二个方程, 得2y =-1, y =-1/2, 再代入第个方程解得x =1+1/2=3/2,⎪⎭⎫ ⎝⎛-21,23方程有唯一解.2) 将第二个方程除以3得35=+y x , 与第一个方程相比较知此方程组为矛盾方程组, 无解,3) 将第2个方程除以2, 可以得到第一个方程, 令y =t 为任意实数, 则x =1+t , 方程组的解集.2. 用Gauss 消元法解下列线性方程组.1)⎪⎩⎪⎨⎧-=-+=++-=-+333693132472321321321x x x x x x x x x2)⎩⎨⎧-=-+=+-223252321321x x x x x x3)⎪⎪⎩⎪⎪⎨⎧=+-=-=--=+54212302433214243241x x x x x x x x x x4)⎪⎩⎪⎨⎧=++=-+=+033803403232132121x x x x x x x x解: 1) 对增广矩阵进行变换:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−−→−-⨯+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----−−−−→−+-⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----0000751010301)2(000075104721)3/1(12115302115304721)3()2(333693131124721123323121r r r r r r r r r则x 3为自由变量, 令x 3=t 为任意实数, 则x 1=10-3t , x 2=5t -7, 方程有无穷多解, 解集为(10-3t , 5t -7, t ).2) 对增广矩阵进行变换:⎥⎦⎤⎢⎣⎡--−−−→−+⨯⎥⎦⎤⎢⎣⎡---−−−→−⨯⎥⎦⎤⎢⎣⎡---−−−−→−+-⨯⎥⎦⎤⎢⎣⎡---121001012121025218/1816802521)3(2123252112221r r r r r则x 3为自由变量, 令x 3=t 为任意实数, 则x 1=-t , x 2=2t -1,解集为(-t , 2t -1, t ).3) 对增广矩阵进行变换:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−−−→−+-⨯+⨯+⨯-⨯⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡-----−−−−→−+-⨯+⨯⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡-----−−−−→−⨯-⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-------−−−→−+⨯+⨯↔⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------−−−−→−+-⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----11000101001001010001)3()32()35()43(34340003235100313201043001)7(461370032351003641043001)12/1()1(613700820120036410430012336410120300112043001)2(50412120300112043001142434443233242324241r r r r r r r r r r r r r r r r r r r r r方程有唯一解x 1=x 2=x 3=x 4=1.4) 此为齐次方程, 对系数矩阵进行变换⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-−−−→−+⨯+⨯⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+-⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-10003000211)6/1(6001301021)3(390130032)4()2(3381340321323312323121r r r r r r r r r r r r r可知方程有唯一零解x 1=x 2=x 3=0.3. 确定下列线性方程组中k 的值满足所要求的解的个数. 1) 无解: 2) 有唯一解:⎩⎨⎧=++=++;486362z y x kz y x⎩⎨⎧-=-=+123214y x y kx3) 有无穷多解:⎪⎩⎪⎨⎧=+-=++=++12524z y x z y x kz y x解:1) 对增广矩阵作变换:⎥⎦⎤⎢⎣⎡--−−−−→−+-⨯⎥⎦⎤⎢⎣⎡143800621)3(486362121k k r r k因此, 要使方程组无解, 须使8-3k =0, 解得k =8/3, 即当k 取值为8/3时, 方程无解. 2) 对增广矩阵作变换:⎥⎥⎦⎤⎢⎢⎣⎡++--−−−−−→−+-⨯⎥⎦⎤⎢⎣⎡--−−−→−↔⎥⎦⎤⎢⎣⎡--14612301232)2(141123212321412121k k r kr k r r k因此, 如要方程组有唯一解, 必须有0123≠+k , 即32-≠k . 3) 对增广矩阵作变换⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−→−+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−−→−+-⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-0440*******1331301110411)1()1(11215121411323121kkk r r k k k r r r r k因此, 如要方程组有无穷多解, 必须4-4k =0, 即当k =1时, 方程组才有无穷多解.4. 证明: 如果对所有的实数x 均有ax 2+bx +c =0, 那么a =b =c =0.证: 既然对所有的实数x 都有ax 2+bx +c =0成立, 那么具体地分别取x =0, x =1, x =2代入上式也成立, 则有⎪⎩⎪⎨⎧=++=++=02400c b a c b a c , 这是关于a ,b ,c 的齐次线性方程组, 对其系数矩阵作变换:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−→−↔↔⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡100320111)4(100124111124111100213221r r r r r r看出此方程只有唯一零解, 因此有a =b =c =0.5. 讨论以下述阶梯矩阵为增广矩阵的线性方程组是否有解; 如有解区分是唯一解还是无穷多解.1)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---0000320003212)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--410030201231 3)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--00004000320040214)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--0000010013201021 解: 1) 方程组有一个自由变元x 2, 因此方程组有无穷多解. 2) 方程组的三个变元均为首项变元, 因此方程组有唯一解. 3) 第三个方程0=4说明此方程无解.4) 方程组的三个变元均为首项变元, 因此方程组有唯一解.6. 对给定方程组的增广矩阵施行行初等变换求解线性方程组..1)⎪⎩⎪⎨⎧=-=+-=+-3284432253y x y x y x 2)⎩⎨⎧=--+=--+302859322207124w z y x w z y x 3)⎪⎩⎪⎨⎧=+-+=--+=+-+222242*********w z y x w z y x w z y x 解: 1) 对增广矩阵进行变换:⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--−−−−→−+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−→−↔⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---781007231032811974190723103281)28/1(74190922803281)3()3(2253443328132814432253322312113r r r r r r r r r方程组无解.2) 对增广矩阵进行变换⎥⎦⎤⎢⎣⎡--−−−−→−+⨯⎥⎥⎦⎤⎢⎢⎣⎡---−−→−⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+-⨯⎥⎥⎦⎤⎢⎢⎣⎡----−−−→−⨯⎥⎦⎤⎢⎣⎡----5452100100960317/4545210021154731422713410021154731)3(302859321154731)4/1(302859322207124122211r r r r r r可以看出y 和w 为自由变元, 则令y =s , w =t , s 与t 为任意常数, 则x =100-3s +96t , z =54+52t . 方程的解集表示为(100-3s +96t , s , 54+52t , t ). 3) 对增广矩阵进行变换()⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡-−−−−−→−+⨯⨯+-⨯⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--−−−−→−+-⨯+-⨯⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡----−−−→−↔⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----0000100021021211)2/1(2/1)2(04002000212121211)4()2(2222411112212121211222242121212111111212232312121r r r r r r r r r r r 可知y 与z 为自由变元, 令y =s , z =t , s 与t 均为任意实数, 则,212121=+-=w t s x , 方程组的解集为⎪⎭⎫ ⎝⎛+-0,,,212121t s t s7. 对给定齐次线性方程组的系数矩阵施行行初等变换求解下列方程组.1) ⎪⎩⎪⎨⎧=-+=+=+-02020z y x yx z y x 2)⎪⎩⎪⎨⎧=+-=+-=+++0202202w z y w y x w z y x解: 1) 对系数矩阵作初等变换.⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−−−→−+-⨯+⨯-⨯⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--−−−−→−+-⨯+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−→−⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+-⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--100010001)3/1()3/2()5/3(350032103101)2(320321011131320230111)1()2(21101211113233321223121r r r r r r r r r r r r r r方程只有零解, x =y =z =0.2) 对系数矩阵作初等变换⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-−−−−−→−+-⨯+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−-⨯⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−+⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−→−↔⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--110000102001)2()2/1(11002/12/1100201)3/1()2/1(3300112002012)1(114011201121112011401121)1(11202021112113233232123221r r r r r r r r r r r r r r因此, w 为自由变元, 令w =t 为任意实数, 则x =-2t , y =0, z =t , 方程组的解集为 (2t , 0, t , t ).8. 设一线性方程组的增广矩阵为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--32223411121α求α的值使得此方程组有唯一解.解: 对增方矩阵求初等变换⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+−−→−+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−+-⨯+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--420034601121126034601121)2(32223411121323121αααr r r r r r因此, 此方程组要有唯一解, 就必须满足α+2≠0, 即α≠-2.9. 设一线性方程组的增广矩阵为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----0410*******β1) 此方程有可能无解吗? 说明你的理由. 2) β取何值时方程组有无穷多解?解: 1) 此方程一定有解, 因为此方程是齐次方程, 至少有零解. 2) 对此增广矩阵做初等变换⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+--−−−→−+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−→−++⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----0500011001216016001100121204103520121323121βββr r r r r r因此, 只有当β+5=0, 即β=-5时,方程才有无穷多解.10. 求λ的值使得下述方程组有非零解.⎩⎨⎧=-+-=+-0)2(0)2(y x y x λλ 解: 对系数矩阵作初等行变换:⎥⎦⎤⎢⎣⎡+---−−−−−→−+-⨯⎥⎦⎤⎢⎣⎡---−−−→−↔⎥⎦⎤⎢⎣⎡---1)2(021)2(1221211222121λλλλλλλr r r r因此, 要使方程有非零解, 必须有(λ-2)2+1=0, 但(λ-2)2+1≥0对λ取任何实数值总是成立, 因此必有(λ-2)2+1≠0, 因此, 无论λ取什么值此方程组都不会有非零解.11. 求出下列电路网络中电流I 1,I 2,I 3的值.解: 根据基尔霍夫定律可得如下方程组:⎪⎩⎪⎨⎧=+=+=+-52384202132321I I I I I I I 对增广矩阵做初等行变换⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−−→−+-⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−−→−-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−+⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-13/1510013/2201013/7001)3()2(13/1510042104301)13/1(151300421043011)5(535042100111)2/1()3(502384200111132331232231r r r r r r r rr r r r最后得I 1=7/13, I 2=22/13, I 3=15/1312. 一城市局部交通流如图所示.(单位: 辆/小时)51) 建立数学模型2) 要控制x 2至多200辆/小时, 并且x 3至多50辆小时是可行的吗? 解: 1} 将上图的四个结点命名为A , B , C , D , 如下图所示:5则每一个结点流入的车流总和与流出的车流总和应当一样, 这样这四个结点可列出四个方程如下:⎪⎪⎩⎪⎪⎨⎧=+=++-=-+=+D x x C x x x Bx x x A x x 3502001503005453243121对增广矩阵进行变换:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---−−−−→−++-⨯+-⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--−−−−→−+-⨯+-⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----−−−−→−+-⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--000000350110002001011050010101)1()1(35011000350110001500111015001101)1()1(35011000200101101500111030000011)1(350110002001011015001101300000111323431232221r r r r r r r r r r r r r可见x 3和x 5为自由变量, 因此令x 3=s , x 5=t , 其中s ,t 为任意正整数(车流量不可能为负值), 则可得x 1=500-s -t , x 2=s +t -200, x 4=350-t .2) 令x 2=200, x 3=s =50, 代入上面的x 2的表达式, 得200=50+t -200, 求出t =350, 则x 1=500-s -t =100, x 4=0, 是可行的.13. 在应用三的货物交换经济模型中, 如果交换系统由下表给出, 试确定农作物的价值x 1, 农具及工具的价值x 2, 织物的价值x 3的比值.313131313131313131CM F C M F解: 根据上表可得关于x 1, x 2,x 3的三个齐次方程如下:⎪⎪⎪⎩⎪⎪⎪⎨⎧=-+=+-=++-032313103132310313132321321321x x x x x x x x x对系数矩阵做行初等变换:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−→−+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−→−-⨯+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+-⨯+⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−→−↔⨯⨯⨯⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡---0001101012000110121)3/1(1330330121)1(221111212133332313131323131313212232312121321r r r r r r r r r r r r r r可见方程有非零解, x 3为自由变量, 令x 3=t 为任意正实数, 则有x 1=x 2=x 3=t , 即三种价值的比值为1:1:1.第二章1. 1. 写出下列方程组的矩阵形式:1) x 1-2x 2+5x 3=-1;2) ⎩⎨⎧=+=-1223231x x x x 3) ⎪⎩⎪⎨⎧=-=+=++002045z x z y z y x 解:1) []15,2,1321=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-x x x ; 2)⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡12110102321x x x ;3) ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-000101120415z y x2. 设⎥⎦⎤⎢⎣⎡=212121A , ⎥⎦⎤⎢⎣⎡--=212234B求: 1) 3A -2B ;2) 若X 满足A T +X T =B T , 求X .. 解: 1)⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡--------=⎥⎦⎤⎢⎣⎡---⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡---⎥⎦⎤⎢⎣⎡=-10110105)4(623)4(64366834244686363632122342212121323B A2)因X 满足A T +X T =B T , 等号两边同时转置, 有 A +X =B ,等号两边同时减去A , 得 X =B -A , 因此有⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡--------=⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡--=-=404113221122122314212121212234A B X3. 计算下列矩阵的乘积:1)[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-213121; 2) []214321-⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡; 3)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎦⎤⎢⎣⎡-103110021212321; 4)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⎥⎦⎤⎢⎣⎡011011120101130213 解:1)[]1211231213121=⨯+⨯+⨯-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-2)[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯-⨯⨯-⨯⨯-⨯⨯-⨯=-⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡8463422124)1(423)1(322)1(221)1(12143213)⎥⎦⎤⎢⎣⎡---==⎥⎦⎤⎢⎣⎡-⨯+⨯+⨯-⨯+⨯+⨯-⨯+⨯+⨯--⨯+⨯+⨯⨯+⨯+⨯⨯+⨯+⨯==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎦⎤⎢⎣⎡-1341410)1(21102021122320112)1(312010312213302111031100212123214)⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡-⨯+-⨯+⨯-⨯+-⨯+⨯=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⎥⎦⎤⎢⎣⎡==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⨯+-⨯+⨯-⨯--⨯+⨯⨯+-⨯+⨯⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⎥⎦⎤⎢⎣⎡83)2(1)2(310)2(2)2(11322113021300)1(11101)1(21001)1(011130213011011120101130213 4. 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=201210003,310120101B A求: 1) (A +B )(A -B );2) A 2-B 2.比较1)和2)的结果, 可得出什么结论? 解: 1)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=-+567063519111110102511330104)201210003310120101)(201210003310120101())((B A B A2)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=-655142418405612009105055041120121000320121000331012010131012010122B A 可得出的结论: 大家知道, 在代数公式上有a 2-b 2=(a +b )(a -b ), 而将此公式中的a 和b 换成矩阵A 与B , 就不一定成立了, 这是因为矩阵乘法一般不满足交换律, 即一般AB ≠BA , 当然也就有A 2-B 2≠(A +B )(A -B ).5. 已知矩阵A ,B ,C , 求矩阵X ,Y 使其满足下列方程:⎩⎨⎧+=+=-T B A Y X CY X )(2解: 将此方程编上号, 用类似解线性方程组一样的办法来解,⎩⎨⎧+=+=-)2()()1(2T B A Y X C Y X将方程(1)的左边和(2)的左边和左边相加, 右边和右边相加, 等号还是成立, 得: 3X =C +(A +B )T 两边同乘1/3, 得TB AC X )(3131++=(3)(2)式等号两边都加上X , 得 Y =(A +B )T -X (4) 将(3)式代入到(4)式, 得CB A B AC B A Y T T T 31)(32)(3131)(-+=+--+=因此⎪⎩⎪⎨⎧-+=++=CB A YC B A X T T T T 3132323131316. 如矩阵AB =BA , 则称A 与B 可交换, 试证:1) 如果B 1, B 2都与A 可交换, 那么B 1+B 2, B 1B 2, 也与A 可交换; 2) 如果B 与A 可交换, 那么B 的k (k >0)次幂B k 也与A 可交换. 证: 1) 因B 1, B 2都与A 可交换, 即AB 1=B 1A , AB 2=B 2A , 则 (B 1+B 2)A =B 1A +B 2A =AB 1+AB 2=A (B 1+B 2) 即B 1+B 2与A 可交换. 而且(B 1B 2)A =B 1(B 2A )=B 1(AB 2)=(B 1A )B 2=(AB 1)B 2=A (B 1B 2), 因此B 1B 2与A 可交换.2)因B 与A 可交换, 即AB =BA , 则用归纳法, 当k =1时, 有B 1=B , 结论显然成立. 假设当k =m 时假设成立, 即AB m =B m A , 则当k =m +1时, 有AB m +1=AB m B =B m AB =B m BA =B m +1A , 结论也成立.7. 如矩阵A =A T , 则称A 为对称矩阵.设A ,B 都是n 阶对称矩阵, 证明AB 是对称矩阵的充分必要条件是AB =BA . 证: 已知A =A T , B =B T ,充分性: 假设AB =BA , 则(AB )T =B T A T =BA =AB , 因此AB 为对称矩阵. 必要性: 如果AB 为对称矩阵, 即(AB )T =AB , 则因 (AB )T =B T A T =BA , 可得BA =AB . 8. 设⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=n a a a A21其中a i ≠a j , 当i ≠j (i , j = 1,2, …, n ). 试证: 与A 可交换的矩阵一定是对角矩阵. 证:假设矩阵B ={b ij }n 与A 可交换, 即有BA =AB , 而BA 相乘得到的矩阵为B 的第j 列所有元素都乘上a j 得到的矩阵, AB 相乘得到的矩阵为B 的第i 行元素都乘上a i 得到的矩阵. 即BA ={a j b ij }n , AB ={a i b ij }n , 但对于任给的i ,j ,i ≠j , 因AB =BA , 因此有a j b ij =a i b ij , 因a i ≠a j , 所以必有b ij =0, 即B 只能是对角矩阵.9. 检验以下两个矩阵是否互为可逆矩阵?⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1000210012100121,1000210032104321B A解: 计算AB 和BA 如下:410000100001000011100012)2(1110013)2(21112)2(111014)2(31213)2(21112)2(11110002100121001211000210032104321I AB =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⨯+-⨯⨯⨯+-⨯+⨯⨯+-⨯⨯⨯+-⨯+⨯⨯+-⨯+⨯⨯+-⨯⨯==⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=41000010000100001110001)2(211100112)2(311)2(21110213)2(41112)2(311)2(21111000210032104321100021********21I AB =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⨯-+⨯⨯⨯+⨯-+⨯⨯-+⨯⨯⨯+⨯-+⨯⨯+⨯-+⨯⨯-+⨯⨯==⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=因此A 与B 确实互为逆矩阵.10. 设A ,B ,C 为n 阶方阵, 且C 非奇异, 满足C -1AC =B , 求证B m =C -1A m C (m 为正整数). 证: 用归纳法, 当m =1时条件已经成立为C -1AC =B , 假设当m =k 时, 命题成立, 即有 B k =C -1A k C , 则当m =k +1时, 有B k +1= B k B =C -1A k CC -1AC = C -1A k (CC -1)AC = C -1A k IAC = C -1A k AC = C -1A k +1C , 命题得证.11. 若n 阶矩阵A 满足A 2-2A -4I =0, 试证A +I 可逆, 并求(A +I )-1. 证: 将A 2-2A -4I =0改写为A 2-2A -3I =I ,先解一元二次方程组x 2-2x -3=0, 根据公式a acb b x 2422,1-±-=其中a =1, b =-2, c =-3, 则⎩⎨⎧-=+±=13212422,1x , 因此可将多项式x 2-2x -3因式分解为x 2-2x -3=(x -3)(x +1), 那么, 根据矩阵相乘相加的性质也就能将A 2-2A -3I 因式分解为 A 2-2A -3I =(A -3I )(A +I )=(A +I )(A -3I ), 因此我们有(A -3I )(A +I )=(A +I )(A -3I )=I , 即A +I 与A -3I 互为逆矩阵, (A +I )-1=A -3I .12. 证明: 如果A =AB , 但B 不是单位矩阵, 则A 必为奇异矩阵.证: 用反证法, 假设A 为可逆, 其逆为A -1, 则对于A =AB 两边同时左乘A -1, 得 A -1A =A -1AB , 即I =B , 这与B 不是单位矩阵相矛盾, 因此A 必为奇异矩阵.13. 判别下列矩阵是否初等矩阵?1) ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-100020001, 2) ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001010100 3) ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡010100201, 4) ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-100410001 解: 1) 是初等矩阵P (2(-2)),2) 是初等矩阵P (1,3), 3) 不是初等矩阵,4) 是初等矩阵P (3(-4), 2).14. 求3阶方阵A 满足⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡333231232221331332123111333231232221131211555a a a a a a a a a a a a a a a a a a a a a A解: 从等式看出A 左乘一矩阵相当于对此矩阵作初等行变换r 3×(-5)+r 1, 因此A 为一相应的初等矩阵, 即⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=-=100010501)1),5(3(P A15. 设A ,B ,C 均为n 阶可逆矩阵, 且ABC =I , 证明BCA =I证: 因B ,C 为可逆矩阵, 则BC 也是可逆矩阵, 且(BC )-1=C -1B -1, 因ABC =I , 对此等式两边右乘(BC )-1, 即ABC (BC )-1=I (BC )-1, 因为BC (BC )-1=I , 因此上式化简为A =(BC )-1, 因此当然有 BCA =BC (BC )-1=I .16. 设A ,B 均为n 阶方阵, 且)(21I B A +=, 证明: A 2=A 的充分必要条件是B 2=I .证: 充分性: 假设B 2=I , 则A IB I B I B B I B A =+=+=++=+=)(21)22(41)2(41)(41222必要性: 如果A 2=A , 则有)2(41)(41)(2122I B B I B I B ++=+=+等式两边乘4得I B B I B ++=+2222,等式两边同时减去2B +I 得 B 2=I 证毕.17. 如果n 阶矩阵A 满足A 2=A , 且A ≠I , 则A 为奇异矩阵.证: 用反证法, 假设A 为可逆, 其逆为A -1, 则上式两边左乘(或者右乘)A -1, 得 AAA -1=AA -1, 即A =I , 但这与A ≠I 相矛盾, 因此A 的逆不存在, 即A 为奇异矩阵.18. 求下列矩阵的逆矩阵:1)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=285421122A ; 2) ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=1111111111111111A 3)),,2,1,0(000000000000121n i a a a a a A i n n=≠⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=-解: 用对[A |I ]进行行初等变换为[I |A -1]的办法来求:1)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−→−↔⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=100285001122010421100285010421001122]|[21r r I A⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−−→−+⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----−−−−→−+-⨯+-⨯11390002196003/13/111)3/1()3(15018180021960010421)5()2(12323121r r r r r r r r ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−→−⨯⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−−−→−+-⨯+9/19/13/11006/16/13/10109/19/23/20019/16/11139001120609/19/23/2001)9/1(321323r r r r r r 因此, 最后得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=-9/19/13/16/16/13/19/19/23/21A 2)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=10001111010011110010111100011111]|[I A⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---------−−−−→−+-⨯+-⨯+-⨯10010220010120200011220000011111)1()1()1(413121r r r r r r ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---------−−−→−↔1001022000112200010120200001111123r r⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------−−−−−→−+⨯+-⨯11002200001122000101202002/102/10101)2/1()1(1242r r r r⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---------−−−−→−+⨯+-⨯11114000001122000101202002/12/1010012/1)1(1343r r r r⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------−−−−→−+⨯+⨯+⨯111140002/12/12/12/102002/12/12/12/100204/14/14/14/100012/12/14/1342414r r r r r r ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------−−−−→−⨯-⨯-⨯4/14/14/14/110004/14/14/14/101004/14/14/14/100104/14/14/14/100014/1)2/1()2/1(432r r r 因此有A A 414/14/14/14/14/14/14/14/14/14/14/14/14/14/14/14/11=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=-3)⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=-10000000000010000001000]|[121n n a a aa I A⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡−−−−→−↔↔↔----01000000100000001000100000012121211n n n n n n a a a a r r r r r r ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡−−−−→−⨯⨯⨯--0/1010000/100100000/10010/1000001/1/1/11211121n n n n n a a a a a r a r a r因此, 最后得⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=--0/10000/10000/1/10001211n n a a a a A19. 解下列矩阵方程, 求出未知矩阵X .1) ⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡12643152X 2) ⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--132321433312120X解: 令⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡=12643152B A , 则要解的方程为AX =B将方程两边左乘上A 的逆A -1, 可得A -1AX =A -1B , 即 X =A -1B 下面求A -1:⎥⎦⎤⎢⎣⎡--−−−−→−+-⨯⎥⎦⎤⎢⎣⎡−−−→−↔⎥⎦⎤⎢⎣⎡=21101031)2(0152103110310152]|[2121r r r r I A⎥⎦⎤⎢⎣⎡--−−−→−-⨯+⨯21105301)1(3212r r r 因此有⎥⎦⎤⎢⎣⎡--=-21531A 因此⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡--==-80232126421531B A X 2) 令⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=132321433312120B A 则矩阵方程为XA =B设A 的逆存在为A -1, 则方程两边右乘A -1, 得XAA -1=BA -1,即X =BA -1 下面求A -1:[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−→−⨯↔⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=10043300112002/102/32/112/1100433010312001120|121r r r I A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-−−−→−⨯+⨯12/302/12/30002/12/11002/102/32/112/13231r r r⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−−−→−+-⨯+⨯12/34/34/100002/12/11002/14/14/701)2/3(2/13212r r r r ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−→−-⨯463100002/12/11002/14/14/701)4(3r⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−−−→−+-⨯+-⨯4631002310107115001)4/7()2/1(1323r r r r因此,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=-46323171151A 最后得⎥⎦⎤⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----⎥⎦⎤⎢⎣⎡-==-47411246323171151323211BA X20. 求矩阵X 满足AX =A +2X , 其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=410011103A解: 将方程两边减去2X , 得AX -2X =A因2X =2IX , 因此上面的方程可以从右边提取公因子X , 得 (A -2I )X =A假设A -2I 可逆, 则方程两边同时左乘(A -2I )-1, 得(A -2I )-1(A -2I )X =(A -2I )-1A , 即X =(A -2I )-1A设B =A -2I , 则X =B -1A , 而⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=210011101200020002410011103B 下面用行初等变换求B 的逆B -1:[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=100210011110001101)1(100210010011001101|21r r I B⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----−−−−→−+-⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--−−−→−-⨯+⨯111100122010112001)1()1(111100011110001101)1(11323232r r r r r r r则⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=-1111221121B最后得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----==-3222342254100111031111221121A B X 验算:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=+1054459341364446844104100111032X A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=10544593413322234225410011103AX21. 利用分块的方法, 求下列矩阵的乘积:1) ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-100110201110021; 2) ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡d d c c b b a a00000010001010001000000解:1) 将乘积分块为[]⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-2|100110201110021I C B A其中[]10,201102,101=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=C B A[][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=+=⎥⎦⎤⎢⎣⎡30111220110210001020110210101|22BI AC I C B A2) 将乘积分块为⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡22222220000001000110001000000dI O cI I bI I O aI d d c c b b a a⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++=⎥⎦⎤⎢⎣⎡+=bd c bd c ac a ac a I bd c I acI aI 010*******)(2222第三章1. 计算下列行列式:1) 4321; 2) 22b b a a ; 3) 7040-解: 1) 26432414321-=-=⨯-⨯=;2) )(2222a b ab b a ab b b a a -=-=;3) 0)4(0707040=-⨯-⨯=-.2. 计算下列三阶行列式:1)241130421--; 2) 320001753-; 3) b a c a c b cb a 解: 1) 将行列式按第一列展开81021342124131241130421=+-=⨯-⨯-=-- 2) 将行列式按第二行展开172353275320001753=⨯-⨯==- 3)3333333c b a abc c b a abc abc abc b a c a c b cb a ---=---++=3. 计算下列行列式:1)000000005544332222211111b a b a b a e d c b a e d c b a ;2)x yy x y x y x D n 0000000000=;3) f e d c b a 0000000000解: 1) 将行列式按第一列展开后, 得到的各子式再按第二列展开, 这样展开后的后三列构成的任何三阶子式都至少包括一行0, 因此后三列任何三阶子式均为0, 整个行列式的值D =0. 2) 将行列式按第一列展开得nn n n n y x y x y x y y x y x y x x D 11)1(0000000)1(0000000++-+=-+=3) 先对第一列展开, 然后对第二列展开, 得abdfbadf fe dbafe dab D -=-=-=-=000004. 利用行列式的性质计算下列行列式1) 2605232112131412-; 2)ef cf bf de cd bd ae ac ab ---;3) 2222222222222222)3()2()1()3()2()1()3()2()1()3()2()1(++++++++++++d d d d c c c c b b b b a a a a解: 下面都将所求行列式的值设为D .1) 因为第1行加到第2行以后, 第2行将和第4行相等, 因此行列式的值D =0; 2) 首先从第1,2,3行分别提取公因子a ,d ,f , 再从第1,2,3列提取公因子b ,c ,e , 得abcdef abcdef adfbce ef cfbfde cd bd ae ac ab 4020200111111111111=-=---=---3) 将第2,3,4列都展开, 并统统减去第1列, 得9644129644129644129644122222++++++++++++=d d d d c c c cb b b b a a a a D 再将第3列减去2倍的第2列, 第4列减去3倍的第2列, 得62126212621262122222=++++=d d c cb b a a D5. 把下列行列式化为上三角形行列式, 并计算其值1) 1502321353140422-----; 2) 2164729541732152-----解:1)121034805350024211203840553004221)2/3(2150232135314042232413121------↔=-----+⨯+⨯+⨯=-----c c r r r r r r 131002050021102042101300520001210024258535034801210024243423242---↔=--+⨯+⨯=-----↔=c c r r r r r r270)27(512270002050021102042)2(43-=-⨯⨯⨯=----+-⨯=r r2)0210311061202251)1()2(12461759243712251216472954173215241312113----+-⨯+-⨯+⨯=------↔=-----r r r r r r c c93000030031102251133000300311022511)2(021061203110225143423232-=--+⨯=--+⨯+-⨯=---↔=r r r r r r r r6. 计算下列n 阶行列式1) 12125431432321-n n n2) a bbba b a解: 1) 设此行列式的值为D , 将第2,3,…,n 列均加于第一列, 则第一列的所有元素均为)1(21321+=++++n n n , 将此公因式提出, 因此有121125411431321)1(21-+=n nn n D再令第n 行减去第n -1行, 第n -1行减去第n -2行, …, 第2行减去第1行, 可得11111111111111111)1(21111011101110321)1(21-----+=--+=n n n n n n n n n n n n D 1)1(21)()1)(1(21)000000111111111)(1(21----+=---++=n n n n n n n n nn n2) 此题和第3题的2)一样, 因此有n n nb a D 1)1(+-+=7. 证明下列行列式1) ))()((111a c c b b a ab ca bc c b a ---=2) nb a n ab a ba b b a b a ba )(222-=证: 1)=----=----+-⨯+-⨯=)()()()(001)1()1(1113221c a b b a c ac a b c a b b a c bc a c a b a c c cc ab ca bc c b a))()(())()((11))((a c c b b a b c c a b a b c c a b a ---=---=----=2) 用归纳法, 设D n 为所求行列式值, 当n =1时,221b a a b ba D -==, 等式成立. 假设当n =k 时假设成立, 即有kk b a k aba b a b b a ba ba D )(222-==当n =k +1时,按第一列展开=+=+221k aba b ab b a b a ba D k=+++=1212k aba b b a ba b bk aa bab ba ba a12222222222)()()()(+-=--=-=-=k kk k k b a b a b a b a D D b D a证毕.8. 求矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=210111302A 的伴随矩阵A *, 并求A -1. 解:31130,32130,12111312111=-==--==--=A A A 11132,42032,22011322212=-=-=-==--=A A A 2112,21002,11011332313-=-=-=-==-=A A A因此得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=221142331332313322212312111*A A A A A A A A A A A 的行列式为5132012||131312121111=⨯+⨯+⨯=++=A a A a A a A 因此有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---==-22114233151||1*1A A A9. 设A 为三阶方阵, A *是A 的伴随矩阵, 且|A |=1/2, 求行列式|(3A )-1-2A *|的值.解: 因11**121||,||1---===A A A A A A A , 以及1131)3(--=A A , 还有2||1||1==-A A ,则27162278||32|32||31||2)3(|13111*1-=⨯-=⎪⎭⎫⎝⎛-=-=-=------A A A A A A10. 设A 为n 阶可逆阵, A 2=|A |I , 证明: A 的伴随矩阵A *=A . 证: 因A 可逆, 则在等式A 2=|A |I 两边乘A -1, 得A =|A |A -1, 即A A A ||11=-, 而因为*1||1A A A =-, 所以有A =A *, 证毕.11. 用克莱姆法则解下列方程组.(1) ⎪⎩⎪⎨⎧=+-=++=++10329253142321321321x x x x x x x x x(2) ⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++=+++24324322256511322121432143214321x x x x x x x x x x x x x x x x解: (1) 方程的系数矩阵A 为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=113215421A , 常数向量T ]102931[=β, 则求A 的逆矩阵:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------−−−−→−+-⨯+-⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-10311700151890001421)3()5(1001130102150014213121r r r r⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−−→−-⨯103117009/19/5210001421)9/1(2r ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−−→−+⨯+-⨯19/79/830009/19/521009/29/10017)2(3212r r r r⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---−−−→−⨯3/127/727/810009/19/521009/29/10013/13r⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−−→−+-⨯3/127/727/81003/227/1127/101009/29/1001)2(23r r 因此得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=-3/127/727/83/227/1127/109/29/11A则方程的解X 为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=-5431029313/127/727/83/227/1127/109/29/11321βA x x x X即x 1=3,x 2=4,x 3=5.(2) 方程的系数矩阵A 为⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=43114312251151132A , 常数向量[]T 2226=β先求A 的逆A -1:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−−→−↔⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡10004311010043120001511320010251110004311010043120010251100015113221r r⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------−−−−→−+-⨯+-⨯+-⨯10102200012007100021111000102511)1()2()2(413121r r r r r r ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------−−−−→−+⨯+-⨯101022000141160000211110003114011)1(3212r r r r ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------−−−−→−-⨯↔014116002/102/1011000021111000311401)2/1(343r r r⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-------−−−−→−+⨯+-⨯+-⨯311150002/102/1011002/102/512010201150016)1()4(332313r r r r r r ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------−−−−→−-⨯5/35/15/15/110002/102/1011002/102/51201020115001)5/1(4r ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-------−−−−→−+⨯+-⨯+-⨯5/35/15/15/1100010/15/110/75/1010010/75/210/295/70010110000011)2()5(342414r r r r r r 因此有⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-------=-5/35/15/15/110/15/110/75/110/75/210/295/711001A则⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-------==⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=-002022265/35/15/15/110/15/110/75/110/75/210/295/7110014321βA x x x x X 即x 1=0, x 2=2, x 3=0, x 4=0.12. 如果齐次线性方程组有非零解, k 应取什么值?⎪⎩⎪⎨⎧=-+=-+=++-0)4(20)6(2022)5(z k x y k x z y x k解: 此方程组的系数矩阵A 为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=k kk A 402062225要使方程组有非零解, 必须有det(A )=0.而k k k k kr r rr k kk A ---+--+⨯+-⨯=---=402242242252)2(402062225)det(2321k kk k r r rr k kk --+---+⨯+-⨯=-----=4022121005)2(2)2(402212225)2(1213)8)(5)(2(80061020122402212201)5)(2(3121----=---+⨯+⨯=-----=k k k kr r rr k k k因此, 只有当k =5或者k =2或者k =8时, 此方程组才有非零解.13. 问λ, μ取何值时, 齐次线性方程组⎪⎩⎪⎨⎧=++=++=++0200321321321x x x x x x x x x μμλ 有非零解?解: 此方程组的系数矩阵A 为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1211111μμλA , 要使方程组有非零解, 必须det(A )=0,而012101111)1()1(1211111)det(3121----+-⨯+-⨯==μλμλλμμλr r rr Aμλμμλμλμλ)1(12111)1(121113-=---=----=列展开按第因此, 只有当λ=1或者μ=0时, 方程组才有非零解.第四章1. 设α1=(1,1,1), α2=(-1,2,1), α3=(2,3,4), 求β=3α1+2α2-α3解: β=3α1+2α2-α3=3(1,1,1)+2(-1,2,1)-(2,3,4)=(3,3,3)+(-2,4,2)-(2,3,4) =(3-2-2, 3+4-3, 3+2-4)=(-1, 4, 1)2. 设3(α1-α)+2(α2+α)=5(α3+α), 求α, 其中α1=(2,5,1,3), α2=(10,1,5,10), α3=(4,1,-1,1) 解: 将上述方程整理: 3α1-3α+2α2+2α=5α3+5α -3α+2α-5α=-3α1-2α2+5α3 (-3+2-5)α=-3α1-2α2+5α3 -6α=-3α1-2α2+5α3 最后得)4,3,2,1()6531023,653521,653125,3103101()65,65,65,310()310,35,31,310()23,21,25,1()1,1,1,4(65)10,5,1,10(31)3,1,5,2(21653121321=-+++-+-+=--+=--+=-+=αααα3. 设R 为全体实数的集合, 并且设}0,,,|),,,({11211=++∈==n n n x x R x x x x x X V 满足, }1,,,|),,,({11212=++∈==n n n x x R x x x x x X V 满足.问V 1,V 2是否向量空间? 为什么?解: (一般的技巧: 凡是对R n 作一个齐次线性方程的约束的集合都是向量子空间, 而作非齐次线性方程的约束的集合则因为它不穿过原点, 就不是向量子空间).V 1是向量空间, 且是R n 的向量子空间, 因为nR V ⊂1, 而任给R k V Y X ∈∈,,1, 设0),,,,(0),,,,(121121=+==++=n n n n y y y y y Y x x x x x X则令),,,(2211n n y x y x y x Y X Z +++=+= ,则因=++++++=+++n n n y x y x y x z z z 221121011=+++++=n n y y x x , 则1V Y X ∈+,因为),,,(21n kx kx kx kX =, 而0)(11=++=++n n x x k kx kx 则1V kX ∈,因此, V 1是R n 的向量子空间.而V 2不是向量空间, 是因为1000≠+++ , 零向量O 不属于V 2, 2V O ∉.4. 试证: 由)1,1,1(),1,1,0(),1,0,0(321===ααα所生成的向量空间就是R 3证: 因为3321),,(R Span ⊂ααα, 只须证),,(3213αααSpan R ⊂, 任给3321),,(R d d d D ∈=, 试求实数x 1,x 2,x 3使。

矩阵方程AXB的解法

矩阵方程AXB的解法

2x2 -3x3 =3,x4 +2x5 -3x6 = -1,x1 +x2 -x3 =2,x4 +x5 -x6 = 0,解 得 x2 = 3-2x1,x3 = 1-x1x5 = 1- 2x4,x6 = 1-x4,于 是 所 求 矩 阵 的 解 为
X = éëêê31--x21xx11
x4 1-2x4 1-x4
方面都有非常广泛的应用.矩阵方程作为代数方程的一类,
在上述应用过程中,涉及核心问题必然是 矩 阵 方 程 解 的 判 定
和求解.文献已对矩 阵 方 AX =B 解 的 判 定 做 了 探 讨,本 文 来探讨矩阵方程 AX=B 的 求 解 方 法. 对 于 AX =B 的 求 解 方法,在我们熟知的线性代数教材中都局限于系数矩阵A 可 逆的情形,系数矩阵A 不可逆的情形基本上没有涉及.本文
例2 设 A =
éëêê341
-1 -3 3
230ùûúú ,B =
éëêê317
9 11 5
777ùûúú ,且 AX
= B,求 X .
解 :对 (A ,B)施 行 初 等 行 变 换 ,得 到 (A ,B)=
éëêê3 4 1
-1 -3 3
2 3 0
3 1 7
9 11 5
777ùûúú ~
éëêê100
给出了系数矩阵 A 不 可 逆 时 的 AX =B 的 两 种 普 适 性 求 解 方法,并通过例子进行了说明.文中出 现 的 数 学 符 号 都 是 标
准 的 ,可 参 见 文 献 .
二 、求 解 方 法 1——— 待 定 法 所谓的待定法,首先设出 未 知 矩 阵 X 中 的 待 定 元 素,然
后转化为线性方 程 组,用 解 方 程 组 的 方 法 求 出 矩 阵 X 的 待 定元素.

不可逆矩阵方程的解法ax=b

不可逆矩阵方程的解法ax=b

不可逆矩阵方程的解法ax=b
不可逆矩阵方程是在线性代数中比较重要的数学问题,一般来说,当矩阵A的行列式为零时 AX=B 就是一个不可逆矩阵方程。

一条不可逆矩阵方程无法求解,原因在于求解不可逆矩阵方程通常需要考虑矩阵A的行列式,若给定的矩阵A的行列式为零,则根本不可能反推矩阵A原有的正确形状,从而无法获得矩阵X的值,也就无法求解不可逆矩阵方程。

因此,如果想求解一个不可逆矩阵方程中的X,就要通过一些方式先让矩阵A可逆,以让其行列式不等于零。

其中有两种思路:
1. 通过行变换,它可以把矩阵A转换为消元矩阵,这样可以获得新的矩阵A的行列式,而且行列式的结果通常不会等于零。

2. 通过左乘一个可逆的矩阵去解AX=B方程,这里的可逆矩阵一般为逆矩阵。

左乘逆矩阵P后,新的方程PX=PB就很容易求解了,因为得到的X矩阵就是需要求解的原矩阵X。

最后,在求解不可逆矩阵方程ax=b时,一定要注意可逆性,确认矩阵A是否可逆,以免出现无法求解的情况。

线性代数-3.矩阵方程及其求解方法

线性代数-3.矩阵方程及其求解方法

1 0 解 ( A B ) 0 0
1 1 1 1 2 3 1 1 0 1 0 0
n 1 0 1 2 n 1 1 0 0 1 n 2 1 0 0 0 1
r1 r2 , r2 r3 , rn1 rn
1 1 P P m 1 A E, 即 〔 A B 〕 初等行变换 〔 E X 〕 . 从而 1 1 P P 1 B X; m
2. XA B,A 可逆,求 X .
1 1 AP P m 1 E, 求解方法II: 类似I,有 1 1 BP P 1 X; m
1 3 设三阶方阵A,B满足A1 BA 6 A BA, A
3 B 2 . 1
1 4
, 求B. 1 7
例5 设 A, B 为同阶方阵且 AB A kBA,若已知矩阵 A,能用什么 方法求矩阵 B ?(其中 k 为常数)
矩阵方程及其求解方法
标准的矩阵方程有三种形式:
AX B,XA B,AXC B,
其中 A , C 均为可逆阵.
1. AX B,A 可逆,求 X .
求解方法I: 因为 A 可逆,故有 X=A-1B.
求解方法II: 由 A 可逆知 A P 1 P m,P i 为初等阵,
1 -1 i 1, 2, , m,故有A1 Pm P1 ,Pi-1 也为初等阵,
1 0 0 2 0 , 1 1 0 1 0
0 1

1 0 1 X 0 2 0 . 1 0 1
2 1 0 ,矩阵 B 满足 ABA* 2 BA * E,求 B. 例4 设 A 1 2 0 0 0 1

用C语言求解N阶线性矩阵方程Ax=b的简单解法

用C语言求解N阶线性矩阵方程Ax=b的简单解法

用C语言求解N阶线性矩阵方程Ax=b的简单解法一、描述问题:题目:求解线性方程组Ax=b,写成函数。

其中,A为n×n的N阶矩阵,x为需要求解的n 元未知数组成的未知矩阵,b为n个常数组成的常数矩阵。

即运行程序时的具体实例为:转化为矩阵形式(为检验程序的可靠性,特意选取初对角线元素为0的矩阵方程组)即为:二、分析问题并找出解决问题的步骤:由高等代数知识可知,解高阶线性方程组有逆矩阵求解法、增广矩阵求解法等,而在计算机C语言中,有高斯列主消元法、LU分解法、雅克比迭代法等解法。

为了与所学的高等代数知识相一致,选择使用“高斯简单迭代消元法”,与高等代数中的“增广矩阵求解法”相一致。

以下简述高斯消元法的原理:算法基本原理:首先,为了能够求解N阶线性方程组(N由用户输入),所以需要定义一个大于N维的数组a[dim+1][dim+1](dim为设定的最大维数,防止计算量溢出),当用户输入的阶数N超过设定值时提示重启程序重新输入。

进而,要判断方程组是否有解,无解提示重启程序重新输入,有解的话要判断是有无数不定解还是只有唯一一组解,在计算中,只有当原方程组有且只有一组解时算法才有意义,而运用高等代数的知识,只有当系数矩阵对应的行列式|A|≠0 时,原方程组才有唯一解,所以输入系数矩阵后要计算该系数矩阵的行列式 |A|(定义了getresult(n)函数计算),当行列式 |A|=0 时同样应提示重启程序重新输入,|A|≠0 时原方程组必然有且仅有唯一一组解。

判断出方程组有且仅有唯一一组解后,开始将系数矩阵和常数矩阵(合并即为增广矩阵)进行初等行变换(以a11 为基元开始,将第j列上j行以下的所有元素化为0),使系数矩阵转化为上三角矩阵。

这里要考虑到一种特殊情况,即交换到第j-1列后,第j行第j列元素a jj=0 ,那此时不能再以a jj 为基元。

当变换到第j列时,从j行j列的元素a jj 以下的各元素中选取第一个不为0的元素,通过第三类初等行变换即交换两行将其交换到a jj 的位置上,然后再进行消元过程。

相关主题
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

矩阵方程xa=b例题解法
在连续时间系统的最优控制问题中,经常遇到连续代数黎卡提矩阵方程(CARE)。

在矩阵维数较大的情况下,求解该方程的解析解比较困难,需要花费大量的时间。

实际中,为得到该方程的近似解,常利用方程解的下界以减小计算复杂度。

通常利用的解的下界有行列式、特征值求和、范数、迹和矩阵约束等类型。

其中,矩阵约束最为常用,但这类约束具有较强的假设条件,比如Q为非奇异矩阵,或矩阵Q奇异但矩阵A在Q的列空间中。

实际求解时这些假设条件通常是不满足的。

在充分利用对称矩阵和半正定矩阵的特征值的性质基础上,运用矩阵不等式和类似于“李方法”的技巧,可以推导出一类新的解的下界。

这类下界的优势在于,在问题解存在的情况下不需要施加上述严格的限制条件。

同时,以这类下界为基础,可设计出求近似解的有效迭代算法。

数值算例表明,在参数选取合理的情况下,这类算法可得到比已有的数值结果更精确的解的下界。

相关文档
最新文档