线性代数实验三向量组的线性相关性和方程组的通解.docVIP

实验三　　向量组的线性相关性与方程组的通解 一、实验目的 1掌握matlab软件分析向量组线性相关的方法； 2 掌握matlab软件求解线性方程组通解的各种方法； 3通过matlab软件进一步理解和认识齐次线性方程组解空间的概念； 二、实验指导 [R,s]=rref(A)----------将矩阵A的行最简行形矩阵赋给R，s是一个行向量，它的元素是行最简行形矩阵的非零行的行数，也称为基准元素； length(s) ---------------计算向量s的长度，即s的维数； null(A,‘r’) ----------计算齐次线性方程组Ax=0的基础解系； x0=A\b-----------求非齐次线性方程组Ax=b的一个特解x0； find(s) -----------计算向量s中非零元素的下标； subs(A,k,n) ---------将A中的所有符号变量k用数值n来替代. D=det(A) ％ 计算矩阵A的行列式D f=factor(D) ％ 对行列式D进行因式分解 X=solve(D) ％ 求解方程“D＝0”的解 三、实验内容 1、求非齐次线性方程组的通解 A=[2 4 -1 4 16;-3 -6 2 -6 -23;3 6 -4 6 19;1 2 5 2 19] b=[-2;7;-23;43] [R,s]=rref([A,b]) %将增广矩阵的行最简形矩阵赋给R [m,n]=size(A) x0=zeros(n,1) %将特解x0初始化为n维零向量； r=length(s) x0(s,:)=R(1:r,end) %将矩阵R的最后一列按基准元素的位置给特解x0赋值 null(A,r) 例1将矩阵化为最简行阶梯形矩阵。 程序运行结果 A=[7 1 -1 10 1;4 8 -2 4 3;12 1 -1 -1 5]; rref(A) 例2 求解线性方程组 程序结果如下： B=[1 3 -2 4 1 7;2 6 0 5 2 5;4 11 8 0 5 3;1 3 2 1 1 -2] rref(B) B = 1 3 -2 4 1 7 2 6 0 5 2 5 4 11 8 0 5 3 1 3 2 1 1 -2 ans = 1.0000 0 0 -9.5000 4.0000 35.5000 0 1.0000 0 4.0000 -1.0000 -11.0000 0 0 1.0000 -0.7500 0 -2.2500 0 0 0 0 0 0 所以原方程组等价于方程组 故方程组的通解为： ，其中 例3 解方程 。 解： ％ 求解符号行列式的方法，编写M文件如下la01.m clear all ％ 清除各种变量 syms x ％ 定义x为符号变量 A=[3 2 1 1;3 2 2-x^2 1;5 1 3 2;7-x^2 1 3 2] ％ 给矩阵A赋值 D=det(A) ％ 计算矩阵A的行列式D f=factor(D) ％ 对行列式D进行因式分解 X=solve(D) ％ 求解方程“D＝0”的解 在matlab命令窗口中输入： la01.m 运行结果为： A = [ 3, 2, 1, 1] [ 3, 2, 2-x^2, 1] [ 5, 1, 3, 2] [ 7-x^2, 1, 3, 2] D = -6+9*x^2-3*x^4 f = -3*(x-1)*(x+1)*(x^2-2) X = -1 1 2^(1/2) -2^(1/2) 实验内容 1 将矩阵化为最简行阶