解线性方程组的迭代法(雅可比、高斯-塞德尔迭代法)
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
解线性方程组的迭代法(雅可比、高斯-塞德尔迭代法)相关推荐
- 计算方法:列主元消去法,LU分解法, 雅可比迭代法,高斯塞德尔迭代法 解线性方程(C++)
Matrix.h包括矩阵类Matrix的定义,Matrix.cpp包括该类成员函数的实现,LinearEqu.h包括线性方程类LinearEqu的定义,继承自Matrix类,其中solve()方法为列 ...
- 解线性方程组的迭代法(高斯-塞德尔迭代法)
MATLAB函数文件 function x=Agui_GS(A,b) %方程Ax=b,x0为初始向量 %e为精度,N为最大迭代次数 n=length(b); N=100; e=1e-4; x0=zer ...
- 高斯--塞德尔迭代法求方程组的解(Python实现)
数值分析题目 求方程组 {5x1+2x2+x3=−12−x1+4x2+2x3=202x1+−3x2+10x3=2\left\{ \begin{array}{c} 5x_1+2x_2 + x_3 = - ...
- 雅可比(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 ...
- MATLAB实现雅可比与高斯塞德尔迭代
概述 用MATLAB编程实现,形成m函数文件.输入A,b矩阵,无返回值,解得x向量直接显示在命令行窗口,同时绘制出x向量的收敛曲线. A = [ 2 − 1 1 1 1 1 1 1 − 2 ] b = ...
- [计算机数值分析]高斯-塞德尔迭代公式解线性方程组
在雅可比迭代公式的基础上,对于收敛的迭代过程,所求出的"新值"常比"老值"更准确些,因此可以用它替代老值作进一步的计算,这样的思想就是著名的高斯-塞德尔迭代公式 ...
- 基于高斯塞德尔方法的超松弛迭代法MATLAB实现
概述 用MATLAB编程实现,形成m函数文件.输入A,b矩阵,与松弛因子 ω \omega ω,无返回值.设置精度要求 ε = 1 0 − 4 \varepsilon=10^{-4} ε=10−4.所 ...
- 方程组求解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 ...
- 分别用雅可比(Jacobi)迭代法和高斯—塞德尔(Gauss—Seidel)迭代法求解线性方程组(转载)
版权声明:本文为博主原创文章,未经博主允许不得转载. https://blog.csdn.net/yyywww666/article/details/42805071 算法介绍(迭代法介绍): 代码C ...
最新文章
- usaco Overfencing 穿越栅栏(BFS)
- C# WinForm使用乐器数字接口
- ORACLE 11G 搭建dataguard详细步骤(所有操作总结)
- 数学习题:求解不定方程a^2 + b^2 = 3025
- 开源 -- 机器学习相关报道
- Spring之Bean作用域
- oracle sql monitor
- 怎么用万用表检测场效应管的好坏?
- UI设计初学者应该如何入门?
- 无法启动计算机上rasman,解决Windows无法启动Remote Access Connection Manager服务
- cad放大_cad快捷键大全amp;鼠标各键用法
- 国徽FLASH SM25QH256M烧录问题总结
- hive 各种 join (left outer join、join、full outer join)
- MAC 启动jar失败 显示 Unable to access jarfile xxx.jar
- 场效应晶体管在电路中的五大作用,你了解哪一种?
- 【莫烦】python基础教程
- Leetcode-545. Boundary of Binary Tree
- 计算机网络:验证性试验
- LeetCode695. 岛屿的最大面积(C++版)
- WP微信机器人插件v2(更新)