[作者]
常用网名: 猪头三
出生日期: 1981.XX.XX
企鹅交流: 643439947
个人网站: 80x86汇编小站
编程生涯: 2001年~至今[共24年]
职业生涯: 22年
开发语言: C/C++、80x86ASM、Object Pascal、Objective-C、C#、R、Python、PHP、Perl、
开发工具: Visual Studio、Delphi、XCode、C++ Builder、Eclipse
技能种类: 逆向 驱动 磁盘 文件 大数据分析
涉及领域: Windows应用软件安全/Windows系统内核安全/Windows系统磁盘数据安全/macOS应用软件安全
项目经历: 股票模型量化/磁盘性能优化/文件系统数据恢复/文件信息采集/敏感文件监测跟踪/网络安全检测
专注研究: 机器学习、股票模型量化、金融分析
[描述]
Example 5: 如果 A A A 是可逆的上三角矩阵, 则其逆矩阵 A − 1 A^{-1} A−1 也必为上三角矩阵. 证明从矩阵方程
A A − 1 = I A\,A^{-1}=I AA−1=I
出发, 考察每一列的结构即可.
[核心解释]
为了说明 A − 1 A^{-1} A−1 是上三角矩阵, 只需证明在方程
A X = I A\,X=I AX=I
中, 未知矩阵 X X X 的每一列向量都在下三行(第 j + 1 j+1 j+1 行至第 n n n 行)为零. 设第 j j j 列向量为 x ( j ) x^{(j)} x(j), 即
A x ( j ) = e ( j ) , j = 1 , 2 , … , n , A\,x^{(j)}=e^{(j)},\quad j=1,2,\dots,n, Ax(j)=e(j),j=1,2,…,n,
其中 e ( j ) e^{(j)} e(j) 是标准基向量, 其第 j j j 个分量为 1, 其余为 0. 以下分三步详细推论:
-
列方程的形式
写出第 j j j 列的线性方程组:∑ k = 1 n a i k x k ( j ) = { 1 , i = j , 0 , i ≠ j . \sum_{k=1}^n a_{ik}\,x_k^{(j)} = \begin{cases} 1, & i=j,\\ 0, & i\neq j. \end{cases} k=1∑naikxk(j)={1,0,i=j,i=j.
由于 A A A 是上三角矩阵, 对角元 a i i ≠ 0 a_{ii}\neq 0 aii=0, 且当 i > k i>k i>k 时 a i k = 0 a_{ik}=0 aik=0.
-
回代求解并保持末尾零
采用自底向上的回代(back substitution):-
对行索引 i i i 从 n n n 递减到 1, 依次求解未知数 x i ( j ) x_i^{(j)} xi(j).
-
当 i > j i>j i>j 时, 右端为 0, 且方程仅涉及 x i ( j ) x_i^{(j)} xi(j) 与已知的更高序未知量 { x i + 1 ( j ) , … , x n ( j ) } \{x_{i+1}^{(j)},\dots,x_n^{(j)}\} {xi+1(j),…,xn(j)}:
a i i x i ( j ) + ∑ k = i + 1 n a i k x k ( j ) = 0. a_{ii}\,x_i^{(j)} + \sum_{k=i+1}^n a_{ik}\,x_k^{(j)} = 0. aiixi(j)+k=i+1∑naikxk(j)=0.
因此可解得
x i ( j ) = − 1 a i i ∑ k = i + 1 n a i k x k ( j ) . x_i^{(j)} = -\frac{1}{a_{ii}}\sum_{k=i+1}^n a_{ik}\,x_k^{(j)}. xi(j)=−aii1k=i+1∑naikxk(j).
但对于初始条件, 所有 x k ( j ) x_k^{(j)} xk(j)( k > i k>i k>i)在此阶段尚未赋非零值, 故必有
x i ( j ) = 0 , ∀ i > j . x_i^{(j)}=0,\quad \forall\,i>j. xi(j)=0,∀i>j.
-
-
对角及其以上分量求值
-
当 i = j i=j i=j 时, 右端为 1, 方程变为
a j j x j ( j ) + ∑ k = j + 1 n a j k x k ( j ) = 1. a_{jj}\,x_j^{(j)} + \sum_{k=j+1}^n a_{jk}\,x_k^{(j)} = 1. ajjxj(j)+k=j+1∑najkxk(j)=1.
因 { x k ( j ) ∣ k > j } = 0 \{x_{k}^{(j)}\mid k>j\}=0 {xk(j)∣k>j}=0, 得
x j ( j ) = 1 a j j ≠ 0. x_j^{(j)} = \frac{1}{a_{jj}} \neq 0. xj(j)=ajj1=0.
-
对于 i < j i<j i<j 的行, 右端仍为 0, 但此时所有更后行分量 { x i + 1 ( j ) , … , x n ( j ) } \{x_{i+1}^{(j)},\dots,x_n^{(j)}\} {xi+1(j),…,xn(j)} 已知, 其中第 j j j 行分量非零, 其它高行分量已知, 则依次可解出每个 x i ( j ) x_i^{(j)} xi(j), 不再影响下三行的零结构.
-
综上, 每个 x ( j ) x^{(j)} x(j) 在分量 j + 1 , … , n j+1,\dots,n j+1,…,n 均为零, 故矩阵 X = [ x ( 1 ) , x ( 2 ) , … , x ( n ) ] X=[x^{(1)},x^{(2)},\dots,x^{(n)}] X=[x(1),x(2),…,x(n)] 为上三角矩阵, 即
A − 1 = X 也是上三角矩阵 . A^{-1}=X\quad\text{也是上三角矩阵}. A−1=X也是上三角矩阵.
[总结]
通过对方程 A X = I A\,X=I AX=I 的列分量逐列回代求解, 证明了可逆的上三角矩阵 A A A 的逆矩阵 A − 1 A^{-1} A−1 保持上三角结构. 此性质在数值代数中至关重要, 保证了用三角分解(如 LU 分解)求解线性方程组时, 矩阵的稀疏性和带状结构不会因取逆而破坏, 从而维持算法的高效与稳定.
3806

被折叠的 条评论
为什么被折叠?



