西北大学数值分析课件第3章线性方程组的直接接法.pptxVIP

;;;求电流方程 ? ; 许多科学和工程技术问题，都归结为求解线性方程组.例如，电学中的网络问题，实验数据的曲线、曲面拟合问题，用差分法解微分方程边值问题，用有限元法解结构力学问题，解非线性方程组等大型稀疏矩阵（即矩阵阶数高且零元素较多）. ;;3.2引言与预备知识;;首先讨论上三角方程组的回代求解：;;;;;;;;;;;;【定理2】 ;function [x,det,index]=Gauss(A,b) %求线性方程组的高斯消去法，其中， %A为方程组的系数矩阵； %b为方程组的右端项； %x为方程组的解； %det为系数矩阵A的行列式的值； %index为指标变量，index=0表示计算失败， index=1表示计算成功. [n,m]=size(A);nb=length(b); %当方程组行与列的维数不相等时， 停止计算，并输出出错信息. if n~=m;%消元过程 for i=k+1:n m=A(i,k)/A(k,k); for j=k+1:n A(i,j)=A(i,j)-m*A(k,j); end b(i)=b(i)-m*b(k); end det=det*A(k,k); end det=det*A(n,n); %回代过程 if abs(A(n,n))1e-10 ;error(The rows and columns of matrix A must be equal!); return; end %当方程组与右端项的维数不匹配时， 停止计算，并输出出错信息. if m~=nb error(The columns of A must be equal the length of b!); return; end %开始计算，先赋初值 index=1;det=1;x=zeros(n,1); for k=1:n-1 ;调用函数Gauss.m解例2.输入 A=[1 4 7;2 5 8;3 6 11]; b=[1 1 1]; [x,det,index]=Gauss(A,b) 得到x = -0.3333 det = -6 index = 1 计算成功. 0.3333 0;Gauss中消元中 ，使得算法数值不稳定。; （2） 用小主元作除数，计算结果不可靠;步骤：