线性代数实践教师班第二讲.pptVIP

用图解法解例6.1(d)第30页，共56页，星期日，2025年，2月5日用解析法解例6.1.1用代入法或消去法等中学方法求解；用矩阵求解的方法用MATLAB符号工具箱求解[x1a,x2a]=solve(s1,s2)其结果当然是相同的： (a)解为x1a=3,x2a=2 (b)解为x1b=[],x2b=[], 表示无解 (c)解为x1c=2*x2-1,x2c=x2， 表示有无数解 (d)解为x1d=[],x2d=[], 表示无解第31页，共56页，星期日，2025年，2月5日求解三元线性方程组的例 →再由两式联立解得： x=1.4286y=-1.8571z=-5.4286用图解方法ezplplot3画此三个方程的平面，用的就是例5.3的程序，得到：第32页，共56页，星期日，2025年，2月5日按例6.2画出的三个平面第33页，共56页，星期日，2025年，2月5日三元一次方程的几何意义这三个方程的几何意义为空间的三个平面，两两消去z，意味着求两者的交线，得到的两个二元一次方程对应于求得的两根交线，最后由交线得到交点，该交点就是方程的解（见图6.3）。当两个平面平行，没有交线，或者两根交线平行，没有交点时，方程组就不相容，因而无解；同样也可能有无穷个解的情况，读者可自行思考想象。第34页，共56页，星期日，2025年，2月5日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,:)第35页，共56页，星期日，2025年，2月5日消元法的要点利用上述初等变换，在保持矩阵等价性和其它各行不变的前提下，把矩阵第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）其中的冒号’：’意味着要对全行的所有元素都做一次乘法和一次加法，相当繁琐。第36页，共56页，星期日，2025年，2月5日消元子程序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第37页，共56页，星期日，2025年，2月5日6.3行阶梯形式具有以下三个特点的矩阵称为行阶梯形式：1。所有非零行都处在全零行的上方；2。行首非零(枢轴)元素的列号比它的上方所有各行的枢轴元素的列号都要大；3。各枢轴元素所处的列中，在枢轴元素下方的所有元素均为零。如果此矩阵还满足以下两个额外特点，4。各行枢轴元素是其列中的唯一非零元素；5。这些枢轴元素都等于1。则称之为简化行阶梯形式。第38页，共56页，星期日，2025年，2月5日例6.3行阶梯方程组的解法用逐次回代法，可以由下而上依次求得 x4=2,x3=2+x4/2=3, x2=-4-8x3=-28, x1=3x2-5x3+2x4=-39。第39页，共56页，星期日，2025年，2月5日消元法生成行阶梯形式步骤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第40页，共56页，星期日，2025年，2月5日第1页，共56页，星期日，2025年，2月5日5.1实验在线性代数中的重要性利用软件