线性方程组的解结构.ppt

《线性方程组的解结构.ppt》由会员分享，可在线阅读，更多相关《线性方程组的解结构.ppt（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、湖湖北北科科技技职职业业学学院院线性方程组的解结构线性方程组的解结构齐次线性方程组的解结构齐次线性方程组的解结构非齐次线性方程组的解结构非齐次线性方程组的解结构湖湖北北科科技技职职业业学学院院齐次线性方程组的解结构齐次线性方程组的解结构例例1.判别方程组判别方程组有无非零解有无非零解,若有若有,写出其通解写出其通解.解解 在在MATLAB中输入该方程组的系数矩阵中输入该方程组的系数矩阵A并将它并将它化为最简行阶梯形矩阵化为最简行阶梯形矩阵,所用命令如下所用命令如下:A=1 2-1;2 5 2;1 4 7;1 3 3;rref(A)运行结果为运行结果为湖湖北北科科技技职职业业学学院院ans=1

2、0 -9 0 1 4 0 0 0 0 0 0由阶梯形矩阵可知由阶梯形矩阵可知R(A)=2 A=1 1 1 1 1;3 2 1 1-3;0 1 2 2 6;5 4 3 3-1;format rat B=null(A,r)%求基础解系求基础解系湖湖北北科科技技职职业业学学院院B=1 1 5 -2 -2 -6 1 0 0 0 1 0 0 0 1 syms k1 k2 k3%定义符号参数定义符号参数 X=k1*B(:,1)+k2*B(:,2)+k3*B(:,3)X=k1+k2+5*k3 -2*k1-2k2-6k3 k1 k2 k3湖湖北北科科技技职职业业学学院院即即为方程组的通解为方程组的通解,其中其

3、中k1,k2,k3为任意实数为任意实数.湖湖北北科科技技职职业业学学院院非齐次线性方程组的解结构非齐次线性方程组的解结构例例3.求解方程组求解方程组解解 在在MATLAB中输入系数矩阵及常数列向量中输入系数矩阵及常数列向量,并检验并检验系数矩阵是否逆系数矩阵是否逆,所用命令及结果如下所用命令及结果如下 A=2 1 1;3 1 2;1-1 0;b=3 3-1;det(A)%检验检验A是否可逆是否可逆ans=2系数矩阵行列式值等于系数矩阵行列式值等于2,是可逆是可逆的的,则可以用矩阵相除来求解则可以用矩阵相除来求解.X=AbX=1 2 -1即是原方程组的解即是原方程组的解.湖湖北北科科技技职职业业

4、学学院院例例4.求解方程组求解方程组解解 先用先用MATLAB函数函数null求出对应的齐次线性方程组求出对应的齐次线性方程组的基础解系的基础解系,再利用其系数矩阵的上、下三角阵求出方再利用其系数矩阵的上、下三角阵求出方程组的一个特解程组的一个特解,这样即可得到该方程组的通解这样即可得到该方程组的通解,程序程序如下如下:A=1 1-3-1;3-1-3 4;1 5-9-8;b=1 4 0;format rat C=null(A,r);%求基础解系求基础解系 L,U=lu(A);%A=LU,L为为上三角阵上三角阵,U为下三角阵为下三角阵 X0=U(Lb)%用用LU求出一个齐次方程的特解求出一个齐次方程的特解湖湖北北科科技技职职业业学学院院 syms k1 k2 X=k1*C(:,1)+k2*C(:,2)+X0运行结果为运行结果为X0=0 0 -8/15 3/5 X=3/2*k1-3/4*k2 3/2*k1+7/4*k2 k1-8/15 k2+3/5即即为该非齐次方程组的通解为该非齐次方程组的通解,其中其中k1,k2为任意实数为任意实数.