数值计算基础(二)线性方程解法篇

2023-10-30

概要:介绍了1.直接法(高斯消去法、列主元消去法、LU分解、平方分解、平方分解改进、追赶法)2.迭代法(雅各比迭代、高斯赛德尔迭代、SOR迭代)求解方程3.迭代法收敛性

1.高斯消去法

用途:解方程
核心
将矩阵直接化为上三角矩阵(注意系数不要化简,计算机无法实现
解释
当高斯消元法进行第一步后,相当于用一个初等矩阵左乘A(1) 。不难看出,这个初等矩阵为
在这里插入图片描述
重复到最后
在这里插入图片描述
在这里插入图片描述
高斯消去法实质上产生了一个将A分解为两个三角形矩阵相乘的因式分解

2.列主元消去法

核心
1.n列矩阵从第1列依次排序到n-1列(排序按列的绝对值进行排序
2.高斯消去形成上三角矩阵
优点
TO DO
例子
在这里插入图片描述

3.矩阵三角分解法

由高斯消元可以得到以下定理 成为LU定理 即秩不为零时,可分解为下三角阵和上三角阵的乘积
这一分解称为Doolittle(杜利脱尔)分解。
在这里插入图片描述

3.1直接三角分法或LU分解或杜利脱尔(Doolittle)分解法

适用范围:ALL(△!=0)
核心
L的主对角线是1 下三角
U上三角
例子
在这里插入图片描述
矩阵计算 卡西欧计算器方法

3.2平方根法和改进的平方根法(针对对称正定矩阵)

适用范围:对称正定矩阵
工程计算中线性方程组的系数矩阵常常具有 对称正定性其各阶顺序主子式及全部特征值均大于 0 。 矩阵的这一特性使它的三角分法具有更简单的形式 从而导出一些特殊的解法 。如平方根法和改进的平方根法 。
定理
在这里插入图片描述
这一分解又称为Cholesky(乔列斯基)分解
表达式
LLTX=Ly

例题
在这里插入图片描述

改进
由于对角线不是1,不好看,我们操作一下
在这里插单位发放完入图片描述
3方程
(1)
在这里插入图片描述

(2)
在这里插入图片描述
(3)
AX=b
推导出
在这里插入图片描述
得到y=L-1b
x=(LT)-1D-1y

3.3追赶法(针对三对角矩阵)

适用范围:这个样子的主对角线三对角矩阵
在这里插入图片描述
此系数矩阵的非零元素集中分布在主对角线及其相邻两条次对角线上,称为三对角矩阵。方程组称为三对角方程组
定理
也就是每一行的主对角线占优+有三条数据均不为0的对角线
在这里插入图片描述
XY确定(方程为AX=d)
在这里插入图片描述
步骤
(1)
在这里插入图片描述
(2)
在这里插入图片描述
(3)
在这里插入图片描述
以下为迭代法解方程的方法

4.Jacobi 迭代法

步骤
将第i个方程转化为Xi的形式

例子
在这里插入图片描述

5.高斯-赛德尔(Gauss-Seidel)迭代法

例子
在这里插入图片描述
在这里插入图片描述

6.SOR迭代法

核心
在高斯赛德尔的基础上进行加权平均,值得注意的是此时括号内包含了第i项(相比于高斯赛德尔方法)
在这里插入图片描述
收敛方式
在这里插入图片描述
一般来说,需要ω的值大于1但不能过大此时收敛速度快
例子

在这里插入图片描述
在这里插入图片描述
之后用Excel或者编写函数迭代即可(有机会写)

3.迭代的收敛性

3.1前置知识

3.1.1范数

列范数(又称为1范数)
列的绝对值加和后取最大的那一个称为列范数
在这里插入图片描述
行范数(又称为无穷范数)
类似于上面的1范数
在这里插入图片描述
2范数
在这里插入图片描述
这里可以看这位老哥的文章(写了编程的相关函数)——求矩阵的1,和2范数

补充一点:最大特征值定义在这里插入图片描述
最大特征值和特征向量求法
(这里不多介绍,我们是选修课,计算机计算估计用不上,暂时没实践,但应该是对的)MATLAB求解
(这里介绍一种在B乎上看到的方法,贼√8香)——原文链接
先求任意俩行向量的垂直向量1,再求这个垂直向量1与其中一个行向量的垂直向量2,没试过,但个人感觉好用
(最传统的方法就是这个线性代数书上的方法,初等行变换嘛,例子的话书上有,上面那个链接也有)

3.1.2迭代公式的矩阵形式

这个拆拆拆分没问题 (甘宁狂喜?你的牌太多了)
在这里插入图片描述
雅克比迭代法的矩阵形式
在这里插入图片描述
高斯赛德尔迭代法的矩阵形式
在这里插入图片描述
SOR迭代法的矩阵形式
在这里插入图片描述
三种迭代方法可统一写成
在这里插入图片描述

3.2迭代法收敛性

核心
B的范数小于1
在这里插入图片描述
误差估计分别为事后估计和事前估计

3.3特殊线性方程组迭代收敛性

前置知识补充
严格对角占优(严格意味着没有等号)
在这里插入图片描述
例子
在这里插入图片描述
行弱对角占优
在这里插入图片描述
例子
在这里插入图片描述
可约矩阵
在这里插入图片描述
三种矩阵的收敛性快速判断
1.对于线性方程组 Ax=b,若A为按行(列) 严格对角占优矩阵或为不可约弱对角占优矩阵 ,则 A 为非奇异矩阵,且雅克比迭代与高斯赛德尔迭代法均收敛
2.对于线性代数方程组 Ax=b ,若 A 为**对称正定矩阵则高斯 赛德尔迭代法收敛。
定理:
3.对于线性代数方程组 Ax=b,若 A 为
对称正定矩阵并且 0< ω <2** 则 SOR 迭代法收敛。
在这里插入图片描述
最后总结一下:
放图————
在这里插入图片描述

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

数值计算基础(二)线性方程解法篇 的相关文章

随机推荐

  • thread ‘main‘ panicked at ‘called `Result::unwrap()` on an `Err` value: Os { code: 2, kind: NotFound

    thread main panicked at called Result unwrap on an Err value Os code 2 kind NotFound message 系统找不到指定的文件 src main rs 31 5
  • 大数据课程I4——Kafka的零拷贝技术

    文章作者邮箱 yugongshiye sina cn 地址 广东惠州 本章节目的 掌握Kafka的零拷贝技术 了解常规的文件传输过程 一 常规的网络传输原理 表面上一个很简单的网络文件输出的过程 在OS底层 会发现数据会被拷贝4次 内核态可
  • 苏移笔试

    1 Linux查看一个目录下文件占用多少空间 du sh path 我们有个系统每天要下载数据 为了检查下载数据是否正常 需要查看下载数据目录大小 因为每天的数据都差不多 如果变化太大就说明不正常 查看目录大小的命令是du 当然也可以查看文
  • n行Python代码系列专栏文章目录

    前往老猿Python博客 https blog csdn net LaoYuanPython 一 专栏简介 本专栏介绍使用少量Python代码 n行 n小于等于10 实现特定的功能或能力 并介绍其背后支撑相关的模块和函数 主要用于介绍Pyt
  • 【游戏引擎】win7配置安装firefly【草稿尝试篇】

    firefly是python的开源服务端游戏引擎 而且它使用了不少软件 数据库mysql 缓存memcached等等 配置firefly其实就是配置完整的一个服务器环境 python2 7安装 安装python2 7 mysql memca
  • 感谢Python

    2012 12 20 不得不说 要不是读了Thinking in Java 要不是受不了作者的的唠叨 不停的推荐Python 我目前所学习到的一切也都无从说起了 感谢Python 到目前为止 Python依然是我的最爱 写代码起来最实用 S
  • 会话技术 Cookie +jsp(实现上次登录时间功能)

    会话技术 1 会话 一次会话中包含多次请求和响应 一次会话 浏览器第一次给服务器资源发送请求 会话建立 直到有一方断开为止 2 功能 在一次会话的范围内的多次请求间 共享数据 3 方式 1 客户端会话技术 Cookie 2 服务器端会话技术
  • 什么是煤矿智能化综采工作面?

    转载 https zhuanlan zhihu com p 354323715 什么是煤矿智能化综采工作面 2021年对于煤矿产业是关键性一年 从国家到集团公司对于煤矿的智能化建设都提出了要求 将人工智能 5G通信 大数据技术引入 基本实现
  • stm32编码器接口测速

    本文章转自B站江科大自化协 目录结构 Encoder c include stm32f10x h Device header void Encoder Init void 第一步开启RCC时钟 开启GPIO和定时器的时钟 RCC APB1P
  • c#窗体程序添加图标详细教程

    一 为什么要给窗体程序添加图标 俗话说 人靠衣装 马靠鞍 一个有个性的窗体程序 不能总是使用微软自带的默认图标 而给程序添加一个自己喜欢的图标 不仅彰显水平 还会增加用户的使用体验 所以 窗体程序加入合适的图标值得程序员的学习 下面就开始教
  • 计算机视觉 图像检索与识别

    目录 一 图像检索相关介绍 1 起因 2 参数介绍 图像纹理 词序无关的文本表述 图像分类 图像特征词典 图像整体描述子 二 图像分类 检索具体内容 图像分类 检索 单词的TF IDF权重 倒排表 Inverted file 图像检索流程
  • QT 添加资源文件

    转自 https blog csdn net dianzishi123 article details 86475565 鼠标点击工程文件目录右键选择 添加新文件 在模板中选择 Qt gt Qt Resource File 资源名称 不要有
  • Shell 基础知识

    Shell 基础知识 number 传递的参数 test sh myvalue1 myvalue2 1是第一个参数 2是第二个参数 echo first 1 second is 2 上一条命令 上一条命令的最后一个参数 Ctrl R 历史搜
  • 无人系统群体智能及其研究进展

    来源 无人机 作者 周兴社 武文亮 西北工业大学 计算机学院 陕西 西安 710129 摘 要 群体智能是人工智能的重要发展方向之一 无人系统群体智能作为人工群体智能的主要形态之一 在许多军用和民用领域都具有广阔且重要的应用前景 同时在基础
  • 服务网格实施周期缩短 50%,丽迅物流基于阿里云 ACK 和 ASM 的云原生应用管理实践

    作者 王夕宁 刘强 华相 公司介绍 丽迅物流是百丽旗下专注于时尚产业 为企业提供专业物流及供应链解决方案的服务商 其产品服务主要包括城市落地配 仓配一体 干线运输及定制化解决方案 通过自研智能化物流管理平台 全面助力企业合作集约化发展 目前
  • Ext_面板_Ext.Panel .

    javascript view plain copy print Ext Panel主要配置表 animCollapse Boolean 设置面板折叠展开是否显示动画 Ext Fx可用默认true 否则false applyTo Mixed
  • orm框架有哪些_Java架构—Spring 核心框架体系结构

    很多人都在用spring开发java项目 但是配置maven依赖的时候并不能明确要配置哪些spring的jar 经常是胡乱添加一堆 编译或运行报错就继续配置jar依赖 导致spring依赖混乱 甚至下一次创建相同类型的工程时也不知道要配置哪
  • nodeJS ejs模板引擎 片段视图+视图助手

    Express 的视图系统还支持片段视图 partials 它就是一个页面的片段 通常是重复的 内容 用于迭代显示 通过它你可以将相对独立的页面块分割出去 而且可以避免显式地使 用 for 循环 让我们看一个例子 在 app js 中新增以
  • WPF Window窗体属性

    XAML的三个顶级元素 Window UserControl 用户控件 Page把窗体以网页形式展现 而一个XAML页面里只能有一个顶级元素 顶级元素只能有一个子元素 在窗体里面设置窗体属性
  • 数值计算基础(二)线性方程解法篇

    概要 介绍了1 直接法 高斯消去法 列主元消去法 LU分解 平方分解 平方分解改进 追赶法 2 迭代法 雅各比迭代 高斯赛德尔迭代 SOR迭代 求解方程3 迭代法收敛性 1 高斯消去法 用途 解方程 核心 将矩阵直接化为上三角矩阵 注意系数