矩陣分析-線性系統-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

相关文章

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

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

Python主动抛出异常的各种用法和场景分析

《Python主动抛出异常的各种用法和场景分析》在Python中,我们不仅可以捕获和处理异常,还可以主动抛出异常,也就是以类的方式自定义错误的类型和提示信息,这在编程中非常有用,下面我将详细解释主动抛... 目录一、为什么要主动抛出异常?二、基本语法:raise关键字基本示例三、raise的多种用法1. 抛

github打不开的问题分析及解决

《github打不开的问题分析及解决》:本文主要介绍github打不开的问题分析及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、找到github.com域名解析的ip地址二、找到github.global.ssl.fastly.net网址解析的ip地址三

Mysql的主从同步/复制的原理分析

《Mysql的主从同步/复制的原理分析》:本文主要介绍Mysql的主从同步/复制的原理分析,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录为什么要主从同步?mysql主从同步架构有哪些?Mysql主从复制的原理/整体流程级联复制架构为什么好?Mysql主从复制注意

java -jar命令运行 jar包时运行外部依赖jar包的场景分析

《java-jar命令运行jar包时运行外部依赖jar包的场景分析》:本文主要介绍java-jar命令运行jar包时运行外部依赖jar包的场景分析,本文给大家介绍的非常详细,对大家的学习或工作... 目录Java -jar命令运行 jar包时如何运行外部依赖jar包场景:解决:方法一、启动参数添加: -Xb

Apache 高级配置实战之从连接保持到日志分析的完整指南

《Apache高级配置实战之从连接保持到日志分析的完整指南》本文带你从连接保持优化开始,一路走到访问控制和日志管理,最后用AWStats来分析网站数据,对Apache配置日志分析相关知识感兴趣的朋友... 目录Apache 高级配置实战:从连接保持到日志分析的完整指南前言 一、Apache 连接保持 - 性

Linux中的more 和 less区别对比分析

《Linux中的more和less区别对比分析》在Linux/Unix系统中,more和less都是用于分页查看文本文件的命令,但less是more的增强版,功能更强大,:本文主要介绍Linu... 目录1. 基础功能对比2. 常用操作对比less 的操作3. 实际使用示例4. 为什么推荐 less?5.

spring-gateway filters添加自定义过滤器实现流程分析(可插拔)

《spring-gatewayfilters添加自定义过滤器实现流程分析(可插拔)》:本文主要介绍spring-gatewayfilters添加自定义过滤器实现流程分析(可插拔),本文通过实例图... 目录需求背景需求拆解设计流程及作用域逻辑处理代码逻辑需求背景公司要求,通过公司网络代理访问的请求需要做请

Java集成Onlyoffice的示例代码及场景分析

《Java集成Onlyoffice的示例代码及场景分析》:本文主要介绍Java集成Onlyoffice的示例代码及场景分析,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要... 需求场景:实现文档的在线编辑,团队协作总结:两个接口 + 前端页面 + 配置项接口1:一个接口,将o

IDEA下"File is read-only"可能原因分析及"找不到或无法加载主类"的问题

《IDEA下Fileisread-only可能原因分析及找不到或无法加载主类的问题》:本文主要介绍IDEA下Fileisread-only可能原因分析及找不到或无法加载主类的问题,具有很好的参... 目录1.File is read-only”可能原因2.“找不到或无法加载主类”问题的解决总结1.File