超松驰迭代法SOR_解线性方程组的迭代法

超松驰迭代法SOR_解线性方程组的迭代法

标签:计算方法实验

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

#define maxn 3

int main()
{
    double a[maxn][maxn + 1], x[maxn] = {0};
    double eps = 1e-9, w = 1.05;
    int n, k, kmax = 100;

    freopen("gauss.txt", "r", stdin);
    scanf("%d", &n);
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n + 1; j++)
            scanf("%lf", &a[i][j]);
    //for(int i = 0; i < n; i++)
    //{
        //for(int j = 0; j < n + 1; j++)  printf("%-15f", a[i][j]);
        //printf("\n");
    //}

    for(k = 0; k < kmax; k++)
    {
         double norm = 0;

         for(int i = 0; i < n; i++)
         {
             double x0 = x[i];
             double sum = 0;
             for(int j = 0; j < n; j++)  if(j != i)  sum += a[i][j] * x[j];  ///
             x[i] = (1 - w) * x0 + (a[i][n] - sum) / a[i][i];  ///w松弛因子
             if(fabs(x[i] - x0) > norm) norm = fabs(x[i] - x0);  //norm计算范数
         }
         printf("\nk = %2d  x =  ", k + 1);
         for(int i = 0; i < n; i++)  printf("%-15f", x[i]);
         if(norm < eps)  break;
    }
    if(k < kmax)
    {
        printf("\n\nk = %d\n", k + 1);
        for(int i = 0; i < n; i++)  printf("x%d = %-15f\n", i + 1, x[i]);
    }
    else  printf("\n\nfailed\n");

    return 0;
}

数据文件
input
实验结果
output

  • 4
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值