实验04向量组的线性相关性及线.ppt

数 学 实 验 空间解析几何与线性代数 (向量组的线性相关性及线性方程组) 成都信息工程学院 计算科学系 实验目的 1、熟悉用MATLAB软件求矩阵的秩及矩阵的行简化阶梯形的命令 2、学会用MATLAB软件求线性方程组 一、向量组的秩及相关性 （1）两个命令 rank（A） 求矩阵A的秩 rref （A）将A化为行简化阶梯形，其中 单位向量对应的列向量即为极大无关组所含向量，且其它列向量的各分量是用极大无关向量组线性表示的组合系数。 [例1] 设向量组A： （1）求A的秩，判断向量组A是否线性相关； （2）求A的一个极大无关组； （3）将其余向量用极大无关组线性表示。 解： A1=[2 1 4 3 -1 1 -6 6 -1 -2 2 -9 1 1 -2 7 2 4 4 9]; A=A1’； r=rank（A） r= 3 A2=rref（A） A2= 1 0 -1 0 4 0 1 -1 0 3 0 0 0 1 -3 0 0 0 0 0 答：（1）r=3，向量组A线性相关； （2）它的一个极大线性无关组是： a1，a2，a4； （3）a3= -a1-a2，a5=4a1+3a2 -3a4 。 二、解线性方程组： 1、Ax=b ，当det（A）≠0时，求解方法 　（1）逆矩阵求解：X=inv（A）*b或 X =A \ b ; 　（2）克莱姆法则求解： 例2. 解线性方程组 解法1： A=[1 -1 1 -2 　　　　　　 2 0 -1 4 　　　　　　 3 2 1 0 　　　　　　 -1 2 -1 2 ] ； b=[2；4；-1；-4]； x=inv（A）*b x= 1 -2 0 1/2 或: x=A \ b x= 1 -2 0 1/2 解法2： D=det（A）； n=size（A）； D1=ones (n(1)，1)； for i=1：n(1) A1=A； A1（：，i）=b； D1（i）=det（A1）； end x=D1/D x= 1 -2 0 1/2 2、一般的线性方程组Ax=b的解法 操作步骤如下： （1）输入系数矩阵A及常数项矩阵b； （2）生成增广矩阵B=[A b]； （3）计算：A的秩r1及B的秩r2； （4）判断：若r1= r2，则转（5）；否则程 序结束； （5）将B化为行简化阶梯形； （6）确定方程组的解。 [例3] 解线性方程组 x1 +5x2 -x3 -x4 = -1 x1 -2x2 +x3+3x4 = 3 3x1+8x2 -x3 +x4 = 1 x1 -9x2 +3x3+7x4 = 7 解： A=[1 5 -1 -1 1 -2 1 3 3 8 -1 1 1 -9 3 7 ] ；