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

实验11线性方程组及向量组的线性相关性 实验目的 掌握齐次线性方程组的通解的求解方法 掌握非齐次线性方程组的特解和通解的求解方法 利用向量组线性相关与无关的几何演示，加深理解线性相关与线性无关 的概念 掌握向量组的秩和最大线性无关组的计算方法 5 .理解线性方程组的解的结构的儿何意义 实验准备 复习线性方程组的解法及其表示； 复习向量的基本运算、向量的线性组合及线性表示 复习向量组的线性相关（无关）、向量组的秩和最大线性无关组的求法等 复习向量空间的基及维数 实验内容 线性方程组的解法 向量的基本运算及线性组合 向量组的线性相关性的判别 4 .向量组的秩和最大线性无关组的求法 线性方程组的解的结构 软件命令 表11-1 Matlab向量操作命令 函数名称 调用格式 说明 syms syms变量名1,变量名2,… 定义符号变量 sym sym(，x1,…) 定义符号变量 rank rank(A) 求矩阵A的秩 rref rref (A) 求矩阵A的行最简形  \ A\b 求人兀=方的解 null null (A, r);注：用于符号计算含义有变 求Ax = 0的基础解 系 format format rat 指定有理数输出格 式 1 insolve x=l i nsol ve (A, B, opt) 解线性方程组 Ax = B 绘图指令plot、surf、patch等 实验示例 求下列线性方程组的通解x + 2y + 3z = 01. 4 求下列线性方程组的通解 x + 2y + 3z = 0 1. 4兀+ 5y + 6z = 0； 7 兀+ 8y + 9z = 0 2x-y + 3z = l 4. 2x +2z = 6 ； 4x + 2y+ 5z = 7 2x-y + 3z = 14x-2y + 5z = 42x-y+ 4z = 0 2x-y + 3z = 1 4x-2y + 5z = 4 2x-y+ 4z = 0 ■ 2. 2x + 4y + 6z = 0； 3. 3兀+ 6y+ 9z = 0 x + 2y + 6z + 7w = 67 x + 9y+ 7z + 14w = 117 x + 2y + 5z + 9w = 79 3兀+ 6y+ 16z + 25w = 225 【原理】线性方程组Ax = b的求解可以分为两类：一类是齐次线性方程组；另 -类是非齐次线性方程组。无论哪一种类型，都需要先判断出该方程组的解的分布情况：唯一解？无解？无穷多解？，然后在进行求解。具体如下(n为未知量 的个数): 【判别原理】:(1)若R(A) = R(A,b) = n.则方程组Ax = b有唯一解; (2)若R(A)=R(A,b)n,则方程组Ax = b有无穷多 解； (3)若R(A)^R(A,b)9则方程组Ax = b无解。 Mat lab 命令：rank (A), rref (A, b)。 【求解方法】原理：高斯消元法。将增广矩阵化成行最简形，然后求解。 【步骤】：Stepl:当人工=方有唯一解时，可采用如下两种方法之一 利用命令：x二A\b； 利用命令：x=l insolve (A, b)； 利用 LU、QR、Cholesky 分解命令:x=U\(L\b): x二Q\(R\b)； x二R\(R \b) Step2:当Ax = b有无穷多解时，可以采用下述方法 ?利用命令：null(A,r)； 第一步，利用z=nul 1 (A,5 r )求出4r = 0的基础解系； 第二步，利用 A\b 或者 rref ([A, b])或者 1 insolve (A, b) 求出Ax = b的一个特解； 第三步，利用Ax = b的解的结构理论构造处通解。 Step3： (3)当Ax = b有无解时，可以采用下述方法求出最小二乘解 命令:x=(A*A)\(A*b)； 命令：x=l insolve (A, B, opt) 【程序】： 【程序 1】：Li nearSy sYesorNo. m % 调用方式：str=LinearSysYesorNo (A, B) %判断线性方程组Ax二b是否有解 %输入：系数矩阵A ,右端矩阵B %输出：解的判别的三种情况：唯一解、无解、无穷多解屮的 一种。 【程序 2】：SolvoLinoarSys. m % 调用方式：x=SolveLi nearSys (A, b) %求解线性方程组Ax二b %输入：A , b %输岀：x %说明：利用髙斯消元法求解线性方程组Ax=b 【程序 3 ］： ExmllDemo02. m %调用前面的函数，具体求解。 【例11.2］参数方程求解 试就参数2的各种情况，讨论下述线性方程组的解的情况： y + z = 1 x+y + Az = A2 【方法一】：根据方程组的特殊性：m二n