用高斯消元法解决线性系统问题
moboyou 2025-04-23 14:15 22 浏览
高斯法解决线性方程组
线性方程组的基本运算
对任何线性方程组进行三种操作可得到一个等价的方程组:
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的逆矩阵。关于逆矩阵的另外一种求法请参见什么是矩阵的逆矩阵 。
- 上一篇:MATLAB矩阵
- 下一篇:越览(44)——Matlab入门学习(3)之矩阵运算
相关推荐
- 电子EI会议!投稿进度查
-
今天为大家推荐一个高性价比的电子类EI会议——IEEE电子与通信工程国际会议(ICECE2024)会议号:IEEE#62199截稿时间:2024年3月25日召开时间与地点:2024年8月15...
- 最“稳重”的滤波算法-中位值滤波算法的思想原理及C代码实现
-
在信号处理和图像处理领域,滤波算法是一类用于去除噪声、平滑信号或提取特定特征的关键技术。中位值滤波算法是一种常用的非线性滤波方法,它通过取一组数据的中位值来有效减小噪声,保留信号的有用特征,所以是最稳...
- 实际工程项目中是怎么用卡尔曼滤波的?
-
就是直接使用呀!个人认为,卡尔曼滤波有三个个关键点,一个是测量,一个是预测,一个是加权测量:通过传感器,获取传感器数据即可!预测:基于模型来进行数据预测;那么问题来了,如何建模?有难有易。加权:主要就...
- 我拿导弹公式算桃花,结果把自己炸成了烟花
-
第一章:学术圈混成“顶流”,全靠学生们把我写成段子最近总有人问我:“老师,您研究导弹飞行轨迹二十年,咋还顺带研究起月老红绳的抛物线了?”我扶了扶眼镜,深沉答道:“同志,导弹和爱情的本质都是动力学问题—...
- 如何更好地理解神经网络的正向传播?我们需要从「矩阵乘法」入手
-
图:pixabay原文来源:medium作者:MattRoss「机器人圈」编译:嗯~阿童木呀、多啦A亮介绍我为什么要写这篇文章呢?主要是因为我在构建神经网络的过程中遇到了一个令人沮丧的bug,最终迫...
- 电力系统EI会议·权威期刊推荐!
-
高录用率EI会议推荐:ICPSG2025(会议号:CFP25J66-PWR)截稿时间:2025年3月15日召开时间与地点:2025年8月18-20日·新加坡论文集上线:会后3个月内提交至S...
- EI论文写作全流程指南
-
推荐期刊《AppliedEnergy》是新能源领域权威EI/SCI双检索期刊,专注能源创新技术应用。刊号:ISSN0306-2619|CN11-2107/TK影响因子:11.2(最新数...
- JMSE投稿遇坑 实验结果被推翻
-
期刊基础信息刊号:ISSN2077-1312全称:JournalofMarineScienceandEngineering影响因子:3.7(最新JCR数据)分区:中科院3区JCRQ2(...
- 斩获国际特等奖!兰理工数学建模团队为百年校庆献礼
-
近日,2019年美国大学生数学建模竞赛(MCM-ICM)成绩正式公布。兰州理工大学数学建模团队再创佳绩,分别获得国际特等奖(OutstandingWinner)1项、一等奖(Meritorious...
- 省气象台开展人员大培训岗位大练兵学习活动
-
5月9日,省气象台组织开展首次基于Matlab编程语言的数值模式解释应用培训,为促进研究性业务发展,积极开展“人员大培训、岗位大练兵”学习活动起到了积极作用。此次培训基于实际业务需求,着眼高原天气特色...
- 嵌入式软件培训
-
培训效果:通过系统性的培训学习,理论与实践相结合,可以胜任相关方向的开发工作。承诺:七大块专业培训,可以任意选择其中感兴趣的内容进行针对性地学习,每期培训2个月,当期没学会,可免费学习一期。本培训内容...
- 轧机支承辊用重载中低速圆柱滚子轴承滚子修形探讨
-
摘 要:探讨了轧机支承辊用重载中低速圆柱滚子轴承滚子修形的理论和方法,确定关键自变量。使用Romax软件在特定载荷工况条件下对轴承进行数值模拟分析,确定关键量的取值范围。关键词:轧机;圆柱滚子轴承;滚...
- 数学建模EI刊,如何避雷?
-
---权威EI会议推荐会议名称:国际应用数学与工程建模大会(ICAMEM)截稿时间:2025年4月20日召开时间/地点:2025年8月15日-17日·新加坡论文集上线:会后2个月内由Sp...
- 制造工艺误差,三维共轭齿面怎样影响,双圆弧驱动的性能?
-
文/扶苏秘史编辑/扶苏秘史在现代工程领域,高效、精确的传动系统对于机械装置的性能和可靠性至关重要,谐波传动作为一种创新的机械传动方式,以其独特的特性在精密机械领域引起了广泛关注。在谐波传动的进一步优化...
- 测绘EI会议——超详细解析
-
【推荐会议】会议名称:国际测绘与地理信息工程大会(ICGGE)会议编号:71035截稿时间:2025年3月20日召开时间/地点:2025年8月15-17日·德国慕尼黑论文集上线:会后2个...
- 一周热门
- 最近发表
- 标签列表
-
- curseforge官网网址 (16)
- 外键约束 oracle (36)
- oracle的row number (32)
- 唯一索引 oracle (34)
- oracle in 表变量 (28)
- oracle导出dmp导出 (28)
- oracle 数据导出导入 (16)
- oracle两个表 (20)
- oracle 数据库 使用 (12)
- 启动oracle的监听服务 (13)
- oracle 数据库 字符集 (20)
- powerdesigner oracle (13)
- oracle修改端口 (15)
- 左连接 oracle (15)
- oracle 标准版 (13)
- oracle 转义字符 (14)
- asp 连接 oracle (12)
- oracle安装补丁 (19)
- matlab三维图 (12)
- matlab归一化 (16)
- matlab求解方程 (13)
- matlab坐标轴刻度设置 (12)
- matlab脚本 (14)
- matlab多项式拟合 (13)
- matlab阶跃函数 (14)