Agui_GS.m:

function x=Agui_GS(A,b)
%方程Ax=b,x0为初始向量
%e为精度,N为最大迭代次数
n=length(b);
N=100;
e=1e-4;
x0=zeros(n,1);
x=x0;
x0=x+2*e;
k=0;
A1=tril(A);
A2=inv(A1);
while norm(x0-x,inf)>e&&k<Nk=k+1;x0=x;x=-A2*(A-A1)*x0+A2*b;format long kdisp(x');
end
if k==N wArning('已达到最大迭代次数');
end

agui_jacobi.m:

function x = agui_jacobi(a,b)
%ax=b,a为系数矩阵,b为右端向量,x0为初始向量n=length(b);        % 精度
N=100;              % 最大迭代次数
e=1e-6;
x0=zeros(n,1);      % 初始向量默认为0向量
x=x0;
x0=x+2*e;
k=0;
D=diag(diag(a));    % 调用 diag 两次将返回一个包含a的对角元素的对角矩阵
L=-tril(a,-1);      % 返回矩阵a的下三角形部分,并取负
U=-triu(a,1);       % 返回矩阵a的上三角形部分,并取负
while norm(x0-x,inf)>e&&k<Nk=k+1;x0=x;x=inv(D)*(L+U)*x+inv(D)*b;  % inv返回逆矩阵kdisp(x');
end
if k==Nwarning('已达到最大迭代次数');
end

解线性方程组的迭代法(雅可比、高斯-塞德尔迭代法)相关推荐

  1. 计算方法:列主元消去法,LU分解法, 雅可比迭代法,高斯塞德尔迭代法 解线性方程(C++)

    Matrix.h包括矩阵类Matrix的定义,Matrix.cpp包括该类成员函数的实现,LinearEqu.h包括线性方程类LinearEqu的定义,继承自Matrix类,其中solve()方法为列 ...

  2. 解线性方程组的迭代法(高斯-塞德尔迭代法)

    MATLAB函数文件 function x=Agui_GS(A,b) %方程Ax=b,x0为初始向量 %e为精度,N为最大迭代次数 n=length(b); N=100; e=1e-4; x0=zer ...

  3. 高斯--塞德尔迭代法求方程组的解(Python实现)

    数值分析题目 求方程组 {5x1+2x2+x3=−12−x1+4x2+2x3=202x1+−3x2+10x3=2\left\{ \begin{array}{c} 5x_1+2x_2 + x_3 = - ...

  4. 雅可比(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 ...

  5. MATLAB实现雅可比与高斯塞德尔迭代

    概述 用MATLAB编程实现,形成m函数文件.输入A,b矩阵,无返回值,解得x向量直接显示在命令行窗口,同时绘制出x向量的收敛曲线. A = [ 2 − 1 1 1 1 1 1 1 − 2 ] b = ...

  6. [计算机数值分析]高斯-塞德尔迭代公式解线性方程组

    在雅可比迭代公式的基础上,对于收敛的迭代过程,所求出的"新值"常比"老值"更准确些,因此可以用它替代老值作进一步的计算,这样的思想就是著名的高斯-塞德尔迭代公式 ...

  7. 基于高斯塞德尔方法的超松弛迭代法MATLAB实现

    概述 用MATLAB编程实现,形成m函数文件.输入A,b矩阵,与松弛因子 ω \omega ω,无返回值.设置精度要求 ε = 1 0 − 4 \varepsilon=10^{-4} ε=10−4.所 ...

  8. 方程组求解matlab实现(朴素高斯求解、LU分解、雅可比迭代方法、高斯-塞德尔方法、连续过松弛(SOR))

    朴素高斯求解 function x= oridinarygauss(a,b) n=length(b); for j=1:n-1for i=j+1:nmult=a(i,j)/a(j,j);for k=j ...

  9. 分别用雅可比(Jacobi)迭代法和高斯—塞德尔(Gauss—Seidel)迭代法求解线性方程组(转载)

    版权声明:本文为博主原创文章,未经博主允许不得转载. https://blog.csdn.net/yyywww666/article/details/42805071 算法介绍(迭代法介绍): 代码C ...

最新文章

  1. usaco Overfencing 穿越栅栏(BFS)
  2. C# WinForm使用乐器数字接口
  3. ORACLE 11G 搭建dataguard详细步骤(所有操作总结)
  4. 数学习题:求解不定方程a^2 + b^2 = 3025
  5. 开源 -- 机器学习相关报道
  6. Spring之Bean作用域
  7. oracle sql monitor
  8. 怎么用万用表检测场效应管的好坏?
  9. UI设计初学者应该如何入门?
  10. 无法启动计算机上rasman,解决Windows无法启动Remote Access Connection Manager服务
  11. cad放大_cad快捷键大全amp;鼠标各键用法
  12. 国徽FLASH SM25QH256M烧录问题总结
  13. hive 各种 join (left outer join、join、full outer join)
  14. MAC 启动jar失败 显示 Unable to access jarfile xxx.jar
  15. 场效应晶体管在电路中的五大作用,你了解哪一种?
  16. 【莫烦】python基础教程
  17. Leetcode-545. Boundary of Binary Tree
  18. 计算机网络:验证性试验
  19. LeetCode695. 岛屿的最大面积(C++版)
  20. WP微信机器人插件v2(更新)

热门文章

  1. 2015年 StoiGdoi 反思总结与未来计划
  2. Biztalk中使用SQL适配器获取数据并用web服务发布的例子
  3. zookeeper源码-分布式锁
  4. Leetcode-952 Largest Component Size by Common Factor(按公因数计算最大组件大小)
  5. LVDS之一_理解SerDes
  6. 如何去除Android应用中的广告
  7. Android 使用AsyncTask 后监听异步加载完毕的动作
  8. 【Python笔记】pygame 游戏框架
  9. Linux NTP时间服务器搭建
  10. 帆软报表登录背后得逻辑