MATLAB科学计算5.pptVIP

3、矩阵转化为列向量 A( : ) 每一列顺序链接成列向量 删除第m行 A(m, : )=[ ] 删除第n列 A( : ,n)=[ ] 4、删除某行或某列 A是矩阵，得到它的对角线元素 5、diag(A) A是向量，得到以它为对角线元素的矩阵 U=triu(A,k) k-对角线位置 k=0, A主对角线 k0，A主对角线上方的第k条对角线 k0，A主对角线下方的第k条对角线 Upper triangular 提取A某条对角线及其以上元素 构造上三角阵U 6、 L=tril(A,k) k-对角线位置 k=0, A主对角线 k0，A主对角线上方的第k条对角线 k0，A主对角线下方的第k条对角线 Lower triangular 提取A某条对角线及其以下元素 构造下三角阵L 7、 函数名 功能 flipud(A) 上下翻转 fliplr(A) 左右翻转 rot90(A,k) 旋转k个90o, k为整数 k0逆时针旋转， k0顺时针旋转 例10 （1）构造一个行向量，元素从1到6，步长为1； （2）以其为对角元素，构造对角阵 例11 求A的第 2 行至第 4 行、 第 2 列至第 3 列 构成的矩阵 例12 交换A的第 1 行和第 3 行 H4= 例13　提取4阶Hilbert矩阵的上三角阵元素 U=triu(A,k) k=0, 1, -1 H4= 例14　提取4阶Hilbert矩阵的下三角阵元素 L=tril(A,k) k=0, 1, -1 已知 编写m文件，求M特征值中模最大者 其中： D是对角阵，其主对角线元素对应A的主对角线元素 -L为严格下三角阵，其主对角线下方元素是A的对应元素 -U为严格上三角阵，其主对角线上方元素是A的对应元素 * 例1 创建矩阵A，并输出维数、第一列、第二行 方法一 x=a+bi 例2 创建复矩阵 例2 创建复矩阵 x=a+bi 方法一 方法二 例3 构造一个从1到10的等比数列,节点数是4 q=logspace(log10(a),log10(b),n) 例4 创建如下7*7矩阵 例5 构造如下方程组的系数矩阵 repmat(A,m,n) spdiags(A,d,m,n) 6.2 数值矩阵的矩阵运算和数组运算 一、矩阵运算 遵循线性代数的矩阵运算法则 1、矩阵加法 C=A + B 一、矩阵运算 当A,B都是矩阵时，要求维数一致 当A为数时， C=A*ones(size(B))+B 2、矩阵减法 C=A - B 一、矩阵运算 当A,B都是矩阵时，要求维数一致 当A为数时， C=A*ones(size(B))-B 3、矩阵乘法 C=A * B 一、矩阵运算 要求A的列数等于B的行数 4、矩阵除法 线代: AB=C, 则B=A-1C MATLAB: B=A\C 或B= inv(A)*C 左除 反斜杠” \” 4、矩阵除法 线代: AB=C, 则A=CB-1 MATLAB: A=C/B 或A=C*inv(B) 右除 斜杠” /” C=mpower(A, n) 或C=A^n 5、方阵乘幂 matrix power A是方阵 n是正整数，输出方阵n个连乘的结果 n是负整数，输出方阵n个连乘后的逆矩阵 eig(A) 返回一个列向量 其元素是A的特征值 6、矩阵特征值和特征向量 eigenvalues and eigenvectors 6、矩阵特征值和特征向量 [v,d] =eig(A) v-以特征向量为列向量的矩阵 d-特征值组成的对角阵 7、方阵的特征多项式 p=poly(A) A:矩阵A必须是方阵； p:A的特征多项式系数向量，按降幂排序。 结合使用roots(p)可求A的特征值。 矩阵运算 数学符号 MATLAB运算符号 行列式 |A| det(A) 转置 AT (或A’) A.’ 共轭 conj(A) 共轭转置 A’ 秩 r(A) rank(A) 逆矩阵 A-1 inv(A)或A^(-1) 6.2 数值矩阵的矩阵运算和数组运算 二、数组运算 把矩阵看做一组数据，在两个矩阵的对应元素之间按运算符号进行数与数的四则运算； 类似于批处理。 1、数组加法 C=A + B 二、数组运算 2、数组减法 C=A - B 3、数组乘法 C=A .* B 当A是数时 运算法