对单位下三角矩阵的意外发现
文章目录
- 背景
- 正文
- 结论
- 证明
背景
今天复习数值分析在对矩阵进行Doolittle分解时发现了一些有意思的事情,请容我慢慢道来…
正文
先介绍矩阵得三角分解即Doolittle分解定理:
定理 :A∈Rn×nA\in R^{n\times n}A∈Rn×n,若detA≠0,Δi≠0(i=1,2,…,n−1)\,detA\ne0,\Delta_i\ne0(i=1,2,\dots,n-1)detA=0,Δi=0(i=1,2,…,n−1),则存在唯一的单位下三角矩阵L\,L\,L和非奇异的上三角矩阵U\,U\,U,使A=LUA = L\,UA=LU在本定理已知条件下可以进行 n - 1 步消去,得到一个非奇异上三角形矩阵,不妨记该矩阵为U\,U\,U,则有Ln−1−1Ln−2−1…L2−1L1−1=UL_{n-1}^{-1}L_{n-2}^{-1}\dots L_{2}^{-1}L_{1}^{-1}=ULn−1−1Ln−2−1…L2−1L1−1=U
所以L=L1L2…Ln−1=[1l211l31l32⋱⋮⋮⋱1ln1ln2⋯ln,n−11]L=L_1L_2\dots L_{n-1}= \begin{bmatrix} 1 \\ l_{21} & 1\\ l_{31} & l_{32} & \ddots \\ \vdots & \vdots & \ddots & 1\\ l_{n1} & l_{n2} & \cdots & l_{n,n-1} & 1 \end{bmatrix} L=L1L2…Ln−1=⎣⎢⎢⎢⎢⎢⎡1l21l31⋮ln11l32⋮ln2⋱⋱⋯1ln,n−11⎦⎥⎥⎥⎥⎥⎤
其中Li=[1⋱1li+1,i1⋮⋱ln,i1]L_i= \begin{bmatrix} 1 \\ &\ddots\\ &&1 \\ &&l_{i+1,i}&1\\ &&\vdots&&\ddots\\ &&l_{n,i}& &&1 \end{bmatrix}Li=⎣⎢⎢⎢⎢⎢⎢⎢⎡1⋱1li+1,i⋮ln,i1⋱1⎦⎥⎥⎥⎥⎥⎥⎥⎤
而Li−1=[1⋱1−li+1,i1⋮⋱−ln,i1]L_i^{-1}= \begin{bmatrix} 1 \\ &\ddots\\ &&1 \\ &&-l_{i+1,i}&1\\ &&\vdots&&\ddots\\ &&-l_{n,i}& &&1 \end{bmatrix}Li−1=⎣⎢⎢⎢⎢⎢⎢⎢⎡1⋱1−li+1,i⋮−ln,i1⋱1⎦⎥⎥⎥⎥⎥⎥⎥⎤
通过观察,我们可以得到结论如下
结论
设Li,Lj∈Rn×n(0<i,j≤n∧i≠j)\,L_i,L_j\in R^{n\times n}(0<i,j\le n\wedge i\ne j)\,Li,Lj∈Rn×n(0<i,j≤n∧i=j),且满足如下形式:
Lk=[1⋱1lk+1,k1⋮⋱ln,k1]L_k= \begin{bmatrix} 1 \\ &\ddots\\ &&1 \\ &&l_{k+1,k}&1\\ &&\vdots&&\ddots\\ &&l_{n,k}& &&1 \end{bmatrix}Lk=⎣⎢⎢⎢⎢⎢⎢⎢⎡1⋱1lk+1,k⋮ln,k1⋱1⎦⎥⎥⎥⎥⎥⎥⎥⎤
则有:
- Li×Lj=Li+Lj−IL_i\times L_j = L_i+L_j-ILi×Lj=Li+Lj−I
- Li−1=2I−LiL_i^{-1}=2I-L_iLi−1=2I−Li
其中I\,I\,I为单位矩阵
证明
我才懒得证呢…
对单位下三角矩阵的意外发现相关推荐
- 评论:26岁成都唐爽发现成果惊动奥巴马--意外发现新材料极可能成下一代电脑芯片...
评论:在这篇文章中,提到马未都老师对他说的一句话,'比时代快一半,你就会很幸福:与时代的节奏相同,你只能庸庸碌碌:比时代快一步,那你可能就会痛苦'.等你快了一步,可能就有如同独孤求败之感,如同莫奈一般 ...
- 小米8计算机快捷键,小米8这4个隐藏功能,使用半年今天才意外发现
原标题:小米8这4个隐藏功能,使用半年今天才意外发现 小米8是今年最出色的小米机型,不仅外观设计好看,机身性能还很强悍.而且相对于同配置的机型,小米8的价格更亲民,当然不仅如此,小米的MIUI系统也是 ...
- 意外发现新效应,量子退火计算能解决世界两大难题:反物质不对称与同位素分离
美国洛斯阿拉莫斯国家实验室(LANL)的一支团队在研究量子退火计算机的运行速度为什么比预期快很多时,意外发现了一个新的效应. 这种效应有可能解决两个著名的科学问题,一是宇宙中物质和反物质分布不平衡的原 ...
- rust全息要啥才能做_绝地求生玩家意外发现红点秘密,网友表示以后终于不用带全息了!...
相信很多玩绝地求生的玩家都曾经纠结过一个问题,那就是红点和全息到底选择用哪个!当然对于很多玩家来说,其实这两种不同的瞄具大家都会使用,但是通常两种用的都不是很好,而对于一些高分段的玩家来说,红点瞄准镜 ...
- 2019年十大科学发现:一个年轻人的意外发现,见证了恐龙灭绝的那一天
来源:大数据实验室 2019年12月20日,国际顶尖学术期刊Science公布了2019年度十大科学突破,其中一项是发表在PNAS的论文,一个年轻人的意外发现,见证了6600万年前小行星撞击地球导致恐 ...
- 博士意外发现秘密:DALL-E 2自创语言,人看不懂但能生成特定图像,可能会被用来搞事!...
丰色 金磊 发自 凹非寺 量子位 | 公众号 QbitAI DALL·E 2,这个AI竟然自己搞了一套秘密语言. 例如这两个非常奇怪的词组: Apoploe vesrreaitais(下文简称A) C ...
- 一场实验室意外爆炸事故,解决了58年量子难题,让科学家意外发现“核电共振”...
晓查 发自 凹非寺 量子位 报道 | 公众号 QbitAI "核磁共振"想必大多数人都听过,可是你听说过"核电共振"吗? 早在1961年,核磁共振的先驱.诺贝 ...
- 意外发现老版本cyberlink decoder存在的问题
三个月前做的多路播放高清H.264视频的SDK,在应用过程中挑机子.大部分机型很稳定,个别机型总是存在着弹出Active Movie窗口现象.以前一直怀疑是自己使用VMR的问题,本周终于在同事的提示下 ...
- 解决spring问题(Unsatisfied dependency)的意外发现
前几天在做一个Dao的单元测试的时候,碰到了一个spring的错误.如下: [quote]org.springframework.beans.factory.UnsatisfiedDependency ...
最新文章
- 关于Linux的一些个人研习感悟
- python导入txt为dataframe-Python提取TXT数据转化为DataFrame
- 学习《FreeRTOS源码详解与应用开发》笔记
- boost::units::conversion_factor相关的测试程序
- 洪水同频率放大的方法_我们应该怎么选择放大器配件?
- RAISERROR (Transact-SQL)
- PyTorch官方权威教程书来了,LeCun力荐!意外的通俗易懂
- Linux下的虚拟Bridge实现
- 接口测试(二)--APP抓包
- 面试官:new Object[5] 一共创建了几个对象?
- 西门子PLC编写抢答系统
- 解决: Couldn‘t find meta-data for provider with authority xxx.fileProvider
- Excel 点击单元格打钩,再点击取消
- mac教程:重置苹果Mac上的NVRAM或PRAM?
- NL驱动表错误导致的性能问题
- 《先知·逸乐》| 《先知·自由》
- [Java并发]の其二
- 2021-11-15 cv2.erode()和cv2.dilate()的原理理解
- 学习Mybatis第三天
- [opencv完整项目详解] 传统图像算法解决路标的检测和识别(改进升级版)