[原创](计算机数学)(Introduction Linear Algebra)(P87): :若A是可逆的上三角矩阵, 则它的逆矩阵A^-1也是上三角矩阵

[作者]
常用网名: 猪头三
出生日期: 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} A1 也必为上三角矩阵. 证明从矩阵方程

A   A − 1 = I A\,A^{-1}=I AA1=I

出发, 考察每一列的结构即可.

[核心解释]

为了说明 A − 1 A^{-1} A1 是上三角矩阵, 只需证明在方程

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. 以下分三步详细推论:

  1. 列方程的形式
    写出第 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=1naikxk(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.

  2. 回代求解并保持末尾零
    采用自底向上的回代(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+1naikxk(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+1naikxk(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.

  3. 对角及其以上分量求值

    • 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+1najkxk(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{也是上三角矩阵}. A1=X也是上三角矩阵.

[总结]

通过对方程 A   X = I A\,X=I AX=I 的列分量逐列回代求解, 证明了可逆的上三角矩阵 A A A 的逆矩阵 A − 1 A^{-1} A1 保持上三角结构. 此性质在数值代数中至关重要, 保证了用三角分解(如 LU 分解)求解线性方程组时, 矩阵的稀疏性和带状结构不会因取逆而破坏, 从而维持算法的高效与稳定.

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

我不是代码教父

我的创作动力离不开你的真诚激励

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值