%---雅克比迭代法-----
%---Jacobi 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   %库函数计算结果
B=zeros(N,N);g=zeros(N,1);
x=zeros(N,1);%迭代初始值
eps=0.001;%相邻解的距离小于该数时,结束迭代
% for i=1:N
%     B(i,:)=-A(i,:)/A(i,i);
%     B(i,i)=0;       %迭代矩阵
%     g(i)=b(i)/A(i,i);
% end
D=diag(diag(A));
E=-tril(A,-1);%下三角
F=-triu(A,1);%上三角
B=inv(D)*(E+F);g=inv(D)*b;

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

雅克比(Jacobi)迭代法解线性方程组(Matlab程序)相关推荐

  1. 牛顿迭代法解线性方程matlab程序,牛顿迭代法matlab程序(解线性方程组)

    <牛顿迭代法matlab程序(解线性方程组)>由会员分享,可在线阅读,更多相关<牛顿迭代法matlab程序(解线性方程组)(4页珍藏版)>请在金锄头文库上搜索. 1.牛顿迭代法 ...

  2. 雅各比迭代法程序matlab,雅克比(Jacobi)迭代法解线性方程组(Matlab程序)

    %---雅克比迭代法----- %---Jacobi iteration method clear;clc; % A=[10,-1,-2;-1,10,-2;-1,-1,5]; % b=[72,83,4 ...

  3. 雅克比 (Jacobi) 迭代法求线性方程组的解及其代码

    简介 实例 可以看到,当迭代次数 k 增大时,迭代值会越来越接近解 x1 = x2 = x3 = 1. 代码 #include<iostream> #include<math.h&g ...

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

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

  5. Jacobi迭代法分量形式matlab,解线性方程组的Jacobi迭代法.pdf

    1 Jacobi 迭代法迭代法 若A为非奇异矩阵且0 ii a 1 2 in 11121 21222 12 n n nnnn aaa aaa A aaa 将原方程组等价变形为 1 n iijj jj ...

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

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

  7. MATLAB Jacobi迭代法 求解线性方程组

    文章目录 前言 一.Jacobi迭代法是什么? 二.对应的编程思想以及公式推导 1.Jacobi迭代法 公式推导 2.Jacobi迭代法求解线性方程组 例子 3.Jacobi迭代法 编程实现 总结 前 ...

  8. 超松弛迭代法解线性方程组c语言,超松弛迭代法解线性方程组.doc

    PAGE PAGE 2 姓名:___________________________ 设计题目:超松弛迭代法解线性方程组 专业: 摘要 本文是在matlab环境下熟悉的运用计算机编程语言并结合超松弛变 ...

  9. Jacobi迭代法分量形式matlab,实验4 Jacobi迭代法和GS迭代

    <数值分析>实验4 一.实验名称:Jacobi迭代法和Gauss-Seidel迭代法 二.实验目的: 熟悉求解线性方程组的Jacobi迭代法和Gauss-Seidel迭代法. 三.实验要求 ...

最新文章

  1. ui动效 unity_Unity - UIWidgets 2. 控件组合
  2. leetcode 79:simplify path
  3. ef1a启动子_常见的真核和原核表达系统的启动子(promoters)
  4. 腾讯云张贤国:北大将我从少年变成了技术研究者
  5. TensorFlow: Could not load requested Qt binding.
  6. 安卓手机中 P 的输入
  7. 云时代架构读后感4--IT架构的本质
  8. 【2017年第3期】专题:面向社会治理和服务的大数据
  9. No module named ‘win32gui‘ 的解决方法(踩坑之旅)
  10. Python--正则表达式
  11. echarts 柱状图 横坐标文字纵向展示(超出显示...)
  12. mysql 循环查询_mysql循环查询(mysql循环语句)
  13. 面试案例(2019)一
  14. oracle漏洞修补,01-oracle漏洞修复
  15. JAVA应用生产问题排查步骤
  16. MACBOOK 连接不上wifi的解决办法
  17. python安装失败0x80240017
  18. 零钱兑换问题c语言编程,leetcode- 零钱兑换 II(背包问题-总结-复盘)
  19. 倍投技巧 - 凯利公式教你如何用正确的方法投资
  20. php根据日期判断星座if函数,js根据日期判断星座的代码

热门文章

  1. 【分享】pushplus入驻集简云平台,实现无代码集成数百款应用
  2. Xcode怎么退回旧版本?
  3. [UE4] LogicDriver状态机盒体颜色切换简单案例
  4. 树莓派基础实验24:超声波测距传感器实验
  5. java native函数库_Java 层调用 Native 层函数的两种方式
  6. 史上最全因果推断合集-12(因果推断在哈啰出行的实践探索)
  7. 计算机网络题库——第3章数据链路层
  8. 服务器添加账号失败是怎么回事啊,outlook添加新账户时失败,该怎么办
  9. MDK5 nRF BLE(蓝牙低功耗)
  10. vuex 设计思路和实现