第4.1讲 线性方程组迭代解法与MATLAB.ppt

调用函数Jacobi_matrix.m解例1 输入 A=[10 3 1;2 -10 3;1 3 10];b=[14 -5 14]; ep=0.005; [x,k,index]=Jacobi_matrix(A,b,ep) 雅可比迭代矩阵描述的MATLAB程序 Jacobi_matrix.m 例1 解方程组 高斯—塞德尔迭代法 第i个方程: 等价方程组: n阶方程组的高斯—塞德尔迭代法 任取初始解向量 计算得迭代序列 高斯-赛德尔(G-S)迭代公式 对于 终止条件: n阶方程组的高斯—塞德尔迭代法 例2 解下面方程组（与例1相同，精确解 ） 解 1） 等价方程组 高斯—塞德尔迭代法 2） 高斯-赛德尔迭代公式 3） 取初始向量 迭代得表 【注】高斯-赛德尔迭代法比雅可比迭代法收敛快。 高斯—塞德尔迭代法 矩阵表示 高斯—塞德尔迭代法的矩阵描述 G-S迭代公式的矩阵形式 【注】 （1）雅可比迭代法和高斯-赛德尔迭代法的分量形式供计算编程使用，矩阵形式供研究迭代序列是否收敛等理论分析使用。 （2）雅可比迭代适合并行计算；不足的是需要存放两个向量空间。 （3）高斯-赛德尔迭代法只需一个向量存储空间。 高斯—塞德尔迭代法的矩阵描述 （4）在某些情况下，G-S迭代法加速收敛，但它是一种典型的串行算法。 高斯-赛德尔迭代的MATLAB程序 Gauss_Seidel.m function [x,k,index]=Gauss_Seidel(A,b,ep,it_max) %解线性方程组的G-S迭代法，其中， %A为方程组的系数矩阵； %b为方程组的右端项； %ep为精度要求，缺省值1e-5； %it_max为最大迭代次数，缺省值100； %x为方程组的解； %k为迭代次数； %index为指标变量，index=0表示迭代失败，%index=1表示收敛到指定要求. [n,m]=size(A);nb=length(b); %当方程组行与列的维数不相等时，停止计算，并输出出错信息. if n~=m 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 if nargin4 it_max=100;end if nargin3 ep=1e-5;end k=0;x=zeros(n,1);y=zeros(n,1);index=1 高斯-赛德尔迭代的MATLAB程序 Gauss_Seidel.m while 1 y=x; for i=1:n z=b(i); for j=1:n if j~=i z=z-A(i,j)*x(j); end end if abs(A(i,i))1e-10|k== it_max index=0;return; end z=z/A(i,i);x(i)=z; end if norm(y-x,inf)ep break; end k=k+1; end 高斯-赛德尔迭代的MATLAB程序 Gauss_Seidel.m 求解实例 调用函数Gauss_Seidel.m 解例1. x = 0.9998 k = 4 index = 1 迭代成功，收敛 0.9998 1.0001 得到 输入 A=[10 3 1;2 -10 3;1 3 10]; b=[14 -5 14];ep=0.005; [x,k,index]=Gauss_Seidel(A,b,ep) 例1 解方程组 例3 分别用雅可比和高斯-赛德尔迭代法解方程组，均取相同初值 1） Jacobi4次达到精度 G-S发散。 2） Jacobi发散， G-S发散. 3） Jacobi 89次达到精度 G-S 8次达到同样的精度。 4） Jacobi发散，而G-S 10次达到精度 雅可比迭代法和高斯-赛德尔迭代法可能同时发散； 也可能同时收敛，但一个快另一个慢； 可能一个收敛而另一个发散。 【注】 结论 ! 谢谢大家 汕头大学新校园 规划效果图 数值计算方法与MATLAB 线性方程组迭代解法与MATLAB 线性方程组求解的工程由来 问题 一种大型输电网络可简化为电路 负载电阻 线路内阻 电源电压V 负载电流 要求：列出各负载上电流的方程 工程问题：输电网络 ?讨论情况 ?设： 及总电流 n=10,求 ?讨论情况 工程问题