线性代数与数值方法--矩阵分解
来源:互联网 发布:怎么成为网络写手 编辑:程序博客网 时间:2024/05/16 18:39
矩阵
正交:正交最早出现于三维空间中的向量分析。 在三维向量空间中, 两个向量的内积如果是零, 那么就说这两个向量是正交的。
正交矩阵:正交矩阵是实数特殊化的酉矩阵,因此总是正规矩阵。如果:AAT=E(E为单位矩阵,AT表示“矩阵A的转置矩阵”。)或ATA=E,则n阶实矩阵A称为正交矩阵, 若A为正交阵,则满足以下条件:
1) AT是正交矩阵
2) AAT=E(E为单位矩阵)
3) A的各行是单位向量且两两正交
4) A的各列是单位向量且两两正交
5) (Ax,Ay)=(x,y) x,y∈R
6) |A| = 1或-1酉矩阵:n阶复方阵U的n个列向量是U空间的一个标准正交基,则U是酉矩阵(Unitary Matrix)。显然酉矩阵是正交矩阵往复数域上的推广。
正规矩阵:正规矩阵是与自己的共轭转置矩阵对易的复系数方块矩阵,A*A=AA*。其中A*是A的共轭转置矩阵。矩阵的正规性是检验矩阵是否可对角化的一个简便方法:任意正规矩阵都可在经过一个酉变换(正交变换)后变为对角矩阵,反过来所有可在经过一个酉变换后变为对角矩阵的矩阵都是正规矩阵。
正规矩阵的性质:属于正规矩阵不同特征值的特征向量两两正交。
定理1:在复数域上,A为正规矩阵的充分必要条件为A有n个两两正交的单位特征向量
定理2:在复数域上,A为正规矩阵的充分必要条件为A酉相似于对角矩阵奇异值分解(SVD)
奇异值:设A为复数域内m*n阶矩阵,A*表示A的共轭转置矩阵,A*·A的n个特征值的算术平方根(即A*·A的特征值的开方)叫作矩阵A的奇异值。记为σi(A) 。
共轭复数:两个实部相等,虚部互为相反数的的复数。
共轭转置:转置,取共轭。奇异值分解,即任意一个MxN的实数矩阵A均可写成:
其中p=min(m,n),矩阵U和V是正交矩阵,特征值分解
如果矩阵C是对称矩阵,那么特可以写成特征值分解的形式:QR因子分解
QR因子分解是一项广泛应用于稳定求解病态最小二乘问题的方法,同时也是一些更复杂方算法的基础,比如SVD及特征值分解。
A=QR
Q是正交矩阵,QQT=I,R是上三角矩阵,且在对角线上的元素为正。
Q可由格拉姆-施密特方法获得。
R=QTA
参考wiki:https://en.wikipedia.org/wiki/QR_decomposition- 格拉姆-施密特正交化
The Gram–Schmidt process then works as follows
参考wiki:https://en.wikipedia.org/wiki/Gram%E2%80%93Schmidt_process#Numerical_stability
6. 乔里斯基分解
参考wiki:https://en.wikipedia.org/wiki/Cholesky_decomposition
- 线性代数与数值方法--矩阵分解
- Matlab与线性代数--矩阵的LU分解
- Matlab与线性代数--矩阵的正交分解
- Matlab与线性代数--矩阵的Cholesky分解
- C++矩阵计算库Eigen3之:线性代数与分解
- 线性代数:矩阵的LU分解
- 数值分析--矩阵QR分解的三种方法
- 数值计算与方法 && 实验6||7 && LU分解求行列式的值和矩阵的逆
- 线性代数 NMF 非负矩阵分解 Matlab
- MIT线性代数--矩阵的LU分解
- 矩阵论与线性代数
- 线性代数与矩阵
- 【线性代数公开课MIT Linear Algebra】 第二十二课 特征分解与矩阵的幂
- Matlab与线性代数--矩阵的奇异值和奇异值分解
- 线性代数与矩阵论简史
- [线性代数] SVD 与 协方差矩阵
- Matlab与线性代数 -- 单位矩阵
- Matlab与线性代数 -- 零矩阵
- DML锁(TX、TM锁)和DDL锁(排他、共享、可中断解析锁)基础
- nyoj 20 吝啬的国度 dfs+建立树
- Dubbo高级篇3
- 下拉选择数据查询思考
- selenium(java)之RemoteWebDriver的基本使用
- 线性代数与数值方法--矩阵分解
- Codeforces Round #373 (Div. 2) ADCE题解
- POJ 3613 Cow Relays k步最短路 二分优化
- C++ Primer(第五版)练习5.5
- 初识docker--Ubuntu16.04
- CODEVS 1091 传染病控制
- MYSQL的23个常用查询语句
- 简单的jsp登录页面跳转功能
- TCP/IP详解:协议(基于UDP的DNS,TFTP,BOOTP)