数学实验课件 第6章6.1.pptx

澡身浴德修业及时6.1线性方程组设有n个未知数m个方程的线性方程组可以写成以向量x为未知元的向量方程(6-1)其中，若b=0，则式（6.1）是齐次线性方程组；若b≠0，则式（6.1）是非齐次线性方程组.定理6.1n元线性方程组（1）无解的充分必要条件是；（2）有唯一解的充分必要条件是；（3）有无穷多解的充分必要条件是.克拉默法则若n元线性方程组系数矩阵A是方阵，且，则方程组有唯一解，且,其中是把系数矩阵A中第j列用方程组右端的常数项b代替后的矩阵.在线性方程组求解的过程中，会用到行列式的介绍、方阵求逆、矩阵的秩和矩阵阶梯化等问题.这些问题都可由MATLAB进行处理.6.1.1行列式的计算在MATLAB中，det(A)可以计算方阵A的行列式.例6.1?计算方阵的行列式.解clearA=[31-1;-513;201];det(A)?运行结果如下：ans=16可得=16.?例6.2计算矩阵的行列式.解clear;symsabcA=[abc;aa+ba+b+c;a2*a+b3*a+2*b+c];det(A)ans=a^3可得.6.1.2矩阵阶梯化在MATLAB中，rref(A)可以求矩阵A的行最简形.例6.3已知矩阵，求A的行最简形．解formatratA=[3102;1-12-1;13-44];rref(A)ans=101/21/401-3/25/40000可知矩阵A的行最简形是6.1.3矩阵的秩在MATLAB中，可以用rank(A)直接求矩阵A的秩，也可以用rref(A)间接求矩阵A的秩.例6.4求例6.3中A的秩.解方法1利用rank法A=[3102;1-12-1;13-44];rank(A)ans=2方法2利用rref法已知例6.3中行最简形B的非零行的行数是2，可知B的秩为2；因为A和B等价，可知A的秩也为2.6.1.4逆矩阵在MATLAB中，可以用inv(A)直接求方阵A的逆矩阵，也可以用rref(A,E)求方阵A的逆矩阵，其中E是与A同阶的单位矩阵.例6.5求例6.1的逆矩阵.解方法1利用inv法formatratA=[31-1;-513;201];inv(A)ans=1/16-1/161/411/165/16-1/4-1/81/81/2方法2利用rref法B=[A,eye(3)];rref(B)ans=1001/16-1/161/401011/165/16-1/4001-1/81/81/2可知矩阵A的逆矩阵存在，且澡身浴德修业及时谢谢