迭代法求解线性方程组的收敛问题总结
本讲之前,先将高斯-赛德尔迭代法和雅克比迭代法以及迭代法求解线性方程组贴出来,毕竟收敛问题研究的是迭代方法的收敛问题。
进入主题:
判断迭代法收敛的办法:
1、首先根据方程组的系数矩阵A的特点判断;
2、可根据迭代矩阵的范数判断;
3、只好根据迭代矩阵的谱半径来判断;
————————————————————————————————————————————————————
下面一一解释:(1、3 很重要!)
————————————————————————————————————————————————————
1、根据方程组系数矩阵A的特点判断;
这个特点其实就是该矩阵是否是严格对角占优矩阵,或者可约不可约问题;
严格对角占优:
也就是说矩阵A的每一行对角元素的绝对值都严格大于同行的其他元素绝对值之和,则称A为按行严格对角占优矩阵。
严格对角占优矩阵有什么好处呢?
定理:若线性方程组的系数矩阵为按行严格对角占优矩阵,则解此方程组的雅克比迭代法和高斯-赛德尔迭代法都收敛。
(定理:严格对角占优矩阵也是非奇异矩阵。证明略!)
定理:若线性方程组的系数矩阵A为对称正定矩阵,则解此线性方程组的高斯-赛德尔迭代法收敛。
可约与不可约问题:
应用:
2、根据迭代矩阵的范数判断:
3、讲讲根据迭代矩阵的谱半径来判断的方法以及实例:
由相关定理可知,迭代矩阵的谱半径小于1,则该迭代法收敛;
由此可见,若要判断迭代法是否收敛,则需要先求得迭代矩阵,下面分别讲解雅克比迭代法以及高斯-赛德尔迭代法的迭代矩阵的求法:
1》 雅克比迭代法:
简言之,雅克比迭代法的迭代矩阵就是上面的B矩阵;
2》 高斯-赛德尔迭代法:
同理可见,高斯-赛德尔迭代法的迭代矩阵为。
之后,便是求迭代矩阵的谱半径了,首先什么是谱半径呢?
简言之,就是特征值模的最大值,因此需要求迭代矩阵的特征值。
据此,拿一道题目练练手:
再来一题:
考题演练:
迭代法求解线性方程组的收敛问题总结相关推荐
- 雅克比(Jacobi)迭代法求解线性方程组
长博文不利于翻阅,于是又将Jacobi迭代法单独出来了. 这篇博文把高斯-赛德尔迭代法和雅克比迭代法都放到一起了,个人觉得看着有点累.(迭代法求解线性方程组),不过还是要看的,因为它引出了迭代法. 进 ...
- 分别用雅可比(Jacobi)迭代法和高斯—塞德尔(Gauss—Seidel)迭代法求解线性方程组(转载)
版权声明:本文为博主原创文章,未经博主允许不得转载. https://blog.csdn.net/yyywww666/article/details/42805071 算法介绍(迭代法介绍): 代码C ...
- MATLAB Jacobi迭代法 求解线性方程组
文章目录 前言 一.Jacobi迭代法是什么? 二.对应的编程思想以及公式推导 1.Jacobi迭代法 公式推导 2.Jacobi迭代法求解线性方程组 例子 3.Jacobi迭代法 编程实现 总结 前 ...
- 【数值分析】Jacobi、Seidel和Sor迭代法求解线性方程组(附matlab代码)
线性方程组迭代解法公式类似非线性方程求根的简单迭代法公式, 有Jacobi迭代法. Seidel迭代法及Sor法等. 题目1-Jacobi迭代法和Seidel迭代法 来源:<数值分析>第5 ...
- 分别用雅可比(Jacobi)迭代法和高斯—塞德尔(Gauss—Seidel)迭代法求解线性方程组
算法介绍(迭代法介绍): 代码C语言实现; # include<stdio.h> # include<math.h> # define N 6 /* *使用雅可比迭代法和高斯- ...
- 【计算方法】雅克比迭代法-高斯赛德尔迭代法求解线性方程组(c语言实现)
计算方法–用雅克比迭代法和高斯赛德尔迭代法求解线性方程组 AX = b 雅克比迭代法的数学基础 x = Gx + b x = (A ^ -1) b A = L + D + U(L是一个下三角阵,他的对 ...
- sor迭代法求解线性方程组(Matlab)
sor迭代法求解线性方程组(Matlab) 函数文件(sor.m) 该函数不含停止的误差判断条件 function [x,n] = sor(A,b,x,w,it_max) % 求线性方程组的sor(s ...
- 松弛迭代法matlab,逐次超松弛迭代法求解线性方程组的MATLAB实现
function [X_reality,n_reality] = SOR(A,b,X_start,w,n_limit,tolerance) %% % A为迭代的系数矩阵 % b为方程组右边的常数项(列 ...
- c语言如何计算出迭代次数,计算方法——C语言实现——迭代法求解线性方程组...
最近在上计算方法这门课,要求是用MATLAB做练习题,但是我觉得C语言也很棒棒啊~ 题目: 和直接法不同,迭代法是一种逐次逼近的方法,将复杂问题简单化,求比较大的方程组时一般都不会用直接法.迭代法有好 ...
最新文章
- 【c语言|Python】求一元二次方程ax^2+bx+c=0的根
- 犹太人的思维习惯 (转载)
- golang 判断 两个slice 是否相等
- 小学生python-小学生学python(三)
- 【LeetCode】103# 二叉树的锯齿形层次遍历
- 网站维护404页面带音乐源码CSS本地化
- 泛型之类型擦除和桥接方法
- 73. 文件上传函数封装
- 图数据库_多款业界流行图数据库对比
- github 不可不知的10个Github功能
- 软考初级信息处理技术员(一)
- 蓝牙驱动卸载后自动安装_外星人的控制中心下载,安装及常见问题处理方法
- VGA接口、DVI接口、HDMI接口
- 上课笔记-机器学习(4)-泰坦尼克号沉没乘客存活分析
- mysql对韵母分组,基于MySQL的中文发音查询的元级实现
- c语言中100m怎么表示,科目三道路驾驶技能和安全文明驾驶常识考试满分分别为100分,成绩分别达到80和90分的为合格。...
- R语言如何绘制韦恩图(6)
- 5个常用的CMD命令
- VSTO C# 操作excel
- ffmpeg源码学习笔记三
热门文章
- Cisco Catalyst 2960系列交换机资料
- 判断两个多边形相交的面积_聊聊3D模型组件—顶点,边,多边形
- 栈和堆存储在计算机RAM中,堆内存和栈内存及C++内存分配
- 太阳花图片_6种漂亮的花,种子已经成熟,别忘采摘种子,来年再播种
- python 百度词典_用 Python 3 写的命令行百度词典
- c语言中菲薄拉次数列,在VB子函数如何写出菲薄拉起数列
- linux下软件编译终止,[2018年最新整理]linux下编译运行程序命令大全.ppt
- numpy 中的三个特别的索引操作 c_, r_, s_
- 测量 XW-HLR26-24G 微波雷达模块的性能测试
- 第十六届智能车竞赛过程中都发生了什么:怎么感到今年更难呢?