线性方程组的直接接法.pptVIP

例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 高斯消去法的MATLAB程序：Gauss.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 index=0;return; end for k=n:-1:1 for j=k+1:n b(k)=b(k)-A(k,j)*x(j); end x(k)=b(k)/A(k,k); end 线性方程组的直接接法 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中消元中 ，使得算法数值不稳定。 例3 解 3位浮点数 3.3 高斯主元素消去法 线性方程组的直接接法 （2） 用小主元作除数，计算结果不可靠 精确解 【注】（1）A 为非奇异，可以用高斯消元法 （3）应避免采用小主元 基本思想：在每轮消元之前，选列主元素（绝对值 最大的元素） 3.3.1 列主元素消去法 线性方程组的直接接法 步骤： 第k步计算如下： （1）选列主元素，即确定使 （2）如果 或者A接近奇异，停止 （3）如果 （4）消元计算： （2）如果 或者A接近奇异，停止 （2）如果 线性方程组的直接接法 第3章 线性方程组的直接接法 理学院应用数学系 立体化教学资源系列——数值分析 线性方程组的直接接法 问题 一种大型输电网络可简化为电路 负载电阻 线路内阻 电源电压V 负载电流 ?列出各负载上电流的方程 ?讨论情况 ?设 n=10,求 及总电流 2.2 输电网络 n=10,求 及总电流 n=10,求 ?讨论情况 分析 根据电路中电流、电压关系,列出 消 和 求电流方程 ? 线性方程组的直接接法 求电流方程 ? 其中 MATLAB计算电流程序 r=1;R=6;v=18;n=10; b1=sparse(1,1,v,n,1); b=full(b1); a1=triu(r*ones(n,n)); a2=diag(R*ones(1,n)); a3=-tril(R*ones(n,n),-1)+ tril(R*ones(n,n),-2); a=a1+a2+a3; I=a\b;I0=sum(I) 线性方程组的直接接法 许多科学和工程技术问题，都归结为求解线性方程组.例如，电学中的网络问题，实验数据的曲线、曲面拟合问题，用差分法解微分方程边值问题，用有限元法解结构力学问题，解非线性方程组等大型稀疏矩阵（即矩阵阶数高且零元素较多）. 线性方程组的直接接法 （2）回代求解 ： （1）消元运算： 线性方程组的直接接法 3.2引言与预备知识 以行初等变换为基础，整理而成的计算机算法. 例1 用消去法解方程组 高斯消去(GaussianElimination) : 线性方程