MATLABSIMULINK实用教程 张化光 刘鑫蕊 孙秋野 第2章数组、矩阵及其运算新.pptVIP

MATLAB/SIMULINK实用教程 张化光 孙秋野 刘鑫蕊 编著 第2章数组、矩阵及其运算 2.1 数组的创建 2.2 矩阵的代数运算 2.2.1 pow2函数 2.2.2 加、减运算 2.2.3 乘法 2.2.4 集合运算 2.2.5 除法运算 2.2.6 矩阵乘方 2.3 矩阵的关系运算 2.3.1 矩阵的比较关系运算 2.3.2 矩阵的逻辑关系运算 2.4 矩阵运算 2.5 符号矩阵运算 2.6 高 维 数 组 2.6.1 高维数组的创建 2.6.2 高维数组的标识 2.6.3 高维数组的操作 2.7 非数和空数组 2.7.1 非数NaN 2.7.2 空数组 2.8 矩 阵 分 解 2.8.1 Cholesky分解 2.8.2 LU分解 2.8.3 QR分解 2.8.4 Schur分解 2.8.5 实Schur分解转化成复Schur分解 2.8.6 特征值分解 2.8.7 奇异值分解 2.8.8 广义奇异值分解 2.8.9 特征值问题的QZ分解 2.8.10 海森伯格形式的分解 2.9 特征值与特征向量 2.9.1 特征值与特征向量的求法 2.9.2 提高特征值的计算精度 2.9.3 复对角矩阵转化为实对角矩阵 2.9.4 正交基 2.10 小 结 求解维数函数：n = ndims(A) 求解大小函数：n = size (A) 求解长度函数：n = length (A) 【例2-58】　维数、大小和长度 2.6.2 高维数组的标识 （1）数组元素对称交换函数flipdim(A,dim) 返回 A按数值dim翻转的结果，当dim=1时，按“行平分面”交换对称位置上的元素，当dim=2时，按“列平分面”交换对称位置上的元素，flipdim(A,1) = flipud(A), flipdim(A,2) =fliplr(A)。 【例2-59】　数组元素对称交换指令flipdim的使用示例。 2.6.3 高维数组的操作 （2）数组的“维序号左移”重组函数shiftdim(X,n) 返回数组X维数移动n后重组的结果。 当n是正数时，维号左移n位；当n是负数时，维号右移n的绝对值位。 【例2-60】　数组的“维序号左移”重组。 （3）广义非共轭转置函数permute（A,order）重组A的维数，使A按照order指定的向量顺序重组。 【例2-61】　广义非共轭转置。 （4）孤维函数B = squeeze(A) 返回一个向量B，它与A数组元素一致。 【例2-62】　“孤维”的撤销和降维。 【例2-63】　赋“空阵”值操作。 MATLAB规定，形如0/0，inf/inf，0*inf（inf代表无穷）等运算都会产生非数NaN（Not-a-Number）。 【例2-64】　非数的产生和性质演示。 【例2-65】　非数元素的寻访。 2.7.1 非数NaN 当数组的某些维数为0时，称该数组为空数组。在复杂的计算机中有时会产生空数组，也可以利用空数组对其他数组进行赋值。但要注意： （1）空数组与全零元素数组不是一回事； （2）空数组有时确实是数组，用which、who、whos可以看出； （3）可以用命令isempty来判断一个数组是否为空数组。 【例2-66】　关于空数组的算例。 2.7.2 空数组 常见的矩阵分解有Cholesky分解、LU分解、QR分解、Schur分解、复Schur分解、特征值分解、奇异值分解、广义奇异值分解、特征值问题的QZ分解和海森伯格形式的分解等。 命令格式： R = chol(X) %如果X为n阶对称正定矩阵，则存在一个实的非奇异上三角阵R， %满足R*R = X；若X非正定，则产生错误信息。 [R,p] = chol(X) %不产生任何错误信息，若X为正定阵，则p=0，R与上相同；若X %非正定，则p为正整数，R是有序的上三角阵。 【例2-67】　Cholesky分解示例。 2.8.1 Cholesky分解 矩阵的三角分解又称LU分解，它的目的是将一个矩阵分解成一个下三角矩阵L和一个上三角矩阵U的乘积，即A=LU。 命令格式： [L,U] = lu(X) %U为上三角阵，L为下三角阵或其变换形式，满足LU=X。 [L,U,P] = lu(X) %U为上三角阵，L为下三角阵，P为单位矩阵的行变换矩阵，满足LU=PX。