雅克比迭代法……

雅克比迭代法……

ID:43470873

大小:149.01 KB

页数:12页

时间:2019-10-04

雅克比迭代法……_第1页
雅克比迭代法……_第2页
雅克比迭代法……_第3页
雅克比迭代法……_第4页
雅克比迭代法……_第5页
资源描述:

《雅克比迭代法……》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、实验课程名称:数值分析实验项目名称运用插值法实验成绩实验者江骏专业班级软件0803组别同组者实验日期年月日第一部分:实验分析与设计一、实验内容描述⑴研究用Jacobi迭代法与Gauss-Seidel迭代法解下列方程组Ax=b的收敛性,通过上机计算,验证分析是否正确,并观察右端项对迭代收敛是否有影响,比较两法的收敛速度⑵松弛因子对超松弛因子迭代法收敛速度的影响,要求对不同的阶数进行迭代⑶观察欧拉显式方法的收敛性⑷观察欧拉隐式方法的收敛性⑸写出实验报告二、实验基本原理与设计①Jacobi迭代法:∑nj=1aijxj=bi(i=1,2,…

2、,n)②Gauss-Seidel迭代法:xi(K+1)=[bi-∑i-1j=1aijxj(k+1)-∑nj=i+1aijxj(k)]/aij③超松弛因子迭代法:Dx(k+1)=Dx(k)+Ex(k+1)+Fx(k)+b-Dx(k)X(k+1)=x(k)+D-1R④三角分解法:利用三角矩阵将方程组化解为两个方程组,从而简化运算。⑤欧拉显式方法:在任意节点tn+1=t0+(n+1)h处,u(tn+1)的近似值由Euler公式给出:un+1=un+hf(tn,un)⑥欧拉隐式方法:改进的Euler公式n+1=yn+hf(xn,yn)yn+

3、1=yn+[f(xn,yn)+f(xn+1,n+1)]三、主要仪器设备及耗材WindowsXpVC++6.0第二部分:实验调试与结果分析一、调试过程①Jacobi迭代法#include#defineN10floatABS(float,float);intmain(void){inti,j,n;floata[N][N],b[N];floatx[N],y[N];floate,total;printf("Pleaseinputthedimension:");scanf("%d",&n);printf("Pleaseinp

4、uttheCoefficientMatrix:");for(i=0;i

5、i!=j)total+=a[i][j]*x[j];}y[i]=(b[i]-total)/a[i][i];}e=0.0f;for(j=0;j0.0001);for(i=0;i

6、elsetotal=x-y;returntotal;}结果:流程图:开始阶数n,系数矩阵,初始向量for(i=0;iintmain(){doublex[3]={0,0,0};doublea[3][3]={6,2,-1,1,4,-2,-3,1,4};doubley[3]={-3,2,4};doubled[3][3],g[3];intround=5,i,j;for(i=0;

7、i<3;++i){g[i]=y[i]/a[i][i];for(j=0;j<3;++j){d[i][j]=i==j?0:-a[i][j]/a[i][i];}}while(round--){for(i=0;i<3;++i){x[i]=g[i];for(j=0;j<3;++j){x[i]+=d[i][j]*x[j];}printf("%lf",x[i]);}printf("");}}结果:流程图:开始系数矩阵增广阵for(i=0;i<3;++i)d[i][j]=i==j?0:-a[i][j]/a[i][i]for(j=0;j<3;++

8、j)x[i]+=d[i][j]*x[j]结束③超松弛因子迭代法#include#includeusingnamespacestd;float*one_array_malloc(intn);float*

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

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

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