关键词:高斯消元法,变治法


文章目录

  • 一、算法针对的问题和应用
  • 二、算法介绍
    • 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】相关推荐

  1. 高斯消元(Gaussian elimination)

    我们将分两个部分来说明高斯消去,第一部分是数学方面,第二个部分是如何用 C++ 来实现. 第一部分 数学知识 什么是高斯消去 高斯消元法(或译:高斯消去法),是线性代数规划中的一个算法,可用来为线性方 ...

  2. Matrix 高斯消元Gaussian elimination 中的complete pivoting和partial pivoting

    首先科普下Pivoting的含义 一般翻译为"主元",在对矩阵做某种算法时,首先进行的部分元素.在线性规划的单纯形法中常见. wiki的解释如下: Pivot element (t ...

  3. 高斯消元法java语言设计_高斯消元法(Gauss Elimination)【超详解模板】

    高斯消元法,是线性代数中的一个算法,可用来求解线性方程组,并可以求出矩阵的秩,以及求出可逆方阵的逆矩阵. 高斯消元法的原理是: 若用初等行变换将增广矩阵 化为 ,则AX = B与CX = D是同解方程 ...

  4. ACM数学模板1 高斯消元 Gaussian Elimination

    声明: 题目来源: https://www.luogu.com.cn/problem/P3389 https://www.luogu.com.cn/problem/P2447 https://ac.n ...

  5. C#,数值计算,高斯消元法与列主元消元法的源代码及数据动态可视化

    高斯消元法! 一.高斯消元法 Gaussian Elimination 高斯消元法(或译:高斯消去法),是线性代数中的一个常用算法,常用于求解线性方程组和矩阵的逆. 本程序的运行效果: 1.高斯消元法 ...

  6. 计算机数值方法:高斯消元法解线性方程组

    题目 用高斯消元法解下列线性方程组(要求按三位小数计算) 高斯消元法简介 高斯消元法(Gaussian elimination)是求解线性方阵组的一种算法,它也可用来求矩阵的秩,以及求可逆方阵的逆矩阵 ...

  7. [POJ 1222] EXTENDED LIGHTS OUT

    题目 http://acm.pku.edu.cn/JudgeOnline/problem?id=1222 描述 给你一个5行6列的矩阵分别表示30个灯,矩阵map[i][j]为1表示灯亮着, 0表示灯 ...

  8. NTL密码算法开源库——模二整数上的矩阵(mat_GF2)

    2021SC@SDUSC 模二整数上的矩阵(mat_GF2) 矩阵运算 高斯消元 矩阵运算 具体代码 #include <NTL/matrix.h> #include <NTL/ve ...

  9. 【直观详解】线性代数的本质

    转载自: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 ...

最新文章

  1. 上传大文件,出现: 413 request Entity too Large错误的解决办法
  2. 心情不能平静,慢慢的调节
  3. matlab中prismastic,魅力值是什么意思
  4. 1344 走格子 (前缀和)
  5. pycharm 打开cfg高亮
  6. 如何让gridview中的checkbox根据数据库情况默认选中?
  7. MySQL主备复制原理、实现及异常处理
  8. linuxsed替换字符串后保存_Numpy运用-文件读写、存储及字符串处理
  9. ubuntu nginx配置负载均衡篇(二)
  10. 大数据入门:各种大数据技术的介绍
  11. ubuntu之间传输文件
  12. 从零实现Word2Vec(上)
  13. python查询手册_Python 手册
  14. java百鸡问题_算法-百鸡问题 Java
  15. 微型计算机及接口技术笔记,2010年自考微型计算机及其接口技术笔记串讲
  16. 如何在Windows 10中打开设置?
  17. win10怎么将计算机放桌面壁纸,Win10专业版无法设置电脑桌面壁纸怎么办?
  18. 深度学习项目训练时突然卡住或崩盘问题
  19. 多元线性回归分析c语言,多元线性回归公式推导及R语言实现
  20. python 百度云不限速版_现在各位是怎么应对百度网盘限速的?

热门文章

  1. 为什么有的测试员路越走越窄?原因在这里
  2. 字节跳动小程序担保支付请求签名算法-Java实现
  3. 电销机器人百度百科_电销机器人好用嘛,其功能以及优缺点介绍
  4. 【自习任我行】任务跟踪3
  5. 【锐捷无线】胖AP模式限速
  6. CFA小课堂:CFA知识三大重点
  7. 模拟实现list和vector反向迭代器
  8. php 操作sqlite
  9. Google PageSpeed Insights 前端性能测试工具
  10. 信安Note_day27