高斯消元法【Gaussian Elimination】
关键词:高斯消元法,变治法
文章目录
- 一、算法针对的问题和应用
- 二、算法介绍
- 1.算法讲解
- 2.一些语言的代码实现
以下是本篇文章正文内容
一、算法针对的问题和应用
1.应用一:解线性方程组
2.应用二:计算一个矩阵的逆
3.应用三:计算矩阵的行列式
二、算法介绍
1.算法讲解
(1)基本思路:
把n个线性方程构成的n元联立方程组变换为一个等价的方程组(方程的解不变),该方程组有着一个上三角形的系数矩阵,该系数矩阵的主对角线下方元素全部为0,然后进行反向替换,从xn,逐个求出原方程组的解。流程见下图:
(2)求解线性方程组的两个阶段。
【前向消失】
伪代码为:
【前向消失的改进】
上述算法仍存在:
问题一:并不总是正确
1>比例因子A[j,i]/A[i,i]的除0操作没能避免。
–>用第i行下面的某行与第i行进行交换,该行的第i个元素不为0
2>比例因子A[j,i]/A[i,i]太大,导致A[j,k]的新值会因为舍入误差而歪曲
–>部分选主元法:每次都去找第i列系数的绝对值最大的行,然后把它作为第i次迭代的起点,保证比例因子的绝对值永远不会大于1
问题二:最内层循环的效率太低
改进算法:
【反向替换】
首先从最后一个方程中可以立即求出xn的值;然后可以把这个值代入到倒数第二个方程来求xn-1,依此类推,直到把最后n-1个变量的已知值代入第一个方程,这样九求出了x1的值。
(3)实现高斯消元的另一种路径——LU分解
A:初始矩阵
L:下三角且主对角线为1的矩阵
U:上三角矩阵
已知:A=LU
实现过程:
Ax=b
LUx=b
令y=Ux
Ly=b =>易求出y
代入,Ux=y =>易求的x
(4)计算一个矩阵的逆
退化矩阵/奇异矩阵:不存在逆的矩阵
非退化矩阵:存在逆的矩阵
判断:当且仅当矩阵的某一行是其他行的一个线性组合时,该矩阵是一个退化矩阵。
用高斯消元法判断是否是退化矩阵:
如果高斯消元法生成的上三角矩阵的主对角线包含0,那该矩阵就是退化矩阵。
求一个非退化的n阶方阵A的逆矩阵:
(5)计算矩阵的行列式
2.一些语言的代码实现
【待补充】
本文是作者在学习数据结构、算法设计与分析课程时,结合自己理解所写的,如有错误或补充,请在评论区留言。
高斯消元法【Gaussian Elimination】相关推荐
- 高斯消元(Gaussian elimination)
我们将分两个部分来说明高斯消去,第一部分是数学方面,第二个部分是如何用 C++ 来实现. 第一部分 数学知识 什么是高斯消去 高斯消元法(或译:高斯消去法),是线性代数规划中的一个算法,可用来为线性方 ...
- Matrix 高斯消元Gaussian elimination 中的complete pivoting和partial pivoting
首先科普下Pivoting的含义 一般翻译为"主元",在对矩阵做某种算法时,首先进行的部分元素.在线性规划的单纯形法中常见. wiki的解释如下: Pivot element (t ...
- 高斯消元法java语言设计_高斯消元法(Gauss Elimination)【超详解模板】
高斯消元法,是线性代数中的一个算法,可用来求解线性方程组,并可以求出矩阵的秩,以及求出可逆方阵的逆矩阵. 高斯消元法的原理是: 若用初等行变换将增广矩阵 化为 ,则AX = B与CX = D是同解方程 ...
- ACM数学模板1 高斯消元 Gaussian Elimination
声明: 题目来源: https://www.luogu.com.cn/problem/P3389 https://www.luogu.com.cn/problem/P2447 https://ac.n ...
- C#,数值计算,高斯消元法与列主元消元法的源代码及数据动态可视化
高斯消元法! 一.高斯消元法 Gaussian Elimination 高斯消元法(或译:高斯消去法),是线性代数中的一个常用算法,常用于求解线性方程组和矩阵的逆. 本程序的运行效果: 1.高斯消元法 ...
- 计算机数值方法:高斯消元法解线性方程组
题目 用高斯消元法解下列线性方程组(要求按三位小数计算) 高斯消元法简介 高斯消元法(Gaussian elimination)是求解线性方阵组的一种算法,它也可用来求矩阵的秩,以及求可逆方阵的逆矩阵 ...
- [POJ 1222] EXTENDED LIGHTS OUT
题目 http://acm.pku.edu.cn/JudgeOnline/problem?id=1222 描述 给你一个5行6列的矩阵分别表示30个灯,矩阵map[i][j]为1表示灯亮着, 0表示灯 ...
- NTL密码算法开源库——模二整数上的矩阵(mat_GF2)
2021SC@SDUSC 模二整数上的矩阵(mat_GF2) 矩阵运算 高斯消元 矩阵运算 具体代码 #include <NTL/matrix.h> #include <NTL/ve ...
- 【直观详解】线性代数的本质
转载自:https://charlesliuyx.github.io/2017/10/06/%E3%80%90%E7%9B%B4%E8%A7%82%E8%AF%A6%E8%A7%A3%E3%80%91 ...
最新文章
- 上传大文件,出现: 413 request Entity too Large错误的解决办法
- 心情不能平静,慢慢的调节
- matlab中prismastic,魅力值是什么意思
- 1344 走格子 (前缀和)
- pycharm 打开cfg高亮
- 如何让gridview中的checkbox根据数据库情况默认选中?
- MySQL主备复制原理、实现及异常处理
- linuxsed替换字符串后保存_Numpy运用-文件读写、存储及字符串处理
- ubuntu nginx配置负载均衡篇(二)
- 大数据入门:各种大数据技术的介绍
- ubuntu之间传输文件
- 从零实现Word2Vec(上)
- python查询手册_Python 手册
- java百鸡问题_算法-百鸡问题 Java
- 微型计算机及接口技术笔记,2010年自考微型计算机及其接口技术笔记串讲
- 如何在Windows 10中打开设置?
- win10怎么将计算机放桌面壁纸,Win10专业版无法设置电脑桌面壁纸怎么办?
- 深度学习项目训练时突然卡住或崩盘问题
- 多元线性回归分析c语言,多元线性回归公式推导及R语言实现
- python 百度云不限速版_现在各位是怎么应对百度网盘限速的?