数值分析3.1 J迭代法&&GS迭代法

最新推荐文章于 2024-04-11 12:20:36 发布

原创

最新推荐文章于 2024-04-11 12:20:36 发布 · 7k 阅读

8 ·

CC 4.0 BY-SA版权

文章标签：

#数值分析

本文探讨了数值分析中的两种迭代方法——J迭代法和GS迭代法。J迭代法通过三维数组来存储每次迭代的结果，而GS迭代法在实现上更为简便。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

J迭代算法：

解向量我直接开了三维的数组，第三个分量表示第k次迭代的结果

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <math.h>
using namespace std;
typedef long long LL;
const int MAXN = 501;
double a[MAXN][MAXN];///系数矩阵
double b[MAXN],x[MAXN][MAXN];///右端项及解
int n,N;///n元方程组，最大迭代N次
double eps;///误差限
int main()
{
    int ok,cur;///cur表示迭代结束所执行的最后一步
    ///freopen("in.txt","r",stdin);
    while(scanf("%d%d%lf",&n,&N,&eps) != EOF){
        ok = 0;
        for(int i = 0;i<n;i++){///输入系数矩阵
            for(int j = 0;j<n;j++){
                scanf("%lf",&a[i][j]);
            }
        }
        for(int i = 0;i<n;i++)scanf("%lf",&b[i]);///输入右端项
        for(int i = 0;i<n;i++)scanf("%lf",&x[i][0]);///输入初始解向量

        double maxs = 0;
        for(int k = 1;k<=N;k++){///最大循环步数
            maxs = 0;
            for(int i = 0;i<n;i++){