线性代数实践_MATLAB幻灯片.pptVIP

求解三元线性方程组的例 → 再由两式联立解得： x=1.4286 y=-1.8571 z=-5.4286 用图解方法ezplplot3画此三个方程的平面， 用的就是例5.3的程序，得到 ： 按例6.2画出的三个平面 三元一次方程的几何意义 这三个方程的几何意义为空间的三个平面，两两消去z，意味着求两者的交线，得到的两个二元一次方程对应于求得的两根交线，最后由交线得到交点，该交点就是方程的解（见图6.3）。当两个平面平行，没有交线，或者两根交线平行，没有交点时，方程组就不相容，因而无解；同样也可能有无穷个解的情况，读者可自行思考想象。 6.2 初等行变换和高斯消元子程序 执行以下三种行运算不会影响方程的相等关系，故不会影响方程组的解。矩阵C的行初等变换的MATLAB表示式为： (1)。行交换： c([i,j],:)=c([j,i],:) (2)。行乘数： c(i,:)=k*c(i,:) (3)。行相加： c(j,:) = c(j,:)+ k*c(i,:) 消元法的要点 利用上述初等变换，在保持矩阵等价性和其它各行不变的前提下，把矩阵第j行中第q列的元素a(j,q)变换为零。 做法：采用第三种变换，并取k=-a(j,q)/a(i,q)。此时第j行的数组将成为: a(j,:) = a(j,:) + -a(j,q)/a(i,q)*a(i,:) （6.9） 而其中的第q个元素的取值将成为： a(j,q) = a(j,q) + -a(j,q)/a(i,q)*a(i,q) = 0 （6.10） 其中的冒号’：’意味着要对全行的所有元素都做一次乘法和一次加法，相当繁琐。 消元子程序gauss.m。 function [B]=gauss(A,i,j,q) % --------------------------------- % A为输入矩阵，B为变换后的输出矩阵 % i为基准行的行号 % j为待变换行的行号 % q为基准元的列号，即A(i,q)为基准元，A(j,q)为待消元 % x = A(i,:); y=A(j,:); % 取出A的第i,j两行命名为x,y， z = y - y(q)/x(q)*x; % 实现（6.2.5）式的运算 A(j,:)=z; % 把结果赋值给A第j行， B=A; % 将A作为输出变元B 6.3 行阶梯形式 具有以下三个特点的矩阵称为行阶梯形式： 1。所有非零行都处在全零行的上方； 2。行首非零(枢轴)元素的列号比它的上方所有各行的枢轴元素的列号都要大； 3。各枢轴元素所处的列中，在枢轴元素下方的所有元素均为零。 如果此矩阵还满足以下两个额外特点， 4。各行枢轴元素是其列中的唯一非零元素； 5。这些枢轴元素都等于1。 则称之为简化行阶梯形式。 例6.3 行阶梯方程组的解法 用逐次回代法，可以由下而上依次求得 x4=2, x3=2+x4/2=3, x2=-4-8x3=-28, x1=3x2-5x3+2x4=-39。 消元法生成行阶梯形式 步骤1 把主对角线下方第一列元素消为零 连续(n-1)次调用gauss子程序来实现。即用 a1=gauss(a,1,2,1) a1=gauss(a1,1,3,1) ... a1=gauss(a1,1,n,1) 步骤2 把主对角线下方第二列元素消为零 a2=gauss(a1,2,3,2) a2=gauss(a2,2,4,2) ... 步骤n-1 消元法生成简化行阶梯形式 与行阶梯形式反向进行 步骤1 把主对角线上方第n列元素消为零 a1=gauss(a1,n,n-1,n) ... a1=gauss(a1,n,1,n) 步骤2 把主对角线上方第n-1列元素消为零 a1=gauss(a1,n-1,n-2,n-1) ... 把各行均除以主对角元素之值 例6.4 用行阶梯求方程组的解 逐步消元得到行阶梯形式： 简化行阶梯形式例 再变为简化行阶梯形式 它对应的方程为： 也就是得出了方程的解。 6.4 MATLAB的行阶梯解 上节所介绍的方法比起一个一个元素地计算，效率是高多了。但仍然很麻烦，头脑中要记各个步骤中的很多细节，而实际上这些步骤仍然是刻板而有规律的，容易机械化。MATLAB已经把‘简化行阶梯形式(reduced row echelon form)’的计算过程集成为一个子程序rref.它的输入变元是线性方程组的增广矩阵，键入U0=rref([A,B]),输出的结果就是增广矩阵C的精简行阶梯形式。 上例的计算程序可以简化为： A=[1,4,7;8,5,2,;3,6,-2]; B=[1;3;5]; C=[A,B]; U0=rref(C) 更高阶的情况举例6.6 ag606