第一篇:线性代数_课件LA1-1B
线性代数讲稿
讲稿编者:使用教材:《线性代数》
教学参考:《线性代数典型题分析解集》张 凯 院
西北工业大学出版社 西工大数学系编 西北工业大学出版社 徐 仲 等编
第一章
n阶行列式 §1.2 排列及其逆序数
1.排列:n个依次排列的元素.
例如, 自然数1,2,3,4构成的不同排列有4!=24种.
1234, 1342, 1423, 1432, 1324, 1243
2134, 2341, 2413, 2431, 2314, 2143
3124, 3241, 3412, 3421, 3214, 3142
4123, 4231, 4312, 4321, 4213, 4132
例1 互异元素p1,p2,,pn构成的不同排列有n!种.
解 在n个元素中选取1个
n种取法
在剩余n1个元素中选取1个
在剩余n2个元素中选取1个
n1种取法 n2种取法
„„„„„„
„„„„
在剩余2个元素中选取1个
2种取法
在剩余1个元素中选取1个
1种取法
------------------
总共n!种取法
2.标准排列:n个不同的自然数从小到大构成的排列.
n个不同的元素按照某种约定次序构成的排列. 3.逆序数:
(1)某两个数(元素)的先后次序与标准次序不同时, 称这两个数(元素)
之间有1个逆序.
(2)排列p1p2pn中逆序的总和称为排列的逆序数, 记作(p1p2pn).
算法:固定i(2,,n), 当
满足pjji时,iipi的“pj”的个数记作(称为p的逆序数),那么(p1p2pn)2n.
271032261
4例2 排列6372451中, .
例3 排列13(2n1)(2n)(2n2)42, 求逆序数.
解
记作p
1p2pnpn1pn2p2n1p2n
20, ,n10
n2221, n3422,„, 2n2(n1)
2[12(n1)]n(n1)4.奇偶性:排列p
(p1p2pn)(p1p2pn)p2pn
奇数时, 称为奇排列; 偶数时, 称为偶排列.
5.对换:
相邻对换:p
一般对换:p
1pipi1pnp1pi1pipn pipjpnp1pjpipn(ij)
1定理1 排列经过1次对换, 其奇偶性改变.
证
先证相邻对换:(1)
(2)
ab:对换后aaba1alabb1bma1albab1bm
2增加1, 不变, 故tbbt11; t11. :对换后不变, 减少1, 故ta所以t与t的奇偶性相反.
2再证一般对换:(1)
(2)
(3)
a1alab1bmbc1cn a1alb1bmabc1cn a1albb1bmac1cn
(1)(2)经过m次相邻对换
(2)(3)经过m1次相邻对换
(1)(3)经过2m1次相邻对换, 所以t与t的奇偶性相
31反.
推论 奇排列标准排列, 对换次数为奇数.
偶排列标准排列, 对换次数为偶数.
§1.3 n阶行列式的定义
1.二阶: a11a21a11a12a22a12a22a32a11a22a12a21
a13a23a33a11a23a32a12a21a33a13a22a31 2.三阶: a21a31a11a22a33a12a23a31a13a21a32
(1)乘积中三个数不同行、不同列:a
行标(第1个下标):标准排列 123
列标(第2个下标):p11p1a2p2a3p3
p2p3是
1,2,3的某个排列(共6种)
(2)正项:123, 231, 312为偶排列
负项:132, 213, 321为奇排列
a11a12a22a32a13a23a3
3于是 a21a312(1)(p1p2p3)a1p1a2p2a3p3,(p1p2p3).
3.n阶:n个数aij(i,j1,2,,n), 称
a1na2nanna11a12a22an2
Da21an1
为n阶行列式, 它表示数值
(1)(p1p2pn)a1p1a2p2anpn, (p1p2pn)
其中, 求和式中共有n!项.
a11a12a22a1na2nanna11a1,n1a2,n1a1na21an例3 计算D1,D2.解 数为
D1中只有一项a11a22ann不显含0, 且列标构成排列的逆序(12n)0, 故D1(1)a11a22anna11a22ann.
D2中只有一项a1na2,n1an1不显含
0, 且列标构成排列的逆序数为
(n21)12(n1)n(n1)2
故D
n(n1)2(1)a1na2,n1an1(1)2a1na2,n1an1.
结论:以主对角线为分界线的上(下)三角行列式的值等于主对角线上元素的乘积.
以副对角线为分界线的上(下)三角行列式的值等于副对角线上元素
的乘积, 并冠以符号(1)
特例:
1n(n1)2.
1212n,2n(n1)(1)212n
nn
定理2
a11Da21an1
(2)
证a12a22an2a1na2nann(1)(q1q2qn)(q1q2qn)aq11aq22aqnn
(p1p2pn)由定义知
D (1)a1p1a2p2anpn6(p1p2pn)
(1)
先证(2)中的项都是(1)中的项:交换乘积次序可得
(3)
①
(q1q2qn)偶数
q1q2qn12n(1)(q1q2qn)aq11aq22aqnn(1)(q1q2qn)a1p1a2p2anpn
偶数次对换
偶数次对换 偶数 12np1p2pn
所以(p ② (q
11p2pn)q2qn)奇数
奇数次对换
奇数次对换 奇数 , 由(3)可得
(p1p2pn)q1q2qn12n12np1p2pn
所以(p
因此(1)
(1)1p2pn)(1)(q1q2qn)(p1p2pn)(q1q2qn)aq11aq22aqnn(1)a1p1a2p2anpn
同理可证(1)中的项都是(2)中的项.
课后作业:习题一
1,2,3
§1.4 行列式的性质
a11a1nanna11an1ann
性质1 设D
证 令b
DΤan1,DΤa1n, 则DΤD.
ijaji(i,j1,2,,n), 则
b11b1nbn1(1)(p1p2pn)b1p1b2p2bnpn
(p1p2pn)
bnn
(1)(p1p2pn)ap11ap22apnnDai1ainajn
(根据Th2)
aj1ajnain
性质2 设ij,Daj1,D1ai1, 则D1D.
证
bikajk,bjkaik(k1,2,,n)
li,j:blkalk(k1,2,,n)
bi1tbinbjn
D1bj1(1)(bipibjpj)
(pipj)
推论1 (1)(1)(bjpjbipi)
t(pjpi)
(1)(1)(aipjajpi)ttqipj,qjpili,j:qlpl
(1)(1)(aiqiajqj)DDt(qiqj)
对调两列得D2D2D.
TT证 因为D对调两列得D, 相当于D对调两行得D
2所以D 2D2DTTD
推论2 D中某两行(列)元素对应相等D0.
证 因为对调此两行(列)后,D的形式不变
所以DDD0
12b23c03
例如, 对于任意的a,b,c, 都有a1a11a1nkainkDanna11.
ka1jkanja1nkD
性质3 kai1an1,an1
ann
证(1)左端(1)[a1p
推论1
推论2 DD1(kaipi)anpn]
(ppp)
1ink(1)(a1p1aipianpn)kD
中某行(列)元素全为0D0. D0. 中某两行(列)元素成比例ij
性质4 若对某个i, 有aa11a1na11bijcij(j1,2,,n), 则
a1na11a1ncinann
ai1an1ainbi1annan1binci1annan1
证 左端(1)(a1p
1aipianpn)
(ppp)
1in(1)(a1p1bipianpn)(1)(a1p1cipianpn)
右端(1)+ 右端(2)[注] 性质4对于列的情形也成立.
ai1ain
性质5
rikraji1aj1ainajn(ij)
aj1ajnaj1ajn [注] 性质5对于列的情形也成立.
1533
例5 计算D20113112.413115331533153
解 D01055211110161011500023(5)0002021911011102115331533
1110111(5)00023(5)002355
0031000112xaa
例6 计算Dxana.
aax111
解 r1(r2rn)Dxan[x(n1)a]a
aax111
[x(n1)a]0xa0
00xa
[x(n1)a](xa)n1
3131
12210030103010n00
例7 计算Dn3n.
n01t2100
解 c1jc0jDnj2,,n0001(2n)
221
§1.5 行列式按行(列)展开
余子式:在n阶行列式中,将元素a所在的行与列上的元素划
ij去,其余
元素按照原来的相对位置构成的n1阶行列式,称为元素a的
ij
余子式,记作M.
ij
代数余子式:元素a的代数余子式Aijij(1)ijMij.
a11a12a22an2a1na2nann
定理3 Da21an1
ai1Ai1ai2Ai2ainAin
(i1,2,,n)(j1,2,,n)a1jA1ja2jA2janjAnj
证
证明第一式, 分以下3步.
a11a1,n1
第1步:Mnnan1,1(1)(p1pn1)a1p1an1,pn1
an1,n1(1pin1)
a11a1,n1an1,n10a1nan1,nann
an1,10(1)(p1pn1pn)a1p1an1,pn1an,pn
(1)pnnpnn(p1pn1pn)a1p1an1,pn1an,pn+
(1)
ann(p1pn1pn)a1p1an1,pn1an,pn a1p1an1,pn1
MannAnn(1)nn(p1pn1n)
(p
D21pn1n)(p1pn1)
annMann(1)nnnna1jD1ai1,j
第2步: D(i,j)00aijai1,janj00
D3D4a1jD1D2ai1,j
(1)(ni)(nj)ai1,jD3000D40anjaij
(1)(ij)aijMijaijAij
第3步:DD(i,1)D(i,2)D(i,n)
ai1Ai1ai2Ai2ainAin
15011712331133121
例8 计算D234.
160010211解 D23116(1)3221471321
20010a51(1)(1)122
(1)272075155
baaccbddd00(1)(2n1)12nb
例9 计算D2n.
c
解 D2n(1)11a0D2(n1)00db0c0D2(n1)0
(2n1)
(1)(2n1)(2n1)adD2(n1)(1)(1)n1(2n1)1bcD2(n1)
(adbc)D2(n1)(adbc)acbdD2
D2adbc
D2n(adbc)n
111200023003n10n1n
例10 计算Dn111.
解 DnnDn1(1)n1(n1)!
(n1)1
n(n1)Dn2(1)(n11)!(1)n1(n1)!
n(n1)Dn2(1)nn!n1(1)n1n!n
n
4n(n1)3D2(1)n!3(1)n!n1(1)n1n!n
D2111221(1)2(1)1
2334n1(1)2(1)(1)(1)Dn(n!)123n
课后作业:习题一(1)(2)(1)(2)(3)(1)(2
1x11x2x2x2n121xn1xn1xn1n121xnxnxnn12 例11 证明Dnx1x12(x1jinixj).n1 14
11(x2xn)x2(x2xn)x2n21(xn1xn)xn1(xn1xn)xn1(xn1xn)n21000 证 D(i)xn(i1)nin,,2(x1xn)x1(x1xn)x1n2
(x1xn)(x2xn) (1)(x D D Dkn1n(x1xn)(x2xn)(xn1xn)Dn1
xn1)(xnxn2)(xnx1)Dn1
(xkxk1)(xkxk2)(xkx1)Dk1(kn,n1,,3)
21x11x2x2x1
n(xnxn1)(xnxn2)(xnx2)(xnx1)
(xn1xn2)(xn1x2)(xn1x1)
„„„„„„
(x3x2)(x3x1)
(xa11am1a1mamm00b11bn100b1nbnn2x1)
例12 证明 D
a11a1mb11b1n a11a1mammam1p1
ammbn1bnn 证 D行倍加1am1pmp1pm
b11b1nbnn D列倍加q102bn1p1q1qnqn00
pm0q1 D前m行“行倍加”后n列“列倍加”(p1pm)(q1qn)D1D2
qn 定理
ai1Aj1ai2Aj2ainAjn04 设ij, 则 a1iA1ja2iA2janiAnj0.
证 只证第一式.ij时, 有
ai1ainajnainainai1Aj1ai2Aj2ainAjnaj1Aj1aj2Aj2ajnAjn Daj1ai1
0ai1
[注]结合定理3与定理4可得
Dai1Aj1ai2Aj2ainAjn0a1iA1ja2iA2janiAnj1241433334122(ij)(ij)(ij)(ij)
D0
例13 D241, 求A11A21A31A41.
12414333341220 解法1 因为D1111 D与D的第1列元素的代数余子式相同
所以将D按第1列展开可得AA21A31A410.
解法2 因为D的第3列元素与D的第1列元素的代数余子式相乘求和
为0,即 3A 所以
§1.7 Cramer法则
考虑线性方程组
a11x1a12x2a1nxnb1a21x1a22x2a2nxnb2 an1x1an2x2annxnbna11a12a22an2a1na2nanna21an1113A213A313A410
A11A21A31A410
D
b1a12a22an20a1na2nanna11b1b2bna13a23an3a1na2nann(j)D(1)b2bn,D(2)a21an1, „„
定理5 若D 证 存在性., 则方程组存在唯一解xjDD(j1,2,,n).17 bib1~ Dbibnai1a11ai1an1aija1jaijanjaina1nr1ainri1ann0(r1ri1)
第1行中元素a的代数余子式为
ij~ Ab1ija11an1a11a1,j1an,j1a1,j1an,j1a1,j1an,j1b1bna1n(1)1(j1)bn
a1nD(j)a1,j1an,j1ann (1)j2(1)j1an1
ann~ 将D按第1行展开可得
biDai1(D(1))aij(D(j))ain(D(n))0
因为D ai10, 所以
aijD(j)D(1)DDjainD(j)D(n)Dbi(i1,2,,n)
故方程组有解 xD(j1,2,,n)
*1 唯一性.设方程组还有解xa11a1,j1an,j1a1,j1an,j1a1,j1an,j1*,x2,,xn, 则
**a1jxjanjxj*a1,j1an,j1*a1n x*jDan1a11
***anna1,j1**(a11x1a1jxja1nxn)(an1x1anjxjannxn)b1bna1,j1an,j1a1nannD(j)a1n an1a11
an,j1ann an1
同理可得 xjDD(j)于是 x*jD xjDxjxj(j1,2,,n)
* 例
x1x22x1x214 解线性方程组3x12x2x1x29, D(1)9, D(2)18(3)x32x40x3x40x35x45x3x41.解 D
D x 127, D(4)9
1, x22, x33, x41
a11x1a12x2a1nxn0a21x1a22x2a2nxn0 齐次方程组
an1x1an2x2annxn0
定理6 若D0, 则齐次方程组只有零解.D0.推论 齐次方程组有非零解 [注] D0齐次方程组有非零解.(定理3.5之推论)
例
x1x2x3015 已知 x1x2x30 有非零解, 求xxx0231.111(2)(1)02 解 D111, 故1或2. 19
a1b1a2a2b2a2anananbn 例16 计算Dna1a1(bi0). 解 采用加边法.10a1a1b1a1a111a1b100a20b20a2a2a2b2a2an0ananananbnt0a1b100a20b20an00bn Dn00
1100bn0
tbb1ana1a21bbbb2n12nb1b2bn
课后作业:习题一
8,9
第二章
矩阵及其运算 §2.1 矩阵
1.方程组由其系数和右端项确定
a11x1a12x2a1nxnb1a21x1a22x2a2nxnb2 am1x1am2x2amnxnbma11a12a22am2a1na2nammb1b2bm
a21am1
2.矩阵
设mn个数a表
ij(i1,2,,m;j1,2,,n)排成m行n列的数
a11a12a22am2a1na2namn
a21am1
用括号将其括起来, 称为mn矩阵, 并用大写字母表示, 即
a11a21Aam1a12a22am2a1na2n, amn
简记为A(aij)mn.(1)方阵
(2)矩阵
(3)矩阵 aij
称为A的i行j列元素
(4)mn
称A为aijR
称A为实矩阵
(5)m1,n1 称A为行aijC
称A为复矩阵
(6)m1,n1 称A为列
零矩阵:所有元素都是0 的矩阵.100010010;对角矩阵 01000
单位矩阵 En200n
3.线性变换与矩阵
设变量y为
1,y2,,ym可由变量x1,x2,,xn表示y1a11x1a12x2a1nxny2a21x1a22x2a2nxn ymam1x1am2x2amnxn
称之为由变量x矩阵
A(aij)mn1,x2,,xn到变量y1,y2,,ym的线性变换, 它与是一一对应关系.
§2.2
矩阵的基本运算
同阶矩阵:指行数相等、列数相等的矩阵.
矩阵相等:设A(a
ij)mn,B(bij)mn, 若 , 称AB.aijbij(i1,2,,m;j1,2,,n))mn,B(bij)mn
1.线性运算:A(a
加法:AB(aijijbij)mna11b11am1bm1a1nb1n amnbmn
数乘:kA(kaij)mnka11kam1)mnka1n kamn
负矩阵:A(1)A(a
减法:AB(abij)mnij
a1nb1namnbmnija11b11am1bm122
算律:设A,B,C为同阶矩阵,(1)
(2)
(3)
(4)
例1 ABBAk,l为常数, 则有
(5)1AA
(6)(kl)Ak(lA)
(kl)AkAlAk(AB)kAkB(AB)CA(BC)AOA
(7)
(8)
08, B552364 A(A)O231设A4
满足2AX
解 1B2X, 求X.
212X(B2A)3121
2.矩阵乘法:
q1q2qn
特殊情形 P1np1p2pn,Qn1
一般情形 PQp1q1p2q2pnqn
A(aij)ms,B(bij)sn
b1jb2jabababi11ji22jissjbsj
cijai1ai2ais
a11ABam1a1samsb11bs1b1nc11bsncm1c1ncmn
[注] A的列数 = B的行数.
AB的行数 = A的行数;AB的列数 = B的列数.
A与B的先后次序不能改变.
例2 3A01130,1B0021110,3AB01260231301
[注] BA无意义.
例3
1A122,1B111
1AB1BA10, BA100 0O [注] AB ;AO, B, 但是BAO.
算律:(1)
(2)
(3)
(4)(AmsBsn)CnlA(BC)Ams(BsnCsn)ABAC(AmsBms)CsnACBCk(AmsBsn)(kA)BA(kB)
则
nEmAmnA,AmnEnA
验证(1)设A(a
ij)ms,B(bij)sn,C(cij)nl,[(AB)C]ijaikbk1k1sc1jaikbknk1cnjst1saikbktctjk1
[A(BC)]ijai1nbc1ttjt1aisnbstctjt1nabcikkttjk1t1s
nt1saikbktctj[(AB)C]ijk1
(i,j)
a11a应用:A21am1a12a22am2a1nx1a2nx, x2amnxnb1b2bbmy1y2yym,,线性方程组的矩阵形式 Axb
线性变换的矩阵形式 yAx
3.方阵的幂:
Ann, k,l为正整数
A1A, Ak1AkA(k1,2,)
算律:(1)AkAlAkl
(2)(Ak)lAkl
101
例4 A20k, 求A(k2,3,).
1101101102 解法1 A22020220
11110210110
A3A2A22020231110k
可以验证:Akk20
1 25
30
1 解法2
1A021101k20010k000k110BC0k
BCCB(BC)C2BkBkCCk1
OAk(BC)BkkBC
k1
112k1k10k01020010110000002k100
000
4.矩阵的转置:
a11a21Aam1a12a22am2(A)(kA)TT2kk01
a1na11a2na, AT12amna1na21a22a2nTam1am2amnATTT
算律:(1)
(3)
验证(4)
A kA
(2)
T(AmnBmn)(AmsBsn)TBT
T
(4)
BA
A(aij)ms, B(bij)snTT
ABC(cij)mn, BAD(dij)nm
左ijcjiaj1ajsb1iaj1b1iajsbsibsi
右
故
ijdijb1iaj1bsib1iaj1bsiajscjiajs
dijcji(i1,2,,n;j1,2,,m),即(AB)TBTAT.
T
对称矩阵:指A满足AnnA26,即aijaji(i,j1,2,,n)
反对称矩阵:指A满足AnnTA,即aijaji(i,j1,2,,n)
5.方阵的行列式:指A(a的
行列
算律:(1)
(3)detATij)nn的元素按照原来的相对位置构成, 记作detA, 或者A.
(2)
(4)
det(lA)ldetA
nkdetAdet(AB)(detA)(detB)
detA(detA)k
[注] 方阵是数表, 而行列式是数值.
6.伴随矩阵:A(aa11a21Aan1a12a22an2*ijAnnBnnBA, 而det(AB)det(BA).)nn, detA中元素a的代数余子式为A.
ijij
*a1nA11a2nA,A*12annA1nA21A22A2nAn1An2Ann
重要性质:AAAA(detA)E
.
7.共轭矩阵:复矩阵A(a
算律:(1)
(3)
§2.3 逆矩阵
定义:对于Annij)mn的共轭矩阵记作A(aij)mn(AB)AB
(2)
(kA)kAT
T记作H(AB)AB
(4)
(A)(A)A , 若有Bnn满足ABBAE1, 则称A为可逆矩阵,且B为A的逆矩阵, 记作A
B.
定理1 若Ann为可逆矩阵, 则A的逆矩阵唯一.
证
设B与C都是A的逆矩阵, 则有
定理2
ABBAE,ACCAE
BBEB(AC)(BA)CECCAnnAnn为可逆矩阵为可逆矩阵1detA0;
A11detAA*.
证
必要性.已知A存在,则有
AA1EdeAtdeAt1deAt01
充分性.已知detA0,则有
AA*AA(detA)EA*AEdeAtdeAt1A*A*
*
由定义知A为可逆矩阵,且A [注]detA0时, 亦称A1detAA.
为非奇异矩阵;
detA0时, 亦称A为奇异矩阵.
推论1 对于A
证
nn, 若有Bnn满足ABE, 则A可逆, 且A1B.
ABEdetAdetB1detA0A可逆 , 则A可逆, 且A1A1A1EA1(AB)(A1A)BEBB
推论2 对于A
算律:
(1)Ann, 若有Bnn满足BAEB.
可逆1A1可逆, 且(AA1)1A.
1
对于A, 取B
(2)A, 有A1BAAE1k.
A1可逆, k0kA可逆, 且(kA)A11.
1
对于kA, 取B
1k, 有(kA)B(kA)(28
1kA)AA1E.
(3)Ann与Bnn都可逆B1ABA1可逆, 且(AB), 有
11B1A1.
对于AB, 取C
(4)(AB)C(AB)(BA1A1)A(BBT)A1E)T.
可逆TAT可逆, 且(A(A1)1(AT1.
1
对于A, 取B
(5)
(6)A)T, 有ABA(AT)T(A1A)TE.
可逆detA与Bnn11detA.
***Ann都可逆(AB)BA1.
1
证(AB)[det(AB)](AB)[(deBt)B01*[(detA)(detB)][B1A1]
][(deAt)A]BA1**
则有
负幂:A可逆, 定义A
例1
例2 设A 解 2El,klAk(A)k(k1,2,), AAA3A21klkl,(A)Ak
(k,l为整数)
4121131111014,A151*1A10550
nn满足A22A4EO,2求(A
E)1.
A2A4EOA2A3EE
应用:
(AE)(A3E)E(AE)1A3E
(1)n阶线性方程组求解
(2)求线性变换的逆变换
(3)矩阵方程求解
设A
Annxb, detA0xAyAnnx1b1
y,detA0xA
mm可逆,Bnn可逆, 且Cmn已知, 则
例3 AXCXAC11
XBCXCB
1AXBCXA1CB
1055设A22131016,2C23 满足AXC2X, 求X.
解
并项:
计算:(A2E)XCX(A2E)1
131223130751011C
1设A11511050111*412111 1
例4 满足A
*XA12X, 求X.
解
并项:
左乘A:
计算:
密码问题:
(A2E)XA1[(detA)E2A]XEdeAt4
X(4E2A)112(2EA)111041110011
a1, b2,c3, „ ,z02126
111
1A1021132 2121,A1
action:1, 3, 20, 9, 15, 14
167加密:A3442043,981A15521443
发出∕接收密码:67, 44, 43, 81, 52, 43
解密:A1671443 4320,A181952154314
明码:1, 3, 20, 9, 15, 14表示action
§2.4 分块矩阵
11A0011A00000000001120112010A111A21310B113
A12 A22
B2B3B4
用若干条横线与纵线将矩阵A划分为若干个小矩阵, 称这些小矩阵
为A的子矩阵, 以子矩阵为其元素的矩阵称为分块矩阵.
特点:同行上的子矩阵有相同的“行数”;
同列上的子矩阵有相同的“列数”.
1.加法:AA11As1A1rAsrmn,BmnB11Bs1B1rBsr
A11B11ABAs1Bs1A1rB1rAsrBsr
要求:A与B同阶, 且分块方式相同.
2.数乘:kAkA11kAs1kA1rkAsrA1tAstmn
B11Bt1B1rBtr
3.乘法:AmlA11As1,Bln
32
CijAi1B1jAitAi1B1jAitBtjBtj
C11ABCs1C1rCsr
要求:A的列划分方式与B的行划分方式相同.
10A11B***100E0A211104201B111B210OE
例1
E B22
B11ABA21B11B21EA21B221121TA11AT1r024110330131T
4.转置:AmnA11As1A1rAsr,ATAs1TAsr
特点:“大转”+“小转”
5.准对角矩阵:设A,A,,A都是方阵, 记
12s
A1Adia(gA1,A2,,As)A2As
性质:(1)
(2)detA(detA1)(detA2)(detAs)A可逆Ai(i1,2,,s)可逆
(3)Ai(i1,2,,s)可逆A1A11A211As
例2 5A000320A11O1OA2
012013
A1A11O15O01A20
OB
例3 设A 解
mm与Bnn都可逆,Cnm,AMC, 求M.
1detM(detA)(detB)0M可逆
X1X3X1X2X3X4M1X1X3X2X4,ACOBX2EmX4OAOBB111OEn
AX1EmAX2OCX1BX3OCXBXE24n1
CA1
MA11BCAO1B
课后作业:习题二 7(1)(3)(5), 8(2)(4), 10~14
第三章
矩阵的初等变换 §3.1 矩阵的秩
1.子式:在A原来的 mn中, 选取k行与k列, 位于交叉处的k个数按照
相对位置构成k阶行列式, 称为A的一个k阶子式, 记作D.
k
对于给定的k, 不同的k阶子式总共有C 2.矩阵的秩:在AmnkmCn个.
k中,若
r
(1)有某个r阶子式D0;
(2)所有的r1阶子式D
称rankO0 Ar1. 0(如果有r1阶子式的话)
r(A)r的秩为r, 记作rankAr, 或者
rankAmnmin{m,n} k0.规定:
性质:(1)
(2)
(3)
(4)
(5)例1 2A21时rank(kA)rankA
TrankAAArankA
中的一个D中所有的D8212124r0rankAr
r10rankAr3123, 求r(A).
解
位于1,2行与1,2列处的一个2阶子式D
计算知, 所有的3阶子式D[注] Amn3222312300
0, 故r(A)2. , 若rankAm, 称A为行满秩矩阵;
若rankAn, 称A为列满秩矩阵.
Ann, 若rankAn, 称A为满秩矩阵(可逆矩阵, 非奇异矩阵);
若rankAn, 称A为降秩矩阵(不可逆矩阵, 奇异矩阵).
§3.2 矩阵的初等变换
1.初等变换
行变换
列变换
① 对调
rirj
cicj
② 数乘(k0)kri kci
③ 倍加 rikrj cikcj
Amn经过初等变换得到Bmn, 记作AmnBmn.
2.等价矩阵:若
有限次AmnBmn, 称
Amn与
Bmn等价, AmnBmn.
(1)自反性:AA
(2)对称性:AmnBmnBmnAmn
(3)传递性:AmnBmn,BmnCmnAmnCmn
定理1 AmnBmnrankArankB.
证
只需证明1次AmnBmnrankArankB.
设rankAr, 仅证行变换之(3)的情形:
记作
irkriijAjkjBj
(1)若rmin{m,n}, 则有
Dr1不含ri(B):D(B)r1Dr10(B)r1(A)
(A)(A)Dr1含ri(B), 不含r:DjjDr1kDr10倍加
D(B)r1含r, 且含r:Di(B)r1Dr10
(A)
故B中所有的r1阶子式D
kranBk, BAranArikrj(B)r10rankBrrankA
于是可得rankArankB.
(2)若rm或者rn, 构造矩阵
O O(m1)(n1)
AA1OOB, B1O(m1)(n1)Orikrj
由(1)可得A
1B1rankA1rankB1
ranAk1ranAkkranBkranAranBk1ranBk
其余情形类似.
例2 2A2131239636418212124, 求r(A).
360140440 解 0A01行61行4040, 故r(A)2.
1行最简形:A00行310010000123041行230000103230223B0
1标准形:A00行与列00H0
定理2 若rankA0行A0mnr(r0), 则
b1irb2irbrir00***B00b1i1b1i2b2i1
:行阶梯形
0行A0[i1]101[i2][ir]00100
***H00
:行最简形
定理3 形. 若rankAmnr(r0),Er则AOOO, 称为A的等价标准 推论1 若A 推论2 nn满秩, 则AEn.
. AmnBmnrankArankB§3.3 解线性方程组的消元法
例如
2x1x23x314x12x25x342x362x12x1x23x314x2x32x2x352x1x23x31x2x353x318(1)(2)(3)(4)(5)(6)(7)(8)(9)
(2)2(1)(3)(1)
x19x21x63
(5)4(6)(5)(6)
解线性方程组的初等变换:(1)互换两个方程的位置(2)用非零数乘某个方程
(3)将某个方程的若干倍加到另一个方程
用矩阵的初等变换表示方程组的求解过程如下:
A2b42200行12011035212行4060313141311010125001
91 6
11行50180
方程组: a11a21am1a12a22am2b
a1na2namnx1b1x2b2 xbnm 或者
Axb
~A
增广矩阵:A
设rankAr, 且A的左上角r阶子式D10行~A0000100000100b1,r1b2,r1br,r100b1nb2nbrn00r0, 则
d1d2drdr10: 行最简形
Axb的同解方程组为
x1b1,r1xr1b1nxnd1xb2,r1xr1b2nxnd22xbxr1brnxndrr,r1r0dr1
(3.4)
若d
若dr10, 0, 则方程组(3.4)无解:rankA~~r1rrankArankA
r1则方程组(3.4)有解:rankAr方程组(3.4)成为
„,xndn
(1)
(2)rn时, x1d1, x2d2,是其唯一解
rn时, 方程组(3.4)成为
x1d1b1,r1xr1b1nxnx2d2b2,r1xr1b2nxnxrdrbr,r1xr1brnxn
一般解为
x1d1b1,r1k1b1nknrxd2b2,r1k1b2nknr2xrdrbr,r1k1brnknrxk1r1knrxn40
其中k
定理4
(1)
(2)1,k2,,knr为任意常数.
Amn,~AAb
~Axb有解rankArankA;
Axb有解时, 若rankAn, 则有唯一解;
若rankAn, 则有无穷多组解.
定理5(1)
(2)
课后作业:习题三
1, 2, 3, 4
例3 求解Ax1~A21100行Amnx0有非零解rankAnAnnx0有非零解detA0;
.
b,341A2146224234120011046232221 0,b422583
242200310
解
141051行803020001052 2
51行1000~ranAkranAk24Axb有无穷多解
x122x2x4x31x4
同解方程组:x1x2一般解:x3x4
22k1k21k1k2k2
(k1,k2为任意常数)
例4 求解Ax~A111b,111A1111111111,1b210
100121111
解 11101行121100110001
1010100
111行1012行1111(1)11
(1)
1
x21x1同解方程组:x3(1)x1x(1)(2)x14x1x2一般解:x3x41(1)kkk
(k为任意常数)
(1)(2)k
(2)1
1(x2x3x4)
同解方程组:xx1x2一般解:x3x41
1k1k2k3k1k2k3
(k1,k2,k3为任意常数)
例5 讨论方程组Ax
解
计算可得
(1)
b何时有唯一解, 无穷多解, 无解?
1A121111,3b44
detA(1)
0且1:根据Cramer法则, 方程组有唯一解.
(2)
0: 1~A100111131行40400011011行1043030101104313
(3)
ranAk2,~rankA3,故方程组无解.
1且0: 1~A11101行2101111110131行404121行1040101010110031行1041201011214
1211212时, 时,~rankA3, rankA2, 故方程组无解. 故方程组有无穷多解.
~rankArankA23, §3.4 初等矩阵
定义
对单位矩阵进行一次初等变换得到的矩阵, 称为初等矩阵.
[注] 对单位矩阵进行一次初等列变换, 相当于对单位矩阵进行一次
同类型的初等行变换.因此, 初等矩阵可分为以下3类:
ErirjE01E10(i)ΔE(i,j)(j)E
1.43
2.EEkrikΔE[i(k)] E
(k0)
3.ErikrjEEcjkciE1Ek1(i)ΔE[i,j(k)] (j)E(i)ΔE[i,j(k)](j)E1Ek1
Amna11a21am1a12a22am21a1nia2njamnm,Amn1,,i,,j,,n
性质1
1jEm(i,j)Aim,1kiEm[i(k)]Ajm,iEm[i,j(k)]Akjjm1
因此可得:对A进行一次初等行变换, 相当于给A左乘一个
同类型的初等矩阵.(定理6的结论之一)
性质2
AEn(i,j)1,,j,,i,,n AEn[i(k)]1,,ki,,j,,n 44
AEn[i,j(k)]1,,i,,jki,,nB3
cjkciΔ
注意:AB3
因此可得:对A进行一次初等列变换, 相当于给A右乘一个
同类型的初等矩阵.(定理6的结论之二)
性质3
定理7 AnndetE(i,j)1,[E(i,j)]1E(i,j)
1detE[i(k)]k0, detE[i,j(k)]1,[E(i(k))]E[i(1k)]
[E(i,j(k))]1E[i,j(k)]
可逆A可以表示为有限个初等矩阵的乘积.
AEn,证
必要性.已知detA0, 则A满秩
P1,,Ps故存在初等矩阵
及Q1,,Qt, 使得,AP1PsQtQ11111PsP1AQ1QtEn1
而P与Q都是初等矩阵.
i1j
充分性.显然成立.
矩阵求逆方法之二(初等行变换法):
deAtnn0AP1P2Ps
(P都是初等矩阵)
i1111PsP2P1EAPsP2P1AE111
PsP2P1111AEEA1
由此可得:对n2n矩阵A
E
施行“初等行变换”,当前n列
(A的位置)成为E时,则后n列(E的位置)为A.
1
例6 1A2121332, 4213213010324求A.
1
解 A1E21100行1003140001001行00101001111231010010341121111010315111001001
213
1226501行100011行4030
100行0012651111143
故A12651431001a.
例7 1aA2a3a1aa21a01aa2, 求A.
1
解 A1aE2a3a***00001
依次作初等行变换
10E00010000100001r4ar3,r3ar2, r2ar1可得
1a0001a0001a
A00 01 46
故 A11a1a1a1.
定理8 设Amn,Bmn, 则AmmB , 使得PAQB
存在可逆矩阵P
证
必要性.已知A等
矩阵Q1和Qnn.
1B, 则存在m阶初等矩阵P,,Ps和n阶初,,Qt, 使得PsP1AQ1QtB, 令
则有PAQ
充分性.已知PAQ为
PP1,,Ps,QQ1,,Qt
BB. , 则由定理7知,P和Q都可以表示
有限个初等矩阵的乘积, 即
故P
第四章
向量组的线性相关性 §4.1 向量及其运算
1.向量:n个数a
1sPP1,,Ps,QQ1,,QtB
P1AQ1QtB, 也就是A. ,a2,,an构成的有序数组,47
记作(a1,a2,,an),称为n维行向量.
ai–– aiR称为向量的第i个分量
–– 称为实向量(下面主要讨论实向量)称为复向量 aiC––
零向量:(0,0,,0),a2,,an)
负向量:()(a 2.线性运算:
相等:若a
加法:i1(a1,a2,,an), (b1,b2,,bn)
bi(i1,2,,n), Δ称.
(a1b1,a2b2,,anbn)
Δ
数乘:k(ka
减法: 3.算律:
(1)
(2)Δ1,ka2,,kan)
()(a1b1,a2b2,,anbn)(a1,a2,,an), (b1,b2,,bn), (c1,c2,,cn)
(5)
(6)
1
()()k(l)(kl)
(3)
(4)
(7)
(8)
k()kk(kl)kl()
4.列向量:n个数a1,a2,,an构成的有序数组,a1a2记作an,或者(a1,a2,,an)T, 称为n维列向量.
00零向量:0
a1a2负向量:()an
5.内积:设实向量
(a1,a2,,an), (b1,b2,,bn),称实数
[,]a1b1a2b2anbn为与的内积.
算律:(a1,a2,,an), (b1,b2,,bn), (c1,c2,,cn)
(1)[,][,]
(2)[k,]k[,]
(k为常数)
(3)[
(4),][,][,]2 时,[,]0. 时, [,]0;(5)[,]证(5)[,][,]
tR, 由[t,t]0可得
[,]2[,]t[,]t20
04[,]4[,][,]0[,][,][,]
6.范数:设实向量, 称实数
性质:(1)
(2)
(3)
(4)
证(3)
[,]为的范数.
k时, 0;时, 0.
k
(kR)
2
[,][,]2[,][,]
证(4)
222
2,()
[,]() 7.夹角:设实向量,, 称
arccos(0)
为与之间的夹角.
正交:若[,]0, 称与正交, 记作.
(1)
(2)
单位化:若
§4.2 向量组的线性相关性
1.线性组合:对n维向量及
1,时, 2;
或时, 有意义, 而无意义.
1, 称0为与同方向的单位向量. ,,m, 若有数组k11,,km使得
k11kmm,,m, 称为,,m的线性组合,或可由例1 1110, 2111411线性表示.,4531,3311
判断可否由解
设4,2,3线性表示? k11k22k33,比较两端的对应分量可得
第二篇:数值线性代数课设课件资料
数值线性代数课程设计报告
姓名:陶英 学号:081410124
任课教师:杨熙
南京航空航天大学
2016 年 6 月 22日
求解线性方程组的三种迭代法及其结果比较
摘要
当今的环境下,数值计算越来越依赖于计算机。大规模科学计算和工程技术中许多问题的解决,最终归结为大型稀疏线性方程组的求解,其求解时间在整个问题求解时间中占有很大的比重,有的甚至达到80%。由于现今科学研究和大型项目中各种复杂的可以对计算精度和计算速度的要求越来越高。因此,作为大规模科学计算基础的线性代数方程组的高效数值求解引起了人们的普遍关注。这种方程组的求解一般采用迭代法。
关于迭代法,是有很多种解决公式的:Jacobi,G-S和超松弛迭代法。这三种方法的原理大致相同,Jacobi需要给定初向量,G-S则需要给定初值,超松弛法是对Guass-Seidel迭代法的加权平均改造。而本文则是对大型稀疏线性方程组迭代求解与三种迭代法(Jacobi,Gauss-Seidel和超松弛迭代法)的收敛速度与精确解的误差比较做出研究。
关键词:Jacobi迭代法;Gauss-Seidel迭代法;SOR迭代法;线性方程组 方法与理论的叙述
1.1迭代法简介
1.Jacobi迭代法:
对于非奇异线性方程组Ax=b,令A=D-L-U,其中
则原方程组可改写为:
其中
给定初始向量:
由(2.2)可以构造迭代公式:
其分量形式为:
2.2)(2.Guass-Seidel迭代法: 类似于Jacobi迭代法,给定初值:
令
则得到Guass-Seidel公式:
其分量形式为:
3.超松弛迭代法(SOR 迭代法):
SOR迭代法是对Guass-Seidel迭代法的加权平均改造,即
为Guass-Seidel迭代解,即
它的分量形式为:
其中ω称为松弛因子,当ω>1时称为超松弛;当ω<1时叫低松弛;ω=1时就是
Guass-Seidel迭代。
上述三种经典迭代法收敛的充分必要条件是迭代矩阵谱半径小于1。
谱半径不易求解,而在一定条件下,通过系数矩阵A的性质可判断迭代法的收敛性。定理1:
若系数矩阵A是严格对角占优或不可约对角占优,则Jacobi迭代法和Gauss-Seidel迭代法均收敛。定理2:
(1)SOR迭代法收敛的必要条件是0 (2)若系数矩阵A严格对角占优或不可约对角占优且0 2.1问题 考虑两点边值问题: d2ydya,0a12 dxdxy(0)0,y(1)1.1a(1e)ax 容易知道它的精确解为:y1/1ex为了将微分方程离散,把[0,1]区间n等分,令h=1/n,xiih,i1,2,...n1,得到差分方程 (h)yi1(2h)yiyi1ah2,从而得到迭代方程组的系数矩阵A。 对=1,a=1/2,n=100,分别用jacobi,G-S,超松弛迭代法分别求线性方程组的解,要求4位有效数字,然后比较与精确解的误差。 对=0.1,=0.01,=0.001,考虑同样问题。 1.方程的表示及存储 由于本题中线性方程组的系数矩阵为三对角矩阵,所以可以采用紧缩方法存储,即 然后在矩阵乘法时对下标处理一下即可。但是考虑到三种迭代方法的一般性,且本题中n=200并不是很大,所以实验中并没有采用紧缩存储,而是采用了直接存储。2.边值条件的处理 由于差分得到的方程组的第一行和最后一行中分别出现了边值y(0)与y(1)作为常数项,因此要在常向量的第一项和最后一项作一些修改: 3.迭代终止条件 首先确定要求的精度tol,我们希望当 则停止迭代。对于迭代格式,若 且,则迭代序列的 第k 次近似解和精确解之间有估计式由题目要求知我们需要有 。,而由上面的迭代估计,只要,即取为,因此最后令迭代终止条件为 即可。而本题中q可近似 4.SOR 迭代中最佳松弛因子的选取 由于SOR 迭代法的效果和其松弛因子w的选取有关,所以有必要选取合适的松弛因子。当选择最佳松弛因子 时,SOR 方法的迭代速度最快。 Matlab实现: 迭代矩阵是n-1阶的,不是n阶; 等号右端向量b的最后一项,不是ah^2,而是ah^2-eps-h 2.2精确解 1ay(1e)ax 1/1ex带入a=1/2,=1 代码: >> clear >> x=linspace(0,1);truy=(1-0.5)/(1-exp(-1/1))*(1-exp(-x./1))+x.*0.5;figure;plot(x,truy,'g','LineWidth',1.5);hold on;Grid 图: 2.3三种迭代法 Jacobi法:代码见附录 Eps=1 结果: 迭代次数k:22273 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) Eps=0.1 结果: 迭代次数k:8753 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) Eps=0.01 结果: 迭代次数k:661 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) G-S迭代法:代码见附录 Eps=1 结果: 迭代次数k:11125 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) Eps=0.1 结果: 迭代次数k:4394 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) Eps=0.01 结果: 迭代次数k:379 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) 超松弛法:代码见附录 Eps=1 w=1.56 结果: 迭代次数k:3503 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) Eps=0.1 w=1.56 结果: 迭代次数k:1369 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果) Eps=0.01 w=1.56 结果: 迭代次数k:131 结果与精确解的比较图(绿色粗线是精确解,黑色细线是迭代结果)分析讨论及心得体会 3.1三种方法的比较 Jacobi、G-S、超松弛法,三者都能够取得对精确解的良好逼近,但是,在相同的精度条件下,三者的收敛速度是不一样的,jacobi 3.2心得体会 这次课程设计,平时感觉挺简单的那些枯燥单调的代码和数学公式,真正到了自己运用的时候却无从下手,但是,解决问题的过程恰是不断学习的过程:数学算法转换为代码的过程要对题目有深入的了解,然后对程序函数定义还要有一定的掌握能力,通过这个的过程让我巩固了自己的数学知识,对数学专业知识和MATLAB的操作有了更深的体会。 课程设计中遇到的问题只凭自己苦思冥想是不能全部解决的,这是同学老师的建议和网络给了我很大的帮助。遇到自己解决不了的问题时,多多向老师同学请教,或许问题就能迎刃而解。 4参考文献 [1]徐树方.数值线性代数.北京:北京大学出版社,1995.[2]马昌凤.现代数值分析.北京:国防工业出版社.2013.[3]刘春凤,米翠兰.实用数值分析教程.北京冶金工业出版社.2006 5附录 源代码 1.Jacobi: function [y,k]=jacobi2(a,eps,h,delta)n=1.0/h;A=ones(n-1);y=zeros(n-1,1);z=zeros(n-1,1);k=0;for i=1:n-1 for j=1:n-1 A(i,j)=0;end end for i=1:n-1 A(i,i)=-(2*eps+h);end for i=1:n-1 for j=1:n-1 if i==j+1 A(i,j)=eps;end if i==j-1 A(i,j)=eps+h;end end end b=zeros(n-1,1);for i=1:n-2 b(i,1)=a*h^2;end b(n-1,1)=a*h^2-eps-h;D=zeros(n-1);for i=1:n-1 D(i,i)=A(i,i);end L=zeros(n-1);for i=1:n-1 for j=1:n-1 if i>j L(i,j)=-A(i,j);end end end U=zeros(n-1);for i=1:n-1 for j=1:n-1 if i end end B=D(L+U);g=Db;while 1 z=B*y+g;if norm(z-y,inf) y=z;k=k+1;end x=linspace(0,1);truy=(1-a)/(1-exp(-1/eps))*(1-exp(-x./eps))+x.*a;figure;plot(100*x,truy,'g','LineWidth',5);hold on;grid hold on;plot(y,'b') 2.G-S: function [y,k]=gs2(a,eps,h,delta)n=1.0/h;A=ones(n-1);y=zeros(n-1,1);z=zeros(n-1,1);k=0;for i=1:n-1 for j=1:n-1 A(i,j)=0;end end for i=1:n-1 A(i,i)=-(2*eps+h);end for i=1:n-1 for j=1:n-1 if i==j+1 A(i,j)=eps;end if i==j-1 A(i,j)=eps+h;end end end b=zeros(n-1,1);for i=1:n-2 b(i,1)=a*h^2;end b(n-1,1)=a*h^2-eps-h;D=zeros(n-1);for i=1:n-1 D(i,i)=A(i,i);end L=zeros(n-1);for i=1:n-1 for j=1:n-1 if i>j L(i,j)=-A(i,j);end end end U=zeros(n-1);for i=1:n-1 for j=1:n-1 if i end end B=D(L+U);g=Db;while 1 z=(D-L)U*y+(D-L)b;if norm(z-y,inf) y=z;k=k+1;end x=linspace(0,1);truy=(1-a)/(1-exp(-1/eps))*(1-exp(-x./eps))+x.*a;figure;plot(100*x,truy,'g','LineWidth',5);hold on;grid hold on;plot(y,'b') 3.SOR: function [y,k]=sor(a,eps,h,delta,w)n=1.0/h;A=ones(n-1);y=zeros(n-1,1);z=zeros(n-1,1);k=0;for i=1:n-1 for j=1:n-1 A(i,j)=0;end end for i=1:n-1 A(i,i)=-(2*eps+h);end for i=1:n-1 for j=1:n-1 if i==j+1 A(i,j)=eps;end if i==j-1 A(i,j)=eps+h;end end end b=zeros(n-1,1);for i=1:n-2 b(i,1)=a*h^2;end b(n-1,1)=a*h^2-eps-h;D=zeros(n-1);for i=1:n-1 D(i,i)=A(i,i);end L=zeros(n-1);for i=1:n-1 for j=1:n-1 if i>j L(i,j)=-A(i,j);end end end U=zeros(n-1);for i=1:n-1 for j=1:n-1 if i end end B=D(L+U);g=Db;Lw=((D-w*L)^-1)*((1-w)*D+w*U);while 1 z=Lw*y+w*(D-w*L)^-1*b;if norm(z-y,inf) y=z;k=k+1;end x=linspace(0,1);truy=(1-a)/(1-exp(-1/eps))*(1-exp(-x./eps))+x.*a;figure;plot(100*x,truy,'g','LineWidth',5);hold on;grid hold on;plot(y,'b') 线性代数学习心得 各位学友好! 首先让我们分析一下线性代数考试卷(本人以1999年上半年和下半年为例) 我个人让为,先做计算题,填空题,然后证明题,选择题等(一定要坚持先易后难的原则,一定要。旁边有某些同志说:“这些都是屁话,我们都知的快快转入正题吧!”) 把选择题第8题拉出来让大家看看 n(n>1)阶实对矩阵A是正定矩阵的充份必要条件是() A.A是正定二次型f(x)=x(A)x的矩阵 B.A是各阶顺序主子式均大于等于零(书本的p231定5.9知,大于零就可以了,明显也是错的) C.二次型f(x)=xTAx的负惯性指数为零 D.存在n阶矩阵C,使得A=CTC(由书本的P230知,存在非奇异N阶矩阵C,使A=CTC)很明显,这个选择是错了) 各位学友在做选择题时要仔细呀! 证明题 先讲1999年下半年 设A,B,C均为n阶矩阵,若ABC=I,这里I为单位矩阵,求证:B为可逆矩阵,且写出的逆矩阵? 证的过程:己知ABC=I,|ABC|=|I|不等于零,|A|*|B|*|C|不等于零,得出|B|不等于零。所以B是可逆矩阵。 求其逆矩阵,ABC=I,两边同时右乘C-1得AB=C-1,接下来左乘以A-1得B=A-1C-1,最后BC=A-1,BCA=I,于是得B-1=CA(不知各位学友有没有更简便的方法谢谢告之) 对这题做后的心得,本人认为一定要记得,a逆阵可逆的充分必要条件是行列式|a|不等零(切记,还有如ab=i,那么a-1=b) 对了还有,在求解逆矩阵,最简单方法是用初等行变换 公式法吗!容易出错,只适合求解比较特殊的 下面这些是相关的证明题 设B矩阵可逆,A矩阵与B矩阵同阶。且满足A2+AB+B2=O,证明A和A+B都是可逆矩阵?(相信大家都能做出) 己知i+ab可逆,试证I+BA也可逆? 接下来看看1999年上半年的 设n阶方阵A与B相似,证明:A和B有相同的特征多项式? 应搞清楚下面的概念 什么是特征多项式呢(1) 什么是特征值呢(2) 什么还有特征向量(3) 什么是相似矩阵(4) λI-A称为A的特征矩阵;|λI-A|称为A的特征多项式;|λI-A|=0称为A的特征矩阵,而由些求出的全部根,即为A的全部特征值。 对每一个求出特征值λ,求出齐次方程组(λI-A)x=o的基础解是&1,&2,&3...&s,则k1&1+k2&2+...ks&s即是A对应于 λ的全部特征向量(其中,k1...ks不全为零) 相似矩阵:设A,B都是n阶方阵,若存在n阶可逆阵p,使得p-1ap=b,则称A相似于B,记为A~B(相拟矩阵有相同的行列式,相同的秩,相同的特征值) 我觉得有这么一题使终我还是一知半解的,拉出来让大家看看: 设A为4阶方阵,A*为A的伴随矩阵,若|A|=3,则|A*|=?,|2A*|=? 这题答案是27,432 怎么算的呢?这个具体我也不太清楚,我是用自己的方法,|A|N-1=|A*|,这个N代表多少阶,如是4阶那么3^3=27,后面那个,切记:把2提出行列式以外,看A是几阶行列式,4阶就提4次,2^4*3^3=432(可能书上不是这样的,我只是根据其习题答案推论出来的) 应注意的问题:区为行列式和矩阵之间的区别,特别是用一个不为零的数K乘以行列式或矩阵,前者只是乘以某一行或列,后者则是每一个元素都要乘! 很容易搞不零清的:线性相关或无关和什么情况下线性方程组有解或无解,还有什么极大无关组,基础解系,特征值,多项式,特征向量,相似矩阵有哪些性质,正交矩阵的充分心要条件,二次型化成标准型。 厦门理工学院继续教育学院20 第 学期期末试卷 线性代数(考试时间:120分钟) 专业 姓名 层次形式 成绩 一、选择题(每小题4分,共16分)1.A,B为三阶方阵,矩阵X满足AXABXBBXAAXBE则().22111(A)X(AB);(B)X(AB)(AB)(C)X(AB)(AB)(D)以上答案都不对.2.11; A、B、C为n阶方阵,且ABC,A、B、C的列向量组分别为1,2,,n;1,2,,n(A); 1,2,,n.若 1,2,,n线性相关,则().1,2,,n线性相关;(B) 1,2,,n线性相关; (C)(A)与(B)都成立;(D)(A)或(B)成立.3.设A,B为三阶矩阵,且r(A3A2E)3,若r(B)2则r(ABB)().(A)1 ;(B)2; (C)3;(D)无法判断. A22334.设三阶矩阵 B22,3,其中,,2,3均为三维行向量,已知A18,2B2,则AB().(A)1 ;(B)2; (C)3;(D)4.二、填空题(每小题4分,共16分) En10ABOB为n阶非零矩阵,5.设A、,且A的阶梯形为1Da1111b1111c1111n00,则矩阵B的秩=.6.已知,则此行列式的所有代数余子式之和i,j1Aij.1 1A0Tx(1,1)7.已知是1a的一个特征向量,则a.8.为已知A是3阶方阵,1,2,3是三维线性无关的向量.若A112,A223,A313,则A的行列式等于.三、计算下列各题(每小题7分,共28分) 01D1110111110111110111110.9.计算n阶行列式 10.若二次型 1f(x1,x2,x3)2x18x2x32ax1x2222正定,求a的取值范围.411.已知(1,1,1),(1,0,1),且A.求A.TTT 2A02 030110B002010000 12.已知矩阵X满足AX2BBA2X,求X. 四、解答下列各题(每小题14分,共28分) 2x13x23x3ax1x2x313x4x2(a2)x3a1x2xax12313.求a使方程组1与1有公共解,并求公共解.14.已知二次型 f(x1,x2,x3)XAXx1x32ax1x22x1x32bx2x3T22的秩为2,Tf(x1,x2,x3)(1,1,1)是A的特征向量.(1)求a,b的值;(2)求经正交变换所得的标准型,并写出相应的正交矩阵.3 五.解答下列各题(每小题4分,共12分) 15.设1,2,,t是线性方程组AxO的基础解系,向量满足AbO.证明1,2,,t,线性无关.16.已知A是n阶方阵且可对角化,问BAAE可否对角化?证明你的结论.2 T17.已知A为n阶矩阵.证明方程组AxO与AAxO的解相同. 线性代数试题 请考生按规定用笔将所有试题的答案涂、写在答题纸上。 说明:在本卷中,AT表示矩阵A的转置矩阵,A*表示矩阵A的伴随矩阵,E是单位矩阵,|A|表示方阵A的行列式,r(A)表示矩阵A的秩。 选择题部分 一、单项选择题(本大题共5小题,每小题2分,共10分) 在每小题列出的四个备选项中只有一个是符合题目要求的,请将其选出并将“答题纸”的相应代码涂黑。错涂、多涂或未涂均无分。1.设行列式A.-3 C.1 2.设4阶矩阵A的元素均为3,则r(A)= A.1 C.3 3.设A为2阶可逆矩阵,若A1B.2 D.4 a1a2b1acabc1,112,则111 b2a2c2a2b2c2B.-1 D.3 13A. 2553C. 21A.r=m时,Ax=0必有非零解 C.r ,则A= 251B.25D.23 53 14.设A为m×n矩阵,A的秩为r,则 B.r=n时,Ax=0必有非零解 D.r 2225.二次型f(xl,x2,x3)=x12x23x38x1x312x2x3的矩阵为 1A.081C.04 08212 1230426 631B.001D.4008212 034026 63═══════════════════════════════════════════════════════════════════════════════ 2 非选择题部分 注意事项: 用黑色字迹的签字笔或钢笔将答案写在答题纸上,不能答在试题卷上。 二、填空题(本大题共10小题,每小题2分,共20分)6.设A为3阶矩阵,且|A|=2,则|2A|=______. 7.设A为2阶矩阵,将A的第1行加到第2行得到B,若B=8.设矩阵A=12,则A=______.34a12a11a12a11,B=,且r(A)=1,则r(B)=______.a21a22a11a21a12a229.设向量α=(1,0,1)T,β=(3,5,1)T,则β-2α=________. 10.设向量α=(3,-4)T,则α的长度||α||=______. 11.若向量αl=(1,k)T,α2=(-1,1)T线性无关,则数k的取值必满足______.12.齐次线性方程组xl+x2+x3=0的基础解系中所含解向量的个数为______. 12210013.已知矩阵A=212与对角矩阵D=010相似,则数a=______ 22100a14.设3阶矩阵A的特征值为-1,0,2,则|A|=______. 22215.已知二次型f(x1,x2,x3)=x1正定,则实数t的取值范围是______. x2tx 3三、计算题(本大题共7小题,每小题9分,共63分) abc16.计算行列式D=2a2a2b2cbac2b.2ccab17.已知向量α=(1,2,k),β=1,,且βαT=3,A=αTβ,求(1)数k的值;(2)A10. 11231231218.已知矩阵A=231,B=00,求矩阵X,使得AX=B.3401019.求向量组α1=(1,0,2,0)T, α2=(-1,-1,-2,0)T, α3=(-3,4,-4,l)T, α4=(-6,14,-6,3)T的秩和一个极大线性无关组,并将向量组中的其余向量由该极大线性无关组线性表出. 2x3yz020.设线性方程组2xyz1,问: xyz1═══════════════════════════════════════════════════════════════════════════════ 3(1)λ取何值时,方程组无解? (2)λ取何值时,方程组有解?此时求出方程组的解. 00121.求矩阵A=010的全部特征值与特征向量. 1002222.用配方法化二次型f(x1,x2,x3)=2x12x24x1x38x2x3为标准形,并写出所用的可逆线性变换. 四、证明题(本题7分) 23.设向量组α1,α2线性无关,且β=clα1+c2α2,证明:当cl+c2≠1时,向量组β-α1,β-α2线性无关. ═══════════════════════════════════════════════════════════════════════════════第三篇:线性代数学习心得
第四篇:线性代数试卷
第五篇:线性代数试卷