矩陣分析-線性系統-2 高斯消元法、高斯-若爾當消元法

2023-10-12 01:18

本文主要是介绍矩陣分析-線性系統-2 高斯消元法、高斯-若爾當消元法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

http://www.cnblogs.com/pegasus/archive/2011/07/31/2123195.html  

1. 高斯消元法

高斯消元法(Gaussian elimination是求解線性方陣組的一種算法,它也可用來求矩陣的秩,以及求可逆方陣的逆矩陣。它通過逐步消除未知數來將原始線性系統轉化為另一個更簡單的等價的系統。它的實質是通過初等行變化(Elementary row operations),將線性方程組的增廣矩陣轉化為行階梯矩陣(row echelon form)。總結起來,如下步驟所示

 

以下面方程組為例,它的執行步驟為

                          image

1)構造增廣矩陣,即系數矩陣A增加上常數向量b(A|b)

                          image

2)通過以交換行、某行乘以非負常數和兩行相加這三種初等變化將原系統轉化為更簡單的三角形式(triangular form)

     注:這裡的初等變化可以通過系數矩陣A乘上初等矩陣E來實現

                         image

3)從而得到簡化的三角方陣組,注意它更容易解

                         image

4)這時可以使用向後替換算法(Algorithm for Back Substitution)求解得

    z=-4/-4=1,  y=4-2z=4-2=2,  x= (1-y-z)/2=(1-2-1)/2=-1

 

總結上面過程,高斯消元法其實就是下面非常簡單的過程

                            原線性方程組       ——>       高斯消元法     ——> 下三角或上三角形式的線性方程組           ——>  前向替換算法求解(對於上三角形式,采用後向替換算法)

image         \begin{matrix}l_{1,1} x_1 &   &             &            &             & = &    b_1 \\l_{2,1} x_1 & + & l_{2,2} x_2 &            &             & = &    b_2 \\     \vdots &   &      \vdots &     \ddots &             &   & \vdots \\l_{m,1} x_1 & + & l_{m,2} x_2 & + \dotsb + & l_{m,m} x_m & = &   b_m  \\\end{matrix}                 image        

2.高斯-若爾當消元法(Gauss-Jordan Elimination

相對於高斯消元法,高斯-若爾當消元法最後的得到線性方程組更容易求解它得到的是簡化行列式。其轉化後的增高矩陣形式如下,因此它可以直接求出方程的解,而無需使用替換算法。但是,此算法的效率較低。

                             image

 

例子如下:

image          解為image

 

3.實際應用中的高斯消元法

前面介紹了最基本的高斯消元法,現在看看應用於實際問題的實用算法。

3.1 誤差

因為實際應用中,我們總是利用計算機來分析線性系統,而計算機中以有限的數來近似無限的實數,因此產生舍入誤差(roundoff error),進而對解線性系統產生很多影響。

 

一個t位(即精度為t)以image為基的浮點數的表達形式為:imageimage對於一個實數x,其浮點近似值image為最接近x的浮點數,必要時進行近似image

例1:對2位以10為基的浮點算法,image

例2:同樣考慮imageimage

 

以下面系統為例,看看在高斯消元中采用浮點算法會產生什麼效果。

                                                                           image

當以精確解法時,通過將第一行乘以m=89/47,並從第二行中減去得到image,進而利用後向替換算法得x=1,y=-1。

當以3位以10為基的浮點算法時,乘子變為image,因為image因此第一步高斯消元後得

image。此時,因為不能將第2行第1列位置變為0,所以不能將其三角化。從而,我們只能接受將這個位置值賦為0,而不管其實際浮點值。因此,3位浮點高斯消元的結果為image後向算法計算結果為image

3.2 部分主元消元(Partial Pivoting)

盡管無法消除近似誤差的影響,可以采用一些技術來盡量減小這類機器誤差。部分主元消元法在高斯消元的每一步,都選擇列上最大值為軸(通過行變換將其移動)。

3.

这篇关于矩陣分析-線性系統-2 高斯消元法、高斯-若爾當消元法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/192188

相关文章

SpringBoot中六种批量更新Mysql的方式效率对比分析

《SpringBoot中六种批量更新Mysql的方式效率对比分析》文章比较了MySQL大数据量批量更新的多种方法,指出REPLACEINTO和ONDUPLICATEKEY效率最高但存在数据风险,MyB... 目录效率比较测试结构数据库初始化测试数据批量修改方案第一种 for第二种 case when第三种

解决1093 - You can‘t specify target table报错问题及原因分析

《解决1093-Youcan‘tspecifytargettable报错问题及原因分析》MySQL1093错误因UPDATE/DELETE语句的FROM子句直接引用目标表或嵌套子查询导致,... 目录报js错原因分析具体原因解决办法方法一:使用临时表方法二:使用JOIN方法三:使用EXISTS示例总结报错原

MySQL中的LENGTH()函数用法详解与实例分析

《MySQL中的LENGTH()函数用法详解与实例分析》MySQLLENGTH()函数用于计算字符串的字节长度,区别于CHAR_LENGTH()的字符长度,适用于多字节字符集(如UTF-8)的数据验证... 目录1. LENGTH()函数的基本语法2. LENGTH()函数的返回值2.1 示例1:计算字符串

Android kotlin中 Channel 和 Flow 的区别和选择使用场景分析

《Androidkotlin中Channel和Flow的区别和选择使用场景分析》Kotlin协程中,Flow是冷数据流,按需触发,适合响应式数据处理;Channel是热数据流,持续发送,支持... 目录一、基本概念界定FlowChannel二、核心特性对比数据生产触发条件生产与消费的关系背压处理机制生命周期

怎样通过分析GC日志来定位Java进程的内存问题

《怎样通过分析GC日志来定位Java进程的内存问题》:本文主要介绍怎样通过分析GC日志来定位Java进程的内存问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、GC 日志基础配置1. 启用详细 GC 日志2. 不同收集器的日志格式二、关键指标与分析维度1.

MySQL中的表连接原理分析

《MySQL中的表连接原理分析》:本文主要介绍MySQL中的表连接原理分析,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、背景2、环境3、表连接原理【1】驱动表和被驱动表【2】内连接【3】外连接【4编程】嵌套循环连接【5】join buffer4、总结1、背景

python中Hash使用场景分析

《python中Hash使用场景分析》Python的hash()函数用于获取对象哈希值,常用于字典和集合,不可变类型可哈希,可变类型不可,常见算法包括除法、乘法、平方取中和随机数哈希,各有优缺点,需根... 目录python中的 Hash除法哈希算法乘法哈希算法平方取中法随机数哈希算法小结在Python中,

Java Stream的distinct去重原理分析

《JavaStream的distinct去重原理分析》Javastream中的distinct方法用于去除流中的重复元素,它返回一个包含过滤后唯一元素的新流,该方法会根据元素的hashcode和eq... 目录一、distinct 的基础用法与核心特性二、distinct 的底层实现原理1. 顺序流中的去重

关于MyISAM和InnoDB对比分析

《关于MyISAM和InnoDB对比分析》:本文主要介绍关于MyISAM和InnoDB对比分析,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录开篇:从交通规则看存储引擎选择理解存储引擎的基本概念技术原理对比1. 事务支持:ACID的守护者2. 锁机制:并发控制的艺

MyBatis Plus 中 update_time 字段自动填充失效的原因分析及解决方案(最新整理)

《MyBatisPlus中update_time字段自动填充失效的原因分析及解决方案(最新整理)》在使用MyBatisPlus时,通常我们会在数据库表中设置create_time和update... 目录前言一、问题现象二、原因分析三、总结:常见原因与解决方法对照表四、推荐写法前言在使用 MyBATis