线性代数实验指导.PDFVIP

线性代数实验指导 实验四 线性方程组求解  实验目的及要求：熟练应用MATLAB 命令进行线性方程组Ax=b 的求解， 包括齐次线性方程组和非齐次线性方程组的求解。 实验设备：数学软件MATLAB 、计算机 实验内容： 1、使用MATLAB 求线性方程组的步骤 (1)对齐次线性方程组Ax=0 来说，设A 是m 行n 列的矩阵 1) 当rank(A)=n 时，方程组有唯一解，且x=0; 2) 当rank(A)n 时，方程组有无穷多解，即存在非零解，使用Matlab 可以求出无穷解的线性表达式，这个表达式由基础解系线性表示，其中基础解系 含向量的个数为n-rank(A) ，求基础解系的函数为null(A,’r’) (2)对非齐次线性方程组Ax=b 来说， 1）当m=n ，且rank(A)=rank(Ab)=n 时，方程有唯一解x=A-1b （Ab 为A 的增广矩阵），使用Matlab 可以立即求出方程组的解为x=inv(A)*b 2)当rank(A)=rank(Ab)n 时，方程组Ax=b 有无穷解，此时用null(A,’r’ ) 求出对应齐次方程组的基础解系，再用pinv(A)*b 求出一个特解，则可得出方程 组的通解。 3 ）当rank(A) ≠rank(Ab) ，方程组Ax=b 无解。 2 、操作过程 （1）在命令行中输入线性方程组的系数矩阵A 和向量b （齐次方程b=0 可 以不用输入） （2 ）在命令行中通过调用相关函数进行运算 例1 求齐次线性方程组： 1 x +2x +3x 0 ⎧1 2 3 ⎪ ⎨2x −x −x 0 1 2 3 ⎪ x −2x −2x 0 ⎩1 2 3 解：在命令行输入： 由此可知该方程组的系数矩阵的秩为3，与变量的个数n=3 相等，因此方程 组只有零解，即x1=x2=x3=0 。 例2 求齐次线性方程组 + − − x x x x 0 ⎧1 2 3 4 ⎪ x + x +x ⎨2 1 2 3 4 0 ⎪ x x x x + + +2 0 ⎩1 2 3 4 解：在命令行输入： 因为系数矩阵的秩r=2n=4 因此方程组有无穷解，继续输入命令 2 得到n-r=4-2=2 个基础解系，故齐次方程组的解为： −1 −0.5 ⎡ ⎤ ⎡ ⎤ ⎢ ⎥ ⎢ ⎥ 1 0 ⎢ ⎥ ⎢ ⎥ x k =+k 1 2 ⎢ ⎥ ⎢ ⎥ − .5 0 1 ⎢ ⎥ ⎢ ⎥ 0 1 ⎣ ⎦ ⎣ ⎦ 例3 解非齐次线性方程组 x 2x 3x 2 ⎧ + + ⎪1 2 3 2x x x 1 ⎨ − − 1 2 3 ⎪ 2 2 1