高斯—赛德尔求解线性方程+C代码

高斯—赛德尔迭代求矩阵特征值公式如下：

举例：

实现代码：

#include "stdio.h"
#include "stdlib.h"
#include "math.h"

#define N 3
#define MAXITER 10000

void main()
{
int i, j, k;
int iter = 0;
int count;
double a[N][N] = {10, -2, -1, -2, 10, -1, -1, -2, 5};
double b[N] = {3, 15, 10};
double new_x[N] = {0};
double old_x[N] = {0};
double e = 0.0000001;
bool flag = false;

//开始迭代
while(iter<MAXITER && !flag)
{
count = 0;
iter++;

for(i=0; i<N; i++)
old_x[i] = new_x[i];

for(i=0; i<N; i++)
new_x[i] = 0;

//高斯—赛德尔计算
for(i=0; i<N; i++)
{

for(j=0; j<N; j++)
{
if(i==j)
continue;

new_x[i] -= a[i][j]*old_x[j];
}

new_x[i] += b[i];

new_x[i] = new_x[i]/a[i][i];
}

//检测是否收敛
for(i=0; i<N; i++)
if(fabs(new_x[i]-old_x[i])>e)
{
count++;
break;
}

if(count==0)
flag = true;

}

for(i=0; i<N; i++)
{
for(j=0; j<N; j++)
printf("%lf ", a[i][j]);
printf("%lf\n", b[i]);
}

for(i=0; i<N; i++)
printf("%lf ", new_x[i]);

printf("\n迭代次数：%d \n", iter);

}

运行结果：

高斯—赛德尔求解线性方程+C代码相关推荐

高斯-赛德尔（Gauss-Seidel）解线性方程组的Matlab实现
高斯-赛德尔(Gauss-Seidel)解线性方程组的Matlab实现代码运行手算例题迭代法解线性方程组的基本思想是构造一串收敛到解的序列,即建立一种从已有近似解计算新的近似解的规则,有不同的 ...
数值计算方法线性方程组的数值解法(4)---向量和矩阵范数(norm) 高斯-赛德尔(Gauss-Seidel)迭代、共轭梯度(Conjugate Gradient)迭代
(范数部分matlab有现成函数,若有需要直接参照matlab_norm) 向量范数设x∈Rn\boldsymbol x\in \boldsymbol R^nx∈Rn则范数||x||满足:∣∣x∣∣ ...
运用雅可比(Jacobi)和高斯-赛德尔(GS公式)求解方程组
实验名称: 运用雅克比和高斯赛德尔公式求解方程组实验目的: a. 比较两种方法的收敛速度 b. 验证收敛条件的正确性实验内容 1.比较两种方法的收敛速度 Jacobi公式原理 xi(k+1)=1a ...
高斯-赛德尔（Gauss-Seidel）迭代法求解线性方程组matlab程序设计
算法步骤例题程序代码 function x=GaussSeidel(A,b,x0,eps,N) n=length(b); if nargin<5N=500; end if nargin< ...
电力系统潮流计算（牛顿-拉夫逊法、高斯-赛德尔法、快速解耦法）【6节点 9节点 14节点 26节点 30节点 57节点】（Matlab代码实现）
【深入浅出强化学习原理入门】高斯-赛德尔（Gauss-Seidel）迭代法
基于模型的动态规划方法中,使用高斯-赛德尔迭代算法求解值函数. 这种迭代算法到底是怎么求解线性方程组的,文章迭代法求解线性方程组给出了非常详细的解释. 但理解该文还需要了解相关的基础知识: 1.范数 ...
三种迭代法解方程组（雅可比Jacobi、高斯-赛德尔Gaisi_saideer、逐次超松弛SOR）
分析用下列迭代法解线性方程组 4 -1 0 -1 0 0 0 -1 4 -1 0 -1 0 5 0 -1 4 -1 0 -1 -2 -1 0 -1 4 -1 0 ...
高斯—赛德尔(Gauss - Seidel)迭代法解线性方程组（Matlab程序）
%---高斯-赛德尔迭代法----- %---Gauss - Seidel iteration method clear;clc; % A=[10,-1,-2;-1,10,-2;-1,-1,5]; % ...
matlab中用高斯-赛德尔（Gauss-Seidel）迭代法解线性方程组
原理矩阵形式: %GS迭代法矩阵形式 clear; clc; a=[4,-1,0,-1,0,0;-1,4,-1,0,-1,0;0,-1,4,-1,0,-1;-1,0,-1,4,-1,0;0,-1,0 ...

高斯—赛德尔求解线性方程+C代码

高斯—赛德尔求解线性方程+C代码相关推荐

最新文章

热门文章