迭代法(雅可比+高斯塞尔德+SOR）

肥仔甜甜

于 2019-05-31 18:10:00 发布

阅读量1.9k

点赞数

分类专栏：计算方法

本文链接：https://blog.youkuaiyun.com/weixin_43242836/article/details/90720925

版权

本文深入探讨了迭代法在数值解法中的应用，重点讲解了雅可比法、高斯塞尔德法以及改进的SOR（松弛法）的原理与实现步骤。通过对这些迭代方法的比较，阐述了它们在解决线性和非线性方程组中的优势和适用场景。

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

#include <iostream>
#include <vector>
#include <math.h>
using namespace std;
#define eps 0.00001
#define maxn 200
typedef vector<vector<double>> Matrix;
void dispRes(vector<double> r);
vector<double> Jacob(Matrix a);
vector<double> Gauss_Seidel(Matrix a);
vector<double> SOR(Matrix a, double w);
bool check(vector<double> x, vector<double> y);
//上面三个方法只做了求解没有验证收敛性
int main()
{
   
    ios::sync_with_stdio("false");
    int n;
    cin >> n;
    Matrix M1(n, vector<double>(n + 1, 0));
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n + 1; j++)
            cin >> M1[i][j];
    dispRes(