matlab第三章课件.pptVIP

第三章 数学运算 3.1 矩阵运算 提取矩阵的对角线元素设A为m×n矩阵，diag(A)函数用于提取矩阵A主对角线元素，产生一个具有min(m,n)个元素的列向量。diag(A)函数还有一种形式diag(A,k)，其功能是提取第k条对角线的元素。(主对角线向上第1条、第2条…,向下第-1条、第-2条…) 构造对角矩阵设V为具有m个元素的向量，diag(V)将产生一个m×m对角矩阵，其主对角线元素即为向量V的元素。diag(V)函数也有另一种形式diag(V,k)，其功能是产生一个n×n(n=m+|k|)对角阵，其第k条对角线的元素即为向量V的元素。 2．三角阵 三角阵又进一步分为上三角阵和下三角阵，所谓上三角阵，即矩阵的对角线以下的元素全为0的一种矩阵，而下三角阵则是对角线以上的元素全为0的一种矩阵。 (2) 下三角矩阵 在MATLAB中，提取矩阵A的下三角矩阵的函数是tril(A)和tril(A,k)，其用法与提取上三角矩阵的函数完全相同。 1．矩阵的转置 转置运算符是单撇号(‘)。 3．矩阵的左右翻转 对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第二列调换，…，依次类推。MATLAB对矩阵A实施左右翻转的函数是fliplr(A)。 4．矩阵的上下翻转 MATLAB对矩阵A实施上下翻转的函数是flipud(A)。 MATLAB提供的部分矩阵求值函数如下表所示 矩阵或向量的范数用来度量矩阵或向量在某种意义下的长度。范数有多种方法定义，其定义不同，范数值也就不同。 附加：矩阵求逆 若方阵A,B满足等式A*B = B*A = I (I为单位矩阵) 则称A为B的逆矩阵，或称B为A的逆矩阵。这时 A，B都称为可逆矩阵(或非奇异矩阵、或满秩矩 阵)，否则称为不可逆矩阵(或奇异矩阵、或降秩 矩阵)。 A=[1 -1 1;5 -4 3;2 1 1]; B=inv(A) B= -1.4000 0.4000 0.2000 0.2000 -0.2000 0.4000 2.6000 -0.6000 0.2000 A*B ans = 1.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 1.0000 1. 矩阵求逆解法 利用求系数矩阵A的逆阵A-1，我们可以得到矩阵求逆解法。对于线性代数方程组Ax=b，等号两侧各左乘A-1，有： A-1Ax=A-1b 由于A-1A=I，故得： x=A-1b A=[1 -1 1;5 -4 3;2 1 1]; b=[2;-3;1]; x=inv(A)*b x = -3.8000 1.4000 7.2000 2. 直接解法 对于线性代数方程组Ax=b，我们也可以运用左除运算符“\”象解一元一次方程那样简单地求解： x=A\b 【例】求下列三阶线性代数方程组的近似解 (1) LU分解 方阵A的特征值λ和其对应的特征向量ν满足下式： A*ν=λ*ν 续表 3.2.3 复数函数 * * 3.1.1 矩阵分析 1．对角阵 只有对角线上有非0元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为1的对角矩阵称为单位矩阵。 对角阵与三角阵 【例】A=[1 2 3；4 5 6]，求其主对角线和第1条对角线的元素 【例】产生一个主对角线元素分别为[1 2 -1 4]的对角阵和第-1条对角线元素分别为[1 2 3]的对角阵。 上三角矩阵 与矩阵A对应的上三角阵B是与A具有相同的行数和列数的一个矩阵，并且B的对角线以上（含对角线）的元素和A对应相等。 求矩阵A的上三角阵的MATLAB函数是triu(A)。另一种形式triu(A,k)，其功能是求矩阵A的第k条对角线以上的元素。例如，提取矩阵A的第2条对角线以上的元素，形成新的矩阵B。 【例】 A=[23 5 71；-45 18 4；92 -4 67]分别求其上三角和下三角矩阵 矩阵的转置与旋转 2．矩阵的旋转 利用函数rot90(A,k)将矩阵A逆时针旋转90o的k倍，当k为1时可省略。 矩阵求值 在MATLAB中，求向量或矩阵的范数