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

相关文章

Go使用pprof进行CPU,内存和阻塞情况分析

《Go使用pprof进行CPU,内存和阻塞情况分析》Go语言提供了强大的pprof工具,用于分析CPU、内存、Goroutine阻塞等性能问题,帮助开发者优化程序,提高运行效率,下面我们就来深入了解下... 目录1. pprof 介绍2. 快速上手:启用 pprof3. CPU Profiling:分析 C

MySQL表锁、页面锁和行锁的作用及其优缺点对比分析

《MySQL表锁、页面锁和行锁的作用及其优缺点对比分析》MySQL中的表锁、页面锁和行锁各有特点,适用于不同的场景,表锁锁定整个表,适用于批量操作和MyISAM存储引擎,页面锁锁定数据页,适用于旧版本... 目录1. 表锁(Table Lock)2. 页面锁(Page Lock)3. 行锁(Row Lock

Springboot中分析SQL性能的两种方式详解

《Springboot中分析SQL性能的两种方式详解》文章介绍了SQL性能分析的两种方式:MyBatis-Plus性能分析插件和p6spy框架,MyBatis-Plus插件配置简单,适用于开发和测试环... 目录SQL性能分析的两种方式:功能介绍实现方式:实现步骤:SQL性能分析的两种方式:功能介绍记录

最长公共子序列问题的深度分析与Java实现方式

《最长公共子序列问题的深度分析与Java实现方式》本文详细介绍了最长公共子序列(LCS)问题,包括其概念、暴力解法、动态规划解法,并提供了Java代码实现,暴力解法虽然简单,但在大数据处理中效率较低,... 目录最长公共子序列问题概述问题理解与示例分析暴力解法思路与示例代码动态规划解法DP 表的构建与意义动

C#使用DeepSeek API实现自然语言处理,文本分类和情感分析

《C#使用DeepSeekAPI实现自然语言处理,文本分类和情感分析》在C#中使用DeepSeekAPI可以实现多种功能,例如自然语言处理、文本分类、情感分析等,本文主要为大家介绍了具体实现步骤,... 目录准备工作文本生成文本分类问答系统代码生成翻译功能文本摘要文本校对图像描述生成总结在C#中使用Deep

Redis主从/哨兵机制原理分析

《Redis主从/哨兵机制原理分析》本文介绍了Redis的主从复制和哨兵机制,主从复制实现了数据的热备份和负载均衡,而哨兵机制可以监控Redis集群,实现自动故障转移,哨兵机制通过监控、下线、选举和故... 目录一、主从复制1.1 什么是主从复制1.2 主从复制的作用1.3 主从复制原理1.3.1 全量复制

Redis主从复制的原理分析

《Redis主从复制的原理分析》Redis主从复制通过将数据镜像到多个从节点,实现高可用性和扩展性,主从复制包括初次全量同步和增量同步两个阶段,为优化复制性能,可以采用AOF持久化、调整复制超时时间、... 目录Redis主从复制的原理主从复制概述配置主从复制数据同步过程复制一致性与延迟故障转移机制监控与维

Redis连接失败:客户端IP不在白名单中的问题分析与解决方案

《Redis连接失败:客户端IP不在白名单中的问题分析与解决方案》在现代分布式系统中,Redis作为一种高性能的内存数据库,被广泛应用于缓存、消息队列、会话存储等场景,然而,在实际使用过程中,我们可能... 目录一、问题背景二、错误分析1. 错误信息解读2. 根本原因三、解决方案1. 将客户端IP添加到Re

Redis主从复制实现原理分析

《Redis主从复制实现原理分析》Redis主从复制通过Sync和CommandPropagate阶段实现数据同步,2.8版本后引入Psync指令,根据复制偏移量进行全量或部分同步,优化了数据传输效率... 目录Redis主DodMIK从复制实现原理实现原理Psync: 2.8版本后总结Redis主从复制实

锐捷和腾达哪个好? 两个品牌路由器对比分析

《锐捷和腾达哪个好?两个品牌路由器对比分析》在选择路由器时,Tenda和锐捷都是备受关注的品牌,各自有独特的产品特点和市场定位,选择哪个品牌的路由器更合适,实际上取决于你的具体需求和使用场景,我们从... 在选购路由器时,锐捷和腾达都是市场上备受关注的品牌,但它们的定位和特点却有所不同。锐捷更偏向企业级和专