Gauss-Jacobi迭代法

Question

我正在尝试编写一个使用 Gauss-Jacobi 迭代法求解方程组 Ax=B 的程序。

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

int main(void) {
    double **a, *b, *x, *f, eps = 1.e-2, c;  
    int n = 3, m = 3, i, j, bool = 1, d = 3;
    /* printf("n=") ; scanf("%d", &n);
       printf("m=") ; scanf("%d", &n) */
   
   
    a =malloc(n * sizeof *a);
    for (i = 0; i < n; i++) 
        a[i] = (double*)malloc(m * sizeof(double));

    b = malloc(m * sizeof *b);
    x = malloc(m * sizeof *x) ;  
    f = malloc(m * sizeof *f) ;
    for (i = 0; i < n; i++) {
        for (j = 0; j < m; j++) { 
            printf("a[%d][%d]=", i, j); 
            scanf("%le", &a[i][j]); 
            if(fabs(a[i][i])<1.e-10) return 0 ; 
        }

        printf("\n") ;
    }
        
    printf("\n") ;
        
    for (i = 0; i < n; i++) {
        for (j = 0; j < m; j++) { 
            printf("a[%d][%d]=%le  ", i, j, a[i][j]); 
        }
         
        printf("\n") ;
    }
    
    for (j = 0; j < m; j++) { 
        printf("x[%d]=", j); 
        scanf("%le", &x[j]); 
    } //intial guess
    
    printf("\n") ;
    
    for (j = 0; j < m; j++) { 
        printf("b[%d]=", j); 
        scanf("%le", &b[j]); 
    }
    
    printf("\n")  ;

    while (1) {
        bool = 0;
        for (i = 0; i < n; i++) {
            c = 0.0;
            for (j = 0; j < m; j++) 
                if (j != i) 
                    c += a[i][j] * x[j];  
            f[i] = (b[i] - c) / a[i][i];
        }
       
        for (i = 0; i < m; i++)  
            if (fabs(f[i] - x[i]) > eps) 
                bool = 1;
       
        if (bool == 1) 
            for (i = 0; i < m; i++) 
                x[i] = f[i];
        else if (bool == 0) 
            break;
    }

    for (j = 0; j < m; j++) 
        printf("%le\n", f[j]);

    return 0;
}

停止循环的条件是所有x的前一近似值减去当前近似值小于epsilon。好像我按照算法做了所有事情，但程序不起作用。我哪里弄错了？

Answer 1

虽然不是最严格的条件，但保证雅可比和高斯-赛德尔方法收敛的通常条件是对角优势，

abs(a[i][i]) > sum( abs(a[i][j]), j=0...n-1, j!=i)

作为迭代前对运行的检查，此测试也很容易实施。

所有这些不等式的相对差距越大，该方法的收敛速度就越快。

Gauss-Jacobi迭代法

Gauss-Jacobi iteration method

c

numerical-methods