方程组的直接解法和迭代法 python_数据与算法总结——基本数值算法2(线性方程组)...

d736e0329b0a5b2418cfd64ea53a515f.png

4 基本数值算法

4.2 线性方程组

4.2.1 线性方程组的特性

解的存在性和唯一性

满足下面条件之一,A非奇异,可逆:

如果b属于A的列向量张成的空间,则称方程组是相容的。

范数需要满足次可加性(三角不等式)。

对于n维矢量x,可以定义p范数,p为大于0的整数:

。所以1范数为相加,2范数为平方和开根,
范数为最大的x。

定义矩阵范数:

根据相应的向量范数,我们可以得到:1范数为列向量求和最大,2范数为行向量求和最大。

矩阵范数:满足这些性质:

,如果

对于任一标量

,有

对任意矢量x,有

问题的敏感性和病态性

:A的条件数

如果A奇异,则

矩阵的条件数刻画了矩阵对于非零矢量最大的延伸和压缩能力,矩阵的条件数越大,说明矩阵越接近奇异。

矩阵条件数的性质:

对于任意矩阵A,

;对于单位阵I,
;对于任意矩阵A和标量
,有
;对于任意对角阵

4.2.2 线性方程组的直接解法

主要是两种:高斯消去法和LU分解

高斯消去法主要的计算量消耗在消元过程,时间复杂度为

在消去过程中,对角线上的元素会作为除数,如果它很小,就会发生上溢从而严重影响求解精度。于是就有了一个操作叫做:列选主元。讲的是一个矩阵计算到

的时候,从
向下看,找到一个最大的家伙,然后把那一行整个搬过来,把这一行搬过去,于是就很好的控制了上溢的问题。用时间换精度。

然后LU分解在线性代数里面已经学过了,这里就不再说了......

线性方程组解的精度

残差向量:

解为
,定义残差向量:
。有

当A为良态时,小的相对残差意味着解的相对误差也小。A如果病态,稳定的算法可以得到小的残差,但解的精度不一定高。

高斯-约当法:把A变成一个对角阵。

乔列斯基分解:如果A是对称正定阵则可以使用。

,L为下三角阵。

4.2.3 线性方程组的迭代解法

迭代求解:直接解法的时间复杂度都是

,迭代运算复杂度不超过
,其中k为迭代步数。

不动点迭代

,定义谱半径
为矩阵M所有特征值的最大值。

如果

,则不动点迭代收敛。

分裂A是实现不动点迭代的基本方法,令

,得到
,得到不动点迭代的形式:

,当
时,收敛

雅可比方法

,D为对角阵,L和U为严格的上三角阵和下三角阵。

,带入迭代式:

高斯-赛德方法

,D为对角阵,L和U为严格的上三角阵和下三角阵。

,带入迭代式:

高斯-赛德方法能够及时利用更新过的分量参与下一步计算,因此收敛速度要比雅可比大约快一倍,但这两种方法收敛通常很慢。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/340562.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

高效的企业测试-集成测试(3/6)

本系列的这一部分将展示如何通过代码级以及系统级集成测试来验证我们的应用程序。 (代码级)集成测试 集成测试一词有时在不同的上下文中使用不同。 根据Wikipedia的定义,我指的是在代码级别上验证多个组件之间相互作用的测试。 通常&#x…

mysql level用法_MYSQL使用方法

1.查询一张表: select * from 表名;2.查询指定字段:select 字段1,字段2,字段3….from 表名;3.where条件查询:select 字段1,字段2,字段3 frome 表名 where 条件表达式…

python程序设计之文件_Python程序设计之文件操作(2)

print(sub_path)if os.path.isdir(sub_path):visitdir(sub_path)path1C:UsersQinHsiuPythonProjectsStringoovisitdir(path1)方法二:使用walk()函数来实现#方法二,通过walk()方法指定遍历目录def visidir2(path):if not os.path.isdir(path):print(error!,endn)retu…

带Prometheus的Spring Boot和测微表第4部分:基础项目

在以前的文章中,我们介绍了Spring Micrometer和InfluxDB。 所以你要问我为什么普罗米修斯。 原因是Prometheus在InfluxDB的拉模型与推模型上进行操作。 这意味着,如果将千分尺与InfluxDB一起使用,则在将结果推送到数据库中时肯定会有一些开…

前端如何实现网络速度测试功能_分析Web前端测试要点,从架构原理上进行分析,希望大家能够掌握...

基于Web前端分析过程,大概有十几个测试要点,我们今天主要来讲解结合前五个要点进行详细解说。前端测试点主要针对前端展开,什么叫前端分析呢?就是我们所有的分析和测试要点所站的视角都是针对客户端或者浏览器来对系统进行分析和测…

mysql数据库表中的类型_MySQL数据库中表类型MyISAM与InnoDB的区别

MyISAM 和 InnoDB 讲解InnoDB和MyISAM是许多人在使用MySQL时最常用的两个表类型,这两个表类型各有优劣,视具体应用而定。基本的差别为:MyISAM类型不支持事务处理等高级处理,而InnoDB类型支持。MyISAM类型的表强调的是性能&#xf…

将Websocket与Spring Framework和Vuejs结合使用

Websocket是客户端和服务器之间的全双工(持久)连接,因此两者可以彼此共享信息,而无需重复建立新的连接。 这消除了从客户端重复轮询以从服务器获取更新的需要。 并非所有浏览器都支持Websocket,因此我们利用SockJS ja…

python函数和模块的使用方法_Python学习06_函数和模块的使用

引入在写有些代码的时候,会发现有些步骤重复了多次,他也不像循环,都是相同的东西在重复,而是指做某件事情的步骤方法,做事的人或对象发生了改变,但是方法却没有改变。要想写出高质量的代码,首先…

tmemo 选择消除行_Divi模块,行和部分加入高级动画选项

一切元素的动画选项每个Divi模块,行和部分都带有高级动画选项,你可以使用这些选项来吸引访问者并使页面更加耀眼。Divi引入一个全新的动画系统,并将这些高级动画选项扩展到每个Divi模块,行和部分!这些新选项已合并到一…

java8 streams_Java 8 Friday:使用Streams API时的10个细微错误

java8 streams在Data Geekery ,我们喜欢Java。 而且,由于我们真的很喜欢jOOQ的流畅的API和查询DSL ,我们对Java 8将为我们的生态系统带来什么感到非常兴奋。 Java 8星期五 每个星期五,我们都会向您展示一些不错的教程风格的Java …

python带参数装饰器 函数名_python 全栈开发,Day11(函数名应用,闭包,装饰器初识,带参数以及带返回值的装饰器)...

一、函数名应用函数名是什么?函数名是函数的名字,本质:变量,特殊的变量。函数名(),执行此函数。python 规范写法1. #后面加一个空格,再写内容,就没有波浪线了。2.一行代码写完,下面一…

python逐行写入excel_快来看看Python如何玩转Excel

来源:ID(innerV)如何用Python来操作Excel文件呢?首先,使用pip 包管理器来安装两个包,安装命令:pip install xlrd pip install xlwt我们来看读取excel的例子,第1行,import 导入xlrd包第4行&#…

python中函数参数_Python函数的参数

本章将介绍函数中参数的用法,更多内容请参考:Python学习指南定义函数的时候,我们把参数的名字和位置确定下来,函数的接口定义就完成了。对于函数的调用者来说,只需要知道传递正确的参数,以及函数就返回什么样的值就够了…

Java面试准备:15个Java面试问题

并非所有的访谈都将重点放在算法和数据结构上—通常,访谈通常只侧重于您声称是专家的语言或技术。在此类访谈中,通常没有任何“陷阱”问题,而是它们要求您利用内存和使用该语言的经验–换句话说,它们测试您对编程语言的了解。 但…

python中字符编码使用_python中字符编码是什么?如何转换字符?

屏幕前我们看到的文字,在计算机中也是这样表示吗?事实上,我们所有的信息电脑是无法直接理解的,同样我们所要表达的信息,计算机也不能看懂。这就需要架起一个理解的桥梁,也就是编码帮助我们相互之间进行翻译…

mysql grant proxy on_MySQL 5.7权限的介绍

MySQL支持的权限如下:ALL或ALL PRIVILEGES 代表指定权限等级的所有权限。ALTER 允许使用ALTER TABLE来改变表的结构,ALTER TABLE同时也需要CREATE和INSERT权限。重命名一个表需要对旧表具有ALTER和DROP权限,对新版具有CREATE和INSERT权限。AL…

python 3.9 发布计划_Python 3.9.0 beta4 发布

Python 3.9.0b4 发布了,Python 3.9 目前仍在开发中,此版本 3.9.0b4 是五个计划的 beta 版本预览中的第四个。Beta 版本预览旨在为广大社区提供测试新功能和错误修复以及准备其项目以支持新功能版本的机会。 Python 3.9 部分主要新功能和变更:…

java junit 异常_JUnit:使用Java 8和Lambda表达式测试异常

java junit 异常在JUnit中,有许多方法可以在测试代码中测试异常,包括try-catch idiom JUnit Rule和catch-exception库。 从Java 8开始,我们还有另一种处理异常的方法:使用lambda表达式。 在这篇简短的博客文章中,我将演…

mysql排插问题_MySQL一次数据插入故障记录

某天突然收到报警,数据库大量事务等待,进到数据库后发线大量的插入操作被阻塞,且都是同一个表的。通过 show engine innodb status 发现插入操作都是在等待索引 idx_create_time(create_time) 的 insert intention lock(跟 gap 锁互斥)&#…

python getopt.getopt_python 参数格式化getopt

命令行参数格式化需要python getopt模块:1.python 命令行参数格式化写道try:options,args getopt.getopt(sys.argv[1:],"hp:i:",["help","ip","port"])except getopt.GetoptError:sys.exit()for name,value in options:i…