大学数学实验1_MATLAB操作基础教程.pptVIP

例1.2-10 产生5阶随机方阵 A，其元素为 [10,90]区间的随机整数，然后判断A的元素是否能被3整除。 生成5阶随机方阵A，输入命令： A=fix((90-10+1)*rand(5)+10) (2) 判断A的元素是否可以被3整除： P=rem(A,3)==0 其中，rem(A,3)是矩阵A的每个元素除以3的余数矩阵。此时，0被扩展为与A同维数的零矩阵，P是进行等于(==)比较的结果矩阵。 例1.2-11 建立矩阵A，然后找出大于4的元素的位置。 建立矩阵A，输入命令： A=[4,-65,-54,0,6;56,0,67,-45,0] (2) 找出大于4的元素的位置： find(A4) %返回索引值 2.4 矩阵分析 2.4.1 对角阵与三角阵 1. 对角阵 只有对角线上有非 0 元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为 1 的对角矩阵称为单位矩阵。 ?V=[1?2?3]V?=?????1?????2?????3?diag(V,3)ans?=?????0?????0?????0?????1?????0?????0?????0?????0?????0?????0?????2?????0?????0?????0?????0?????0?????0?????3?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0?????0 ?B=[1?2?3;4?5?6]B?=?????1?????2?????3?????4?????5?????6 ?diag(B)ans?=?????1?????5 ?diag(B,2)ans?=?????3 ?diag(B,1)ans?=?????2?????6 提取矩阵的对角线元素 命令：diag(A),diag(A,k) 构造对角阵V为向量 命令：diag(V),diag(V,k)， diag(1:3)*A %将A的每一行增大行数倍数 ans = 1 2 3 8 10 12 21 24 27 2. 三角阵 三角阵又进一步分为上三角阵和下三角阵，所谓上三角阵，即矩阵的对角线以下的元素全为 0 的一种矩阵，而下三角阵则是对角线以上的元素全为0的一种矩阵。 triu(A) ：求A上三角阵 triu(A,k):求矩阵A的第k条对角线以上的元素 tril(A):求A下三角阵 tril(A,k)：求矩阵A的第k条对角线以下的元素 a’ 求矩阵a的转置 det(a) 求矩阵a的行列式 b=eig(a) 求矩阵a的特征值,构成向量b [v,d]=eig(a) a的特征值构成对角阵d，特征向量构成v的列向量 inv(a)或a ^ (-1) 求矩阵a的逆矩阵 rank(a) 求矩阵a的秩 trace(a) 求矩阵a的迹（对角线元素之和） rref(a) 求矩阵a的行最简形 null(A,’r’) 求系数矩阵为A的齐次方程的基础解系 pinv(A)*b 求非齐次方程的特解 rot90(A,k) 将矩阵A旋转90o的k倍 fliplr(A) 将A左右翻转 flipud(A) 将A上下翻转 3.其它命令 例2.3-4 用求逆矩阵的方法解线性方程组。 Ax=b 其解为： x=A-1b 命令为： a=[2,-3,1;8,3,2;45,1,-9]; b=[4;2;17]; x=inv(a)*b 例2.3-5 用求特征值的方法解方程： 3x5-7x4+5x2+2x-18=0 p=[3,-7,0,5,2,-18]; A=