第一篇:Matlab 与线性代数教案
Matlab 与线性代数
一、Matlab 入门:
1.启动、退出、运行: 2.窗口介绍: 3.基本符号: =:赋值符号
[ ]:数组定义符号 , 区分列 函数参数分隔符;区分行 取消运行显示 % 注释标记
: 具有多种应用功能
4.matlab的变量(区分大小写): 预定义变量: ans
pi 相关命令: format(显示格式 rat long short)
who whos clear
5.M 文件(纯文本文件,扩展名为.m)建立 修改 保存 运行
二、Matlab 与线性代数的基本运算
1.矩阵的输入
数字矩阵:A=[1 2 3;3 2 1]
或 A=[1, 2, 3;3, 2, 1] 或 A=[1 2 3 2 1]
符号矩阵(显示出来元素之间有逗号): 定义符号变量 sym syms
用法:(1).sym(‘[a,b,c;b,c,a]’)或 sym(‘[a b c;b c a]’)
(2).syms a b c
A=[a b c;b c a]
2.产生特殊矩阵的函数:
zeros(m,n)zeros(n)
ones(m,n)ones(n)eye(n)
magic(n)rand(m,n)randn(n)% 产生(0,1)区间均匀分布的随机矩阵
3.相关命令:
round(A)% 表示对矩阵A中所有元素进行四舍五入 length(A)% 返回A的长度(列数)size(A)% 返回A的尺寸,行数 列数 A(i,j)% 引用矩阵A的第i行第j列元素
4.矩阵的基本运算
(1).+-*.*
(2).转置 A’
(3).方阵的幂:A^3
5.求向量组的极大无关组
A[1,2,3 ]
(1).U=rref(A)% U为A的行最简形
(2).[U,s]=rref(A)% U为A的行最简形, s为首非零元所在列组成的向量
(3).rrefmovie(A)% 返回A的行最简形,且给出每一步化简过程
6.求线性方程组的解
情形1。Ax=b,其中A为n阶可逆阵
法1: x=inv(A)*b 或 x=A^(-1)*b
法2: U=rref([A,b])% 返回值U为矩阵的行最简形,最后一列即为解x。
情形2。Ax=0, 其中A 为m*n 矩阵,R(A)=r 法1:U=rref(A), 选定自由变量,得到一组基础解系 法2:z=null(A) % z的列向量为Ax=0的一组标准正交基。 情形3。Ax=b, 其中A 为m*n 矩阵, 求通解 U=rref([A,b])从最后一列找特解,前n列找导出组的基础解系,然后按格式写 出Ax=b的通解。(或先写出以U为增广矩阵的同解方程组也可。) 6x13x22x33x44x554x12x2x32x43x54 例子: .4x2x3x2xx0234512xx7x3x2x112345(4).(5).(6).(7).方阵行列式 det(A)方阵的秩 rank(A)方阵的逆 inv(A)或 A^(-1)矩阵的除法 左除 右除/ AB=C 则 A=C/B B=AC 输入:A=[6 3 2 3 4;4 2 1 2 3;4 2 3 2 1;2 1 7 3 2]; b=[5 4 0 1]’; U=rref([A,b])10得到:U001/2000010000103/417/203/22 603/20取x2,x5为自由变量,令x20,x50得Ax=b的特解*2 60 1/23/410x210分别令和得导出组的基础解系为:10,21 x50107/2013x24x5x112x3x5或:导出组Ax=0的同解方程组:,x2,x5为自由变量,分别令x47x521/23/410x21,x50和x20,x51得导出组的基础解系为:10,21。 07/2017.求矩阵的特征值与特征向量 (1).d=eig(A)% d为矩阵A的特征值构成的向量 (2).[V,D]=eig(A)% D为A 的特征值构成的对角阵,V 的列为A的单位特征向 量,与D中的特征值对应,满足:AVDV8.Schmidt 正交化方法 B=orth(A)% B的列向量为A的列空间的一组标准正交基,换句话说,B的列是 A的列向量的正交标准化,满足B*Beye(rank(A))。 9.用正交变换化二次型为标准形 先写出所给二次型的矩阵A,则A为实对称矩阵,[V,D]=eig(A)% D 为A的特征值构成的对角阵,V的列向量为A的正交单位特征 向量,次序与D的元素对应。满足VAVDVT1'1,即AVVD。 AV。 《数学实验》在线习题3 Matlab程序设计部分 一.分析向量 组 a1[1T2a23],T[a31T2,0],a4[121]T,a5[246]T的线性相关性,找出它们的最大无关组,并将其余向理表示成最大无关组的线性组合。 解,a1=[1 2 3]'; a2=[-1-2 0]';a3=[0 0 1]';a4=[1-2-1]';a5=[2 4 6]';A=[a1,a2,a3,a4,a5];[R,S]=rref(A)r=length(S) R = 1.0000 0 0.3333 0 2.0000 0 1.0000 0.3333 0 0 0 0 0 1.0000 0 S = 4 r = 线性相关 a1,a2,a3,a4,a5 最大无关组是a1,a2,a4 其余向量的线性组合是a3=1/3a1+1/3a2 a5=2a1 二.计算行列式 x13D4x23x33x43x12y1x22y2x32y3x42y4x1y12x2y22x3y32x4y42y13y23y3323的值。其中1解,syms x1 x2 x3 x4 y1 y2 y3 y4 xxxy43 x42357,y1y2y3y44567。 D=[x1^3 x1^2*y1 x1*y1^2 y1^3;x2^3 x2^2*y2 x2*y2^2 y2^3;x3^3 x3^2*y3 x3*y3^2 y3^3;x4^3 x4^2*y4 x4*y4^2 y4^3];d=det(D)x1=2;x2=3;x3=5;x4=7;y1=4;y2=5;y3=6;y4=7;eval(d) d = ans = 153664 三.已知向量a1,1,0,b1,0,1,求向量a与b的夹角的度数。解,a=[1-1 0];b=[-1 0-1]; x=a.*b;x1=sum(x,2);x2=norm(a);x3=norm(b);y=x1/(x2*x3)y1=acos(y)y = -0.5000 y1 = 2.0944 四.已知线性方程组 clear 2x1x23x32x409xx14x2x112343x12x25x34x414x15x27x310x42,求系数矩阵的秩和方程组的通解。 a=[2-1 3 2;9-1 14 2;3 2 5-4;4 5 7-10];b=[0 1 1 2]';[R,s]=rref([a,b]);[m,n]=size(a);x0=zeros(n,1);r=length(s);x0(s,:)=R(1:r,end); disp('非齐次线性方程组的特解为:') x0 disp('对应的线性方程组的基础解系为:')x=null(a,'r') 非齐次线性方程组的特解为: x0 = 0.1429 0.2857 0 0 对应的齐次线性方程组的基础解系: x = -1.5714 0-0.1429 2.0000 1.0000 0 0 1.0000 则方程组的通解为: x1x2x412x2x322x3xxx0234五.求齐次方程组1的通解。 clear a=[-1 1 0 1;0 2 1 0;2 3-1-1];b=[1 2 0]'; [R,s]=rref([a,b]);[m,n]=size(a);x0=zeros(n,1);r=length(s); x0(s,:)=R(1:r,end); disp('非齐次线性方程组的特解为:')x0 disp('对应的线性方程组的基础解系为:')x=null(a,'r') 非齐次线性方程组的特解为: x0 = -0.4286 0.5714 0.8571 0 对应的齐次线性方程组的基础解系: x = 0.8571-0.1429 0.2857 1.0000 232A36112115,求正交矩阵P及对角形矩阵B,使P1APB。六.clear a=[2 3-2;3 6 11;-2 11 5];[v,d]=eig(a)v = -0.3684 0.9280 0.0562 0.6512 0.2144 0.7280-0.6635-0.3047 0.6833 d = -6.9057 0 0 0 3.3500 0 0 0 16.5556 七.求下列向量的秩和最大无关组,并将其余向量用该最大无关组线性表出: 11,2,1,324,1,5,631,3,4,7a1=[1 2 1 3]';a2=[4-1-5-6]';a3=[1-3-4-7]';A=[a1,a2,a3];[R,S]=rref(A)r=length(S)R = 1.0000 0-1.2222 0 1.0000 0.5556 0 0 0 0 0 0 S = r = 最大线性无关组为:a1 a2 a3=-1.2222a1+0.5556a2 八.判断方程组否有解,如果有,求其通解: x1x23x3x413x1x23x34x44x5x9x8x02341 clear a=[1 2-3-1;3-1-3 4;1 5-9-8];b=[1 4 0]';[R,s]=rref([a,b]);[m,n]=size(a);x0=zeros(n,1);r=length(s);x0(s,:)=R(1:r,end); disp('非齐次线性方程组的特解为:')x0 disp('对应的线性方程组的基础解系为:')x=null(a,'r') 非齐次线性方程组的特解为: x0 = 1.5000 0 0.1667 0 对应的线性方程组的基础解系为: x = -2.5000 0 -1.1667 1.0000 a112,a2021,求两向量的点积(数量积)和叉积(向九.已知向量1量积),以及它们之间的夹角的大小。 a1=[1 1 2]';a2=[0 2 1]'; TTy1=norm(a1);y2=norm(a2);y3=dot(a1,a2);y=y3/(y2*y3);c=acos(y)c*180/pi c = 1.1071 ans = 63.4349 十.计算行列式: 1x1y1D1x1y21x1y31x1y41x2y11x2y21x2y31x2y41x3y11x3y21x3y31x3y41x4y11x4y21x4y31x4y4 的值。其中x1x2x3x42357,y1y2syms x1 x2 x3 x4 y1 y2 y3 y4 D=[1+x1*y1 1+x1*y2 1+x1*y3 1+x1*y4;1+x2*y1 1+x2*y2 1+x2*y3 1+x2*y4;1+x3*y1 1+x3*y2 1+x3*y3 1+x3*y4;1+x4*y1 1+x4*y2 1+x4*y3 1+x4*y4]; x1=2;x2=3;x3=5;x4=7;y1=4;y2=5;y3=6;y4=7;d=det(D);eval(d) ans = 0 十一.y3y44567。 a1122,a20215,a32051,分析向量组1TTTTa43386的线性相关性,找出它们的最大无关组,并将其余向量表示成最大无关组的线性组合。 a1=[1 1 2 2]';a2=[0 2 1 5]';a3=[2 0 5-1]';a4=[3 3 8 6]';A=[a1,a2,a3,a4];[R,S]=rref(A)r=length(S)R = 0 0 0 0 0 0 0 0 0 0 S = r = 最大线性无关组为:a1 a2 a3; a4=a1+a2+a3 十二.求解五阶方程组 注:在系数矩阵中没有数据的地方,矩阵元素均为零。 a=[4 1 0 0 0;1 4 1 0 0;0 1 4 1 0;0 0 1 4 1;0 0 0 1 4];b=[2 1 1 1 2]';inv(a)*b 41x12141x12141x31141x41142 x5 ans = 0.4808 0.0769 0.2115 0.0769 0.4808 第一章 线性方程组的消元法与矩阵的初等变换 教学目标与要求 1.了解线性方程组的基本概念 2.掌握矩阵的三种初等变换 教学重点 运用矩阵的初等变换解一般的线性方程组 教学难点 矩阵的初等变换 §1.1 线性方程组的基本概念 一、基本概念 定义:m个方程n个未知数的线性方程组为如下形式: a11x1a12x2a1nxnb1axaxaxb2112222nn(1)am1x1am2x2amnxnbm称(1)为非齐次线性方程组;当b1b2bm0时则称为齐次线性方程组。方程组(1) a12a22am2a1na2n为系amna11a21TA的一个解为:x(c1,c2,,cn)(或称为解向量);此时称am1a11a12a1na21a22a2n数矩阵,称Bam1am2amn 二、线性方程组的消元法 b1b2为增广矩阵。bm2x1x23x31例1:解线性方程组4x12x25x34 2x2x6312x1x23x312x1x23x312x1x23x31解:4x2x32,x2x35,x2x35; xx54xx23x18232332x1x23x312x1x2192x118x19 x2x35,x21,x21,x21 x6x6x6x63333从上面可以看出,整个消元过程和回代过程都只与x1,x2,x3的系数有关,且仅用了以下3种变换:①交换两行;②某行乘k倍;③某行乘k倍加至另一行(即初等行变换)。 故我们隐去x1,x2,x3,,得到一个数字阵(即矩阵B),对B进行初等行变换: 213121312131B425404120115 2026011504121213121019213011501150101 003180016001620018100901010101 0016001612131009其中0115称为行阶梯形矩阵,0101称为行最简形矩阵。 003180016 三、小结 例1告诉我们求解一般的线性方程组的基本方法:对其增广矩阵B进行3种初等行变换,把它变为行阶梯形矩阵,再最终变成行最简形矩阵,然后从中读出所需的解。 四、一般解和通解 x12x2x32x41例2:解方程组2x14x2x3x45 x2x2xx42341解: 212112121121211B241150033300333 12214003330000012121120120011100111 0000000000即x12x2x42x122x2x4,亦即一般解为,其中x2,x4为自由未知量。 x3x41x31x4x122c1c2xc21令x2c1,x4c2,得方程组的通解为 x31c2x4c2注意:自由未知量的取法并不唯一。 a11x1a12x2a1nxn0axaxax02112222nn2、定理:在齐次线性方程组中,若mn(即方程 am1x1am2x2amnxn0的个数小于未知数的个数),则它必有非零解。 五、习题 P11 T1(2) T2 §1.2 矩阵的初等变换 一、矩阵及其初等变换 1、定义:称由mn个数aij(i1,2,,m;j1,2,,n)排成的m行n列的数表 a11a21Aam1a12a22am2a1na2n为矩阵,简记为A(aij)mn。amn 二、矩阵的初等行(列)变换 ①交换两行(列); ②某行(列)乘k倍; ③某行(列)乘k倍加至另一行(列)。 三、矩阵的标准形 定理:任意一个mn的矩阵A,总可以经过初等变换(包括行变换和列变换)化为如10下的标准形:F00000001000Er0100即AmnFO00000000O O其中1的个数r就是行阶梯形矩阵中非零行的行数。 四、习题 P18 T1(4)(5) T2(1) T3 P19 总复习题:T3 T4 第二章 行列式 教学目标与要求 1.会用对角线法则计算二阶行列式和三阶行列式 2.理解排列、逆序数的概念,掌握n阶行列式的定义及其重要性质 3.理解并会灵活运用行列式的展开公式,掌握范德蒙德行列式的结论 4.掌握克拉默法则及其应用 教学重点 1.n阶行列式的重要性质 2.n阶行列式展开公式的运用以及范德蒙德行列式的结论 3.克拉默法则的运用 教学难点 1.n阶行列式的重要性质及其展开公式 2.克拉默法则的运用 §2.1 二阶和三阶行列式 一、二阶行列式 a11x1a12x2b1a11a12 1、引例:对于线性方程组(1),其系数矩阵为A a21x1a22x2b2a21a22 用消元法解得 (a11a22a12a21)x1b1a22b2a12(2) (a11a22a12a21)x2b2a11b1a21a12a11a22a12a21称为二阶行列式,记DAdetA a12a11b1,D2 a22a21b22、定义:Da11a21a22a11a12b1Dx1D1那么(2)可以表示为,其中D,D1aab2DxD212222从而x1 二、三阶行列式 D1D,x22。DDa11x1a12x2a13x3b1a11a12axaxaxb1、定义:对于三元线性方程组211a222222332,记Aa21axaxaxba331a32311322333a11称DAdetAa21a13 a23,a33a12a22a32a13a23a11a22a33a12a23a31a13a21a32 a33a 31a11a23a32a12a21a33a13a22a31 为三阶行列式。 a112、三对角线法则(记忆):Da21a12a22a32a13a11a23a21a33a31a12a22 a32a31 三、习题 P25 T1(2)(3)(5) T2 T3 §2.2 n阶行列式的定义和性质 一、排列与逆序数 1.定义1:由1,2,,n组成的一个有序数组称为一个n级排列。(n级排列共有n!个)定义2:在一个排列中,如果一对数的前后位置与大小顺序相反,即前面的数大于后面的数,那么它们就称为一个逆序,一个排列中逆序的总数称为这个排列的逆序数,记作。)402107(奇排列)例:(25431;)1412108(偶排列) (5243。 定理:对换改变排列的奇偶性;在全部n级排列中,奇、偶排列的个数相等,各有 二、n阶行列式的定义 n!个。21.定义:n阶矩阵A(aij)nna11a21am1a12a22am2a1na2n,则n阶行列式定义如下: amna11 DAa12a1np1p2pna21an1a22a2nan2ann(1)(p1p2pn)a1p1a2p2anpn 这里,表示对1,2,,n这n个数的所有排列p1p2pn求和。即n阶行列式是指n!项取自不同行不同列的n个元素乘积的代数和。 2、例:(常用结论) a11(1) a11a22ann0a11a22ann0n(n1)2a12a1na11000 a22a2na210annan1a22an2ann1(2)2(1)12n n3、n阶行列式的等价定义 定理:D12(1)ai1j1ai2j2ainjn;其中1为行标排列i1i2in的逆序数,2为列标排列j1j2jn的逆序数。 三、行列式的性质 设n阶矩阵A(aij)nn的行列式为DA,则D有如下性质: T①AA; ②交换两行(列),则D变号; ③提公因子:某行(列)所有元素的公因子可以提到D的外面。 特别地,若某行(列)为0,则D0;若某两行(列)成比例,则D0。④拆和:若D中某行(列)的元皆为两项之和,则D等于两个行列式的和。⑤某行(列)乘k倍加至另一行(列),则D不变。 123例:②如211111211234234;③如33932113 ***123123④如456123333; 112112112111111111111⑤如23340120120120 45345012000 注意:计算行列式的常用方法:(1)利用定义; (2)利用性质把行列式化为上(下)三角形行列式,从而算得行列式的值;(3)利用展开公式(下一节)。 四、习题 P36 T1 T4 T5(3)(4)(8) T6(1) §2.3 行列式的展开公式 一、余子式与代数余子式 1、定义:在n阶行列式det(aij)中,划去元aij所在的第i行和第j列的元后,剩下的元按原来的顺序所构成的n1阶行列式称为aij的余子式,记作Mij;又记Aij(1)ijMij,称Aij为aij的代数余子式。 142.如:***中,a111的余子式为M11412,代数余子式为 23411234A11(1)11M11M11,a214的余子式为M21412,代数余子式为 341A21(1)21M21M21,二、展开公式 定理:n阶行列式等于它的任一行(列)的各元素与其对应的代数余子式的乘积之和。即可按第i行展开 Dai1Ai1ai2Ai2ainAin(i1,2,,n) 或可按第j列展开 Da1jA1ja2jA2janjAnj(j1,2,,n) 14如:3221433214431A112A123A134A141A114A213A312A41 21 2、讲解P42例2和例3 三、范德蒙德行列式 1x1Dnx12x1n1 1x22x2n1x21x32x311xn2xn1ijn(xjxi) n1n1x3xn推论:行列式某行(列)的元素与另一行(列)的对应元素的代数余子式乘积之和等于零。即 ai1Aj1ai2Aj2ainAjn(ij)或 a1iA1ja2iA2janiAnj(ij) 11例证:如3222433314441A112A123A134A14a21A11a22A12a23A13a24A140 21四、习题 P46 T2(3)(4)(5) §2.4 克拉默法则 一、克拉默法则 定理1:含有n个未知数x1,x2,,xn与n个方程的线性方程组 a11x1a12x2a1nxnb1axaxaxb2112222nn 2 (1) an1x1an2x2annxnbn 称(1)为非齐次线性方程组;当b1b2bn0时称为齐次线性方程组。 如果线性方程组(1)的系数行列式DA0(这里A(aij)nn),那么(1)有唯一解,且解为xjDjD(j1,2,,n),其中Dj(j1,2,,n)是把D中第j列元素用方程组右端的常数项替代后所得到的n阶行列式。 推论: (1)如果线性方程组(1)无解或至少有两个不同的解,那么它的系数行列式D0。 (2)如果齐次线性方程组的系数行列式D0,那么它只有零解;如果齐次线性方程组有非零解,那么它的系数行列式D0。 注意:用克拉默法则解线性方程组的两个条件:①方程个数等于未知数个数;②系数行列式不等于零。克拉默法则的意义主要在于建立了线性方程组的解和已知的系数以及常数项之间的关系。它主要适用于理论推导。 二、习题 P50 T2 T3 ; P51 总复习题:T1 T2 T3 T6 第三章 矩阵 教学目标与要求 1.理解矩阵的概念,掌握矩阵的3种运算(加法、数乘、乘法),以及它们的运算律 2.熟记几种特殊矩阵(单位阵、对角阵、数量矩阵、三角阵、转置矩阵、对称和反对称阵)及其性质,掌握方阵行列式的性质 3.掌握伴随矩阵和逆矩阵的定义及其性质,熟悉逆矩阵的运算规律 4.了解分块矩阵的运算律,以及常用结论 5.理解初等矩阵与初等变换之间的关系,掌握初等变换求逆矩阵的方法 6.掌握矩阵的秩的概念及其性质,会用初等变换求矩阵的秩 教学重点 1.矩阵乘法的运算律和方阵行列式的性质 2.逆矩阵和伴随矩阵的运算性质,以及初等变换法求逆矩阵 3.矩阵的秩的性质,以及初等变换法求矩阵的秩 教学难点 1.逆矩阵的概念,以及求逆的方法 2.矩阵的秩的概念,以及求秩的方法 §3.1 矩阵的概念及其运算 一、矩阵的概念 1、定义:称由mn个数aij(i1,2,,m;j1,2,,n)排成的m行n列的数表 a11a21Aam1a12a22am2a1na2n为矩阵,简记为A(aij)mnAmn。amn矩阵的相等:AmnBmnaijbij(i1,2,,m;j1,2,,n) b1b2行矩阵(行向量):A(a1,a2,,an);列矩阵(列向量):A bn 二、矩阵的运算 1、矩阵的加法 定义1:设A(aij)mn,B(bij)mn,则AB(aijbij)mn 注意:两个矩阵是同型矩阵时才能进行加法运算。 矩阵的加法满足下列运算律(设A,B,C都是mn矩阵):(1)交换律:ABBA; (2)结合律:(AB)CA(BC)(3)负矩阵A(A)0,规定减法运算:ABA(B) 2、矩阵的数乘 a11a21定义2:数与矩阵A的乘积记作A或A,规定为Aam1a12a1na22a2nam2amn; 矩阵的数乘满足下列运算律(设A,B都是mn矩阵,,为数):(1)()A(A); (2)()AAA;(3)(AB)AB; (4)1AA;(5)A00或A0 3、矩阵的乘法 定义3:设A(aij)ms,B(bij)sn,那么矩阵A与矩阵B的乘积是一个mn矩阵C(cij)mn,其中 cijai1b1jai2b2jaisbsjaikbkj(i1,2,,m;j1,2,,n) k1s记为CmnAmsBsn(A的列数等于B的行数)。 例1:求矩阵A4242与B36的乘积AB与BA。12 解:AB41632242 1612368 BA424002AB 361200例1说明:矩阵的乘法不满足交换律,即一般地ABBA。若ABBA,则称方阵A与B可交换。矩阵的乘法满足下列运算律: (1)结合律:(AB)CA(BC) (2)(AB)(A)BA(B)(3)分配律:A(BC)ABAC,(BC)ABACA 例2:举例说明下列命题是错误的(1)若A0,则A0; 2(2)若AA,则A0或AE; 2(3)若AXAY,且A0,则XY。 11101010 解:(1)A(2)A(3)AX11;00;00,Y01。 三、方阵的幂及方阵多项式 1、定义:设A是n阶方阵,则A1A,A2AA,,Ak1AkA klklklkl方阵的幂满足的运算律:(1)AAA;(2)(A)A 2、方阵多项式 设f(x)a0xma1xm1am1xam(a00)为m次多项式,A为n阶方阵,则 称f(A)为方阵A的多项式。f(A)a0Ama1Am1am1AamE仍为一个n阶方阵,四、习题 P61 T2(3)(4)(5)(8) T3 T4 T6 §3.2 特殊矩阵与方阵行列式 一、特殊矩阵 1、单位矩阵 10En01000010,性质:EAAEA 01nn0 2、对角矩阵 020diag(1,2,,n) 0nmm 性质:[diag(1,2,,n)]mdiag(1,m2,,n),m为正整数。 3、数量矩阵 00EE00 4、三角矩阵 00,性质:EAAEA a12a1na11a22a2na21或0annan 1性质:Aa11a22ann 5、转置矩阵 a110A000a220 an2ann如果A(aij)mn,则AT(aij)nm。 性质:(1)(A)A; (2)(AB)AB; (3)(A)A; (4)穿脱原理:(AB)BA 6、对称矩阵和反对称矩阵 TT设A(aij)nn,如果AA,则称A为对称矩阵;如果AA,则称A为反对称TTTTTTTTTT矩阵。 二、方阵行列式 性质:①ABABBA(A,B都是n阶方阵) n ②AA n ③kAknA 三、伴随矩阵 定义:n阶行列式A的各个元素的代数余子式Aij所构成的如下矩阵 A11A12A1n称为A的伴随矩阵。 A21An1A22An2 A2nAnnn1* 例1:试证:(1)AAAAAE; (2)当A0时,AA 证明:(1)因为 a11a21*故AAan1A,ijai1Aj1ai2Aj2ainAjn(i,j1,2,,n) 0,ija12a1nA11A21An1A00a22a2nA12A22An20A0AE an2annA1nA2nAnn00A同理可得A*AAE。 (2)对A*AAE两边取行列式,得AAAE * 即 AAAEA,所以当A0时,AAnnn1。 四、习题 P69 T1 T2 T6 T7 T8(2) §3.3 逆矩阵 一、逆矩阵 1、定义:对于n阶方阵A,如果有一个n阶方阵B,使 ABBAE 则称A是可逆的,并把矩阵B称为A的逆矩阵,记为BA。 2、可逆的判定定理 定理:方阵A可逆A0;当A可逆时,A11 A,其中A为A的伴随矩阵。 AE。证明:必要性.因为A可逆,即存在A,使AA111 1故AAAAE1,所以A0 充分性.由§3.3的例1可知 AAAAAE;因为A0,故有 A11AAAE AA1A。 A按照逆矩阵的定义,即有 A1注意:当A0时,称A为非奇异矩阵,否则称A为奇异矩阵。可见,可逆矩阵就是非奇异矩阵。同时,定理也提供了一种求逆矩阵的方法——伴随矩阵法(公式法)。 13、推论:若ABE(或BAE),则BA。 证明:ABABE1,故A0,从而A存在,于是 1BEB(A1A)BA1(AB)A1EA1 二、逆矩阵的运算律 方阵的逆矩阵满足下列运算律: ①若n阶方阵A可逆,则A也可逆,且(A)②若A可逆,数0,则A可逆,且A1111A; 1A1; 1③若A,B均为n阶可逆方阵,则AB也可逆,且(AB)④若A可逆,且ABAC,则BC; ⑤若A可逆,则A也可逆,且(A)T; B1A1(穿脱原理) T1(A1)T; ⑥若A可逆,则A也可逆,且(A*)1(A1)*; ⑦若A可逆,则(A*)T(AT)*; 1⑧若A可逆,则AA1* ⑨若A,B均为n阶可逆方阵,则(AB)*B*A*(穿脱原理) 证明: ①因为AA1E,由推论可知,(A1)1A ②因为A1A1AA1E,由推论可知,A11A1 1③(AB)(B1A1)A(BB1)A1AEA1AA1E,由推论有,(AB)11④因为A可逆,则AABAAC,即EBEC,故BC B1A1 ⑤AT(A1)T(A1A)TETE,由推论有,(A)⑥因为A可逆,故A1T1(A1)T 1*AA1A,且AAE,从而(A*)1A; AAAA 1又A(A)(A)A11*1*A1E,即(A1)*AA1E1A A 所以(A)*1(A1)*。 T*TT11T⑦因为(A*)T(AA1)TA(A1)T,(A)A(A)A(A) 所以(A)(A) 111⑧因为AAE1,即AA1,所以A*TT*11A A⑨由ABAB0可知,AB也可逆。又(AB)(AB)*ABE,所以(AB)*AB(AB)1ABB1A1BB1AA1B*A* ab1例 1、问Acd满足什么条件时可逆,并求A。 解:Aadbc,Acdb,当Aadbc0时,A可逆; a且 A 11db adbcca例 2、设A是三阶方阵,且A解:(3A)118A*11*,求(3A)18A 271112A18AA1A1A1 333(1)A1(1)3A11 3327A1 例 3、解矩阵方程2571913X411 解:X25171935719113411124111 三、习题 P75 T2 T3(3) T6 T7 T9 23 §3.4 分块矩阵和初等矩阵 一、分块矩阵 设AnnOA1OB1,BnnOA2O,其中Ai与Bi(i1,2)是同阶的子方块,则 B2O A2B2O 1A21A2 OA1B1①ABOA1k③AOkA1B1; ②ABOA2B2OA11O1; ④AkOA21O⑤AA; ⑥A12A2A1O1AO1 二、初等矩阵 1、定义:由n阶单位阵E经过一次初等变换得到的矩阵称为n阶初等矩阵。 2、三种初等变换对应三种初等矩阵 (1)交换第i行和第j行; 对应En(i,j)(2)第i行乘k倍; 对应En(i(k))(3)第j行乘k倍加至第i行; 对应En(i,j(k)) 24例 1、将A13化为标准形。 解:A2413131310B 1324020101则 0100113101/22110AB 0112即 E2(1,2(3))E2(2())E2(2,1(2))E2(1,2)AB 3、初等变换与初等矩阵的关系 定理1:设A是一个mn矩阵,对A施行一次初等行变换,相当于对A左乘一个相应的m阶初等矩阵;对A施行一次初等列变换,相当于对A右乘一个相应的n阶初等矩阵。 三、初等变换求逆矩阵 定理2:对任意一个mn矩阵A,总存在有限个m阶初等矩阵P1,P2,,Ps和n阶初等矩阵Ps1,Ps2,,Pk,使得P1PsAPs1PkO ErOFmn Omn定理3:对于n阶可逆矩阵A,总存在有限个n阶初等矩阵P1,,Ps,Ps1,,Pk,使得P1PsAPs1PkEnn 定理4:设A为可逆矩阵,则有限个初等矩阵P1,P2,,Pk,使得AP1P2Pk 推论:mn矩阵A与B等价存在m阶可逆矩阵P和n阶可逆矩阵Q,使 PAQB,记为AB。(等价关系具有反身性、对称性、传递性) 因此,由定理3可知,方阵A可逆AE 由定理4可知,方阵A可逆AP,2,,k为初等矩阵)1P2Pk(Pi,i 1由推论可知,AB存在可逆矩阵P,Q,使PAQB1、求逆方法的推导: 111由定理4的AP1P2Pk,得 PkP2P1AE (1)1111(1)式两端分别右乘A,得 PkP2P1EA (2) 1上述两式表明,用一样的初等行变换将A变成E的同时,会将E变成A。 2、求逆矩阵的基本方法 初等变换法:(A|E)初等行变换(E|A1)或(3、解矩阵方程AXB或XAB(A可逆) 初等变换法:(A|B)初等行变换(E|A1B)或()(四、习题 P91 T1 T2(1)(2) T3 1AE)初等列变换(1)EAAB初等列变换E)BA1§3.5 矩阵的秩 一、k阶子式的概念 2m,n}),其交叉处的k个元素定义:在mn矩阵A中,任取k行k列(1kmin{按原来的位置构成的一个k阶行列式,称为矩阵A的一个k阶子式。 11111111例:A1234,1,0等都是A的一个2阶子式。 12000000kk可知,mn矩阵A的k阶子式共有Cm个。Cn 二、矩阵的秩 定义:矩阵A的非零子式的最高阶数,称为矩阵A的秩,记为R(A)。若R(A)r,则A中至少有一个r阶子式不为0,且所有r1阶子式都为0。 三、矩阵秩的性质 m,n} ① 1R(A)min{② R(A)R(A) ③ R(A)rA的行阶梯形含r个非零行A的标准形FO④ 若A~B则R(A)R(B)(矩阵的初等变换不改变矩阵的秩) ⑤ 若P,Q可逆,则R(PAQ)R(A) ⑥ max{A,B}R(A,B)R(A)R(B); 特别地,当B为列向量b时,有R(A)R(A,b)R(A) 1⑦ R(AB)R(A)R(B) ⑧ R(AB)min{R(A),R(B)} ⑨ 若AmnBnsO,则R(A)R(B)n 例 1、设A为n阶矩阵A的伴随矩阵,证明 *TErO OR(A)nn,R(A*)1,R(A)n1 0,R(A)n1 证明: **(1)当R(A)n时,则A可逆,即A0;由AAAE知AAn10。故A*可逆,从而R(A)n (2)若R(A)n1,则AAAE0。故R(A)R(A)n,R(A)nR(A)1。又由R(A)n1知矩阵A中至少有一个n1阶子式不为零,也就是说A中至少有一个元素不为零。所以R(A)1,从而有R(A)1。 *(3)若R(A)n1,则A的任意一个n1阶子式都为零。故A0,即R(A)0。 ********21113例 2、求A42232的秩 21561211132111321113解:422320045400454 215610045200006 故R(A)3 12例 3、已知矩阵A1212a32314的秩为3,求a的值 01153554a3112a311200112a200112a2解:A 0111a20111a201152a200063a0a31120111a2 因为R(A)3,所以63a0,即a2 00112a200063a0 四、习题 P96 T2 T3(2) T7 T8 P97 总复习题:T1 T2 T3 T4 T5 第四章 线性方程组理论 教学目标与要求 1.掌握齐次和非齐次线性方程组解的判定定理和解的结构定理 2.理解向量组的线性相关与线性无关的概念,以及它们的判定方法 3.掌握向量组的秩和最大无关组的概念,会求向量组的秩 4.理解基础解系的概念,会求齐次与非齐次线性方程组的通解 教学重点 1.齐次与非齐次线性方程组解的判定定理以及通解的求法 2.向量组线性相关与线性无关的判定方法 3.向量组的最大无关组的求法和秩的求法 教学难点 1.齐次与非齐次线性方程组解的判定方法 2.向量组秩的概念及其求法 3.基础解系的概念及其求法 §4.1 线性方程组有解的条件 一、线性方程组解的判定 1、非齐次线性方程组 定理1:对于非齐次线性方程组Amnxb(1),则 ① 有唯一解R(A)R(A,b)n ② 有无穷多解R(A)R(A,b)n ③ 无解R(A)R(A,b) 2、齐次线性方程组 定理2:对于齐次线性方程组Amnx0(2),则 ① 仅有零解R(A)n ② 有非零解R(A)n 推论:当mn时,Annx0有非零解R(A)nA0 定理3:矩阵方程AXB有解R(A)R(A,B) 二、线性方程组的解法 x12x23x30例 1、求下列线性方程组的通解2x15x23x30 x8x041301090123012解:253001300130 1008023800980100810900108/3 0130018/90018/9x18x4x18 x8/382x2x4,令x41,得通解为:k(kR)x8/933 1x84x3x49 例 2、问取何值时,下列线性方程组(1)有唯一解;(2)无解;(3)有无穷多解?并在有无穷多解时求其通解。 x1x22x3 1x1(21)x23x31 x1x2(3)x32122解:A213011(1)(1)3001由克拉默法则知,当0,1,1时,方程组有唯一解。 当0时,B002101310101310021000031003100因R(A)2,R(B)3,R(A)R(B),所以方程组无解。 1121112当1时,B133110210 11230004因R(A)2,R(B)3,R(A)R(B),所以方程组无解。 11211121当1时,B1131001011010010114100200000因R(A)R(B)23,所以方程组有无穷多解。 即xxx11k112x0,令x2k,得其通解为:x2k(kR)3x30 三、习题 P106 T1 T2 T3(2)T4 T5 T6 T7 312105 2 §4.2 向量组的线性相关性 一、n维向量及其线性运算 1.定义:由n个数a1,a2,,an组成的有序数组称为n维向量。称n1矩阵 a1a2a为n维列向量;其转置aTa1,a2,,an称为n维行向量。其中ai称为a的第ian个分量(i1,2,,n)。 2.运算 ①n维向量的相等;②零向量;③负向量;④加法;⑤数乘 二、向量组的线性组合 1.向量组 定义:由若干个同维的列向量(或行向量)所组成的集合,称为一个向量组。 2.向量组与矩阵 a1ja2j(j1,2,,n)为矩阵A的列设A(aij)mn,则A1,2,,n,其中jamj12向量组;或A,其中iai1,ai2,,ain(i1,2,,m)为矩阵A的行向量组。 m3.向量组与线性方程组 一个线性方程组Amnxb可以写成:x11x22xnnb 4.向量组的线性组合 定义:设向量组A:1,2,,m,对于数k1,k2,,km,我们称k11k22kmm为向量组A的一个线性组合,k1,k2,,km称为这个线性组合的系数。 5.线性表示 给定向量组A:1,2,,m和向量b,若存在一组数1,2,,m,使得 b1122mm 则称向量b是向量组A的线性组合,也称向量b可以由向量组A线性表示。 例:任何一个n维向量aa1,a2,,an都可以由n维单位向量组: Te1(1,0,0,,0)T,e2(0,1,0,,0)T,,en(0,0,,0,1)T 线性表示。即aa1e1a2e2anen。 显然,向量b能由向量组A线性表示,也就线性方程组:x11x22xnnb有解。 6.定理1:向量b能由向量组A:1,2,,m线性表示的充要条件是R(A)R(A,b),其中A(1,2,,m)。 三、向量组的线性相关与线性无关 设齐次线性方程组Amnx0,写成向量形式:x11x22xnn0。若它有非零解,即存在一组不全为零的数k1,k2,,kn,使得k11k22knn0。因此,我们引入如下概念。 1.线性相关与线性无关 定义:设有n维向量组A:1,2,,m,如果存在一组不全为零的数k1,k2,,km使 k11k22knn0 则称向量组A线性相关;否则称它线性无关。 注意:(特殊情形) ① 只有一个向量a的向量组线性相关a0 ② 两个向量a,b的向量组线性相关ab(即两向量共线:对应分量成比例)③ 三个向量线性相关:几何意义是三个向量共面。 ④ 含有零向量的向量组一定线性相关。 定理2:向量组1,2,,m(m2)线性相关的充要条件是其中至少有一个向量可由其余m1个向量线性表示。 定理3:设向量组A:1,2,,m构成矩阵A(1,2,,m),则向量组A线性相关的充要条件是R(A)m;向量组A线性无关的充要条件是R(A)m。 推论1:当向量的个数等于向量的维数时,向量组A线性相关的充要条件是A0;向量组A线性无关的充要条件是A0。 推论2:m(mn)个n维向量组成的向量组一定线性相关。推论3:任一个n维向量组中线性无关的向量最多有n个。 定理4: (1)设向量组A:1,2,,m线性无关,而向量组B:1,2,,m,b线性相关,则向量b必能由向量组A线性表示,且表示法是唯一的。 (2)若向量组1,2,,r线性相关,则向量组1,2,,r,r1,,n(nr)必线性相关;反之,若向量组1,2,,r,r1,,n(nr)线性无关,则向量组1,2,,r必线性无关。(部分相关,整体相关;整体无关,部分无关。) (3)若m个n维向量1,2,,m线性相关,同时去掉其第i个分量(1in)得到的m个n1维向量也线性相关;反之,若m个n1维向量1,2,,m线性无关,同时增加其第i个分量(1in)得到的m个n维向量也线性无关。 四、习题 P116 T1(3)(4)T2 T3 T4(1)(2)T5 T6 T7 T8 T9(1)(3) §4.3 向量组的秩 一、向量组的等价 定义1:设有向量组A:1,2,,m;向量组B:1,2,,s,若向量组A中的每一个向量都能由向量组B线性表示,则称向量组A能由向量组B线性表示。如果向量组A和向量组B能相互线性表示,则称这两个向量组等价。 命题1:若A,B为有限个列向量组成的向量组,则向量组B能由向量组A线性表示的充要条件是矩阵方程BAX有解。 命题2:若矩阵A经过初等行(列)变换变成B,则矩阵A的列(行)向量组与矩阵B的列(行)向量组等价。 定理1:设向量组A:1,2,,m和向量组B:1,2,,s均为列向量组成的向量组,则向量组B能由向量组A线性表示的充要条件为R(A)R(A,B) 推论:向量组A:1,2,,m和向量组B:1,2,,s等价的充要条件是 R(A)R(B)R(A,B) 其中A和B是向量组A和向量组B所构成的矩阵。 讲教材P118例1 二、向量组的秩 1.最大无关组 定义2设向量组A0:1,2,,r是向量组A:1,2,,m(mr)的一个部分组,若(1)向量组A0:1,2,,r线性无关; (2)A中的任意向量均可由向量组A0:1,2,,r线性表示; 则称A0:1,2,,r为A的一个最大线性无关向量组(简称最大无关组)。 显然,最大无关组一般不唯一;任意向量组都与它的最大无关组等价。 2.最大无关组的求法 定理:矩阵的初等行变换不改变(部分或全部)列向量之间的线性关系; 矩阵的初等列变换不改变(部分或全部)行向量之间的线性关系。 注意:上述定理提供了求向量组最大无关组的方法 定理2:设向量组B:1,2,,r可由向量组A:1,2,,s线性表示,(1)若向量组B线性无关,则rs;(2)若rs,则向量组B线性相关。 推论1:两个等价的线性无关的向量组必含有相同个数的向量。推论2:两个等价的向量组的最大无关组含有相同个数的向量。推论3:一个向量组的任意两个最大无关组所含向量个数相等。 3.向量组的秩 定义3:向量组的最大无关组所含向量的个数,称为该向量组的秩。 定理2':若向量组B能由向量组A线性表示,则向量组B的秩不大于向量组A的秩。 三、矩阵的秩与向量组的秩的关系 定理3:对矩阵A(aij)mn,则 R(A)A的行秩A的列秩。即矩阵的秩等于它的行向量组的秩也等于它的列向量组的秩。 四、矩阵的秩的性质 性质1:R(AB)R(A)R(B) 性质2:R(AB)min{R(A),R(B)} 性质3:若P,Q可逆,则R(PAQ)R(PA)R(AQ)R(A) 五、习题 P124 T1 T2 T3 T9 §4.4 线性方程组解的结构 一、齐次线性方程组解的结构 1.解的性质 对于齐次线性方程组 Amnx0 (1)性质1:若1,2都是Ax0的解,则12也是Ax0的解。性质2:若是Ax0的解,则k也是Ax0的解。 2.解的结构 定义1:设1,2,,k是Ax0的非零解,且满足 (1)1,2,,k线性无关; (2)Ax0的任一个解都可由1,2,,k线性表示,即c11c22ckk 则称1,2,,k是齐次线性方程组Ax0的基础解系;且Ax0的通解可表示为如下形式:c11c22ckk(c1,c2,,ck为任意常数)。 定理1:若n元齐次线性方程组Ax0的系数矩阵A的秩R(A)rn,则Ax0的基础解系恰含有nr个线性无关的解向量。 讲教材P128 例1和例2 二、非齐次线性方程组解的结构 1.解的性质 对于非齐次线性方程组 Amnxb (2)性质1:若1,2都是Axb的解,则12是Ax0的解。 性质2:若是Ax0的解,是Axb的解,则是Axb的解。 2.解的结构 *定理2:设是非齐次线性方程组Axb的一个解,1,2,,nr是对应的导出组Ax0的基础解系,则Axb的通解为 *k11k22knrnr 其中k1,k2,,knr为任意常数。 讲教材P132 例3和例4 三、习题 P134 T1 T2(1)T3 T4 T5 T6 T7 T8 P141 总复习题:T1 T2 T4 T5 T6至T13 第五章 特征值和特征向量 矩阵的对角化 教学目标与要求 1.理解内积和正交向量组的概念,掌握施密特正交化方法和正交矩阵的性质 2.理解特征值与特征向量的定义,掌握它们的性质及其求法 3.理解相似矩阵的定义,掌握相似矩阵的性质 4.掌握矩阵可对角化的条件,熟悉实对称矩阵的对角化方法 教学重点 1.施密特正交化方法的运用 2.特征值与特征向量的求法 3.实对称矩阵的对角化方法 教学难点 1.施密特正交化方法 2.特征值与特征向量的性质及其求法 3.实对称矩阵的对角化方法 §5.1 预备知识 一、向量的内积 定义1:设有n维向量xx1,x2,,xn,yy1,y2,,yn,令 TTx,yx1y1x2y2xnyn,称x,y为向量x与y的内积。 内积的性质: (1)x,yy,x (2)x,yx,y (3)xy,zx,zy,z (4)x,x0,当且仅当x0时等号成立 定义2:令xx,x22x12x2xn,称为n维向量x的长度(或范数)。当x1时,称x为单位向量。 向量的长度具有以下性质: (1)非负性:x0 (2)齐次性: 定义3:当x0,y0时,称arccosxx (3)三角不等式:xyxy (4)柯西不等式:x,yxy x,yxy为n维向量x与y的夹角。 定义4:当x,y0时,称向量x与y正交。 定义5:若一个向量组中任意两个向量都正交,则称此向量组为正交向量组。若正交向量组中的每一个向量都是单位向量,则称此向量组为规范正交向量组或标准正交向量组。 定理1:若n维向量1,2,,r是一组两两正交的非零向量,则1,2,,r线性无关。 二、施密特正交化方法 施密特正交化方法是将一组线性无关的向量1,2,,r,化为一组与之等价的正交向量组1,2,,r的方法。令 2,1;; 1,11,,r,r1。rrr11r221,12,2r1,r1r111; 22 讲教材P147 例2和例3 三、正交矩阵 定义6:如果方阵A满足AAAAE(即Acos例如:En,sinAT),则称A为正交矩阵。 01/21/2sin,2/61/61/6都是正交阵。cos1/31/31/3TT1 定理2:A为正交矩阵A的行(列)向量组为规范正交向量组。即 1,ijATAEiTj(i,j1,2,,n)(其中A(1,2,,n)) 0,ij 定理3:设A,B都是n阶正交方阵,则 (1)A1;(2)A,A,AB也是正交方阵。 定义7:若P为正交矩阵,则线性变换yPx称为正交变换。 四、习题 P149 T1(2)T2(2)T3 T4 T5 §5.2 特征值和特征向量 T 1一、特征值与特征向量的概念 定义1:设A是n阶方阵,如果存在数和非零列向量x,使得Axx,称为方阵A的特征值,非零列向量x称为A的属于特征值的特征向量。 特征方程:Axx(AE)x0 或者(EA)x0 (AE)x0有非零解AE0特征矩阵:(AE)或者(EA) EA0 a11特征多项式:AEa12an2a1na2n() a21an1a22annnn1aaan1an0[a0(1)n] 二、求n阶方阵A的特征值与特征向量的步骤 (1)求出特征方程()AE0的全部根1,2,...,n,即是A的特征值;(2)对于每个特征值i求解线性方程组AiEx0,得出的基础解系就是A的属于特征值i的特征向量;基础解系的线性组合就是A的属于特征值i的全部特征向量。 讲教材P152 例3和例4 三、特征值与特征向量的性质 性质1:设A是n阶方阵,则A与A有相同的特征值。性质2:设是方阵A的特征值,k,mN,则(1)是方阵A的特征值; (2)f()a0a1am是f(A)a0Ea1AamA的特征值。 性质3:设n阶方阵A(aij)nn的n个特征值为1,2,...,n,则(1) mmkkTaii1i1nnii,其中 ai1niitr(A)称为A的迹; (2)iA i1n 证明: 由特征值的定义可得 a11 a12a1na2n ()AEa21an1a22an2ann (a11)(a22)(ann) (1)nn(1)n1(a11a22ann)n1 由题设可知 ()AE(1)(2)(n) (1)nn(1)n1(12n)n1(12n)比较多项式同次幂的系数可得 a11a22ann12n,A(0)12n 推论:A0 0是A的特征值;A可逆A0A不含零特征值。 讲教材P154 例5和例6 性质4:1,2,,m是方阵A的互异特征值,其对应的特征向量依次为 p1,p2,,pm,则向量组p1,p2,,pm线性无关。 四、习题 P157 T1 T2 T3 T4 §5.3 相似矩阵 一、相似矩阵的概念 定义1:设A,B都是n阶方阵,若存在可逆矩阵P,使PAPB,则称矩阵A与B相似,记为A~B,可逆矩阵P称为相似变换矩阵。 相似矩阵的基本性质: 1、(1)反身性:对任意方阵A,都有A~A (2)对称性:若A~B,则B~A (3)传递性:若A~B,B~C,则A~C 2、定理1:若A~B,则 ① A与B有相同的特征多项式和特征值; ② AB; ③ R(A)R(B); mm④ A与B也相似(m为正整数); 1⑤ tr(A)tr(B) 二、矩阵可对角化的条件 定义:n阶方阵A可以相似于一个对角矩阵,则称A可对角化。 定理2:n阶方阵A可对角化A有n个线性无关的特征向量。 推论:n阶方阵A有n个互异的特征值A可对角化。 定理3:n阶方阵A可对角化A的每个k重特征值对应有k个线性无关的特征向量(或R(AE)nk)。即A的几何重数nR(AE)等于代数重数k。 讲教材P160 例1和例2 三、小结 n阶方阵A对角化的步骤: (1)解特征方程AE0,求出A的全部特征值1,2,...,s,其中i是ni重特征值(i1,2,,s),sni1in。 (2)对每个i,解齐次线性方程组AiEx0,得基础解系i1,i2,...,ini;(3)令P(11,12,,1n1,21,22,,2n2,,s1,s2,,sns),则PAP,其中diag(1,,1,2,,2,,s,,s),这里i的个数为ni个(i1,2,,s)。 四、习题 P162 T1 T2 T3 T4 T5 T6 §5.4 实对称矩阵的相似矩阵 1一、实对称矩阵的特征值性质 定理1:实对称矩阵的特征值都是实数。 定理2:实对称矩阵A的属于不同特征值的特征向量相互正交。 定理3:设是n阶实对称矩阵A的r重特征值,则R(AE)nr,即对应特征值恰有r个线性无关的特征向量。 二、实对称矩阵的相似理论 定理4:任意实对称矩阵A都与对角矩阵相似。即实对称阵一定可以对角化。 1T定理5:设A是n阶实对称矩阵,则存在正交矩阵P,使PAPPAP。其中diag(1,2,,n),且1,2,...,n是A的n个特征值。 三、实对称矩阵对角化方法 n阶实对称矩阵A对角化的步骤: (1)解特征方程AE0,求出A的全部特征值1,2,...,s,其中i是ni重特征值(i1,2,,s),sni1in。 (2)对每个i,解齐次线性方程组AiEx0,得基础解系i1,i2,...,ini;(3)利用施密特正交化方法将i1,i2,...,ini正交化,得正交向量组i1,i2,...,ini,再单位化得规范正交向量组i1,i2,...,ini(i1,2,,s); (4)令P(11,12,,1n1,21,22,,2n2,,s1,s2,,sns),则P为正交矩阵,且P1APPTAP,其中diag(1,,1,2,,2,,s,,s),这里i的个数为。ni个(i1,2,,s) 讲教材P164 例1和例2 四、习题 P167 T1 T2 T4 P167 总复习题:T1 T2 T3 T4 T5 T6; T8 T9 T10 T11 T12 T13 T14 T15 T16 第六章 特征值和特征向量 矩阵的对角化 教学目标与要求 1.理解二次型及其秩的相关概念,了解矩阵的合同关系 2.掌握二次型的标准形,以及用配方法、正交变换法和初等变换法化二次型为标准型 3.理解惯性定理和二次型的规范形,掌握二次型正定的判别方法 教学重点 1.用正交变换法化二次型为标准型 2.二次型正定的判别方法 教学难点 1.用正交变换法化二次型为标准型 2.二次型正定的判别方法 §6.1 二次型及其矩阵表示 一、二次型及其矩阵表示 定义1:含有n个变量的二次齐次函数: 22f(x1,x2,...,xn)a11x12a22x2annxn 2a12x1x22a13x1x32an1,nxn1xn称为二次型。当aij全为实数时,f称为实二次型。 为了便于用矩阵讨论二次型,令aijaji,则二次型为: f(x1,x2,...,xn)a11x12a12x1x2a1nx1xn2 a21x2x1a22x2a2nx2xn.................................................2 an1xnx1an2xnx2annxn a11a21记 Aan1a12a22an2i,j1anijxixj a1nx1a2nx2x,,xannnT则二次型f(x1,x2,,xn)xAx,其中A为对称矩阵。 由此可见,对称矩阵A与二次型f是一一对应关系,故称对称矩阵A为二次型f的矩阵,也称二次型f为对称矩阵A的二次型,R(A)也称为二次型f的秩。 讲教材P173 例1和例2 二、线性变换 x1c11y1c12y2c1nynxcycycy22112222nn 定义2:称为由变量x1,x2,,xn到变量y1,y2,,yn.................................................xncn1y1cn2y2cnnyn的一个线性变量替换,简称线性变换。 c11c21其中,矩阵Ccn1c1nc22c2n称为线性变换的矩阵。cn2cnnc12x1y1x2y2记x,y,则线性变换可用矩阵形式表示:xCy。 xynn若C0,则称线性变换xCy为非退化的(或满秩变换);否则,称为退化的(或降秩变换)。若C是正交矩阵,则称线性变换xCy为正交变换。因此,我们有 f(x)xTAx(Cy)TA(Cy)yTCTACyyTBy,其中BCTAC,而且 BT(CTAC)TCTATCCTACB 三、矩阵的合同 1.定义3:设A,B为两个n阶方阵,如果存在n阶可逆矩阵C,使得CACB,则 TB。称矩阵A与B合同,记为:A~B(合同)定理:若A~,则AB(等价),且R(A)R(B)。 2.合同的性质 A ① 反身性:对任意方阵A,都有A~B,则B~A ② 对称性:若A~C B,B~C,则A~③ 传递性:若A~3.定理:任何一个实对称矩阵A都合同于一个对角阵(是以A的n个特征根为对角元的对角阵),即存在可逆矩阵C,使得CAC。 四、习题 P175 T1 T3 T4 §6.2 二次型的标准形 T一、二次型的标准形 222定义:形如d1x1的二次型称为二次型的标准形。d2x2dnxn 二、化二次型为标准形 (1)配方法 对任意一个二次型fxTAx,都可用配方法找到满秩变换xCy,将f化为标准形。步骤:若f中含变量项xi的平方项,则先将所有含xi的项合并在一起配成完全平方,依次类推直到都配成完全平方项;若f中不含任何平方项,则令x1y1y2,x2y1y2,xkyk,使f中出现平方项,再按照前面的思路进行配方。 (2)正交变换法 定理:任给二次型f(x)xTAx,总存在正交矩阵Q,使QTAQQ1AQ,其中diag(1,2,,n),1,2,,n是A的全部特征值。 22即存在正交变换xQy使f化为标准形:(其中1,2,,n1x122x2nxn是对称矩阵A的全部特征根) 讲书上P176 例1 (3)初等变换法 由于任意对称阵A都存在可逆矩阵C,使CAC为对角阵;由于C是可逆阵,故可表 TTTT示一系列初等矩阵的乘积。设CP1P2PS,则CPsP2P1,因此 TCTACPsTP2TP1AP1P2Ps ① T CP1P2PSEP1P2PS ② ①式表示对实对称矩阵A施行初等列变换的同时也施行相应的行变换,将A化为对角阵;②表示单位阵E在相同的初等列变换下就化为C。即(三、习题 P181 T1 T3 T4 §6.3 惯性定理和二次型的正定性 A)合同变换()EC 一、惯性定理和规范形 定理1:设实二次型fxTAx的秩为r,有两个实满秩线性变换xCy及xPz,222使得 fk1y1kpy2,2,,r) (1)pkp1yp1kryr(ki0,i12222及 f1z1qzqq1zq,2,,r)1rzr(i0,i1则pq;且称p为二次型f的正惯性指数,rp为二次型f的负惯性指数。 对二次型f的标准形(1)式再作满秩线性变换 (y1,,yr,yr1,,yn)Tdiag(11,,1,,1)(t1,,tr,tr1,,tn)T k1kr2222则有ft1tptp1tr,称之为二次型f的规范形。 惯性定理的等价表述:任意一个秩为r的实二次型f都可以经过满秩线性变换化为规范形,且其规范形是唯一的。即规范形中正项的个数p与负项的个数rp都是唯一确定的。 定理2:实对称阵A与B合同A与B的正负惯性指数相同 A与B的规范形相同R(A)R(B),且A与B的正惯性指数相同 二、二次型的正定性 定义1:设实二次型f(x)f(x1,x2,,xn)xTAx,若对任意x0,都有f(x)0,则称f为正定二次型,并称其对称矩阵A为正定矩阵。三、二次型正定的判别方法 定理3:设A是n阶实对称矩阵,则 fxTAx正定(或A正定)A的n个特征值全为正; f的标准形的n个系数全为正f的正惯性指数pn; 存在可逆矩阵P,使APTPA与单位矩阵合同; A的各阶顺序主子式全为正,即 a11a1na11a120 a110,0,,a21a22an1ann讲教材P184 例3 四、习题 P185 T1(1)(3) T2(3) T3 T4 T5 T6 P186 总复习题: T4 T5 T6 T7 ; T9 T12 T13 线性代数教案第一章 第一章 行列式(12学时) 教学时数:12学时 教学目的与要求:理解并掌握行列式的概念和性质,行列式按行(列)展开定理,行列式的计算,克莱姆法则解方程组。 教学重点:行列式的性质,行列式按行(列)展开,克莱姆法则解方程组。教学难点:行列式按行按列展开。本章主要阅读文献资料: 1.吴赣昌主编,《线性代数》(第4版),中国人民大学出版社,2008年2月。2.戴斌祥主编,《线性代数》,北京邮电大学出版社,2005年10月。3.陈维新主编,《线性代数》(第二版),科学出版社,2010年8月。 4.赵树嫄主编,《线性代数学习与考试指导》,中国人民大学出版社,2008年5月。 教学内容: 第一节 二阶与三阶行列式 一.二阶行列式 引入新课: 我们从二元方程组的解的公式,引出二阶行列式的概念。 在线性代数中,将含两个未知量两个方程式的线性方程组的一般形式写为 (1) 用加减消元法容易求出未知量x1,x2的值,当 时,有 (2)这就是二元方程组的解的公式。但这个公式不好记,为了便于记这个公式,于是引进二阶行列式的概念。 (一)定义:我们称记号 为二阶行列式,它表示两项的代数和: 即定义 (3) 二阶行列式所表示的两项的代数和,可用下面的对角线法则记忆:从左上角到右下角两个元素相乘取正号,从右上角到左下角两个元素相乘取负号,即 - + 由于公式(3)的行列式中的元素就是二元方程组中未知量的系数,所以又称它为二元方程组的系数行列式,并用字母D表示,即有 如果将D中第一列的元素a11,a21 换成常数项b1,b2,则可得到另一个行列式,用字母D1表示,于是有 按二阶行列式的定义,它等于两项的代数和:,这就是公式(2)中x1 的表达式的分子。同理将D中第二列的元素a a b2,12,22 换成常数项b1,可得到另一个行列式,用字母D2表示,于是有 按二阶行列式的定义,它等于两项的代数和:a11b2-b1a21,这就是公式(2)中x2的表达式的分子。 于是二元方程组的解的公式又可写为 其中D≠0 例1 计算51=5×2-(-1)×3=13 32例2 设D231 问:(1)当λ为何值时D=0(2)当λ为何值时D≠0 解:D231=23 (1)当λ=0或3时,D=0(1)当λ≠0且λ≠3时,D≠0 二.三阶行列式 含有三个未知量三个方程式的线性方程组的一般形式为 (1) 还是用加减消元法,即可求得方程组(1)的解的公式,当 时,有 (2) 这就是三元方程组的解的公式。这个公式更不好记,为了便于记它,于是引进三阶行列式的概念。 (二)定义: 我们称记号 为三阶行列式。三阶行列式所表示的6项的代数和,也用对角线法则来记忆:从左上角到右下角三个元素相乘取正号,从右上角到左下角三个元素取负号,即 (3) 由于公式(3)的行列式中的元素是三元方程组中未知量的系数,所以称它为三元方程组的系数行列式,也用字母D来表示,即有 同理将D中第一列、第二列、第三列的元素分别换成常数项得到另外三个三阶行列式,分别记为 于是有 就可以 按照三阶行列式的定义,它们都表示6项的代数和;并且分别是公式(2)中x1,x2,x3 的表达式的分子,而系数行列式D是它们的分母。 123例3 405 106解:原式=-58 例4 实数a,b满足什么条件时 ab0ba00 101ab0解:ba0a2b2 a,b为实数,若要a2b20,则a,b需同时等于零。 a10例5 1a0>0的充分必要条件是什么? 411a10a10解:1a0=a21,即a>1时,1a0>0,411411a10所以1a0>0的充分必要条件a>1 411作业:课本35页,1,2,3,4,5 第二章MATLAB软件 2.1 MATLAB语言及其特点 人机界面友好 强大而简易的作图功能 功能丰富,可扩展性强 超强的数值运算功能 实用的程序接口 2.2 MATLAB语言的工作环境 2.2.1 MATLAB启动 2.2.2命令窗(Command Window)中的菜单与功能 MATLAB命令窗口是用户和MATLAB交互的主要窗口。命令窗口菜单条下有6个下拉子菜单: File、Edit、Window、debug、desktop和Help。File子菜单: New及其子菜单:允许用户打开一个新的文件。 Open....:从指定路径和文件名打开一个已经存在的文件。Close Command Window:关闭命令窗口。 Import Data....:从指定路径和相应文件中获取数据,在工作空间中生成一变量。Save Workspace As....:保存在指定路径下的相应文件中。Set Path....:设置MATLAB的搜索路径。Preferences...:对系统性能参数进行设置。 Page Setup:设置页面布局、页眉、页脚中所用的字体等 Edit子菜单: Undo、Redo、Cut、Copy和Paste等:分别用于撤销、重复上一次操作、剪切、复制和粘贴。 Paste Special:将存放在缓冲区内的内容按照一定的格式存放在剪贴板变量中。Select All:用于选定当前窗口中所有内容。Delete:删除预先选定的内容。 Find:当前目录中的M-文件中查找相应内容。 Clear Command Window:清除命令窗口中显示内容,但不会删除工作空间中内容。(命令:clc) Clear Command History:清除历史记录。Debug子菜单: Open M-File „:打开一个要调试的M-文件 Step:逐个语句运行 Step In:进入下一个语句 Step Out:退出调试 Clear breakpoints in all file:清除所有文件断点 Stop if errors/warning„:遇到错误或警告断点停止运行 2.2.3 工作空间 允许用户查看当前工作空间的内容,它的作用与命令“whos”相同(“whos”的作用是:键入后可在命令窗口中查看当前工作空间的内容),不同的是用图形化的表示方法来显示。而且,通过它可以对工作空间中的变量进行删除、保存、修改等操作。用鼠标双击工作空间中所选变量,则进入数组编辑器(Array Editor),用户可对变量的维数、内容等进行修改。 2.2.4命令历史窗口(Command History)主要显示曾经在Command Window窗口执行过的命令。2.2.5当前路径窗口(Current Directory)主要显示当前工作路径,包括M文件的打开路径等。当前路径窗口允许用户对MATLAB的路径进行查看和修改。通常启动MATLAB的默认当前路径是MATLABwork,不改变当前目录,用户自己的工作空间和文件都将保存到该目录。2.2.6 MATLAB搜索路径 MATLAB是采用路径搜索的方法来查找文件系统中的M文件。如果在命令窗口中输入如下命令: >>test MATLAB对这一命令的搜索顺序为: (1)检查“test”是否是工作空间中的变量;(2)检查“test”是否为内部函数;(3)检查当前目录中是否有test.m、test.mex或test.dll文件; (4)检查MATLAB搜索路径上是否存在test.m、test.mex或test.dll文件;(5)如不满足上述任何一条件,则返回出错信息。如果在搜索路径中存在两个或多个同名函数,则只能发现搜索路径中的第一个函数,而其它同名函数不被执行。 为了运行已有程序,用户可事先设置好MATLAB的搜索路径,然后再运行程序。方法一:在当前路径窗口(Current Directory)将文件所在的路径设置为当前路径; 方法二:单击File菜单中的Set Path„项,弹出路径设置对话框如图1.11所示。2.2.7 MATLAB帮助系统 帮助浏览器 help命令 (1)在线帮助总览 >> help % 在命令窗口直接输入help,显示在线帮助总览(2)显示基本函数信息 >> help elfun(3)显示具体函数的帮助信息 >> help round %显示具体函数的详细信息,本例为round函数 注意:MATLAB对字母的大小写是敏感的,变量A与变量a表示两个不同的变量。MATLAB所有的命令和函数都必须用小写。lookfor命令 模糊查询 在线帮助页 2.3 MATLAB基本运算 2.3.1矩阵运算 第3章 MATLAB绘图 数据可视化是MATLAB的强大功能之一,数据可视化不仅仅是二维的,还可以在三维空间展示数据,而数据或者图形的可视化也是进行数据处理或者图形图像处理的第一步骤。 MATLAB绘图都是在MATLAB图形窗体中进行的,而所有图形数据可视化的工作也都以图形窗体为主。 MATLAB的图形窗体主要包括如下几个部分: * 菜单栏:图形窗一般包括一个菜单栏,利用这个菜单栏可完成对窗体中各种对象的基本操作。* 工具条:图形窗体的工具条用来完成对图形对象的一般性操作。* 绘图区域:图形窗体的绘图区域是面积最大的一部分,在这个区域中可以绘制各种曲线,显示图形图像文件,以及完成对图形图像或者曲线的编辑。 数据可视化的步骤: * 准备绘图数据。* 创建图形窗体。 * 使用绘图函数绘制图形或者曲线。* 设置图形属性第二篇:matlab线性代数例题
第三篇:线性代数教案
第四篇:线性代数教案第一章
第五篇:MATLAB教案