计算方法jacobi迭代法与高斯seidel迭代法实验.doc

计算方法jacobi迭代法与高斯seidel迭代法实验.doc

ID:58648684

大小:113.50 KB

页数:6页

时间:2020-10-16

计算方法jacobi迭代法与高斯seidel迭代法实验.doc_第1页
计算方法jacobi迭代法与高斯seidel迭代法实验.doc_第2页
计算方法jacobi迭代法与高斯seidel迭代法实验.doc_第3页
计算方法jacobi迭代法与高斯seidel迭代法实验.doc_第4页
计算方法jacobi迭代法与高斯seidel迭代法实验.doc_第5页
资源描述:

《计算方法jacobi迭代法与高斯seidel迭代法实验.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、实验名称:实验3jacobi迭代法与Gauss-seidel迭代法实验题目:给定线性方程组Ax=b如下实验目的:掌握用jacobi迭代法与Gauss-seidel迭代法求解线性方程组的基本步骤。基础理论:jacobi迭代法基本思路是方程组Ax=b等价于x=Bx+f,然后通过迭代算出方程组的解;Gauss-seidel迭代法主要是通过对A的分解,构造迭代公式,进而迭代算出方程组的解。实验环境:操作系统:WindowsXP;实验平台:matlab实验过程:方法一:jacobi迭代法程序:n=4;A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];b=[6,2

2、5,-11,15];x0=[0,0,0,0];x0=zeros(n,1);x=x0;epsilon=input('精度=');N=input('最大迭代次数N=');fprintf('%d:',0);fori=1:nfprintf('%f',x0(i));end%以下是迭代过程fork=1:N%这是第k步迭代,迭代前的向量在x0[]中,迭代后的在x[]中normal=0;fori=1:nx(i)=b(i);forj=1:nifj~=ix(i)=x(i)-A(i,j)*x0(j);endendx(i)=x(i)/A(i,i);temp=abs(x(i)-x0(i));%求范数与迭

3、代在同一个循环中iftemp>normalnormal=temp;%这里用的是无穷范数endend%第i步迭代结束fprintf('%d:',k);fori=1:nx0(i)=x(i);%为下一次迭代准备初值fprintf('%f',x(i));%输出迭代过程endifnormal

4、8];b=[6,25,-11,15];x=[0,0,0,0];epsilon=input('精度=');N=input('最大迭代次数N=');fprintf('%d:',0);fori=1:nfprintf('%f',x(i));end%以下是迭代过程fork=1:N%这是第k步迭代,迭代前的向量在x0[]中,迭代后的在x[]中normal=0;fori=1:nt=x(i);x(i)=b(i);forj=1:nifj~=ix(i)=x(i)-A(i,j)*x(j);endendx(i)=x(i)/A(i,i);temp=abs(x(i)-t);%求范数与迭代在同一个循环中if

5、temp>normalnormal=temp;%这里用的是无穷范数endend%第i步迭代结束fprintf('%d:',k);fori=1:nfprintf('%f',x(i));%输出迭代过程endifnormal

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

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

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