雅可比迭代法和高斯-赛德尔迭代法(共4页).doc

《雅可比迭代法和高斯-赛德尔迭代法(共4页).doc》由会员分享，可在线阅读，更多相关《雅可比迭代法和高斯-赛德尔迭代法(共4页).doc（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上题目：使用雅可比迭代法和高斯-赛德尔迭代法 求解线性方程组，C语言实现本人编写了一个特定的解决方案程序，并且利用其对课本156页例1进行验证，主要代码和运行结果如下：1. 代码：# include# include# define N 3/*使用雅可比迭代法和高斯-赛德尔迭代法 求解线性方程组 */ main()float NF2(float *x,float *y);float ANN,bN,sum=0;float xN,yN=0,x0N;int i,j,n=0;/输入系数矩阵 printf(please input the Coefficient Matrix:)

2、;for(i=0;iN;i+)for(j=0;jN;j+)scanf(%f,&Aij);/输入常数矩阵printf(please input the Constant Matrix:);for(i=0;iN;i+) scanf(%f,&bi); /输入解的初值 printf (Please input the initial vector: ); for(i=0;iN;i+) scanf(%f,&x0i); /输出系数矩阵printf(show your the Coefficient Matrix：n); for(i=0;iN;i+)for(j=0;jN;j+)printf(%3.1f ,A

3、ij);printf(n); /输出成数矩阵printf(show your the Constant Matrix：n);for(i=0;iN;i+)printf(%3.1fn,bi); /输出解的迭代初值printf(show your the Initial Vector：n);for(i=0;iN;i+) printf(%3.1fn,x0i);/*利用雅可比迭代法求解线性方程组 */for(i=0;iN;i+)xi=x0i;for(n=0;n+)/计算下一个值 for(i=0;iN;i+)sum=0;for(j=0;j0.01)for(i=0;iN;i+)xi=yi; elsebrea

4、k; printf(经过%d次雅可比迭代解出方程组的解：n,n+1);for(i=0;iN;i+)printf(%f ,yi);/*利用高斯-赛德尔迭代法求解线性方程组 */for(i=0;iN;i+)xi=x0i;yi=0;for(n=0;n+)/计算下一个值 for(i=0;iN;i+)sum=0;for(j=0;ji;j+)sum=sum+Aij*yj;for(j=i+1;j0.01)for(i=0;iN;i+)xi=yi; elsebreak; printf(n经过%d次高斯-赛德尔迭代解出方程组的解：n,n+1);for(i=0;iN;i+)printf(%f ,yi);/求两个向量差的二范数函数 float NF2(float *x,float *y)int i;float z,sum1=0;for(i=0;iN;i+)sum1=sum1+pow(yi-xi,2);z=sqrt(sum1);return z;2. 运行结果专心-专注-专业