【数值计算方法】Gauss消元法及其Python/C实现

文章目录

  • 一、基础理论
    • 1. 线性方程组
    • 2. Gauss消元法的详细步骤
    • 3. 注意事项
  • 二、具体计算过程
    • 1. 用Gauss 消元法求A的LU分解,并由此求解方程组 Ax =b
      • a. 将A进行LU分解。
      • b. 使用LU分解求解方程组Ax=b
  • 三、代码实现
    • 1. Python代码实现
    • 2. C语言代码实现

  Gauss消元法,也称为高斯消元法或高斯-约当消元法,是一种用于求解线性方程组的数值方法。它是由德国数学家卡尔·弗里德里希·高斯在18世纪末发展起来的。

  Gauss消元法的基本思想是通过一系列的行变换将线性方程组转化为一个上三角形的方程组,然后通过回代过程求解方程组的解。

一、基础理论

1. 线性方程组

  线性方程组是由一组线性方程组成的方程集合。每个线性方程都可以表示为形如 “a₁x₁ + a₂x₂ + … + aₙxₙ = b” 的形式,其中 a₁, a₂, …, aₙ 是已知的常数,x₁, x₂, …, xₙ 是未知的变量,b 是已知的常数。方程中的每一项都是变量的一次幂与常数的乘积,且没有乘法运算符连接变量。
  线性方程组可以包含多个线性方程,这些方程共同描述了一组变量的关系。解线性方程组就是找到满足所有方程的变量值,使得所有方程都成立。解线性方程组的目标是找到一组变量的值,使得方程组中的每个方程都得到满足。
  线性方程组的解可以有多个或者没有解。如果存在至少一个满足所有方程的变量值组合,那么方程组有解。如果不存在这样的变量值组合,那么方程组无解。
  解线性方程组的方法包括高斯消元法、矩阵法、克莱姆法则等。这些方法可以用于求解不同规模和形式的线性方程组。线性方程组在数学、物理、工程等领域中广泛应用,用于描述和解决各种实际问题。

2. Gauss消元法的详细步骤

  1. 将线性方程组写成增广矩阵的形式,即将系数矩阵和常数向量合并在一起。
  2. 选取第一个未知数的系数不为零的方程作为主元方程,如果没有这样的方程,则交换两行或者两列,使得主元系数不为零。
  3. 将主元方程的系数除以主元系数,使主元系数变为1。
  4. 用主元方程的系数乘以其他方程的主元系数,并将得到的结果从相应的方程中减去,以消除其他方程中的主元系数。
  5. 重复步骤2到步骤4,直到所有的未知数的系数都变为上三角形矩阵的形式。
  6. 进行回代过程,从最后一行开始,依次求解每个未知数的值。回代的过程是通过将已知的未知数代入到方程中,求解出未知数的值。

3. 注意事项

&emps;&emps;Gauss消元法的优点是可以精确地求解线性方程组,适用于任意个数的未知数和方程,在数值计算和科学工程领域有广泛的应用,然而,它也有一些限制和注意事项:

  1. 如果方程组的系数矩阵是奇异的(即行列式为零),则无法使用Gauss消元法求解。

  2. 在进行消元过程中,需要注意避免除以零的情况,如果遇到主元系数为零的情况,需要进行行交换或列交换。

  3. 如果方程组的系数矩阵很大,消元的计算量会很大,可能需要较长的计算时间。

二、具体计算过程

1. 用Gauss 消元法求A的LU分解,并由此求解方程组 Ax =b

A = [ [ 1 , 2 , 1 , − 2 ] , [ 2 , 5 , 3 , − 2 ] , [ − 2 , − 2 , 3 , 5 ] , [ 1 , 3 , 2 , 3 ] ] A=[ [1, 2, 1, -2], [2, 5, 3, -2], [-2, -2, 3, 5], [1, 3, 2, 3] ] A=[[1,2,1,2],[2,5,3,2],[2,2,3,5],[1,3,2,3]]

b = [ 2 , 8 , 4 , 9 ] b=[ 2, 8, 4, 9 ] b=[2,8,4,9]

a. 将A进行LU分解。

  1. 选取第一个未知数的系数不为零的方程作为主元方程,即第1行第1列元素不为零,因此选择第1行为主元方程。

  2. 将主元方程的系数除以主元系数,即第1行的所有元素除以1,得到:

1   2   1  -2
2   5   3  -2
-2 -2   3   5
1   3   2   3
  1. 用主元方程的系数乘以其他方程的主元系数,并将得到的结果从相应的方程中减去,以消除其他方程中的主元系数。对第2行、第3行和第4行进行消元操作:
1   2   1  -2
0   1   1   2
0   4   4   1
0   1   1   5
  1. 选择第二个未知数的系数不为零的方程作为主元方程,即第2行第2列元素不为零,因此选择第2行为主元方程。

  2. 将主元方程的系数除以主元系数,即第2行的所有元素除以1,得到:

1   2   1  -2
0   1   1   2
0   4   4   1
0   1   1   5
  1. 用主元方程的系数乘以其他方程的主元系数,并将得到的结果从相应的方程中减去,以消除其他方程中的主元系数。对第3行和第4行进行消元操作:
1   2   1  -2
0   1   1   2
0   0   0  -7
0   0   0   3

现在,我们得到了上三角形矩阵U和下三角形矩阵L:

U = 
1   2   1  -2
0   1   1   2
0   0   0  -7
0   0   0   3L = 
1   0   0   0
2   1   0   0
-2  -4  1   0
1   -1  -1  1

b. 使用LU分解求解方程组Ax=b

  1. 首先,根据LU分解,我们可以得到Ly=b,其中y是一个新的未知向量。
1   0   0   0   |  y1  =  2
2   1   0   0   |  y2  =  8
-2  -4  1   0   |  y3  =  4
1   -1  -1  1   |  y4  =  9

通过前向代入法,我们可以求解出y的值:

y1 = 2
y2 = 8 - 2y1 = 8 - 2(2) = 4
y3 = 4 - 2y1 + 4y2 = 4 - 2(2) + 4(4) = 18
y4 = 9 - y1 + y2 - y3 = 9 - 2 + 4 - 18 = -7
  1. 然后,根据LU分解,我们可以得到Ux=y,其中x是我们要求解的未知向量。
1   2   1  -2   |  x1  =  y1
0   1   1   2   |  x2  =  y2
0   0   0  -7   |  x3  =  y3
0   0   0   3   |  x4  =  y4

通过回代法,我们可以求解出x的值:

x1 = y1 = 2
x2 = y2 - x1 = 4 - 2 = 2
x3 = y3 / (-7) = 18 / (-7)-2.571
x4 = y4 / 3 = (-7) / 3-2.333

  因此,方程组Ax=b的解为x = [2, 2, -2.571, -2.333]。

三、代码实现

1. Python代码实现

import numpy as npA = np.array([[1, 2, 1, -2],[2, 5, 3, -2],[-2, -2, 3, 5],[1, 3, 2, 3]])b = np.array([2, 8, 4, 9])def gauss_elimination(A, b):n = len(A)for i in range(n-1):for j in range(i+1, n):factor = A[j, i] / A[i, i]A[j, i:] -= factor * A[i, i:]b[j] -= factor * b[i]return A, bdef back_substitution(U, y):n = len(U)x = np.zeros(n)x[-1] = y[-1] / U[-1, -1]for i in range(n-2, -1, -1):x[i] = (y[i] - np.dot(U[i, i+1:], x[i+1:])) / U[i, i]return xdef solve_linear_equations(A, b):U, y = gauss_elimination(A, b)x = back_substitution(U, y)return xx = solve_linear_equations(A, b)
print("Solution x:", x)

2. C语言代码实现

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/125652.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

基于单片机设计的防煤气泄漏装置

一、前言 煤气泄漏是一个严重的安全隐患,可能导致火灾、爆炸以及对人体健康的威胁。为了提高家庭和工业环境中煤气泄漏的检测和预防能力,设计了一种基于单片机的防煤气泄漏装置。 单片机选择STC89C52作为主控芯片。为了检测煤气泄漏,采用了…

cocosCreator 调用wxAPI 及后台授权设置、获取用户昵称和头像

版本: 3.8.0 语言: TypeScript 环境: Mac 官方文档: 微信官方文档 - 开放能力 微信 API 小游戏环境 在cocosCreator的3.x版本项目开发中,TypeScript最终会被转换为JavaScript语言。 JavaScript的运行时调用的API…

Linux学习笔记1-入门

前言:之前的基于单片机的闭环控制步进电机项目其实已经完成了,但很多时间都花在调试和生产上,实在没时间去做总结笔记,现在又开始做新项目了,从单片机到了Linux,想用这个平台来督促自己继续学习&#xff0c…

联发科MT6893(天玑1200)_MTK5G芯片规格参数性能_安卓手机主板方案

联发科天玑1200集成MediaTek 5G调制解调器,通过包含6大维度、72个场景测试的德国莱茵TV Rheinland认证,支持高性能5G连接,带给用户全场景的高品质5G连网体验。 进入5G时代,AI多媒体成为主流应用,天玑1200以强劲的平台…

【高效开发工具系列】你真的会使用Mac吗?

💝💝💝欢迎来到我的博客,很高兴能够在这里和您见面!希望您在这里可以感受到一份轻松愉快的氛围,不仅可以获得有趣的内容和知识,也可以畅所欲言、分享您的想法和见解。 推荐:kwan 的首页,持续学…

ROS中/odom,/map,/base_link几个坐标系的含义:从ROS的REP105严谨介绍

这两日看lego-loam,其中的坐标系定义很混乱。 在lego-loam的一个issue中,作者指出他坐标系定义是完全按照ROS的REP103和105的建议的。 然后又看到了关于/odom, /map, /base_link 几个坐标系的含义的说明,因此整理一下。 关于ROS的REP 关于R…

uniapp实现瀑布流

首先我们要先了解什么是瀑布流: 瀑布流(Waterfall Flow)是一种常见的网页布局方式,也被称为瀑布式布局或砌砖式布局。它通常用于展示图片、博客文章、商品等多个不同大小和高度的元素。 瀑布流布局的特点是每个元素按照从上到下…

现代的简洁,诠释轻奢的精致!福州中宅装饰,福州装修

轻奢风是一种生活新时尚 优雅、低调、舒适、简单,不断地推陈出新 站在时尚的前沿,引领潮流 中宅装饰集团轻奢风格产品 追求高品质生活细节 以设计精致的空间构造营造出 一种优雅、时尚生活氛围 将低调奢华之美注入现代家居设计中 客厅|The Sitt…

世界电信日 | 人大金仓助力中国移动租赁核算系统升级上线

世界电信日 5月17日恰逢第五十四个世界电信日,运营商作为新型基础设施建设以及维护网信安全的主力军,掌握关键核心技术,实现科技自立自强刻不容缓。 作为数据库领域国家队,人大金仓坚持原始创新,低难度、低成本、低风…

早安心语微语早读,保持一颗平常心,坐看云起落花开谢得之淡然

1、保持一颗平常心,坐看云起落花开谢得之淡然,失之坦然,让生命中每一天都充满着阳光和希望! 2、每个人都一样,都有一段独行的日子,或长或短,这都是无可回避的。不必总觉得生命空空荡荡&#xf…

react动态插入样式

在开发组件过程中,偶尔需要动态的插入css,比如在在iframe中渲染组件后,iframe中是没有样式的,所以需要手动插入样式。 插入样式 通常是在useLayoutEffect中动态创建style标签 useLayoutEffect(() > {if (!ref.current) {cons…

ArcGIS计算土地现状容积率

本文讲解在ArcGIS中,基于建筑数据和地籍边界数据,计算土地容积率。 一、容积率介绍 容积率(Plot Ratio/Floor Area Ratio/Volume Fraction)是指一个小区的地上建筑总面积与净用地面积的比率。又称建筑面积毛密度。 二、数据分析 (1)建筑数据(dwg) (2)地籍边界数据…

C语言 写一个简易音乐播放器

#include <stdio.h> #include <stdlib.h> #include <stdbool.h> #include <math.h>#define SAMPLE_RATE 44100 // 采样率 #define AMPLITUDE 32767 // 振幅 #define NO_SAMPLES 44100 // 样本数// 声明一个结构体用于表示音符 typedef struct {double …

[USACO23OPEN] Field Day S题解

远古的回忆。 把变换一个字符视为边权为 1 1 1 的边&#xff0c;即求最长路。 最长路不好搞&#xff0c;考虑转补集最短路&#xff08;容易感性理解&#xff09;&#xff0c;BFS 即可。 #include<bits/stdc.h> #define int long long using namespace std;const int …

利用AI Chat 将电子书自动截屏并保存成pdf文件

电子书如果要下载下来&#xff0c;无非就两种类型的方法&#xff0c;一种是从内部破解&#xff0c;通常是某些极客将软件破解成免费版&#xff0c;但是风险也大。另一种是从外部破解&#xff0c;就是截屏保存&#xff0c;然后将所有图片拼成pdf文件。 如果要将整本电子书截屏保…

LeetCode刷题---简单组(六)

文章目录 &#x1f352;题目一 69. x 的平方根&#x1f352;解法一&#x1f352;解法二&#x1f352;题目二 70. 爬楼梯&#x1f352;解法一 &#x1f352;题目一 69. x 的平方根 &#x1f352;解法一 class Solution(object):def mySqrt(self, x):""":type x:…

制作一个简单的C语言词法分析程序

1.分析组成 C语言的程序中&#xff0c;有很单词多符号和保留字。一些单词符号还有对应的左线性文法。所以我们需要先做出一个单词字符表&#xff0c;给出对应的识别码&#xff0c;然后跟据对应的表格来写出程序 2.程序设计 程序主要有循环判断构成。不需推理即可产生的符号我…

SSM培训报名管理系统开发mysql数据库web结构java编程计算机网页源码eclipse项目

一、源码特点 SSM 培训报名管理系统是一套完善的信息系统&#xff0c;结合SSM框架完成本系统&#xff0c;对理解JSP java编程开发语言有帮助系统采用SSM框架&#xff08;MVC模式开发&#xff09;&#xff0c;系统具有完整的源代码和数据库&#xff0c;系统主 要采用B/S模式开…

Windows2008系统怎么隐藏或打开文件后缀

打开服务器的控制面板-选择小图标-文件夹选项 在文件夹选项那边点击查看-隐藏一直文件类型的扩展名 选择勾选&#xff08;隐藏一直文件类型的扩展名&#xff09;-下图示文件后缀不显示 选择不勾选&#xff08;隐藏一直文件类型的扩展名&#xff09;-下图示文件后缀显示

Go-ansible

文章目录 安装包AdhocPlaybookExecuteDefaultExecute自定义执行器 Optionsansible adhoc and ansible-playbook common options Stdout CallbackResultsTransformersDefaultJSONManage JSON output 示例 Go ansible是一个用于从Golang应用程序运行ansible剧本或ansible命令的包…