经常遇到求解包含个未知数、由个方程构成的线性方程组.docVIP

在科技、工程、医学和经济等各个领域中，经常遇到求解包含个未知数、由个方程构成的线性方程组 （3.1） 的问题.线性方程组（3.1）的解法一般有直接法和迭代法.迭代法是求解大型稀疏矩阵方程组，尤其是由微分方程离散后得到的大型方程组的重要方法.在这章中我们要学习线性方程组的直接法，特别是适合用数学软件在计算机上求解的方法. 3.1 方程组的逆矩阵解法及其MATLAB程序 3.1.1 逆矩阵、行列式及其MATLAB命令 当矩阵为方阵时，的行列式通常表示为或.当时，可逆，且可以用表 3-1中的MATLAB命令求的逆矩阵和行列式. 表 3-1 命 令 能 NA=inv(A) 输入矩阵，运行后输出的逆矩阵； HA=det(A) 输入矩阵，运行后输出的行列式的值. 3.1.2 方程组的逆矩阵解法及其MATLAB程序 在MATLAB中引进了矩阵除法的概念，它包括左除和右除.设两矩阵为阶方阵，且，皆可逆，则矩阵方程，,的解都可以用表 3-2中的命令求出. 表 3-2 命 令 能 X=A\C X= inv(A)*C 输入矩阵和C，运行后输出矩阵方程的解 Y=D/B Y=D*inv(B) 输入矩阵和，运行后输出矩阵方程的解 Z=A\ F/B Z= inv(A)*F* inv(B) 输入矩阵，和，运行后输出矩阵方程的解 3.1.3 线性方程组有解的判定条件及其MATLAB程序 解非齐次线性方程组，将其增广矩阵化成行阶梯形矩阵，便可判断其是否有解．若有解，化成行最简形矩阵，便可写出其通解. 判定线性方程组是否有解的MATLAB程序 输入的量：系数矩阵和常数向量； 输出的量：RA和RB分别是系数矩阵和增广矩阵的秩, n是方程组中未知量的个数，另外还输出有关方程组解的信息. 说明：(1)求矩阵的秩用MATLAB命令：RA=rank(A). (2)如果有唯一解，则用表 3-2方法求解; 如果有无穷多解，则将其增广矩阵化成行最简形矩阵，便可写出其通解. 根据线性方程组有解的判定条件，现提供判定线性方程组是否有解的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==n disp(请注意：因为RA=RB=n，所以此方程组有唯一解.) else disp(请注意：因为RA=RBn，所以此方程组有无穷多解.) end end 例3.1.4 判断下列线性方程组解的情况.如果有唯一解，则用表 3-2方法求解. (1) (2) (3) (4) 解 在MATLAB工作窗口输入程序 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 = 4,RB =4,n =4 在MATLAB工作窗口输入X=A\b, 运行后输出结果为 X =(0 0 0 0)’. (2) 在MATLAB工作窗口输入程序 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 (3)　在MATLAB工作窗口输入程序 A=[4 2 -1;3 -1 2;11 3 0]; b=[2;10;8]; [RA,RB,n]=jiepb(A,B) 运行后输出结果 请注意：因为RA~=RB，所以此方程组无解. RA =2,RB =3,n =3 (4)在MATLAB工作窗口输入程序 A=[2 1 -1 1;4 2 -2 1;2 1 -1 -1]; b=[1; 2; 1]; [RA,RB,n]=jiepb(A,b) 运行后输出结果 请注意：因为RA=RBn，所以此方程组有无穷多解. RA =2,RB =2,n =3 3.2 三角形方程组的解法及其MATLAB程序 用初等行变换解线性方程组是解线性方程组的一般方法.从例3.1.6可见，当方程的个数和未知量的个数很大时，用初等行变换解线性方程组的工作量相当的大，如果在计算过程中某一步发生错误，那么就会得出错误的结果.人们非常想用计算机的程序解线性方程组，而高斯（Gauss）消元法以及各种变形的高斯消元法使得用初等行变换解线性方程组