当前位置: 首页 > news >正文

线代——求逆矩阵的快捷方法

通常,求逆矩阵有两种方法:
方法一:
在这里插入图片描述
方法二:

在这里插入图片描述
但是,对于特殊矩阵,如:

1、二阶矩阵

A = [ a b c d ] A = \begin{bmatrix} a & b\\ c & d \end{bmatrix} A=[acbd],其逆矩阵 A − 1 = 1 a d − b c [ d − b − c a ] A^{-1}=\frac{1}{ad-bc}\begin{bmatrix} d & -b\\ -c & a \end{bmatrix} A1=adbc1[dcba]

2、分块矩阵


分块矩阵在主对角位置,直接对分块矩阵取逆矩阵:

A = [ X Y ] A = \begin{bmatrix} X & \\ & Y \end{bmatrix} A=[XY],其逆矩阵 A − 1 = [ X − 1 Y − 1 ] A^{-1}=\begin{bmatrix} X^{-1} & \\ & Y^{-1} \end{bmatrix} A1=[X1Y1]


分块矩阵在副对角位置,副对调,再取逆

A = [ X Y ] A = \begin{bmatrix} & X \\ Y & \end{bmatrix} A=[YX],其逆矩阵 A − 1 = [ Y − 1 X − 1 ] A^{-1}=\begin{bmatrix} & Y^{-1}\\ X^{-1} & \end{bmatrix} A1=[X1Y1]


分块矩阵为右上三角形状,首先主对角直接取逆,然后再对右上角子矩阵左乘其行,右乘其列,再添符号

A = [ X W Y ] A = \begin{bmatrix} X & W\\ & Y \end{bmatrix} A=[XWY],其逆矩阵 A − 1 = [ X − 1 − X − 1 W Y − 1 Y − 1 ] A^{-1}=\begin{bmatrix} X^{-1} & -X^{-1}WY^{-1}\\ & Y^{-1} \end{bmatrix} A1=[X1X1WY1Y1]

同理,对于左下三角形状,首先主对角直接取逆,然后再对左下角子矩阵左乘其行,右乘其列,再添符号

A = [ X W Y ] A = \begin{bmatrix} X & \\ W & Y \end{bmatrix} A=[XWY],其逆矩阵 A − 1 = [ X − 1 − Y − 1 W X − 1 Y − 1 ] A^{-1}=\begin{bmatrix} X^{-1} & \\ -Y^{-1}WX^{-1} & Y^{-1} \end{bmatrix} A1=[X1Y1WX1Y1]

它们相同之处,都是分块三角矩阵占据主对角位置。


分块矩阵为左上三角形状,首先副对调,再取逆,然后将左上角子矩阵换到右下角,最后再对该子矩阵左乘其行,右乘其列,再添符号

A = [ W X Y ] A = \begin{bmatrix} W & X \\ Y & \end{bmatrix} A=[WYX],其逆矩阵 A − 1 = [ Y − 1 X − 1 − X − 1 W Y − 1 ] A^{-1}=\begin{bmatrix} & Y^{-1}\\ X^{-1} & -X^{-1}WY^{-1} \end{bmatrix} A1=[X1Y1X1WY1]

同理,对于右下三角形状,首先副对调,再取逆,然后将右下角子矩阵换到左上角,最后再对该子矩阵左乘其行,右乘其列,再添符号

A = [ X Y W ] A = \begin{bmatrix} & X \\ Y & W \end{bmatrix} A=[YXW],其逆矩阵 A − 1 = [ − Y − 1 W X − 1 Y − 1 X − 1 ] A^{-1}=\begin{bmatrix} -Y^{-1}WX^{-1} & Y^{-1}\\ X^{-1} & \end{bmatrix} A1=[Y1WX1X1Y1]

它们相同之处,都是分块三角矩阵占据副对角位置。


综上,对于形状是上、下三角的分块矩阵求逆,如果分块子矩阵占据主对角位置,不需要对调位置;如果分块子矩阵占据副对角位置,都需要对调位置。

相关文章:

  • 分布式存储之 etcd 的集群管理
  • std::unique_ptr(基础和仿写)
  • UG旋钮盖模具设计与加工
  • 【JAVA高级】——封装JDBC中的DaoUtils工具类(Object类型方法)
  • MYSQL事务原理分析
  • ORB-SLAM3算法学习—Frame构造—基于SAD滑窗的双目特征匹配
  • 【最佳实践】gorm 联表查询 joins
  • http请求走私漏洞原理,利用,检测,防护
  • simulink中比scope模块还好用的平替出图工具?
  • 微信小程序和H5之间互相跳转、互相传值
  • web前端期末大作业 html+css家乡旅游主题网页设计 湖北武汉家乡介绍网页设计实例
  • PS-HDR图像编辑与应用
  • 2022亚太杯C题思路代码分析
  • 涨知识!Python 的异常信息还能这样展现
  • JDBC编程
  • __proto__ 和 prototype的关系
  • C# 免费离线人脸识别 2.0 Demo
  • CentOS6 编译安装 redis-3.2.3
  • Electron入门介绍
  • js ES6 求数组的交集,并集,还有差集
  • Node项目之评分系统(二)- 数据库设计
  • php的插入排序,通过双层for循环
  • socket.io+express实现聊天室的思考(三)
  • Yii源码解读-服务定位器(Service Locator)
  • 搭建gitbook 和 访问权限认证
  • 两列自适应布局方案整理
  • 前端每日实战:61# 视频演示如何用纯 CSS 创作一只咖啡壶
  • 如何实现 font-size 的响应式
  • 什么软件可以剪辑音乐?
  • 它承受着该等级不该有的简单, leetcode 564 寻找最近的回文数
  • 微服务入门【系列视频课程】
  • 以太坊客户端Geth命令参数详解
  • 白色的风信子
  • ### Error querying database. Cause: com.mysql.jdbc.exceptions.jdbc4.CommunicationsException
  • #android不同版本废弃api,新api。
  • #单片机(TB6600驱动42步进电机)
  • #调用传感器数据_Flink使用函数之监控传感器温度上升提醒
  • #我与Java虚拟机的故事#连载04:一本让自己没面子的书
  • (1综述)从零开始的嵌入式图像图像处理(PI+QT+OpenCV)实战演练
  • (BFS)hdoj2377-Bus Pass
  • (C语言)输入一个序列,判断是否为奇偶交叉数
  • (附源码)ssm失物招领系统 毕业设计 182317
  • (一)spring cloud微服务分布式云架构 - Spring Cloud简介
  • (转)IIS6 ASP 0251超过响应缓冲区限制错误的解决方法
  • (转)visual stdio 书签功能介绍
  • (轉貼) 2008 Altera 亞洲創新大賽 台灣學生成果傲視全球 [照片花絮] (SOC) (News)
  • *1 计算机基础和操作系统基础及几大协议
  • .apk 成为历史!
  • .NET 使用 ILRepack 合并多个程序集(替代 ILMerge),避免引入额外的依赖
  • .NET正则基础之——正则委托
  • :中兴通讯为何成功
  • [ 云计算 | AWS ] AI 编程助手新势力 Amazon CodeWhisperer:优势功能及实用技巧
  • [20180312]进程管理其中的SQL Server进程占用内存远远大于SQL server内部统计出来的内存...
  • [AIGC 大数据基础]hive浅谈
  • [CakePHP] 在Controller中使用Helper