第三章解线性方程组的直接方法的matlab程序.docVIP

J第三 J 第三章 解线性方程组的直接方法 高等教育出版社 教育电子音像出版社 ：任玉杰 第三章 解线性方程组的直接方法的MATLAB程序 在科技、工程、医学和经济等各个领域中，经常遇到求解包含〃个未知数、由加个 方程构成的线性方程组 ⑷內+%2兀2+_ +知心=5 色內+。22兀2+…+。2朋严為 小八 、细內+知2兀2十…+ ^加兀=5 的问题?线性方程组(3.1)的解法一般有直接法和迭代法.迭代法是求解大型稀疏矩 阵方程组，尤其是由微分方程离散后得到的大型方程组的重要方法?在这章中我们要学习 线性方程组的直接法，特别是适合用数学软件在计算机上求解的方法. 3. 1方程组的逆矩阵解法及其MATLAB程序 3.1.1逆矩阵、行列式及其MATLAB命令 当矩阵A为方阵时，A的行列式通常表示为det(A)或|A| ?当det(A)^O时，A可 逆，且可以用表3-1中的MATLAB命令求A的逆矩阵和行列式. 表3-1 命令 功 能 NA=inv(A) 输入矩阵4,运行后输出A的逆矩阵A1； HA=det(A) 输入矩阵A ,运行后输出A的行列式det(A)的值. 3.1.2方程组的逆矩阵解法及其MATLAB程序 在MATLAB中引进了矩阵除法的概念,它包插左除和右除.设两矩阵A, 〃为兄阶方阵, 且A, 〃皆可逆，则矩阵方程AX=C, YB = D, AZB = F的解都可以用表3-2中 的命令求出. 表3-2 命令 相同功能的命令 功 能 X=A\C X= inv(A)*C 输入矩阵人和C,运行后输出矩阵方程 AX =C 的解 X = A C Y=D/B Y=D*i nv(B) 输入矩阵B和D,运行后输出矩阵方程 YB = D 的解 Y = DB1 Z=A\ F/B Z= inv(A)*F* inv(B) 输入矩阵A, B和F,运行后输出矩阵 方程AZB二F的解Z = A ]FB } 3.1.3线性方程组有解的判定条件及其MATLAB程序 解非齐次线性方程组，将其增广矩阵化成行阶梯形矩阵，便可判断其是否有解.若 有解，化成行最简形矩阵，便可写出其通解. 判定线性方程组A ”呦X=b是否有解的MATLAB程序 根据线性方程组有解的判泄条件，现提供判泄线性方程组是否有解的MATLAB 根据线性方程组有解的判泄条件，现提供判泄线性方程组是否有解的MATLAB程序 function [RA,RB,n]=jiepb(A,b) B=[A b];n=length(b); RA=rank(A); RB=rank(B);zhica=RB-RA; if zhica0, disp(请注意：因为RA~=RB,所以此方程组无解.) return end if RA==RB if RA= dispC请注意：因为RA=RB=n,所以此方程组有唯一解.) else disp(请注意：因为RA=RBn,所以此方程组有无穷多解.) end end 例3.1.4判断下列线性方程组解的情况?如果有唯一解，贝9用表3-2方法求解. 2%| + 3x 2%| + 3x2 一兀3 + 5勺=0, (]) 3%| + 兀2 + 2兀3 一 7兀4 = 0, 4兀[+ x2 - 3巾 + 6巾-0, x{ 一 2x2 + 4无3 一 7无4 - 0; 4“ + 2x2 一兀3 =2, ⑶ 3兀| -兀2 + 2*3 = 10, 11 %)+ 3 兀2 = 8; 解 在MATLAB I作窗口输入程序 3兀]+ 4x2 一 5兀3 + 7兀4 = 0, 2%1 一 3兀2 + 3小 一 2x4 = 0, 4兀]+ 1 1 兀2 - 13心 + 16*4 = 0, 7兀[一 2x2 + 心 + 3*4 = 0； 2尢 + y — z + vv = 1, 4x + 2y - 2z + w = 2, 2x+y-z-w = l. ? A=[2 3 -1 5;3 1 2 -7;4 1 ?3 6;1 -2 4 ?7]; b=[ 0; 0; 0; 0]; [RA,RB,n]=jiepb(A,b) 运行后输出结果为 请注意：因为RA=RB=n,所以此方程组有唯一解. RA = 4,RB =4,n =4 在MATLAB工作窗口输入 ?X=A\b, 运行后输出结果为X=(0 0 0 0)? 在MATLAB I作窗口输入程序 ? A=[3 4 -5 7;2 -3 3 -2;4 11 -13 16;7 -2 1 3];b=[ 0; 0; 0; 0]; [RA,RB,n]=jiepb(A,b) 运行后输出结果 请注意：因为RA=RBn,所以此方程组有无穷多解. RA =2,RB =2,n =4 在MATLAB工作窗口输入程序 ? A=[4 2 -1 ;3 -1 2;11 3 0];b=[2;10