Linear Algebra lecture4 note
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的更多相关文章
- Linear Algebra lecture1 note
Professor: Gilbert Strang Text: Introduction to Linear Algebra http://web.mit.edu/18.06 Lecture 1 ...
- Linear Algebra lecture9 note
Linear independence Spanning a space Basis and dimension 以上概念都是针对a bunch of vectors, 不是矩阵里的概念 Supp ...
- Linear Algebra lecture10 note
Four fundamental subspaces( for matrix A) if A is m by n matrix: Column space C(A) in Rm (列空间在m维实 ...
- Linear Algebra lecture8 note
Compute solution of AX=b (X=Xp+Xn) rank r r=m solutions exist r=n solutions unique example: 若想方程有解 ...
- Linear Algebra lecture7 note
Computing the nullspace (Ax=0) Pivot variables-free variables Special solutions: rref( A)=R rank o ...
- Linear Algebra lecture6 note
Vector spaces and subspaces Column space of A solving Ax=b Null space of A Vector space requiremen ...
- Linear Algebra Lecture5 note
Section 2.7 PA=LU and Section 3.1 Vector Spaces and Subspaces Transpose(转置) example: 特殊情况,对称 ...
- Linear Algebra lecture3 note
Matrix multiplication(4 ways!) Inverse of A Gauss-Jordan / find inverse of A Matrix multiplication ...
- 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 ...
随机推荐
- as3 中文转拼音
private static const PinYin:Object = {"a":"\u554a\u963f\u9515","ai":&q ...
- 基本矩阵运算的Java实现
一: 矩阵的加法与减法 规则:矩阵的加法与减法要求两个矩阵的行列完全相等,方可以完成两个矩阵的之间的运算. 举例说明如下 二:矩阵的乘法 规则:矩阵的乘法要求两个矩阵符合A(mx k), B( ...
- a标签
a链接是一种触发行为元素,行内元素 属性:href/target/class/id/title href: <a href="www.baidu.com"></a ...
- Google Authentication 机制原理
Google Authenticator,谷歌身份认证器,Google公司推出的一款动态口令工具,旨在解决大家Google账户遭到恶意攻击的问题.该工具主要基于TOTP(Time-Based One- ...
- 硬件抽象层——HAL
(一)为什么要在android中加入HAL Linux系统中Linux驱动有两种类型的代码:访问硬件寄存器的代码——调用的Linux内核的标准函数进行的标准操作 业务逻辑代码——有些企业或个人并不想 ...
- Spring中的自动装配
src\dayday\Person.java package dayday;/** * Created by I am master on 2016/11/28. */public class Per ...
- mac上spacemacs体验小记
project: blog target: note-of-spacemacs-on-mac.md date: 2016-01-04 status: publish tags: - emacs - s ...
- JAVA Arrays.binarySearch
转自:http://blog.csdn.net/somebodydie/article/details/8229343 package com.jaky; import java.util.*; pu ...
- 虚机centos和本机Windows之间文件的拷贝无法用xftp时用FileZilla也行
步骤如下: 1.如果Centos没有安装ssh,则需要先安装: 2.查看虚拟机中IP地址,命令如下: ifconfig 3.在windows中安装ftp软件 FileZilla启动软件如图: 6 这 ...
- system_call中断处理过程分析
本文所有的分析内容都是基于Linux3.18.6内核,鉴于对应不同内核版本,系统调用的实现不相同.若需要分析其他版本内核的系统调用的实现过程,请谨慎参考. system_call函数的功能是用来响应外 ...