高斯-赛德尔(Gauss-Seidel)迭代法求解线性方程组matlab程序设计
算法步骤
例题
程序代码
function x=GaussSeidel(A,b,x0,eps,N)
n=length(b);
if nargin<5N=500;
end
if nargin<4eps=1e-6;
end
if nargin<3x0=zeros(n,1);
end
k=0;
x=zeros(n,1);
while k<Nx(1)=(b(1)-A(1,2:n)*x0(2:n))/A(1,1);for i=2:n-1x(i)=(b(i)-A(i,1:i-1)*x(1:i-1)-A(i,i+1:n)*x0(i+1:n))/A(i,i);endx(n)=(b(n)-A(n,1:n-1)*x(1:n-1))/A(n,n);if norm(x-x0,inf)<epsbreak;endx0=x;k=k+1;
end
if k==Nwarning('已达到迭代次数上限');
end
fprintf('k= %d',k)
end
运行结果
>>A=[0.76,-0.01,-0.14,-0.16;-0.01,0.88,-0.03,0.05;-0.14,-0.03,1.01,-0.12;-0.16,0.05,-0.12,0.72];
>> b=[0.68,1.18,0.12,0.74]';
>> x=GaussSeidel(A,b)
k= 8
x =
1.2762
1.2981
0.4890
1.3027
高斯-赛德尔(Gauss-Seidel)迭代法求解线性方程组matlab程序设计相关推荐
- 高斯—赛德尔(Gauss - Seidel)迭代法解线性方程组(Matlab程序)
%---高斯-赛德尔迭代法----- %---Gauss - Seidel iteration method clear;clc; % A=[10,-1,-2;-1,10,-2;-1,-1,5]; % ...
- 雅可比(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 ...
- 高斯-赛德尔(Gauss-Seidel)解线性方程组的Matlab实现
高斯-赛德尔(Gauss-Seidel)解线性方程组的Matlab实现 代码 运行 手算例题 迭代法解线性方程组的基本思想是构造一串收敛到解的序列,即建立一种从已有近似解计算新的近似解的规则,有不同的 ...
- 高斯—赛德尔求解线性方程+C代码
高斯-赛德尔迭代求矩阵特征值公式如下: 举例: 实现代码: #include "stdio.h" #include "stdlib.h" #include &q ...
- 数值计算方法 线性方程组的数值解法(4)---向量和矩阵范数(norm) 高斯-赛德尔(Gauss-Seidel)迭代、共轭梯度(Conjugate Gradient)迭代
(范数部分matlab有现成函数,若有需要直接参照matlab_norm) 向量范数 设x∈Rn\boldsymbol x\in \boldsymbol R^nx∈Rn则范数||x||满足:∣∣x∣∣ ...
- 数值计算大作业:Jacobi与Gauss -Seidel迭代求解线性方程组(Matlab实现)
作为研究生的入门课,数值计算的大作业算是所有研究生开学的重要编程作业. 我把Jacobi与Gauss -Seidel迭代求解线性方程组的数值计算作业在MATLAB中编程实现.具体的程序详细标注后放在文 ...
- sor迭代法求解线性方程组(Matlab)
sor迭代法求解线性方程组(Matlab) 函数文件(sor.m) 该函数不含停止的误差判断条件 function [x,n] = sor(A,b,x,w,it_max) % 求线性方程组的sor(s ...
- 【深入浅出强化学习原理入门】高斯-赛德尔(Gauss-Seidel)迭代法
基于模型的动态规划方法中,使用高斯-赛德尔迭代算法求解值函数. 这种迭代算法到底是怎么求解线性方程组的,文章迭代法求解线性方程组给出了非常详细的解释. 但理解该文还需要了解相关的基础知识: 1.范数 ...
- 三种迭代法解方程组(雅可比Jacobi、高斯-赛德尔Gaisi_saideer、逐次超松弛SOR)
分析用下列迭代法解线性方程组 4 -1 0 -1 0 0 0 -1 4 -1 0 -1 0 5 0 -1 4 -1 0 -1 -2 -1 0 -1 4 -1 0 ...
最新文章
- Navicat for Oracle
- Reallusion Character Creator 3中文版
- C# 学习笔记(7) 进程
- 报错 findMergedAnnotation activemq
- IT从业者(程序员等)之强身健体-武术
- CodeForces - 817D Imbalanced Array(单调栈)
- 不依赖任何系统API,用c语言实现gbk/utf8/unicode编码转换
- 同等质量下那种图片格式小_最实用的Window小工具合集,总有一款适合你!
- 河南淅川高考成绩查询2021,2021年河南高考成绩查询网址,河南高考成绩查询系统时间安排...
- C语言正交表测试用例,用正交表设计测试用例
- 西威变频器avo下载调试资料_超能士变频器报输出缺相维修
- 制定自动化测试实施计划
- 负压电源设计 TPS54160 负压 Level Shifting Control for an Inverting Buck-boost
- 苏宁易购,淘宝网,京东商城,百万级价格数据海量抓取
- MySQL(二) DDL、DML、DQL、DTL概述
- Android TextView 文字设置不同字体大小和颜色样式
- 文件操作细致详解(下)
- 微信商户平台配置JSAPI支付目录,vue项目,好多坑呐
- JAVA毕业设计公立医院绩效考核系统计算机源码+lw文档+系统+调试部署+数据库
- SVG绘制圆环进度条