精品PPT课件MATLAB从入门到精通教程第章数值计算.pptVIP

第5章 数值计算 5.1 特殊矩阵 5.2 矩阵分析 5.3 矩阵分解与线性方程组求解 5.4 数据处理与多项式计算 5.1 特殊矩阵5.1.1 对角阵与三角阵 1. 矩阵的对角元素 (1)提取矩阵的对角线元素 设A为m×n矩阵，diag(A)函数用于提取矩阵A主对角线元素产生一个具有min(m,n)个元素的列向量。 diag(A)函数还有更进一步的形式diag(A,k)，其功能是提取第k条对角线的元素。(与主对角线平行，往上第1条，往下第-1条，以此类推) (2)构造对角矩阵 设V为具有m个元素的向量，diag(V)将产生一个m×m对角矩阵，其主对角线元素即为向量V的元素。 diag(V)函数也有更进一步的形式diag(V,k)，其功能是产生一个n×n(n=m+|k|)对角阵，其第k条对角线的元素即为向量V的元素。 例5.1 先建立5×5矩阵A，然后将A的第1行元素乘以1，第2行乘以2，…，第5行乘以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,2,3,4,5]); D*A 2. 矩阵的三角阵 (1)下三角矩阵 求矩阵A的下三角阵的MATLAB函数是tril(A)。 tril(A)函数也有更进一步的一种形式tril(A,k)，其功能是求矩阵A的第k条对角线以下的元素。 (2)上三角矩阵 在MATLAB中，提取矩阵A的上三角矩阵的函数是triu(A)和triu(A,k)，其用法与提取下三角矩阵的函数tril(A)和tril(A,k)完全相同。 5.2 矩阵分析5.2.1 矩阵结构变换 1. 矩阵的转置 转置运算符是单撇号()。 2. 矩阵的旋转 矩阵的旋转利用函数rot90(A,k)，功能是将矩阵A逆时针旋转90o的k倍，当k为1时可省略。 3. 矩阵的左右翻转 对矩阵实施左右翻转是将原矩阵的第1列和最后1列调换，第2列和倒数第2列调换，…。 对矩阵A实施左右翻转的函数是fliplr(A)。 4. 矩阵的上下翻转 对矩阵实施上下翻转是将原矩阵的第1行和最后1行调换，第2行和倒数第2行调换，…。 对矩阵A实施上下翻转的函数是flipud(A)。 5.2.2 矩阵的逆与伪逆 1. 矩阵的逆 求一个矩阵的逆非常容易。求方阵A的逆可调用函数inv(A)。 例5.4 用求逆矩阵的方法解线性方程组。 命令如下 Ax=b A=[1,2,3;1,4,9;1,8,27]; b=[5,–2,6]; x=inv(A)*b 一般情况下，用左除比求矩阵的逆的方法更有效，即x=A\b。 2. 矩阵的伪逆 A不是一个方阵，或者A是一个非满秩的方阵。找到矩阵B，与A’同型，使得 ABA=A BAB=B 矩阵B为矩阵A的伪逆，记作A+ MATLAB中，求一个矩阵伪逆的函数是pinv(A)。 例5.5 求A的伪逆，并将结果送B。 命令如下： A=[3,1,1,1;1,3,1,1;1,1,3,1]; B=pinv(A) 例5.6 求矩阵A的伪逆。 在MATLAB命令窗口，输入命令： A=[0,0,0;0,1,0;0,0,1]; pinv(A) 5.2.3 方阵的行列式 求方阵A所对应的行列式的值的函数是det(A)。 D=[2,2,-1,1;4,3,-1,2;8,5,-3,4;3,3,-2,2]; det(D) 5.2.4 矩阵的秩 MATLAB中，求矩阵秩的函数是rank(A)。例如，求矩阵D的秩，命令是： D=[2,2,-1,1;4,3,-1,2;8,5,-3,4;3,3,-2,2]; r=rank(D) r = 4 说明D是一个满秩矩阵。 5.2.5 向量和矩阵的范数 1. 计算向量3种常用范数的函数 (1)norm(V)或norm(V,2) 计算向量V的2—范数 (2)norm(V,1) 计算向量V的1—范数 (3)norm(V,inf) 计算向量V的∞—范数 例5.8 已知V，求V的3种范数。 命令如下： V=[-1,1/2,1]; v1=norm(V,1) %求V的1—范数 v2=norm(V) %求V的2—范数 vinf=norm(V,inf) %求∞—范数 2. 矩阵的范数及其计算函数 MATLAB中提供了求3种矩阵范数的函数，其函数调用格式与求向量的范数的函数完全相同 例5.9 求矩阵A的三种范数。 命令如下