%算法 p188

%例题 p209

clc;

clear;

A=[1 0.4 0.4

0.4 1 0.8

0.4 0.8 1];

b=[ 1

2

3];

%迭代四次

height=size(A,1);

D=diag(diag(A));

L=tril(A,-1);%得到L

U=triu(A,1);%得到U

% x=[1 ;1 ;1]%初始向量x

x=ones(height,1);

k=0;

for i=1:1:length(A)

for j=1:1:length(A)

sum=0;

if j~=i

sum=sum+A(i,j)*x(j);

end

x(i)=(b(i)-sum)/A(i,i);

k=k+1;

end

% %迭代4次

% %迭代矩阵B

% B=inv(D)*(L+U)

% f=inv(D)*b

% x(1)=B*x(1)+f

% x(2)=B*x(2)+f

% x(3)=B*x(3)+f

end

x1=A\b

x

%上面的三行代码等于下面注释掉的多行代码

% d=diag(A)

% D=zeros(height,height);

% for i=1:1:height

% D(i,i)=d(i)

% end

% L=zeros(height,height);

% U=zeros(height,height);

% for i=1:1:height

% for j=1:1:height

% if i

% L(i,j)=A(i,j)

% else if i>j

% U(i,j)=A(i,j)

% end

% end

% end

% end

%A=D+L+U

%k表示迭代次数

雅克比迭代法的例子matlab,matlab-雅克比迭代法相关推荐

  1. 用MATLAB实现雅克比迭代、高斯-赛德尔迭代以及超松弛迭代

    文章目录 前言 一.解线性方程组的方法 二.解线性方程组的迭代法及其代码实现 1. 迭代法的收敛性 2. 基本参数设置 3. 雅克比(Jacobi)迭代 4. 高斯-塞德尔(Gauss-Seidel) ...

  2. Jacobi迭代法分量形式matlab,mtalab中jacobi迭代法

    一.实验目的及题目 1.1 实验目的: (1)学会用高斯列主元消去法,LU 分解法,Jacobi 迭代法和 Gauss-Seidel 迭代法解线性 方程组. (2)学会用 Matlab 编写..... ...

  3. MATLAB数学实验——Jacobi迭代法Gauss-Seidel迭代法

    MATLAB数学实验--Jacobi迭代法&Gauss-Seidel迭代法 一.迭代算法的数学知识 线性方程组的数值求解方法,有经典的Jacobi和Gauss-Seidel迭代方法. 二者通过 ...

  4. 用matlab求解jacobi,用jacobi迭代法求解线性方程,求助matlab大师

    问题描述: 用jacobi迭代法求解线性方程,求助matlab大师 用jacobi迭代法求解线性方程 9x1-x2+x3=10 -x1+10x2-2x3=7 -2x1+x2+10x3=6,设迭代初值为 ...

  5. matlab循环语句从小到大排列,三个例子讲解MATLAB三种循环

    三个例子讲解MATLAB三种循环 FOR循环 在for和end语句之间的{commands}按数组中的每一列执行一次.在每一次迭代中,x被指定为数组的下一列,即在第n次循环中,x=array(:, n ...

  6. matlab加速迭代法方程求根,【源码】迭代法求根的matlab算法

    [源码]迭代法求根的matlab算法 [源码]迭代法求根的matlab算法 本篇是在课程学习中自己编程实现的迭代法计算非线性方程或者超越方程近似根的算法,写一下,后边便于复习和期末课程设计引用. % ...

  7. 雅可比(Jacobi)迭代法解线性方程组的Matlab实现

    雅可比(Jacobi)迭代法解线性方程组的Matlab实现 代码 运行 代码 迭代法解线性方程组的基本思想是构造一串收敛到解的序列,即建立一种从已有近似解计算新的近似解的规则,有不同的计算规则得到不同 ...

  8. 【源码】牛顿迭代法求根的matlab实现

    牛顿迭代法求根的matlab实现 本篇是在课程学习中自己编程实现的牛顿迭代法计算非线性方程或者超越方程近似根的算法,写一下,后边便于复习和期末课程设计引用. 牛顿迭代法本质上是一种特殊的不动点迭代,只 ...

  9. python实现gauss-seidel迭代公式_python实现高斯(Gauss)迭代法的例子

    python实现高斯(Gauss)迭代法的例子 我就废话不多说了,直接上代码大家一起看吧! #Gauss迭代法 输入系数矩阵mx.值矩阵mr.迭代次数n(以list模拟矩阵 行优先) def Gaus ...

  10. 向量合并 matlab,MATLAB追加向量

    如果有两个行向量 r1 和 r2 这两个行向量中各有 n 和 m 个元素,现在创建行向量 r 并将n和m个元素都放在行向量 r 中,通过附加这些载体,编写:r = [r1,r2] 通过追加这两个向量, ...

最新文章

  1. 测试各种型号LED的反向光敏电流与光伏信号
  2. 石墨变为石墨烯的方法
  3. 固定资产拆分比例怎么计算_持有的基金要被份额拆分?对收益有影响吗?
  4. 第二部分面向对像基础第五章
  5. 显示数量_SOLIDWORKS 标注螺纹不显示数量?Why
  6. react native text换行_基于React+Koa实现React SSR服务端渲染
  7. Kali Linux Web 渗透测试秘籍 第一章 配置 Kali Linux
  8. C++之访问命名空间/名字空间元素的三种方式
  9. 利用WinSock进行无连接的通信
  10. Perl导入代码文件
  11. 半边数据结构(The_Half-Edge_Data_Structure)
  12. python最小值最大化和最大值最小化_OJ 21658::Monthly Expense(二分搜索+最小化最大值)...
  13. yacc语法分析minipascal_用Yacc实现语法分析器-4-编译原理
  14. 黑客帝国 数字雨 屏幕保护程序 linux ubuntu
  15. DynamipsGUI笔记
  16. WebGL中实现三维融合视频投影技术
  17. asan c/c++内存检测
  18. ES6 --promise了解
  19. excel锁定单元格不能修改_简单的Excel技巧,却能让造价人变得如此逆天!
  20. 解读蔚来Q3财报:亏损额度收窄 蔚来汽车幸运的软着陆?

热门文章

  1. c语言中按位异或运算,^按位异或运算符
  2. RHCE(firewalld——查看及使用)
  3. 风场可视化:随机重置
  4. Android实现银行卡号扫描识别
  5. LLDB(二):基础命令详解
  6. 宝鲲财经:六点简单让你做单盈利
  7. 云时代的腾讯运维团队转型:ECUG 10周年大会演讲
  8. 我可以用计算机控制挖掘机,幽默段子:我能用计算机控制挖掘机炒菜……
  9. 022WXS脚本—WXS的特点
  10. isNotBlank与isNotEmpty