雅克比(Jacobi)迭代法解线性方程组(Matlab程序)
%---雅克比迭代法-----
%---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程序)相关推荐
- 牛顿迭代法解线性方程matlab程序,牛顿迭代法matlab程序(解线性方程组)
<牛顿迭代法matlab程序(解线性方程组)>由会员分享,可在线阅读,更多相关<牛顿迭代法matlab程序(解线性方程组)(4页珍藏版)>请在金锄头文库上搜索. 1.牛顿迭代法 ...
- 雅各比迭代法程序matlab,雅克比(Jacobi)迭代法解线性方程组(Matlab程序)
%---雅克比迭代法----- %---Jacobi iteration method clear;clc; % A=[10,-1,-2;-1,10,-2;-1,-1,5]; % b=[72,83,4 ...
- 雅克比 (Jacobi) 迭代法求线性方程组的解及其代码
简介 实例 可以看到,当迭代次数 k 增大时,迭代值会越来越接近解 x1 = x2 = x3 = 1. 代码 #include<iostream> #include<math.h&g ...
- 雅可比(Jacobi)迭代法解线性方程组的Matlab实现
雅可比(Jacobi)迭代法解线性方程组的Matlab实现 代码 运行 代码 迭代法解线性方程组的基本思想是构造一串收敛到解的序列,即建立一种从已有近似解计算新的近似解的规则,有不同的计算规则得到不同 ...
- 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 ...
- Jacobi迭代法分量形式matlab,mtalab中jacobi迭代法
一.实验目的及题目 1.1 实验目的: (1)学会用高斯列主元消去法,LU 分解法,Jacobi 迭代法和 Gauss-Seidel 迭代法解线性 方程组. (2)学会用 Matlab 编写..... ...
- MATLAB Jacobi迭代法 求解线性方程组
文章目录 前言 一.Jacobi迭代法是什么? 二.对应的编程思想以及公式推导 1.Jacobi迭代法 公式推导 2.Jacobi迭代法求解线性方程组 例子 3.Jacobi迭代法 编程实现 总结 前 ...
- 超松弛迭代法解线性方程组c语言,超松弛迭代法解线性方程组.doc
PAGE PAGE 2 姓名:___________________________ 设计题目:超松弛迭代法解线性方程组 专业: 摘要 本文是在matlab环境下熟悉的运用计算机编程语言并结合超松弛变 ...
- Jacobi迭代法分量形式matlab,实验4 Jacobi迭代法和GS迭代
<数值分析>实验4 一.实验名称:Jacobi迭代法和Gauss-Seidel迭代法 二.实验目的: 熟悉求解线性方程组的Jacobi迭代法和Gauss-Seidel迭代法. 三.实验要求 ...
最新文章
- ui动效 unity_Unity - UIWidgets 2. 控件组合
- leetcode 79:simplify path
- ef1a启动子_常见的真核和原核表达系统的启动子(promoters)
- 腾讯云张贤国:北大将我从少年变成了技术研究者
- TensorFlow: Could not load requested Qt binding.
- 安卓手机中 P 的输入
- 云时代架构读后感4--IT架构的本质
- 【2017年第3期】专题:面向社会治理和服务的大数据
- No module named ‘win32gui‘ 的解决方法(踩坑之旅)
- Python--正则表达式
- echarts 柱状图 横坐标文字纵向展示(超出显示...)
- mysql 循环查询_mysql循环查询(mysql循环语句)
- 面试案例(2019)一
- oracle漏洞修补,01-oracle漏洞修复
- JAVA应用生产问题排查步骤
- MACBOOK 连接不上wifi的解决办法
- python安装失败0x80240017
- 零钱兑换问题c语言编程,leetcode- 零钱兑换 II(背包问题-总结-复盘)
- 倍投技巧 - 凯利公式教你如何用正确的方法投资
- php根据日期判断星座if函数,js根据日期判断星座的代码