用高斯消元法解决线性系统问题
liebian365 2024-12-28 23:51 14 浏览 0 评论
高斯法解决线性方程组
线性方程组的基本运算
对任何线性方程组进行三种操作可得到一个等价的方程组:
1. 将任意两个方程交换
2. 将系统中任何方程的所有项乘以任何不等于零的数
3. 将任意两个方程相加/相减(左右同时)
矩阵行的运算
- 交换两行
- 将一行的倍数添加到另一行
- 将一行乘以一个非零常数
从上面可以看出方程组的变换与矩阵的行变换是一致的，因此可以用矩阵变换解方程组。
行阶梯形矩阵遵循以下规则:
- 如果一行不都是零，那么第一个非零数字，称为主元。
- 对于连续两个以1开头的行，下面一行的1在上面一行的1的右边。
- 任何只有0的行都位于矩阵的底部
阶梯矩阵形式:
最简形的阶梯矩阵:
.
通过将系统的增广矩阵改写为行阶梯形来求解下列线性方程组
解：
系统的增广矩阵如下：
步骤1:在第一列中使用行操作使其生成一个主元1，但本例已有，不用这一步了。
将第(2)行加第（1）行乘-2行，在第(3)行加第（1）行乘-5。
步骤2:在第2列中使用行操作或它们的组合生成一个1(如果没有的话)，本例已有。
将- 1乘以行(2)加到行(3)
上面的矩阵是行阶梯形。相应的线性系统为：
可以将z带回上一个方程得出y, 然后求出z。
最后得出解：
前面谈到最简阶梯形矩阵，我们注意到在主元的1上下都是0。求矩阵的最简阶梯形的方法称为高斯法。
我们继续对上面最后一个增广矩阵做行变换。
将第二行加上第三行乘以6：
接着将第一行加上第三行：
最后将第一行减去第二行：
将增广矩阵改写为最简阶梯形的优点是，无需进一步计算就能给出给定方程组的解，如下所示:
总结一下高斯消元法转换为最简形的阶梯矩阵的方法是：
- 构造一个需要的增广矩阵。
- 互换行，使第一行是的首位是1，如果没有一般可以将一个合适的放在首行。
- 将首行的第一列数a去除第一行的全部元素，使首行第一个数变成1.
- 将首行乘以一个系数消掉首行下面第一列的所有元素，使其它行的首项都是0.
- 重复3-4步，使得其它非零行的首位是1，直到形成一个阶梯矩阵。
- 最后利用行运算把所得的阶梯矩阵变成最简阶梯矩阵。
上面的高斯法也可以用来求矩阵A的逆矩阵，其方法就是：
上述式子就是把增广矩阵AlI经过一系列高斯法的行变换，使得AlI变为IlC, C就是A的逆矩阵。关于逆矩阵的另外一种求法请参见什么是矩阵的逆矩阵 。
- 上一篇：基于 PID 控制的离散控制系统实现
- 下一篇：数学建模最常用的10种数学算法
相关推荐
- 4万多吨豪华游轮遇险 竟是因为这个原因……
-
（观察者网讯）4.7万吨豪华游轮搁浅，竟是因为油量太低？据观察者网此前报道，挪威游轮“维京天空”号上周六（23日）在挪威近海发生引擎故障搁浅。船上载有1300多人，其中28人受伤住院。经过数天的调...
- “菜鸟黑客”必用兵器之“渗透测试篇二”
-
"菜鸟黑客"必用兵器之"渗透测试篇二"上篇文章主要针对伙伴们对"渗透测试"应该如何学习？"渗透测试"的基本流程？本篇文章继续上次的分享，接着介绍一下黑客们常用的渗透测试工具有哪些？以及用实验环境让大家...
- 科幻春晚丨《震动羽翼说“Hello”》两万年星间飞行，探测器对地球的最终告白
-
作者|藤井太洋译者|祝力新【编者按】2021年科幻春晚的最后一篇小说，来自大家喜爱的日本科幻作家藤井太洋。小说将视角放在一颗太空探测器上，延续了他一贯的浪漫风格。...
- 麦子陪你做作业(二):KEGG通路数据库的正确打开姿势
-
作者：麦子KEGG是通路数据库中最庞大的，涵盖基因组网络信息，主要注释基因的功能和调控关系。当我们选到了合适的候选分子，单变量研究也已做完，接着研究机制的时便可使用到它。你需要了解你的分子目前已有哪些...
- 知存科技王绍迪:突破存储墙瓶颈，详解存算一体架构优势
-
智东西（公众号：zhidxcom）编辑|韦世玮智东西6月5日消息，近日，在落幕不久的GTIC2021嵌入式AI创新峰会上，知存科技CEO王绍迪博士以《存算一体AI芯片：AIoT设备的算力新选择》...
- 每日新闻播报(September 14)_每日新闻播报英文
-
AnOscarstatuestandscoveredwithplasticduringpreparationsleadinguptothe87thAcademyAward...
- 香港新巴城巴开放实时到站数据 供科技界研发使用
-
中新网3月22日电据香港《明报》报道，香港特区政府致力推动智慧城市，鼓励公私营机构开放数据，以便科技界研发使用。香港运输署21日与新巴及城巴（两巴）公司签署谅解备忘录，两巴将于2019年第3季度，开...
- 5款不容错过的APP: Red Bull Alert，Flipagram，WifiMapper
-
本周有不少非常出色的app推出，鸵鸟电台做了一个小合集。亮相本周榜单的有WifiMapper's安卓版的app，其中包含了RedBull的一款新型闹钟，还有一款可爱的怪物主题益智游戏。一起来看看我...
- Qt动画效果展示_qt显示图片
-
今天在这篇博文中，主要实践Qt动画，做一个实例来讲解Qt动画使用，其界面如下图所示(由于没有录制为gif动画图片，所以请各位下载查看效果)：该程序使用应用程序单窗口，主窗口继承于QMainWindow...
- 如何从0到1设计实现一门自己的脚本语言
-
作者：dong...
- 三年级语文上册 仿写句子 需要的直接下载打印吧
-
描写秋天的好句好段1.秋天来了，山野变成了美丽的图画。苹果露出红红的脸庞，梨树挂起金黄的灯笼，高粱举起了燃烧的火把。大雁在天空一会儿写“人”字，一会儿写“一”字。2.花园里，菊花争奇斗艳，红的似火，粉...
- C++|那些一看就很简洁、优雅、经典的小代码段
-
目录0等概率随机洗牌：1大小写转换2字符串复制...
- 二年级上册语文必考句子仿写，家长打印，孩子照着练
-
二年级上册语文必考句子仿写，家长打印，孩子照着练。具体如下：...
欢迎 你 发表评论:
- 一周热门
- 最近发表
- 标签列表
-
- wireshark怎么抓包 (75)
- qt sleep (64)
- cs1.6指令代码大全 (55)
- factory-method (60)
- sqlite3_bind_blob (52)
- hibernate update (63)
- c++ base64 (70)
- nc 命令 (52)
- wm_close (51)
- epollin (51)
- sqlca.sqlcode (57)
- lua ipairs (60)
- tv_usec (64)
- 命令行进入文件夹 (53)
- postgresql array (57)
- statfs函数 (57)
- .project文件 (54)
- lua require (56)
- for_each (67)
- c#工厂模式 (57)
- wxsqlite3 (66)
- dmesg -c (58)
- fopen参数 (53)
- tar -zxvf -c (55)
- 速递查询 (52)