高斯消元详解

算法概述

高斯消元法是一个用来求解线性方程组的算法

那么什么是线性方程组呢?

线性:每个未知数次数都为1次方程组:多个方程,多个未知数。

(a1x1+a2x2+..+anxn=bn)x为一次的

当x是平方的时候就不是线性

简而言之就是有多个未知数,并且每个未知数的次数均为一次,这样多个未知数组成的方程组为一

性方程组。或者我们也可以说是多元一次方程组

问题引入

给定一个线性方程组,对其求解

2x+y-z=8

-3x-y+2z=11

-2x+y+2z=-1

加减消元法和代入消元法

 高斯消元的目的是,先进行加减消元,然后我们可以先求得一个未知数的值,然后可以逐层往回代,也就是代入消元法,依次可以得到定第2个、第3个未知数的值,以至第n个未知数。最终解出方程组中各个未知数。

加减消元

对第一列进行操作(绝对值最大的一行,交换到第一行)

寻找最大值的原因是因为观察是否已经全是0了。

第二行减去若干倍的第一行,用第三行减去若干倍的第一行

高斯消元的代码详解

package gaosixiaoyuan;import java.util.*;
import java.util.BitSet;public class chapter1 {static final int N = 2800;static BitSet[] a = new BitSet[N];static int n, m, x, ans;public static void main(String[] args) {// TODO Auto-generated method stubScanner scanner = new Scanner(System.in);n = scanner.nextInt();for (int i = 1; i <= n; i++) {a[i] = new BitSet(n + 1);m = scanner.nextInt();if (m % 2 != 0) {a[i].set(n + 1);a[i].set(i);}while (m-- > 0) {x = scanner.nextInt();a[i].set(x);}}gauss();}public static void gauss() {int cnt = 0;for (int i = 1; i <= n; i++) {int maxx = cnt + 1;for (int j = i + 1; j <= n; j++) {if (a[j].get(i) && !a[maxx].get(i)) maxx = j;}BitSet temp = a[cnt + 1];a[cnt + 1] = a[maxx];a[maxx] = temp;if (!a[i].get(i)) continue;cnt++;for (int j = 1; j <= n; j++) {if (a[j].get(i) && i != j) a[j].xor(a[i]);}}if (cnt < n) {for (int i = 1; i <= n; i++) {if (!a[i].get(i) && a[i].get(n + 1)) System.out.println("no solution");}}}
}

 综上所述,可以大致分为三种情况:

1.高斯消元完成后,若存在系数全为0、常数不为0的行,则方程组无解。

2.若系数不全为0的行恰好有n个,则主元有n个,方程组有唯一解。

3.系数不为0的行<n个,则有无数个解。

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

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

相关文章

0基础学习Mybatis系列数据库操作框架——多环境配置

大纲 配置代码参考资料 在实际开发中&#xff0c;我们往往会将开发环境分成&#xff1a;开发、测试、线上等环境。这些环境的数据源不一样&#xff0c;比如开发环境就不能访问线上环境&#xff0c;否则极容易出现线上数据污染等问题。Mybatis通过多环境配置分开定义来解决这个问…

Shell脚本之基本语法

目录 一、变量定义 变量命名规则&#xff1a; 变量的赋值&#xff1a; 只读变量&#xff1a; 删除变量&#xff1a; 二、变量的类型 自定义变量&#xff1a; 环境变量&#xff1a; 位置参数&#xff1a; 预定义变量&#xff1a; 三、键盘输入 四、数值运算 为什么…

数据结构—堆

什么是堆 堆是一种特殊的树形结构&#xff0c;其中每个节点都有一个值。堆可以分为两种类型&#xff1a;最大堆和最小堆。在最大堆中&#xff0c;每个节点的值都大于等于其子节点的值&#xff1b;而在最小堆中&#xff0c;每个节点的值都小于等于其子节点的值。这种特性使得堆…

RPA自动化小红书自动化写文以及发文!

1、视频演示 RPA自动化小红书自动写作发文 2、核心功能点 采集笔记&#xff1a;采集小红书上点赞量大于1000的爆款笔记 下载素材&#xff1a;下载爆款笔记的主图 爆款改写&#xff1a;根据爆款笔记的标题仿写新的标题以及新的文案 自动发布&#xff1a;将爆款笔记发布到小红…

希尔排序算法(Java实现)

1.希尔排序&#xff08;Shell Sort&#xff09; &#xff08;1&#xff09;算法思想 先追求表中元素部分有序&#xff0c;再逐渐逼近全局有序。先将待排序表分割成若干形如 L [ i , i d , i 2 d , . . . , i k d ] L[i,id,i2d,...,ikd] L[i,id,i2d,...,ikd]的子表&#xff…

【docker】基础背景 Windows安装docker(WSL2)

docker解释 docker作用&#xff1a; 使用具体例子解释&#xff0c;在没有docker之前&#xff0c;整体项目由程序员在本地编写完成之后&#xff0c;再把执行文件丢给运维人员。但是有问题&#xff1a;因为运维同学的具体环境与编码人员的环境可能不一致&#xff0c;所以可能在环…

Three.js——scene场景、几何体位置旋转缩放、正射投影相机、透视投影相机

个人简介 &#x1f440;个人主页&#xff1a; 前端杂货铺 &#x1f64b;‍♂️学习方向&#xff1a; 主攻前端方向&#xff0c;正逐渐往全干发展 &#x1f4c3;个人状态&#xff1a; 研发工程师&#xff0c;现效力于中国工业软件事业 &#x1f680;人生格言&#xff1a; 积跬步…

第N6周:使用Word2vec实现文本分类

import torch import torch.nn as nn import torchvision from torchvision import transforms,datasets import os,PIL,pathlib,warnings #忽略警告信息 warnings.filterwarnings("ignore") # win10系统 device torch.device("cuda"if torch.cuda.is_ava…

DD3L内存芯片介绍

在数字科技迅猛发展的今天&#xff0c;内存芯片作为硬件的核心组件之一&#xff0c;扮演着至关重要的角色。而DD3L内存芯片以其卓越的性能和独特的设计&#xff0c;成为众多高端电子设备的不二选择。那么&#xff0c;DD3L内存芯片究竟如何应用&#xff0c;又是如何释放数字世界…

导入预览以及解决导入量大引发超时等问题

1、首选解决预览问题 由于使用的是vue3&#xff0c;页面与数据都是交互响应式的&#xff0c;所以可以通过组件或者原生的文件上传&#xff0c;获取到excel的sheet&#xff0c;从而来计算条数&#xff0c;页码&#xff0c;页数&#xff0c;手动实现分页逻辑&#xff0c;也就是把…

flex:1的作用是什么?

占满剩余的高度 <div classfather><div classson1></div><div classson2></div> </div>当给father添加display:flex之后&#xff0c;假设给son2添加flex:1&#xff0c;那么son2将会占满除son1之外的高度

【DevOps工具篇】Keycloak中设置LDAP认证

【DevOps工具篇】Keycloak中设置LDAP认证 目录 【DevOps工具篇】Keycloak中设置LDAP认证本次使用的环境服务器配置LDAP目录结构使用存储在LDAP中的用户进行登录Keycloak配置步骤功能测试从LDAP向Keycloak批量添加用户Keycloak配置步骤功能测试推荐超级课程: Docker快速入门到精…

基于springboot+vue+Mysql的招生管理系统

开发语言&#xff1a;Java框架&#xff1a;springbootJDK版本&#xff1a;JDK1.8服务器&#xff1a;tomcat7数据库&#xff1a;mysql 5.7&#xff08;一定要5.7版本&#xff09;数据库工具&#xff1a;Navicat11开发软件&#xff1a;eclipse/myeclipse/ideaMaven包&#xff1a;…

震惊!!原来阻塞队列消息队列这样理解会更简单!!!

震惊!!原来阻塞队列&&消息队列这样理解会更简单!!! 一:阻塞队列二:消息队列2.1:生产者消费者模型2.1.1:解耦合:2.1.2:削峰填谷: 三:消息队列代码3.1.13.1.2:3.1.3:生产慢,消费快,消费阻塞3.1.3:生产快,消费慢,生产阻塞 二级目录二级目录 一:阻塞队列 阻塞队列:先进先出…

gitcode 配置 SSH 公钥

在 gitcode 上配置SSH公钥后&#xff0c;可以通过SSH协议安全地访问远程仓库&#xff0c;无需每次都输入用户名和密码。以下是配置SSH公钥的步骤&#xff1a; 5分钟解决方案 用 OpenSSH公钥生成器 生成 公钥和私钥&#xff0c;私钥文件&#xff08;id_rsa&#xff09;下载&am…

解决Maven Clean过程因内存溢出导致的“Process terminated”问题

正文&#xff1a; 在今天的开发过程中&#xff0c;我遇到了一个意料之外的问题&#xff0c;当我尝试使用 Maven 进行项目清理&#xff08;maven clean&#xff09;时&#xff0c;命令行反馈了一个错误信息&#xff1a;“Process terminated”。经过一番排查&#xff0c;发现问…

为什么资讯网站选择高防IP防护攻击

为什么资讯网站选择高防IP防护攻击&#xff1f;在当今信息爆炸的时代&#xff0c;各种资讯网站扮演着重要的角色&#xff0c;为人们提供丰富的信息资源和资讯服务。然而&#xff0c;随着网络环境的不断变化和网络安全威胁的增加&#xff0c;资讯网站面临着越来越严重的网络攻击…

【Leetcode】top 100 图论

基础知识补充 1.图分为有向图和无向图&#xff0c;有权图和无权图&#xff1b; 2.图的表示方法&#xff1a;邻接矩阵适合表示稠密图&#xff0c;邻接表适合表示稀疏图&#xff1b; 邻接矩阵&#xff1a; 邻接表&#xff1a; 基础操作补充 1.邻接矩阵&#xff1a; class GraphAd…

JavaScript class继承

如果用新的class关键字来编写Student&#xff0c;可以这样写&#xff1a; class Student { constructor(name) { this.name name; } hello() {alert(Hello, this.name !); }} class的定义包含了构造函数constructor和定义在原型对象上的函数hello()&#xff08;注意没有f…

Open3D(C++) 鲁棒损失函数优化的ICP算法

目录 一、损失函数1、关于2、损失函数3、Open3D实现二、代码实现三、结果展示1、配准前1、配准后本文由CSDN点云侠原创,