实验3-向量组的相关性及方程组的通解..docVIP

3 向量组的相关性及方程组的通解 分析向量组线性相关性的方法； 求解线性方程组通解的各种方法； 例3.0求向量组，，，，的秩，并判断其线性相关性. 解 A=[1 -2 2 3;-2 4 -1 3;-1 2 0 3;0 6 2 3;2 -6 3 4]; k=rank(A) 结果为 k = 3 由于秩为3小于向量组所含向量个数，因此向量组线性相关. 例3.1求非齐次线性方程组 的通解。 解：在MATLAB命令窗口，输入以下命令： A=[2,4,-1,4,16;-3,-6,2,-6,-23;3,6,-4,6,19;1,2,5,2,19]; % 输入系数矩阵A b=[-2;7;-23;43]; % 输入常数列向量b [R,s]=rref([A,b]) % 把增广矩阵的最简行阶梯矩阵赋给R % 而R的所有基准元素在矩阵中的列号构成了行向量s 计算结果为： R = 1 2 0 2 9 3 0 0 1 0 2 8 0 0 0 0 0 0 0 0 0 0 0 0 s = 1 3 程序la06.m给出非齐次方程组的通解。 % 求齐次线性方程组的通解 clear A=[2,4,-1,4,16;-3,-6,2,-6,-23;3,6,-4,6,19;1,2,5,2,19]; % 输入系数矩阵A b=[-2;7;-23;43]; % 输入常数列向量b [R,s]=rref([A,b]); % 把增广矩阵的最简行阶梯矩阵赋给R % 而R的所有基准元素在矩阵中的列号构成了行向量s [m,n]=size(A); % 矩阵A的行数、列数赋给了变量m、n x0=zeros(n,1); % 将特解x0初始化为N维零列向量 r=length(s); % 矩阵A的秩赋给变量r x0(s,:)=R(1:r,end); % 将矩阵R的最后一列按基准元素的位置给特解x0赋值 disp(非齐次线性方程组的特解为：) x0 % 显示特解x0 disp(对应齐次线性方程组的基础解系为：) x=null(A,r) % 得到齐次线性方程组Ax＝0的基础解系x 在MATLAB命令窗口中输入： la06 运算结果为： 非齐次线性方程组的特解为： x0 = 3 0 8 0 0 对应齐次线性方程组的基础解系为： x = -2 -2 -9 1 0 0 0 0 -2 0 1 0 0 0 1 则方程组的通解为： 齐次线性方程组的特解还可以用 Matlab的矩阵左除运算来求得， 直接在MATLAB命令窗口输入以下命令： 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]; x0=A\b % 用矩阵左除运算求得方程组特解x0 x=null(A,r) % 得到齐次线性方程组Ax＝0的基础解系x 运算结果为： Warning: Rank deficient, rank = 2 tol = 4.3099e-014. x0 = 0 0 7.3333 0 0.3333 x = -2 -2 -9 1 0 0 0 0 -2 0 1 0 0 0 1 方程组的通解为： 例3.2 已知向量组 ，，，，， 求出它的最大无关组，并用该最大无关组来线性表示其它向量。 解：用笔计算的过程为： 编写Matlab程序la0