《 线性代数及其应用 (原书第4版)》—— 2.4 分块矩阵

简介:

本节书摘来自华章出版社《 线性代数及其应用 (原书第4版)》一书中的第2章,第2.4节,作者:(美)戴维C. 雷(David C. Lay)马里兰大学帕克学院 著刘深泉 张万芹 陈玉珍 包乐娥 陆 博 译,更多章节内容可以访问云栖社区“华章计算机”公众号查看

2.4 分块矩阵

我们既可以把矩阵看作一个数的矩形表,也可以把它看作一组列向量,后面这种看法起了很重要的作用,因而,我们想考虑A的其他分块,把它用水平线和竖直线分成几块,如下面例1所示. 分块矩阵也出现在线性代数的现代应用中,因为这些记号简化了许多讨论,并使矩阵计算中许多本质的结构显露出来. 本节也给出复习矩阵代数和可逆矩阵定理的机会.
例1 矩阵
screenshot
也可写成screenshot分块矩阵
screenshot
的形状,它的元素是分块(或子矩阵)
screenshot

例2 当某一矩阵A出现在物理问题的数学模型中时,例如,电子网络、传输系统或大公司等,会很自然地把A看作一个分块矩阵. 例如,若一个微型计算机电路板主要由3块超大规模的集成电路芯片组成,如图2-9所示,
screenshot
那么这电路板的矩阵可以写成一般形式
screenshot
A的“对角”线上的子矩阵,即 screenshot 是有关3块超大规模集成电路本身的矩阵,而其他子矩阵则与这三块芯片之间的相互联系有关.
加法与标量乘法
若矩阵A与B有相同维数且被同样地分块,则自然矩阵的和A + B也被同样地分块. 这时A + B的每一块恰好是A和B对应分块的(矩阵)和. 分块矩阵乘以一个数也可以逐块计算.
分块矩阵的乘法
分块矩阵也可用通常的行列法则进行,就如每一块都是数一样,只要A的列的分法与B的行的分法一致.
例3 设
screenshot
A的5列被分成3列一组和2列一组. B的5行按同样方法分块——被分成3行一组和2行一组. 我们称A和B的分块是与分块乘法相一致的. AB的乘积可以被写成
screenshot
重要的是,在AB的表达式中的小乘积,每一项应把来自A的子矩阵写在左边,因矩阵乘法是不可交换的. 例如
screenshot

因此AB的上面一块是
screenshot

分块矩阵乘法的行列法则给出了两个矩阵乘积的最一般观点. 下面每一个有关矩阵乘积的观点已经使用简单的矩阵分块的思想讨论过:(1)使用A的列来给出 Ax的定义;(2)AB的列的定义;(3)计算AB的行列法则;(4)A 的行与矩阵B的乘积作为AB的行. 在下面的定理10仍然应用分块的思想给出AB第5种观点.
下面的例子为定理10做准备. 符号screenshot表示A的第k列,screenshot表示B的第k行.
例4 设screenshotscreenshot. 证明
screenshot
解 上面的每一项都是外积(见2.1节习题27和28),由计算矩阵乘积的行列法则,有
screenshot
于是 screenshot
这个矩阵恰好就是AB. 注意AB的(1, 1)元素是3个外积的(1, 1)元素之和,AB的(1, 2)元素是3个外积的(1, 2)元素之和,等等.
定理10 (AB的列行展开)
若A是 screenshot矩阵,B是 screenshot矩阵,则
screenshot (1)
证 对每个行指标i 和列指标 j,乘积screenshot 的 (i,j)元素是screenshot 中元素screenshotscreenshot中元素 screenshot的积,因此在(1)的和中,(i,j) 元素为
screenshot
而根据行列法则,该和恰好是AB的 (i,j)元素.
分块矩阵的逆
下例说明分块矩阵的逆的求法.
例5 形如 screenshot的矩阵称为分块上三角矩阵,设 screenshotscreenshot 矩阵, screenshotscreenshot 矩阵,且A为可逆矩阵. 求 screenshot的表达式.
解 用B表示screenshot 且把它分块使
screenshot (2)
这个矩阵方程包含了4个有关未知子矩阵 的方程,计算(2)式左边的乘积得
screenshot (3)
screenshot (4)
screenshot (5)
screenshot (6)
方程(6)本身并不能说明screenshot 可逆,因我们还不知道 screenshot,但应用可逆矩阵定理,及screenshot 是方阵的事实,可以断定 screenshot为可逆且 screenshot. 现在我们利用(5)式求得
screenshot
因此(3)式简化为
screenshot
这说明screenshot 是可逆的,且 screenshot,最后由(4),
screenshot
于是
screenshot
分块对角矩阵是一个分块矩阵,除了主对角线上各分块外,其余全是零分块. 这样的一个矩阵是可逆的当且仅当主对角线上各分块都是可逆的. 见习题13和14.
数值计算的注解

  1. 当矩阵太大时,不适于存储在高速计算机内存中,分块矩阵允许计算机一次处理两到三块子矩阵,例如,最近关于线性规划的工作中,一个研究团队把矩阵分为837行和51列以简化问题. 这个问题的解在Cray超计算机上大约需要4分钟 .1
  2. 某些高速计算机,特别是具有向量传输技术的计算机,当把矩阵分块后再进行矩阵运算更有效 .1
  3. 高性能数值计算的线性代数专业软件LAPACK,广泛使用分块矩阵进行计算.
    下面的习题给出了运用矩阵代数的实践,表明了应用中的典型计算.

练习题

  1. 证明screenshot 可逆且求出它的逆.
  2. 计算screenshot ,其中X 分块为screenshot .
    习题2.4

screenshot
screenshot
screenshot
screenshot
screenshot
screenshot

相关文章
|
消息中间件 监控 Linux
手把手教你 centos 7 安装RabbitMQ
手把手教你 centos 7 安装RabbitMQ
1675 1
|
4月前
|
编译器 程序员 Linux
实用工具:VS Code 配置 Markdown 编译器全指南
本文介绍如何在VS Code中配置高效Markdown写作环境,通过启用内置预览、安装Markdown All in One与Markdown Preview Enhanced插件,并配置Princexml实现PDF等格式导出,全面提升编辑、预览与输出效率,适合程序员、学生及内容创作者使用。
2150 2
实用工具:VS Code 配置 Markdown 编译器全指南
|
1月前
|
人工智能 运维 安全
OpenClaw 帮我睡后全自动完成了老板交代的任务
安装 OpenClaw 后,我发现它最大的价值是 7x24 小时远程全自动完成任务。本文分享我的几个实际使用场景和 7 条实战经验。
851 4
|
3月前
|
人工智能 前端开发 开发工具
从 ReAct 到 Ralph Loop:AI Agent 的持续迭代范式
Ralph Loop是一种AI编程范式,通过Stop Hook拦截机制强制模型持续迭代,直至满足可验证的完成条件(如测试通过、输出特定标记),解决AI“半途而废”问题。它依赖外部状态(Git、文件系统)实现自我参照循环,支持无人值守的AFK开发,已集成于Claude Code等主流工具。
从 ReAct 到 Ralph Loop:AI Agent 的持续迭代范式
|
存储 安全 Linux
【实战指南】7个设置/获取接口了解Linux时间管理
本文系统介绍了Linux时间管理中的7个关键设置/获取接口,涵盖时间获取(如`time`、`gettimeofday`、`clock_gettime`)、时间设置(如`stime`、`settimeofday`、`clock_settime`)以及时间转换和格式化等内容。文章详细解析了绝对时间和相对时间的概念,包括GMT、UTC及本地时间的区别,并通过实例测试展示了各接口的使用方法与特性。此外,还探讨了时区设置对时间计算的影响,强调在实际开发中推荐使用UTC作为基准时间以避免时区变化带来的问题。总结部分结合项目经验,提醒开发者注意时间服务的重要性及潜在风险,例如时间跳跃可能引发的应用故障。
847 113
【实战指南】7个设置/获取接口了解Linux时间管理
|
搜索推荐 NoSQL Java
微服务架构设计与实践:用Spring Cloud实现抖音的推荐系统
本文基于Spring Cloud实现了一个简化的抖音推荐系统,涵盖用户行为管理、视频资源管理、个性化推荐和实时数据处理四大核心功能。通过Eureka进行服务注册与发现,使用Feign实现服务间调用,并借助Redis缓存用户画像,Kafka传递用户行为数据。文章详细介绍了项目搭建、服务创建及配置过程,包括用户服务、视频服务、推荐服务和数据处理服务的开发步骤。最后,通过业务测试验证了系统的功能,并引入Resilience4j实现服务降级,确保系统在部分服务故障时仍能正常运行。此示例旨在帮助读者理解微服务架构的设计思路与实践方法。
877 17
|
10月前
|
编解码 搜索推荐
mp4视频压缩小软件,怎么让mp4视频压缩,小丸视频压缩,小丸b站压制视频,b站视频压制小丸绿色免费版
小丸工具箱是一款傻瓜型视频压制工具,支持H264+AAC视频编码,内核采用x264、neroaac、mp4box等开源软件。它集成了强大的编码压缩技术(如x264和x265)、多格式兼容(MP4、MKV、AVI等)、音频处理(AAC编码)以及批量处理功能,同时提供丰富的参数自定义选项。用户可通过简单步骤完成视频导入、参数设置、压缩处理及输出查看,适合不同水平的视频处理需求。
1500 11
|
12月前
|
Java 开发者 微服务
Spring Cloud OpenFeign详解与实践
总结起来说,Spring Cloud OpenFeign提供了一种简单易懂且高效的方式去实现微服务之间通信.它隐藏了许多复杂性,并且允许开发者以声明式方式编写HTTP客户端代码.如果你正在开发基于Spring Cloud 的微服务架构系统,Spring Cloud Open Feign是一个非常好用且强大工具.
793 33
|
关系型数据库 MySQL Linux
Docker安装Mysql5.7,解决无法访问DockerHub问题
当 Docker Hub 无法访问时,可以通过配置国内镜像加速来解决应用安装失败和镜像拉取超时的问题。本文介绍了如何在 CentOS 上一键配置国内镜像加速,并成功拉取 MySQL 5.7 镜像。
2604 3
Docker安装Mysql5.7,解决无法访问DockerHub问题