判断并求解矩阵的逆——高斯-若尔单消元法

这篇博文主要讲解关于线性代数中矩阵的逆的那些事。

没看过前几篇博客的朋友先去看一看。

首先,方阵(n*n)才有单位矩阵的概念。所以下面提到的都会是矩阵都指方阵。

单位矩阵

矩阵的逆与单位矩阵的概念息息相关。

如单位长度,单位1,等等这些前缀是单位的概念们一样,矩阵里的单位矩阵同样特殊。
简而言之,任何矩阵 A 乘以单位矩阵 I 等于它自身。即
A I = A A I =A AI=A

下面的便是一个三乘三的单位矩阵:

1 0 0 0 1 0 0 0 1 \begin{matrix} 1 & 0 & 0\\ 0 & 1 & 0\\ 0 & 0 & 1 \end{matrix} \ 100010001 

有兴趣的小伙伴可以试试矩阵乘法,将单位矩阵与其他同样大小的方阵相乘,看是不是会得到同一个方阵。

总结,单位矩阵即为正对角线上元素为1,其余元素为0的方阵。


矩阵的逆

矩阵的逆的概念为,如果
A B = I ( 单位矩阵 ) AB=I(单位矩阵) AB=I(单位矩阵)
那么,矩阵B是矩阵A的逆。说实话,看到这个概念第一个想到的是倒数(ab=1),虽然书上没写。

有一个特殊的点是矩阵与其逆满足交换律,矩阵乘逆等于单位矩阵,逆乘矩阵也等于单位矩阵。


判断矩阵是否有逆

1 2 3 6 \begin{matrix} 1 & 2 \\ 3 & 6 \end{matrix} 1326
如何快速知道这个矩阵是否有逆呢?

法一 瞪眼 观察法

我们知道,矩阵右乘即为矩阵的列的线性组合

那么,在一些特殊情况下,我们很明显就能看出这个矩阵的列的线性组合能否得到单位矩阵。

比如上面的矩阵,很明显,因为两列成倍数,或者说在同一条直线上(这一点值得细品),他们的线性组合不能覆盖整个空间,而单位矩阵也不在他们所在的直线上,所以没有逆。

这个方法写在这里与其说是用来判断逆,不如说是为了更好帮助你理解线性组合这一概念,并引出了矩阵可能在图像上的这种思路,我想我们后面一定会用到的。

法二 利用单位矩阵

A B = I ( 单位矩阵 ) AB=I(单位矩阵) AB=I(单位矩阵)
首先,矩阵有逆说明存在矩阵B,能使A满足这个关系式。

那么,我们可以推出这个结论:

如果存在 A X = 0 , X 不为非零矩阵 , 那么矩阵 A 没有逆 如果存在AX=0,X不为非零矩阵,那么矩阵A没有逆 如果存在AX=0,X不为非零矩阵,那么矩阵A没有逆

推导逻辑如下:

A X = 0 , 两边同时乘 B → A B X = 0 → I X = 0 → X = 0 AX=0, \\两边同时乘B\rightarrow ABX=0\rightarrow IX=0\rightarrow X=0 AX=0,两边同时乘BABX=0IX=0X=0

可以看到,矩阵A存在逆和AX=0这两个命题矛盾,因此如果能找到矩阵X使AX=0便能证明矩阵A没有逆。


求解矩阵的逆——高斯-若尔单消元法

1 3 2 7 \begin{matrix} 1 & 3 \\ 2 & 7 \end{matrix} 1237

矩阵A如上。

已知矩阵A有逆,要求矩阵的逆,我们第一次想到的应该是如下方法:

[ 1 3 2 7 ] [ a b c d ] = [ 1 0 0 1 ] (1) \left[ \begin{matrix} 1 & 3 \\ 2 & 7 \end{matrix} \right] \left[ \begin{matrix} a & b \\ c & d \end{matrix} \right] = \left[\begin{matrix} 1 & 0 \\ 0 & 1 \end{matrix} \right] \tag{1} [1237][acbd]=[1001](1)


简言之,解方程,根据这个等式能知道如下方程,原理讲解见矩阵乘法

a + 3 c = 1 2 a + 7 c = 0 b + 3 d = 0 2 b + 7 d = 1 a+3c=1\\2a+7c=0\\b+3d=0\\2b+7d=1 a+3c=12a+7c=0b+3d=02b+7d=1

然而,这种解法的确有些繁琐。于是,我们便引入上篇文章所提的矩阵消元法

先将方程写为如下增广炬阵:

∣ 1 3 1 0 2 7 0 1 ∣ \left| \begin{array}{lc|cc} {1}&{3}&{1}&{0}\\ {2}&{7}&{0}&{1} \end{array} \right| 12371001

进行消元处理后得到如下矩阵:

∣ 1 3 1 0 0 1 − 2 1 ∣ \left| \begin{array}{lc|cc} {1}&{3}&{1}&{0}\\ {0}&{1}&{-2}&{1} \end{array} \right| 10311201

这时,高斯停下了脚步,在这里开始进行消元,回到解方程组的步骤。


然而,后来出现了个名为若尔丹的数学家,说,这里还可以再进一步:

若尔单在高斯停下脚步上的矩阵继续操作,对第一行也进行消元,得到如下矩阵:

∣ 1 0 7 − 3 0 1 − 2 1 ∣ \left| \begin{array}{lc|cc} {1}&{0}&{7}&{-3}\\ {0}&{1}&{-2}&{1} \end{array} \right| 10017231

这时,右边的矩阵就变成了我们所求的矩阵A的逆。

脑袋是不是没转过来?有兴趣的小伙伴可以将这个矩阵乘乘上原矩阵A看是否等于单位矩阵看看。


经验证,我们发现右边的矩阵的确是单位矩阵,下面我们来分析为什么。

根据上篇博文对消元矩阵的解释,消元操作等同于矩阵乘一个矩阵E1,再进行一次消元矩阵将左边部分变为单位矩阵,写成等式如下:

A E 1 E 2 = I A E_1E_2=I AE1E2=I

同样的,因为增广炬阵右边进行了同样的操作,可得:

I E 1 E 2 = E 1 E 2 IE_1E_2=E_1E_2 IE1E2=E1E2

此时,观察这两个等式,我想,不难发现,E1E2即为矩阵的逆。这就是高斯-若尔单消元法。

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

如若内容造成侵权/违法违规/事实不符,请联系一条长河网进行投诉反馈,一经查实,立即删除!

相关文章

【AI人工智能】文心智能体,陪爸妈去旅游,国庆假期不容错过,旅游搭子首选

文章目录 背景创作灵感陪爸妈去旅游简介角色与目标思考路径个性化开场白调优 智能体体验总结和感受 背景 文心智能体平台,开启新一轮活动,超级创造营持续百日活动。 在AI 浪潮席卷的今天,如雨后春笋般丛生的 AI 应用,昭告着时代风…

MySQL优化实战 解决CPU100%

问题表象 在24年初有一个日经问题困扰着我们,每到正点03分DB的CPU开始打满,持续1分钟又恢复正常水平。但由于日常业务交付压力较大且权限限制没有登录DB主机的权限,大家也就得过且过一直没有去认真排查。直到某天我来兴趣了也有时间了&#…

基于51单片机的家用防火防盗控制系统设计

本设计基于51单片机的家用防火防盗控制系统,该系统通过模块间的协同作用实现了对烟雾与天然气浓度的监测、温度监测、人体红外监测、通信传输、声光报警等功能。利用按键模块设置报警的阈值,将处理后的信息与阈值进行对比。判断气体浓度和温度是否超过阈…

酒店智能门锁SDK接口pro[V10] 门锁校验C#-SAAS本地化-未来之窗行业应用跨平台架构

一、代码 int 酒店标识_int Convert.ToInt32(酒店标识);StringBuilder 锁号2024 new StringBuilder(8);//信息 "未知返回值:" bufCard_原始;GetGuestLockNoByCardDataStr_原始(酒店标识_int, bufCard_原始.ToString(), 锁号2024);StringBuilder 退…

C++语言学习(4): identifier 的概念

1. 什么是 identifier identifier 中文意思是标识符,在 cppreference 中明确提到,identifier 是任意长度的数字、下划线、大写字母、小写字母、unicode 字符 的序列: An identifier is an arbitrarily long sequence of digits, underscores…

nginx打包部署前端vue项目全过程【保姆级教程】

🤹‍♀️潜意识起点:个人主页 🎙座右铭:得之坦然,失之淡然。 💎擅长领域:前端 是的,我需要您的: 🧡点赞❤️关注💙收藏💛 是我持…

Python字符串string方法大全及使用方法[1]以及FastAPI框架文件上传的处理-client使用postman

一、Python字符串string方法大全及使用方法[1] 1. Python字符串string方法大全及意义解释 #将字符串的第一个字符转换为大写 capitalize() #返回一个指定的宽度 width 居中的字符串,fillchar 为填充的字符,默认为空格。 center(width, fillchar) #返…

磁编码器磁铁要求和安装要求

总结来说, 磁铁需要是径向两极充磁、牌号N35、直径10mm、高度2.5mm的烧结钕铁硼磁铁。 磁铁的固定套必须是非导磁材料,比如铜、铝、塑料。 磁铁要距离电机轴至少2mm以上。 磁铁距离磁编码芯片0.5~3mm,最好1到2mm。 使用磁编码器的伺服&#…

css设置文本样式属性

目录 1.font-size:字体大小 案例:通过font-size属性设置字体的大小 1.代码 2.效果 2.font-family:字体的展现形式 案例:使用font-family属性设置字体的风格 1.代码实现 2.效果 3. font-weight:字体的粗细 案例:使用font-weight定义…

基于SSM+Vue技术的定制式音乐资讯平台

文未可获取一份本项目的java源码和数据库参考。 一、选题的背景与意义: 随着个人计算机的普及和互联网技术的日渐成熟,网络正逐渐成为人们获取信息及消费的主要渠道。然而在当前这个信息时代,网络中的信息种类和数量呈现爆炸性增长的趋势&a…

爱拼才会赢,甲骨文公司智算中心标配英伟达GPU10万颗

【科技明说 | 科技热点关注】 之前有有外媒消息,甲骨文宣布推出了多款智算集群,可通过甲骨文云基础设施提供AI训练服务,其中最顶级的一款配备了超过10万块的NVIDIA Blackwell GPU。 它一共使用了多达131072万块B200 GPU加速卡&…

单链表及其代码实现

目录 前言单链表1.1 单链表的定义1.2单链表代码实现1.2.1 头文件1.2.2 函数实现文件1.2.3 测试文件1.2.4 野指针问题 总结 前言 本文介绍单链表,主要是创销、增删改查代码实现。 注:文章中函数命名采取STL库。 单链表 1.1 单链表的定义 单链表是链线…

北京市大兴区启动乐享生活 寻味大兴 美食嘉年华 系列促销费活动

北京市大兴区启动乐享生活 寻味大兴 系列促销费活动 区商务局副局长 兰莉 致开幕辞 区餐饮行业协会会长 董志明 介绍活动内容 2024年9月30日,由大兴区商务局主办、大兴区餐饮行业协会承办,并得到高米店街道和大兴绿地缤纷城大力支持的“乐享生活 寻味大…

OceanBase—02(入门篇——对于单副本单节点,由1个observer扩容为3个observer集群)——之前的记录,当初有的问题未解决,目前新版未尝试

OceanBase—02(入门篇——对于单副本单节点,由1个observer扩容为3个observer集群)——之前的记录,有的问题未解决,新版未尝试 1、前言—安装单副本单节点集群1.1 docker安装OB 2、查看现有集群情况2.1 进入容器&#x…

SOMEIP_ETS_147: SD_Send_triggerEventUINT8_Eventgroup_2

测试目的: 验证DUT在Tester订阅事件组后,能够响应Tester触发的triggerEventUINT8方法,并将TestEventUINT8事件发送到订阅请求中端点选项指定的IP地址和端口。 描述 本测试用例旨在确保DUT能够正确处理事件组的订阅请求,并且在T…

VSOMEIP代码阅读整理(1) - 网卡状态监听

一. 概述 ​ 在routing进程所使用的配置文件中,存在如下配置项目: {"unicast" : "192.168.56.101",..."service-discovery" :{"enable" : "true","multicast" : "224.244.224.245&q…

在2核2G服务器安装部署MySQL数据库可以稳定运行吗?

阿里云2核2G服务器可以安装MySQL数据库吗?当然可以,并且可以稳定运行MySQL数据库,目前阿里云服务器网aliyunfuwuqi.com使用的就是阿里云2核2G服务器,在云服务器上安装MySQL数据库,可以稳定运行。 目前阿腾云用于运行M…

C++系列-继承补充

🌈个人主页:羽晨同学 💫个人格言:“成为自己未来的主人~” 继承和友元 友元关系不能继承,父亲的朋友不能是你的朋友 比如在这个例子当中: class Student; class Person { public:friend void Display(const Per…

厦门网站设计的用户体验优化策略

厦门网站设计的用户体验优化策略 在信息化快速发展的今天,网站作为企业与用户沟通的重要桥梁,用户体验(UX)的优化显得尤为重要。尤其是在交通便利、旅游资源丰富的厦门,吸引了大量企业进驻。在这样竞争激烈的环境中&am…