算法步骤


例题


程序代码

function x=GaussSeidel(A,b,x0,eps,N)
n=length(b);
if nargin<5N=500;
end
if nargin<4eps=1e-6;
end
if nargin<3x0=zeros(n,1);
end
k=0;
x=zeros(n,1);
while k<Nx(1)=(b(1)-A(1,2:n)*x0(2:n))/A(1,1);for i=2:n-1x(i)=(b(i)-A(i,1:i-1)*x(1:i-1)-A(i,i+1:n)*x0(i+1:n))/A(i,i);endx(n)=(b(n)-A(n,1:n-1)*x(1:n-1))/A(n,n);if norm(x-x0,inf)<epsbreak;endx0=x;k=k+1;
end
if k==Nwarning('已达到迭代次数上限');
end
fprintf('k= %d',k)
end

运行结果

>>A=[0.76,-0.01,-0.14,-0.16;-0.01,0.88,-0.03,0.05;-0.14,-0.03,1.01,-0.12;-0.16,0.05,-0.12,0.72];

>> b=[0.68,1.18,0.12,0.74]';

>> x=GaussSeidel(A,b)

k= 8

x =

1.2762

1.2981

0.4890

1.3027

高斯-赛德尔(Gauss-Seidel)迭代法求解线性方程组matlab程序设计相关推荐

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

    %---高斯-赛德尔迭代法----- %---Gauss - Seidel iteration method clear;clc; % A=[10,-1,-2;-1,10,-2;-1,-1,5]; % ...

  2. 雅可比(Jacobi)及高斯-塞德尔(Gauss_Seidel)迭代法求解线性方程组的matlab现实

    雅可比迭代法的实现代码: function X=Jacobi(A,B,P,delta,max1) %Input -A is a X*N nosingular matrix % -B is a N*1 ...

  3. 高斯-赛德尔(Gauss-Seidel)解线性方程组的Matlab实现

    高斯-赛德尔(Gauss-Seidel)解线性方程组的Matlab实现 代码 运行 手算例题 迭代法解线性方程组的基本思想是构造一串收敛到解的序列,即建立一种从已有近似解计算新的近似解的规则,有不同的 ...

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

    高斯-赛德尔迭代求矩阵特征值公式如下: 举例: 实现代码: #include "stdio.h" #include "stdlib.h" #include &q ...

  5. 数值计算方法 线性方程组的数值解法(4)---向量和矩阵范数(norm) 高斯-赛德尔(Gauss-Seidel)迭代、共轭梯度(Conjugate Gradient)迭代

    (范数部分matlab有现成函数,若有需要直接参照matlab_norm) 向量范数 设x∈Rn\boldsymbol x\in \boldsymbol R^nx∈Rn则范数||x||满足:∣∣x∣∣ ...

  6. 数值计算大作业:Jacobi与Gauss -Seidel迭代求解线性方程组(Matlab实现)

    作为研究生的入门课,数值计算的大作业算是所有研究生开学的重要编程作业. 我把Jacobi与Gauss -Seidel迭代求解线性方程组的数值计算作业在MATLAB中编程实现.具体的程序详细标注后放在文 ...

  7. sor迭代法求解线性方程组(Matlab)

    sor迭代法求解线性方程组(Matlab) 函数文件(sor.m) 该函数不含停止的误差判断条件 function [x,n] = sor(A,b,x,w,it_max) % 求线性方程组的sor(s ...

  8. 【深入浅出强化学习原理入门】高斯-赛德尔(Gauss-Seidel)迭代法

    基于模型的动态规划方法中,使用高斯-赛德尔迭代算法求解值函数. 这种迭代算法到底是怎么求解线性方程组的,文章迭代法求解线性方程组给出了非常详细的解释. 但理解该文还需要了解相关的基础知识: 1.范数 ...

  9. 三种迭代法解方程组(雅可比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 ...

最新文章

  1. Navicat for Oracle
  2. Reallusion Character Creator 3中文版
  3. C# 学习笔记(7) 进程
  4. 报错 findMergedAnnotation activemq
  5. IT从业者(程序员等)之强身健体-武术
  6. CodeForces - 817D Imbalanced Array(单调栈)
  7. 不依赖任何系统API,用c语言实现gbk/utf8/unicode编码转换
  8. 同等质量下那种图片格式小_最实用的Window小工具合集,总有一款适合你!
  9. 河南淅川高考成绩查询2021,2021年河南高考成绩查询网址,河南高考成绩查询系统时间安排...
  10. C语言正交表测试用例,用正交表设计测试用例
  11. 西威变频器avo下载调试资料_超能士变频器报输出缺相维修
  12. 制定自动化测试实施计划
  13. 负压电源设计 TPS54160 负压 Level Shifting Control for an Inverting Buck-boost
  14. 苏宁易购,淘宝网,京东商城,百万级价格数据海量抓取
  15. MySQL(二) DDL、DML、DQL、DTL概述
  16. Android TextView 文字设置不同字体大小和颜色样式
  17. 文件操作细致详解(下)
  18. 微信商户平台配置JSAPI支付目录,vue项目,好多坑呐
  19. JAVA毕业设计公立医院绩效考核系统计算机源码+lw文档+系统+调试部署+数据库
  20. SVG绘制圆环进度条

热门文章

  1. Maven安装和依赖管理详解
  2. Filter - 通过拦截器与动态代理实现敏感词汇过滤
  3. python中的row函数均方_row函数的使用方法
  4. 【PDF操作】简单的PDF文档签名开源工具
  5. Matplotlib系列(六):路径、面片和集合
  6. 测试代理ip是否可用的方法
  7. 40 - xml文档与字典之间的互相转换
  8. 推荐系统学习笔记——特征工程
  9. jquery插件chosen 选择框无法重置
  10. Navicat primium 英文字体出现框框