第3讲线性方程组.ppt

第1页，共34页，星期日，2025年，2月5日打开编辑器，创建ff.m文件A=[21-3-1;3107;-124-2;10-15];B=[2354;1-186;-346-7;2354];C=A-BD=A+BE=A*BF=B‘G=A^2H=det(A)K=inv(A)saveffAB在命令窗口键入clearloadffABAB第2页，共34页，星期日，2025年，2月5日二、在命令空间调用1中的A,B矩阵进行如下操作：(1)将A，B扩展为4×8阶的矩阵C.C=[AB](2)提取C中的1,2,4行；3,5,7列构成的新矩阵D.D=C([124],[357])(3)提取C中的3,5列构成新矩阵.E=C(:,[35])(4)建立与A同阶的单位阵,1矩阵,零矩阵.F=eye(size(A))(5)提取A矩阵中的2行3列的元素.?G=A(2,3)返回第3页，共34页，星期日，2025年，2月5日实验三线性方程组实验目的1、熟悉用MATLAB软件求矩阵的秩及矩阵的行简化阶梯形的命令2、学会用MATLAB软件求解线性方程组3、学会用MATLAB软件求矩阵的特征值和特征向量第4页，共34页，星期日，2025年，2月5日一、向量组的秩及相关性rank(A)%求矩阵A的秩rref(A)%将A化为行简化阶梯形，其中单位向量对应的列向量即为极大无关组所含向量，且其它列向量的各分量是用极大无关向量组线性表示的组合系数。第5页，共34页，星期日，2025年，2月5日[例1]设向量组A：（1）求A的秩；（2）求A的一个极大无关组；（3）将其余向量用极大无关组线性表示。第6页，共34页，星期日，2025年，2月5日A1=[2143-11-66-1-22-911-272449];A=A1;r=rank(A)r=3第7页，共34页，星期日，2025年，2月5日A2=rref(A)A2=10-10401-1030001-300000答：（1）r=3，向量组A线性相关；（2）它的一个极大线性无关组是：a1，a2，a4；（3）a3=-a1-a2，a5=4a1+3a2-3a4。第8页，共34页，星期日，2025年，2月5日二、求线性方程组的唯一解或特解Ax=b，当det（A）≠0时，求解方法1.逆矩阵求解：X=inv（A）*b或X=A\b;2.克莱姆法则求解：.第9页，共34页，星期日，2025年，2月5日[例2]解线性方程组解法1：A=[1-11-220-143210-12-12]；第10页，共34页，星期日，2025年，2月5日b=[2;4;-1;-4];r=rank(A)x=inv(A)*br=x=41-201/2或:x=A\bx=