MATLAB实现高斯赛德尔迭代.doc

MATLAB实现高斯赛德尔迭代.doc

ID:48359488

大小:69.50 KB

页数:2页

时间:2019-11-26

MATLAB实现高斯赛德尔迭代.doc_第1页
MATLAB实现高斯赛德尔迭代.doc_第2页
资源描述:

《MATLAB实现高斯赛德尔迭代.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、%高斯赛德尔迭代clearall;closeall;clc;ticformatlongedisp('请输入参数');K=input('维数K=');ticA=100rand(K);%A元素是0-100fori=1:KA(i,i)=sum(abs(A(i,:)))+25rand(1);%对角占优的量为0~25endb=zeros(K,1);fori=1:K;x=0;forr=1:K;x=x+A(i,r);endb(i,1)=x;end%产生b矩阵,b中的元素为A中对应行的和,目的是使方程解全为1jd=input('控制精度jd=');

2、I=eye(K);%单位阵L=-tril(A,-1);%下三角U=-triu(A,1);%上三角D=diag(diag(A));%对角矩阵B=I-(D-L)A;%迭代矩阵%或B=(D-L)U;f=(D-L)b;%fx0=zeros(K,1);%初始迭代矩阵y=Bx0+f;%迭代公式s=1;%迭代次数whilenorm(y-x0)>=jd&&s<3000;%迭代条件x0=y;y=Bx0+f;%迭代公式s=s+1;%迭代次数+1endystoct=1:K;yy=abs(y'-1)/1;plot(t,yy);title('绝对误差图

3、')

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。