Gauss-Seidel迭代法.doc

Gauss-Seidel迭代法廖绮绮04120049计算机04级系2班实验目的理解和掌握Gauss-Seidel迭代法的方法来解线性方程组,了解Matlab,加深和巩固知识。。由于我是采用用文件读入数据的方法,、=b。。执行(RUN),用户验证一下是否有错。B=10-1-2-110-2-1-15b=728342 然后用户一直按Enter键即可一步一步的得出结果。实验内容用Gauss-Seidel迭代法来演示求解线性方程组的过程。用Malab软件实现。主要是因为Matlab强大的计算力。方法源码:%Createdbycolent12/7%Copyright(C)2007-.\nAuthors:廖绮绮(ID:04120049).%Thanks鞠颖老师的先前程序%Gauss-Seidel迭代法function[]=gauss_Seidel() %定义方程%datainput %B=[10-1-2;-110-2;-1-15]%b=[728342]'B=importdata('') %用文件读入数据b=importdata('')pause; %等待用户响应dim=max(size(B));%矩阵A的行、列数的最大值x=zeros(dim,1);%生成dim维全0矩阵maxCount=6;%最大容许迭代次数forj=2:maxCount x(1)=(b(1)-B(1,2:dim)*x(2:dim))/B(1,1);fori=2:dim-1x(i)=(b(i)-B(i,1:i-1)*x(1:i-1)-B(i,i+1:dim)*x(i+1:dim))/B(i,i);endx(dim)=(b(dim)-B(dim,1:dim-1)*x(1:dim-1))/B(dim,dim);count=j-1x=xpause;end 用数组B存放AX=b中的A 用数组b存放AX=b中的b 用数组x存放迭代解 用count显示当前迭代次数 用dim存放矩阵的行列数的最大数理论分析及测试结果设Bx=b其中为非奇矩阵且(i=1,2,…n),本算法用高斯-塞德尔迭代法解Bx=b,数组x(n)开始存放x(0)=(0,0,…..),后存放,maxCount为最大