用MATLAB线性代数的基本运算.docxVIP

用MATLAB线性代数的基本运算 第五章线性代数的基本运算 本章学习的主要目的： 1复习线性代数中有关队列式、矩阵、矩阵初等变换、向量的线性有关性、线性方程组的求解、相像矩阵及二次型的有关知识. 2学会用MatLab软件进行队列式的计算、矩阵的基本运算、矩阵初等变换、向量的线性有关性的鉴别、线性方程组的求解、二次型化标准形的运算. 5.1队列式 5.1.1n阶队列式定义 由n2个元素aij(i,j1,2,,n)组成的记号 a11 a12 a1n D=a21 a22 a2n an1 an2 ann 111 称为n阶队列式.其值是所有取自不同行不同 列的n个元素的乘积 a1p1a2p2anpn的代数和各项 , 的符号由n级排列p1p2 pn决定,即 D= (1)(p1p2pn)a1p1a2p2 anpn, p1p2 pn 其中 表示对所有n级排列求和,(p1,p2,,pn) p1p2 pn 是排列p1p2pn的逆序数. 队列式的性质 (1) 队列式与它的转置队列式相等. 交换队列式的两行(列),队列式变号. 若队列式有两行(列)完全相同,则此队列式为零. 队列式的某一行(列)中所有的元素都乘以同一数k,等于用数k乘此队列式. 若队列式有两行(列)元素成比率,则此队列式为零. 若队列式的某一列(行)的元素是两数的和,则此队列式等于对应两个队列式之和.即 112 a11 a12 a1i a1i a1n a11 a12 a1i a1n a11 a12 a1i a1n a21 a22 a2i a2i a2n a21 a22 a2i a2n a21 a22 a2i a2n an1 an2 ani ani ann an1 an2 ani ann an1 an2 ani ann 若队列式的某一行(列)的各元素乘以同一数加到另一行(列)对应的元素上去,队列式不变. 队列式等于它的任一行(列)的各元素与其 对应的代数余子式乘积之和,即 n D,i k ,n), D aijAkj (i 1,2, 0,ik j1 或 n D,j k D aijAik ( j 1,2, ,) 0,j k n i1 (9) 设A,B是n阶方阵,则A AT ,kA knA,ABAB, (10)若A是n阶可逆矩阵,则A0,A 1 1 A (11) 设1,2,,n是n阶方阵A 的特点值,则 A n i, i1 (12) 设A*是n 阶方阵A 的陪伴矩阵,则 A* An1 n2 (13)几种特殊队列式的计算: 113 a11 0 0 0 a22 0 a11a22ann , 0 0 ann a11 0 0 a21 a22 0 ann, a11a22 an1 an2 ann  a11 a12 a1n 0 a22 a2n a11a22ann 0 0 ann a11 a12 a1n n(n 1) a21 a22 0 (1) 2 a1na2n1an1 an1 0 0 5.1.3MatLab计算队列式的命令 det(var)%计算方阵 var的队列式 1 3 2 2 例1计算队列式 3 4 0 9的值 262 383 在MatLab命令窗口输入: A=[1,-3,2,2;-3,4,0,9;2,-2,6,2;3,-3,8,3] det(A) 履行结果: A=1 -3 2 2 -3 4 0 9 2 -2 6 2 3 -3 8 3 ans=-50 a 1 0 0 例2 计算队列式 1 b 1 0 的值,其中a,b,c,d是 0 1 c 1 0 0 1 d 参数. 114 在MatLab命令窗口输入: symsabcd A=[a,1,0,0;-1,b,1,0;0,-1,c,1;0,0,-1,d] det(A) 履行结果: A=[ a, 1, 0, 0] [-1, b, 1, 0] [ 0,-1, c, 1] [ 0, 0,-1, d] ans=a*b*c*d+a*b+a*d+c*d+1 1 1 1 1 例3求方程 1 2 2 x 0 的根. 1 4 4 x2 1 8 8 x3 先求队列式的值 在MatLab命令窗口输入: symsx A=[1,1,1,1;1,-2,2,x;1,4,4,x*x;1,-8,8,x^3] y=det(A) 履行结果: A= [ 1, 1, 1, 1] [ 1, -2, 2, x] [ 1, 4, 4,x^2] [ 1, -8, 8,x^3] y=-12*x^3+48*x+12*x^2-48 115 求3次方程的根. 首先经过函数的图形确定根的大概范围, 在MatLab命令窗口输入: gridon ezplot(y) -12x3+48x+12x2-48 3000 2000 1000 0 -1000 -2000 -6-4-2 024 6 图1 察看图1,可知3 个根大概在