[问题2014A04] 解答
[问题2014A04] 解答
(1) 由条件可得 \(AB+BA=0\), 即 \(AB=-BA\), 因此 \[AB=A^2B=A(AB)=A(-BA)=-(AB)A=-(-BA)A=BA^2=BA,\] 从而 \(AB=BA=0\).
(2) 由条件可得 \(0=B(AB)^kA=(BA)^{k+1}\), 因此 \[(I_n-BA)\Big(I_n+BA+\cdots+(BA)^k\Big)=I_n,\] 从而 \(I_n-BA\) 可逆.
(3) 我们给出此小题的三种解法.
解法一(凑因子法)
凑因子法即为将 \(A-BD^{-1}C\) 的逆阵给凑出来, 方法的关键就是不断地变形, 凑出 \(A-BD^{-1}C\) 这个因子. 我们将分成若干个步骤对这一典型例题加以说明.
先设 \(H=(D-CA^{-1}B)^{-1}\), 则 \[(D-CA^{-1}B)H=I_n. \cdots(1)\] (1) 式是我们的出发点, 接下来就开始变形了. 我们的目标是凑出 \(A-BD^{-1}C\), 所以需要的是 \(D^{-1}\), 而不是 \(D\), 于是 (1) 式两边同时左乘 \(D^{-1}\) 可得 \[(I_n-D^{-1}CA^{-1}B)H=D^{-1}. \cdots(2)\] 为了凑出 \(A-BD^{-1}C\), 在 (2) 式两边同时左乘 \(B\) 右乘 \(C\) 可得 \[BHC-BD^{-1}CA^{-1}BHC=BD^{-1}C. \cdots(3)\] (3) 式左边提出公因子 \(A^{-1}BHC\), 右边的 \(BD^{-1}C\) 移到左边, 并且两边同时加上 \(A\) 以凑出 \(A-BD^{-1}C\), 可得 \[(A-BD^{-1}C)A^{-1}BHC+(A-BD^{-1}C)=A. \cdots(4)\] 将 (4) 式左边的公因子 \(A-BD^{-1}C\) 提出, 并将两边同时右乘 \(A^{-1}\) 可得 \[(A-BD^{-1}C)\Big(I_n+A^{-1}BHC\Big)A^{-1}=I_n. \cdots(5)\] 由 (5) 式即得 \[(A-BD^{-1}C)^{-1}=A^{-1}+A^{-1}BHCA^{-1}=A^{-1}+A^{-1}B(D-CA^{-1}B)^{-1}CA^{-1}.\,\,\Box\]
解法二(利用已证结论)
上课时我证明过以下结论:
若 \(I_n-AB\) 可逆, 则 \(I_n-BA\) 也可逆, 且 \((I_n-BA)^{-1}=I_n+B(I_n-AB)^{-1}A\).
当时我用了凑因子法和幂级数展开+验证法这两种方法去证明上述结论, 而且这个结论也是本小题的特例. 由降阶公式易证 \(|A-BD^{-1}C|\neq 0\), 因此 \(A-BD^{-1}C\) 非异. 我们进行如下的变形:
\[(A-BD^{-1}C)^{-1}=\Big(A(I_n-A^{-1}BD^{-1}C)\Big)^{-1}=(I_n-A^{-1}BD^{-1}C)^{-1}A^{-1}.\] 将 \(A^{-1}B\) 与 \(D^{-1}C\) 分别看成两个整体, 利用上述结论可得
\[(A-BD^{-1}C)^{-1}=\Big(I_n+A^{-1}B(I_n-D^{-1}CA^{-1}B)^{-1}D^{-1}C\Big)A^{-1}\]
\[=\Big(I_n+A^{-1}B(D-CA^{-1}B)^{-1}C\Big)A^{-1}=A^{-1}+A^{-1}B(D-CA^{-1}B)^{-1}CA^{-1}.\,\,\Box\]
解法三(分块初等变换法)
按照课本上降阶公式的证法, 分块矩阵 \(\begin{bmatrix} A & B \\ C & D \end{bmatrix}\) 可以通过分块初等变换变为分块对角阵 \(\begin{bmatrix} A & 0 \\ 0 & D-CA^{-1}B \end{bmatrix}\) 以及 \(\begin{bmatrix} A-BD^{-1}C & 0 \\ 0 & D \end{bmatrix}\). 用分块初等阵的乘法去改写上述过程即有
\[\begin{bmatrix} I_n & 0 \\ -CA^{-1} & I_n \end{bmatrix}\begin{bmatrix} A & B \\ C & D \end{bmatrix}\begin{bmatrix} I_n & -A^{-1}B \\ 0 & I_n \end{bmatrix}=\begin{bmatrix} A & 0 \\ 0 & D-CA^{-1}B \end{bmatrix},\]
\[\begin{bmatrix} I_n & -BD^{-1} \\ 0 & I_n \end{bmatrix}\begin{bmatrix} A & B \\ C & D \end{bmatrix}\begin{bmatrix} I_n & 0 \\ -D^{-1}C & I_n \end{bmatrix}=\begin{bmatrix} A-BD^{-1}C & 0 \\ 0 & D \end{bmatrix}.\]
因此我们有
\[\begin{bmatrix} (A-BD^{-1}C)^{-1} & 0 \\ 0 & D^{-1} \end{bmatrix}=\begin{bmatrix} I_n & 0 \\ D^{-1}C & I_n \end{bmatrix}\begin{bmatrix} A & B \\ C & D \end{bmatrix}^{-1}\begin{bmatrix} I_n & BD^{-1} \\ 0 & I_n \end{bmatrix}\]
\[=\begin{bmatrix} I_n & 0 \\ D^{-1}C & I_n \end{bmatrix}\begin{bmatrix} I_n & -A^{-1}B \\ 0 & I_n \end{bmatrix}\begin{bmatrix} A^{-1} & 0 \\ 0 & (D-CA^{-1}B)^{-1} \end{bmatrix}\begin{bmatrix} I_n & 0 \\ -CA^{-1} & I_n \end{bmatrix}\begin{bmatrix} I_n & BD^{-1} \\ 0 & I_n \end{bmatrix}\]
\[=\begin{bmatrix} A^{-1}+A^{-1}B(D-CA^{-1}B)^{-1}CA^{-1} & 0 \\ 0 & D^{-1} \end{bmatrix},\]
从而 \((A-BD^{-1}C)^{-1}=A^{-1}+A^{-1}B(D-CA^{-1}B)^{-1}CA^{-1}\). \(\,\,\Box\)
[问题2014A04] 解答的更多相关文章
- [问题2014A11] 解答
[问题2014A11] 解答 我们需要利用以下关于幂等阵判定的结论,它是复旦高代书第 142 页的例 3.6.4: 结论 设 \(A\) 为 \(n\) 阶方阵, 则 \(A^2=A\) 当且仅当 ...
- 精选30道Java笔试题解答
转自:http://www.cnblogs.com/lanxuezaipiao/p/3371224.html 都 是一些非常非常基础的题,是我最近参加各大IT公司笔试后靠记忆记下来的,经过整理献给与我 ...
- 精通Web Analytics 2.0 (8) 第六章:使用定性数据解答”为什么“的谜团
精通Web Analytics 2.0 : 用户中心科学与在线统计艺术 第六章:使用定性数据解答"为什么"的谜团 当我走进一家超市,我不希望员工会认出我或重新为我布置商店. 然而, ...
- 【字符编码】Java字符编码详细解答及问题探讨
一.前言 继上一篇写完字节编码内容后,现在分析在Java中各字符编码的问题,并且由这个问题,也引出了一个更有意思的问题,笔者也还没有找到这个问题的答案.也希望各位园友指点指点. 二.Java字符编码 ...
- spring-stutrs求解答
这里贴上applicationContext里的代码: <?xml version="1.0" encoding="UTF-8"?> <bea ...
- JavaScript Bind()趣味解答 包懂~~
首先声明一下,这个解答是从Segmentfault看到的,挺有意思就记录下来.我放到最下面: bind() https://developer.mozilla.org/zh-CN/docs/Web/J ...
- CMMI4级实践中的5个经典问题及解答
这五个问题相当经典而且比较深,需要做过CMMI4.5级的朋友才能看懂这些问题.这5个问题是一位正在实践CMMI4级的朋友提出来的,而解答则是我的个人见解. 五个疑问是: A.流程,子流程部分不明白 ...
- 海边直播目标2017全国初中数学竞赛班课堂测试题解答-The Final
1. 设函数 $f(x) = 2^x(ax^2 + bx + c)$ 满足等式 $f(x+1) - f(x) = 2^x\cdot x^2$, 求 $f(1)$. 解答: 由 $f(x) = 2^x( ...
- 知乎大牛的关于JS解答
很多疑惑一扫而空.... http://www.zhihu.com/question/35905242?sort=created JS的单线程,浏览器的多进程,与CPU,OS的对位. 互联网移动的起起 ...
随机推荐
- 关于UGUI Image Sliced模式的一个BUG。
Unity4.6.2f1 在Android/IOS平台下,Image选择Sliced模式,并且对Sprite设置好Border后,会发现并没有按照预计的 情况进行拉伸. 搜了一下是因为Sprite的G ...
- mongodb备忘
1.远程拷贝数据库 db.copyDatabase(fromdb, todb, fromhost, [dbuser, dbpassword]) 2.数据库备份/恢复(导出/导入) mongoexpor ...
- Java NIO 网络编程基础
Java NIO提供了一套网络api,可以用来处理连接数很多的情况.他的基本思想就是用一个线程来处理多个channel. 123456789101112131415161718192021222324 ...
- response实现文件下载
package cn.itcast.response; import java.io.FileInputStream; import java.io.IOException; import java. ...
- 【iCore3 双核心板_ uC/OS-III】例程十一:任务消息队列
实验指导书及代码包下载: http://pan.baidu.com/s/1pLQYiE3 iCore3 购买链接: https://item.taobao.com/item.htm?id=524229 ...
- HDU 5919 Sequence II(主席树+逆序思想)
Sequence II Time Limit: 9000/4500 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others) To ...
- PHP 标准库 SPL 之数据结构栈(SplStack)简单实践
PHP 5.3.0 版本及以上的堆栈描述可以使用标准库 SPL 中的 SplStack class,SplStack 类继承双链表 ( SplDoublyLinkedList ) 实现栈. 代码: & ...
- MySQL查询优化:查询慢原因和解决技巧
在开发的朋友特别是和mysql有接触的朋友会碰到有时mysql查询很慢,当然我指的是大数据量百万千万级了,不是几十条了,下面我们来看看解决查询慢的办法. MySQL查询优化:查询慢原因和解决方法 会经 ...
- iOS cocospods Updating local specs repositories
pod install --verbose --no-repo-update (在安装的时候) pod update --verbose --no-repo-update (在更新库的时候) 如果长时 ...
- DIA
Dia Tutorial: http://www.seanet.com/~hgg9140/comp/diatut/all/all.html Dia Tutorial