数值分析雅克比迭代法高斯—赛德尔迭代法董安葳

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

数值分析第八次实验报告姓名:董安葳学号:5123119题目:雅克比迭代法、高斯—赛德尔迭代法实验方法:本人通过学习雅克比迭代法和高斯—赛德尔迭代法,理解算法的本质,上机编写代码并调试,最终得到两个程序。实验过程:1、雅克比迭代:functionx=yakebi(A,b,p,wucha)%A是系数矩阵,b是常数向量,p是范数类型,wucha是限制的精度[m,n]=size(A);D=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);ifdet(D)==0error('ÎÞ½â')endB=(D^-1)*(L+U);f=(D^-1)*b;x=zeros(n,1);c=50;%设定迭代次数为50次forii=1:cx1=(B*x+f);shijiwucha=norm(x1-x,p);if(shijiwuchawucha)break;elsex=x1;endend2、高斯—赛德尔迭代functionx=gaosisaideer(A,b,p,wucha)%A是系数矩阵,b是常数向量,p是范数类型,wucha是限制的精度[m,n]=size(A);D=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);ifdet(D)==0error('ÎÞ½â')endB=((D-L)^-1)*U;f=((D-L)^-1)*b;x=zeros(n,1);c=50;%设定迭代次数为50次forii=1:cx1=((D-L)^-1)*(U*x+b);shijiwucha=norm(x1-x,p);if(shijiwuchawucha)break;elsex=x1;endend实验截图:

1 / 2
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功