高斯—赛德尔迭代法解线性方程组(C语言)

原创 2012年03月30日 18:45:09

/*-------------高斯—赛德尔迭代法解线性方程组---------

参考教材:《数值分析》李乃成,梅立泉,科学出版社
    《计算方法教程》第二版 凌永祥,陈明逵*/

#include<stdio.h>
#include<math.h>

int main(void)
{
 double A[5][5] = {{28,-3,0,0,0},
      {-3,38,-10,0,-5},
     {0,-10,25,-15,0},
     {0,0,-15,45,0},
     {0,-5,0,0,30}};
 double b[5] = {10,0,0,0,0};
 double x[5] = {0,0,0,0,0}; //第k+1次迭代的结果
 double xx[5] = {0,0,0,0,0}; //第k次迭代的结果
 int size = 5;
 int Max = 100;  //最大迭代次数
 double residual = 0.0;  //
 double sum = 0.0; 
 double dis = 0.0;
 double dif = 1.0;  //相邻迭代的结果差
 double eps = 1.0e-3; //迭代精度

 for(int k=1;(k<Max)&&(dif>eps);k++)
 {
  dif = 0.0;

  for(int i=0;i<size;i++)
  {
   
   for(int j=0;j<size;j++)
   {
    //与Jacobi迭代法唯一的区别就是sum的求法
    if(j<i)
    {
     sum +=A[i][j]*x[j];//本次的x[j]
    }
    if(j>i)
    {
     sum +=A[i][j]*xx[j];//上一次的x[j]
    }

   }
   x[i] = (b[i]-sum)/A[i][i];
   sum=0.0;

  } 
  residual=0.0;
  //计算相邻迭代的结果差
  for(int m=0;m<size;m++)
  {
   dis=fabs(x[m]-xx[m]);
   if(dis>residual)
    residual=dis;
  }
  dif=residual;

  //打印第k次的结果
  printf("\n第%d次迭代的结果:\n",k);
  for(i=0;i<size;i++)
  {
   printf("%12.8f ",x[i]);
   xx[i]=x[i];
  }
  printf("\n与上次计算结果的距离(无穷范数):%12.8f \n",dif);
 }
 printf("\n迭代计算的结果为:\n");
 for(k=0;k<size;k++)
 {
  printf("%12.8f ",xx[k]);
 }
 printf("\n");
 return 0;
}

高斯—赛德尔(Gauss - Seidel)迭代法解线性方程组(Matlab程序)

%---高斯—赛德尔迭代法----- %---Gauss - Seidel iteration method clear;clc; % A=[10,-1,-2;-1,10,-2;-1,-1,5]...
  • zhangchao3322218
  • zhangchao3322218
  • 2012年03月30日 18:37
  • 15903

基于matlab的Guass-Seidel(高斯--赛德尔) 迭代法求解线性方程组

Guass-Seidel(高斯–赛德尔) 迭代法(简称 G−S 迭代)是对 Jacobi 迭代的一种改进. 了解G-S迭代法之前先了解什么是Jacobi迭代?链接如下: http://blog.c...
  • zengxyuyu
  • zengxyuyu
  • 2016年11月06日 20:49
  • 5491

分别用雅可比(Jacobi)迭代法和高斯—塞德尔(Gauss—Seidel)迭代法求解线性方程组

算法介绍(迭代法介绍): 代码C语言实现; # include # include # define N 6 /* *使用雅可比迭代法和高斯-赛德尔迭代法 求解线性方程组 ...
  • yyywww666
  • yyywww666
  • 2015年01月17日 14:26
  • 4418

Gauss–Seidel方法和Jacobi 方法

Gauss–Seidelmethod 对应于形如Ax = b的方程(A为对称正定矩阵或者Diagonally dominant),可求解如下:         Jacobi method 另一种方法是...
  • wangxiaojun911
  • wangxiaojun911
  • 2011年10月20日 12:57
  • 8397

计算方法上机程序——高斯赛德尔迭代法

  • 2008年12月05日 10:15
  • 1.14MB
  • 下载

高斯赛德尔迭代法的C++程序实现

  • 2008年12月23日 12:31
  • 2KB
  • 下载

分别用雅可比迭代法与赛德尔迭代法求解线性方程组Ax=b

  • 2013年03月11日 16:46
  • 607KB
  • 下载

高斯—赛德尔迭代法解线性方程组(C语言)

/*-------------高斯—赛德尔迭代法解线性方程组--------- 参考教材:《数值分析》李乃成,梅立泉,科学出版社     《计算方法教程》第二版 凌永祥,陈明逵*/ #inclu...
  • zhangchao3322218
  • zhangchao3322218
  • 2012年03月30日 18:45
  • 15211

高斯赛德尔迭代求解矩阵的解

  • 2016年01月20日 16:43
  • 121KB
  • 下载

java实现高斯赛德尔算法解线性方程组

java实现高斯赛德尔算法解线性方程组
  • jacke121
  • jacke121
  • 2017年02月18日 11:35
  • 653
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:高斯—赛德尔迭代法解线性方程组(C语言)
举报原因:
原因补充:

(最多只允许输入30个字)