实验五线性方程组地迭代法实验.docx
《实验五线性方程组地迭代法实验.docx》由会员分享,可在线阅读,更多相关《实验五线性方程组地迭代法实验.docx(8页珍藏版)》请在三一文库上搜索。
1、实用标准文档柬霜工大母计算方法实验报告学 院: 信息学院专 业: 计算机科学与技术指导教师:郭卫斌班级学号:10101438 计102姓 名: 闻翰计算机科学与工程系实验五线性方程组的迭代法实验实验目的(1)深入理解线性方程组的迭代法的设计思想,学会利用系数矩阵的性质以保证迭 代过程的收敛性,以及解决某些实际的线性方程组求解问题。(2)熟悉Matlab编程环境,利用Matlab解决具体的方程求根问题。实验要求建立Jacobi迭代公式、Gauss-Seidel迭代公式和超松弛迭代公式,用 Matlab 软件实现线性方程组求解的 Jacobi迭代法、Gauss-Seidel迭代法和超松弛迭 代法,
2、并用实例在计算机上计算。实验内容1 .实验题目迭代求解下列线性方程组,取(1 )分别利用Jacobi迭代和 Gauss-Seidelx 0,0,0,0,0,0 T ,要求精度(2)分别取1、度为 10 54101001.05、1.114101001400110 5100410、1.25010141001014X1X2X3X4X5X6050626用超松弛法求解上面的方程组,要求精2 .设计思想1 .Jacobi 迭代:Jacobi迭代的设计思想是将所给线性方程组逐步对角化,将一般形式的线性方程组 的求解归结为对角方程组求解过程的重复。2 .Gauss-Seidel 迭代:Gauss-Seidel
3、迭代的设计思想是将一般形式的线性方程组的求解过程归结为下三角 方程组求解过程的重复。3 .超松弛迭代:基于Gauss-Seidel迭代,对i=1 , 2,反复执行计算迭代公式,即为超松弛迭代。4 .对应程序1.Jacobi 迭代:function x,k=Jacobimethod(A,b,x0,N,emg)%靠线性方程组白左端矩阵,b是右端向量,x0是迭代初始值% N表示迭代次数上限,emg表示控制精度,k表示迭代次数,x是解 n=length(A);x1=zeros(n,1);x2=zeros(n,1);x1=x0;k=0;r=max(abs(b-A*x1);while remgfor i=
4、1:nsum=0;for j=1:nif i=jsum=sum+A(i,j)*x1(j);endendx2(i)=(b(i)-sum)/A(i,i);endr=max(abs(x2-x1);x1=x2;k=k+1;if kNdisp(迭代失败,返回);return;endendx=x1;2 .Gauss-Seidel 迭代:function x,k=Gaussmethod(A,b,x0,N,emg)%靠线性方程组白左端矩阵,b是右端向量,x0是迭代初始值% N表示迭代次数上限,emg表示控制精度,k表示迭代次数,x是解 n=length(A);x1=zeros(n,1);x2=zeros(n,
5、1);x1=x0;r=max(abs(b-A*x1);k=0;while remgfor i=1:nsum=0;for j=1:nif jisum=sum+A(i,j)*x1(j);elseif jNdisp(迭代失败,返回);return;endendx=x1;3 .超松弛(SOR)代:function x,k=SORmethod(A,b,x0,N,emg,w)%幅线性方程组白左端矩阵,b是右端向量,x0是迭代初始值% N表示迭代次数上限,emgg示控制精度,k表示迭代次数,x是解 %WS小松弛因子n=length(A);x1=zeros(n,1);x2=zeros(n,1);x1=x0;r
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 线性方程组 迭代法
链接地址:https://www.31doc.com/p-14096412.html