实验三线性方程组的求解.doc

实验三 线性方程组的求解 1、 求齐次方程组的基础解系和解空间的维数 对于齐次线性方程组，可用null求其一个基础解系，用size求得基础解系的维数。 例2 求解下列方程组 相应的MATLAB代码为 A=[1,1,-1,-1;2,-5,3,2;7,-7,3,1]; A=sym(A); B=null(A) B = [ 3/4, -1/4] [ 1, 0] [ 0, 1] [ 7/4, -5/4] size(B,2) ans = 2 故方程组的基础解系为：,且基础解系的维数是2. 2、 非齐次线性方程组的直接求解 对于非齐次线性方程组，当时，有唯一解，此时MATLAB中可以直接使用 “\” 解此类线性方程组。 例1 求解下列方程组 相应的MATLAB代码为 A=[2,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6];b=[8,9,-5,0]; x=A\b x = 3.0000 -4.0000 -1.0000 1.0000 所以方程组有唯一解：。 3、 不定线性方程组的通解 当线性方程组的方程个数少于未知数个数时，称这个方程组为不定方程组。这样的方程组通常有无穷组解，MATLAB在求这类方程时只给出其中的一组解，也就不再给出错误信息。此时可以利用命令求系数矩阵的基础解系，引入参数即可得出原方程组的通解。 求解不定方程组 相应的MATLAB代码为 A=[1,2,3,-1;3,2,1,-1;1,-2,-5,1]; A=sym(A) ; b=[2,4,0]; x=A\b Warning: System is rank deficient. Solution is not unique. In sym.mldivide at 41 x = 1 0 0 -1 B=null(A) B = [ 1, 0] [ 0, 1] [ 1, 0] [ 4, 2] syms c1 c2 X=c1*B(:,1)+c2*B(:,2)+A\b Warning: System is rank deficient. Solution is not unique. In sym.mldivide at 41 X = c1+1 c2 c1 4*c1+2*c2-1 4、 含参数的线性方程组解的讨论 例4 取何值时，下列线性方程组有唯一解、无解或有无穷多解？并在无穷多解时求出其解。 解：（1）求该方程组的系数行列式，并求其根 syms a; A=[a,1,1;1,a,1;1,1,a]; p=det(A) p = a^3-3*a+2 solve(p,a) ans = -2 1 1 （2）当时，线性方程组有唯一解 （3）当时，分别求系数矩阵与增广矩阵的秩 A=[-2,1,1;1,-2,1;1,1,-2]; b=[1,-2,4]; rank(A) ans = 2 rank([A,b]) ans = 3 由于，所以当时，线性方程组无解； （4）当时，求其增广矩阵的最简形，然后写出其通解 A=[1,1,1;1,1,1;1,1,1]; b=[1,1,1]; rref([A,b]) ans = 1 1 1 1 0 0 0 0 0 0 0 0 由于，所以当时，线性方程组有无穷多解； 此时，其所有解为 。