高斯—赛德尔(Gauss - Seidel)迭代法解线性方程组(Matlab程序)

原创 2012年03月30日 18:37:28

%---高斯—赛德尔迭代法-----
%---Gauss - Seidel iteration method
clear;clc;
% A=[10,-1,-2;-1,10,-2;-1,-1,5];
% b=[72,83,42]';
A=[ 28,-3,0,0,0;
    -3,38,-10,0,-5;
    -10,0,25,-15,0;
    0,0,-15,45,0;
    0,-5,0,0,30];
b=[10,0,0,0,0]';
N=length(b);    %解向量的维数
fprintf('库函数计算结果:');
x=inv(A)*b      %库函数计算结果

x=zeros(N,1);%迭代初始值
%-----(A=D-E-F)------
D=diag(diag(A));
E=-tril(A,-1);%下三角
F=-triu(A,1);%上三角
B=inv(D-E)*F;g=inv(D-E)*b;
eps=0.001;%相邻解的距离小于该数时,结束迭代

%--------开始迭代-------
for k=1:100 %最大迭代次数为100
    fprintf('第%d次迭代:',k);
    y=B*x+g;
    fprintf('\n与上次计算结果的距离(2范数):%f \n',norm(x-y)^2);
    if norm(x-y)<eps  
        break;
    end
    x=y
end
x


 

高斯—赛德尔迭代法解线性方程组(C语言)

/*-------------高斯—赛德尔迭代法解线性方程组--------- 参考教材:《数值分析》李乃成,梅立泉,科学出版社     《计算方法教程》第二版 凌永祥,陈明逵*/ #inclu...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

基于matlab的Guass-Seidel(高斯--赛德尔) 迭代法求解线性方程组

Guass-Seidel(高斯–赛德尔) 迭代法(简称 G−S 迭代)是对 Jacobi 迭代的一种改进. 了解G-S迭代法之前先了解什么是Jacobi迭代?链接如下: http://blog.c...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

Jacobi迭代法与Gauss-Seidel迭代法

今天刚好有朋友和我讨论泊松图像融合算法,我说我过去文章里给出的是最原始、最直观的实现算法。对于理解泊松融合的原理比较有帮助,但是效率可能并不理想。印象中,泊松融合是有一个以矩阵为基础的快速算法的。但是...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

牛顿迭代法 && 高斯牛顿法

相关知识 多元函数的泰勒展开: 原文:http://blog.sina.com.cn/s/blog_6819cb9b0100mkkx.html 雅克比矩阵和海森矩阵: 原文:http://jaco...
  • lgfun
  • lgfun
  • 2016年11月21日 17:32
  • 4625

数值分析 Gauss-Seidel迭代法求解线性方程组 MATLAB程序实现

Gauss-Seidel迭代法 参考数值分析第四版 颜庆津著 P39 运行输入 运行结果 函数内容(保存为gauss.m文件,在MATLAB中运行) %Gauss-Seidel迭代法求解线性方程组。迭...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:高斯—赛德尔(Gauss - Seidel)迭代法解线性方程组(Matlab程序)
举报原因:
原因补充:

(最多只允许输入30个字)