线性方程组及向量组的线性相关性.docVIP

实验11 线性方程组及向量组的线性相关性 实验目的 1．掌握齐次线性方程组的通解的求解方法 2．掌握非齐次线性方程组的特解和通解的求解方法 3．利用向量组线性相关与无关的几何演示，加深理解线性相关与线性无关的概念 4．掌握向量组的秩和最大线性无关组的计算方法 5．理解线性方程组的解的结构的几何意义 实验准备 1．复习线性方程组的解法及其表示； 2．复习向量的基本运算、向量的线性组合及线性表示 3．复习向量组的线性相关（无关）、向量组的秩和最大线性无关组的求法等 4．复习向量空间的基及维数 实验内容 1．线性方程组的解法 2．向量的基本运算及线性组合 3．向量组的线性相关性的判别 4．向量组的秩和最大线性无关组的求法 5．线性方程组的解的结构 软件命令 表11-1 Matlab向量操作命令 函数名称 调用格式 说 明 syms syms 变量名1，变量名2,… 定义符号变量 sym sym(x,…) 定义符号变量 rank rank(A) 求矩阵A的秩 rref rref(A) 求矩阵A的行最简形 \ A\b 求的解 null null(A,r)；注：用于符号计算含义有变 求的基础解系 format format rat 指定有理数输出格式 linsolve x=linsolve(A,B,opt) 解线性方程组 绘图指令plot、surf、patch等 实验示例 【例11.1】 线性方程组的解法 求下列线性方程组的通解 1．； 2．；3．；　 4．； 5．。 【原理】线性方程组的求解可以分为两类：一类是齐次线性方程组；另一类是非齐次线性方程组。无论哪一种类型，都需要先判断出该方程组的解的分布情况：唯一解？无解？无穷多解？，然后在进行求解。具体如下（n为未知量的个数）： 【判别原理】：（1）若，则方程组有唯一解； （2）若，则方程组有无穷多解； （3）若，则方程组无解。 Matlab命令：rank(A),rref(A,b)。 【求解方法】原理：高斯消元法。将增广矩阵化成行最简形，然后求解。 【步骤】：Step1:当有唯一解时，可采用如下两种方法之一 利用命令：x=A\b； 利用命令：x=linsolve(A,b)； 利用LU、QR、Cholesky分解命令：x=U\(L\b)；x=Q\(R\b)；x=R\(R\b) Step2: 当有无穷多解时，可以采用下述方法 利用命令：null(A,r)； 第一步，利用z=null(A,r)求出的基础解系；第二步，利用A\b或者rref([A,b])或者linsolve(A,b)求出的一个特解；第三步，利用的解的结构理论构造处通解。 Step3：（3）当有无解时，可以采用下述方法求出最小二乘解 命令：x=(A*A)\(A*b)； 命令：x=linsolve(A,B,opt) 【程序】： 【程序1】：LinearSysYesorNo.m % 调用方式： str=LinearSysYesorNo(A,B) % 判断线性方程组 Ax=b 是否有解 % 输入：系数矩阵 A ，右端矩阵 B % 输出：解的判别的三种情况：唯一解、无解、无穷多解 中的一种。 【程序2】：SolveLinearSys.m % 调用方式： x=SolveLinearSys(A,b) % 求解线性方程组 Ax=b % 输入：A , b % 输出：x % 说明：利用高斯消元法求解线性方程组 Ax=b 【程序3】：Exm11Demo02.m % 调用前面的函数，具体求解 。 【例11.2】参数方程求解 试就参数的各种情况，讨论下述线性方程组的解的情况： 。 【方法一】：根据方程组的特殊性：m=n，先利用Cramer法则确定唯一解；然后分别讨论无解和无穷多解的情况。 syms a; A=[a 1 1;1 a 1;1 1 a]; B=[1;a;a^2]; det(A) factor(det(A)) solve(det(A)) 因此，由Cramer法则，当，即时，有唯一解； 下面讨论其否定情况，即时原方程组的解的情况： 当时，将代入上面A和B中的a，利用【例12.1】中的函数计算： A=[-2 1 1;1 -2 1;1 1 -2]; B=[1;-2;4]; LinearSysYesorNo(A,B) SolveLinearSys(A,B) 可得此时，方程组无解，最小二乘解为[1;2;0]； 当时，将代入上面A和B中的a，利用【例12.1】中的函数计算： A=[1 1 1;1 1 1;1 1 1]; B=[1;1;1]; LinearSysYesorNo(A,B) SolveLinearSys(A,B) 可得此时，方程组无穷多解，解为[1-c1-c2;c1;c2