第七章线性代数中的数值计算问题.pptVIP

第1页，共33页，星期日，2025年，2月5日要解决的问题矩阵求逆行列式的求法及其应用线性方程组的求解第2页，共33页，星期日，2025年，2月5日一、对角阵与三角阵

1．对角阵

只有对角线上有非0元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为1的对角矩阵称为单位矩阵。7.1特殊矩阵第3页，共33页，星期日，2025年，2月5日提取矩阵的对角线元素

设A为m×n矩阵，diag(A)函数用于提取矩阵A主对角线元素，产生一个具有min(m,n)个元素的列向量。

构造对角矩阵

设V为具有m个元素的向量，diag(V)将产生一个m×m对角矩阵，其主对角线元素即为向量V的元素。

第4页，共33页，星期日，2025年，2月5日例7-1先建立5×5矩阵A，然后将A的第一行元素乘以1，第二行乘以2，…，第五行乘以5。

A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;...

11,18,25,2,19];

D=diag(1:5);

D*A%用D左乘A，对A的每行乘以一个指定常数

第5页，共33页，星期日，2025年，2月5日2．三角阵

三角阵又进一步分为上三角阵和下三角阵，所谓上三角阵，即矩阵的对角线以下的元素全为0的一种矩阵，而下三角阵则是对角线以上的元素全为0的一种矩阵。第6页，共33页，星期日，2025年，2月5日(1)上三角矩阵

求矩阵A的上三角阵的MATLAB函数是triu(A)。

triu(A)函数有另一种形式triu(A,k)，其功能是求矩阵A的第k条对角线以上的元素。例如，提取矩阵A的第2条对角线以上的元素，形成新的矩阵B。(2)下三角矩阵

在MATLAB中，提取矩阵A的下三角矩阵的函数是tril(A)和tril(A,k)，其用法与提取上三角矩阵的函数triu(A)和triu(A,k)完全相同。第7页，共33页，星期日，2025年，2月5日一、矩阵的转置与旋转

1．矩阵的转置

转置运算符是单撇号(‘)。

2、矩阵的左右翻转

对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第二列调换，…，依次类推。MATLAB对矩阵A实施左右翻转的函数是fliplr(A)。7.2矩阵分析第8页，共33页，星期日，2025年，2月5日3、矩阵的上下翻转

MATLAB对矩阵A实施上下翻转的函数是flipud(A)。第9页，共33页，星期日，2025年，2月5日1．矩阵的逆

对于一个方阵A，如果存在一个与其同阶的方阵B，使得：

A·B=B·A=I(I为单位矩阵)

则称B为A的逆矩阵，当然，A也是B的逆矩阵。

求一个矩阵的逆是一件非常烦琐的工作，容易出错，但在MATLAB中，求一个矩阵的逆非常容易。求方阵A的逆矩阵可调用函数inv(A)。

二、矩阵的逆与伪逆第10页，共33页，星期日，2025年，2月5日例7-2用求逆矩阵的方法解线性方程组。设则原线性方程组可简写为：其解为：第11页，共33页，星期日，2025年，2月5日第12页，共33页，星期日，2025年，2月5日逆矩阵的求解方法1：A^-1方法2：inv(A)方法3：A\eye(4)方法4：U=rref([A,eye(4)]);U(:,5:8)第13页，共33页，星期日，2025年，2月5日rref函数：求矩阵的行最简阶梯形矩阵matlab把“最简行阶梯形式（reducedrowechelomform）”的计算过程集成为一个子程序rref，它的输入变元可以是线性方程组的系数矩阵，也可以是其增广矩阵，输出变元是它的最简行阶梯形式。第14页，共33页，星期日，2025年，2月5日用matlab语言表达行变换：1)将矩阵的第i,j两行进行交换的语句为：a([i,j],:)=a([j,i],:)2)将矩阵的第j行乘以常数k的语句为：a(j,:)=k*a(j,:)3)将矩阵的第i行乘以常数k加到第j行的语句为：a(j,:)=a(j,:)+k*a(i,:)第15页，共33页，星期日，2025年，2月5日当方阵为奇异矩阵（行列式等于零）时，要求其逆矩阵，matlab将给出警告信息：MatrixisclosetosingularorbadlyscaledResultsmaybeinaccurate.第16页，共33页，星期日，2025年，2