Inverse of AB,A^(A的转置)

Product of elimination matrices  A=LU (no row exchanges)

 


Inverse of AB,A^(A的转置):

 


Product of elimination matrices  A=LU (no row exchanges)

E32E31E21A=U (no row exchanges)    EA=U

A=E21`E31`E32`U

L表示下三角矩阵,lower triangle

D表示对角矩阵,diagonal triangle

A=LU

L=E21`E31`E32`

对比EA=U  &   A=LU,哪种形式更好一些?

example:

设E31为单位矩阵,

A=LU的形式更加简洁一些,if no row exchanges,  multipliers go directly into L

思考:How many operations on n * n Matrix(A)? 也就是多少数字改变了?

假设n=100,

row1 不改变,后续行消元,99*100次运算,约看作100*100次;

对第二行进行消元,98*99次运算,约看作99*99次;

继续…

所以总运算次数为

从积分角度考虑可得最终结果

若再考虑等式右侧的b,需要的运算量为

 


Permutation

考虑单位矩阵的行变换,所有情况一一列出(共6种)

这6个矩阵好像构成了一个群,不论是乘还是逆运算,所得结果均在此群中,且满足条件

若是4维矩阵,则满足条件的、构成群的矩阵有4!=24个

若是n维矩阵,则满足条件的、构成群的矩阵有n!个

Linear Algebra lecture4 note的更多相关文章

  1. Linear Algebra lecture1 note

    Professor: Gilbert Strang Text: Introduction to Linear Algebra http://web.mit.edu/18.06   Lecture 1 ...

  2. Linear Algebra lecture9 note

    Linear independence Spanning a space Basis and dimension 以上概念都是针对a bunch of vectors, 不是矩阵里的概念   Supp ...

  3. Linear Algebra lecture10 note

    Four fundamental subspaces( for matrix A)   if A is m by n matrix: Column space  C(A) in Rm (列空间在m维实 ...

  4. Linear Algebra lecture8 note

    Compute solution of AX=b (X=Xp+Xn) rank r r=m solutions exist r=n solutions unique   example: 若想方程有解 ...

  5. Linear Algebra lecture7 note

    Computing the nullspace (Ax=0) Pivot variables-free variables Special solutions: rref( A)=R   rank o ...

  6. Linear Algebra lecture6 note

    Vector spaces and subspaces Column space of A solving Ax=b Null space of A   Vector space requiremen ...

  7. Linear Algebra Lecture5 note

    Section 2.7     PA=LU and Section 3.1   Vector Spaces and Subspaces   Transpose(转置) example: 特殊情况,对称 ...

  8. Linear Algebra lecture3 note

    Matrix multiplication(4 ways!) Inverse of A Gauss-Jordan / find inverse of A   Matrix multiplication ...

  9. Codeforces Gym101502 B.Linear Algebra Test-STL(map)

    B. Linear Algebra Test   time limit per test 3.0 s memory limit per test 256 MB input standard input ...

随机推荐

  1. ldap + kerberos + google authentication 实现两步验证

    第一步:ldap + kerberos 整合  ,参考之前的文章 第二步:google authentication 安装配置,参考之前的文章 第三步:整合 ldap + kerberos + goo ...

  2. popen使用不当引起产生僵尸进程

    FILE * popen(const char * command, const char * type)popen函数会通过fork产生子进程,然后从子进程中调用/bin/sh -c执行参数comm ...

  3. 按月统计tcode和report使用次数的工具

    执行report,输入要查询的日期和user, 工具会按照使用次数从高到低列出输入日期所在的月份内所有该user 曾经使用过的tcode 和report list: REPORT zusertcode ...

  4. Java并发编程学习笔记(三)——对象的组合

    重要概念: 1.在设计线程安全类的过程中,需要包含以下三个基本要素: (1)找出构成对象状态的所有变量. (2)找出约束状态变量的不变性条件. (3)建立对象状态的并发访问管理策略. 2.

  5. UIkit框架之uiUIapplication

    1.继承链:uiresponder:NSObject 2.使用 sharedApplication方法来存取对象 3.这个类可以遵守协议是UIApplicationDelegate 4.获取app的实 ...

  6. [转]使用Scrapy建立一个网站抓取器

    英文原文:Build a Website Crawler based upon Scrapy 标签: Scrapy Python 209人收藏此文章, 我要收藏renwofei423 推荐于 11个月 ...

  7. H3C汇聚层交换机认证在线人数展示系统之需求说明和功能点说明

    一.需求 (一)每五分钟查询一次交换机的连接情况: (二)每2.5分钟更新每栋楼的连接情况. 二.功能点 序号 功能点说明 待定 完成 未完成 完成时间 预计用时(min) 实际用时(min) 备注 ...

  8. Protobuf - 使用scons编译proto文件

    使用protobuf过程中,需要先对消息结构进行定义,文件以.proto格式结尾.然后要使用google提供的protoc命令行,把.proto文件转成对应的代码文件. protoc --proto_ ...

  9. (转)appendChild()、insertBefore()是移动element节点!

    原文地址 appendChild().insertBefore()是移动element节点,看书的时候注意过,也可以做一个简单的例子测试一下: <div id="div1"& ...

  10. 关于Windows7的安装纠结过程

    关于Windows7的安装过程 背景交代 因为自己有两个笔记本,一个是伴我读过四年大学生活的老华硕笔记本,一个是姐夫不用的上网本.自己的构想是将老笔记本装上Ubuntu Server当作平时开发pyt ...