BOOST的LU分解求解线性方程以及求逆
来源:互联网 发布:服装商品数据分析公式 编辑:程序博客网 时间:2024/06/04 22:57
原文网址:http://blog.sina.com.cn/s/blog_644dcdfc0100is9l.html
BOOST的LU分解求解线性方程以及求逆
1.LU分解
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <boost/numeric/ublas/operation.hpp>
#include <boost/numeric/ublas/triangular.hpp>
#include <boost/numeric/ublas/lu.hpp>
int _tmain(int argc, _TCHAR* argv[])
{
}
m=[1,2,1,1;2,1,0,0;1,0,1,0;1,0,0,1];
lu_factorize(m)的结果如下:
m=[1,2,1,1;2,-3,-2,-2;1,2/3,4/3,1/3;1,2/3,1/4,5/4];
其中L,U阵为
L=[1,0,0,0;2,1,0,0;1,2/3,1,0;1,2/3,1/4,1]; 即L=tril(m,-1)+eye(1);m的下三角阵与单位矩阵之和
U=[1,2,1,1;0,-3,-2,-2;0,0,4/3,1/3;0,0,0,4/5];即U=tril(m);m的上三角阵
可得m=L*U;
lu_factorize(m,P)的结果
P=[1,1,2,3];
m=[2,1,0,0;1/2,3/2,1,1;1/2,-1/3,4/3,1/3;1/2,-1/3,1/4,5/4];
L=[1,0,0,0;1/2,1,0,0;1/2,-1/3,1,0;1/2,-1/3,1/4,1];
U=[2,1,0,0;0,3/2,1,1;0,0,4/3,1/3;0,0,0,5/4];
Pm=LU,P=[1,1,2,3]
Pm = A,A的0行为m的1行,A的1行为新的1行(m的第0行),A的2行为m的2行,A的3的行为m的3行。
- BOOST的LU分解求解线性方程以及求逆
- 基于LU分解的矩阵求逆
- LU分解与求解
- Matlab实现——求矩阵的逆(LU分解)
- LU分解求线性方程组的解
- 矩阵LU分解求逆(学习笔记)
- 求解线性方程组之LU分解
- 矩阵的LU分解求解线性方程组(C++实现)
- LU分解求线型方程
- MATLAB的LU分解
- A的LU分解
- 矩阵LU的分解
- LU分解的实现
- OpenCV中LU分解实现矩阵求逆invert(DECOMP_LU)
- 数值计算与方法 && 实验6||7 && LU分解求行列式的值和矩阵的逆
- 三对角阵的LU分解和三对角方程组的求解(C语言)
- 工作规划(1) 基于LU分解的千万阶稀疏矩阵求解器
- 带状对角矩阵的LU分解及回代求解算法实现
- 新劳动法年假计算工具
- UVA 489 - Hangman Judge
- JS 中如何判断 undefined和null
- boost序列化代码
- myeclipse 6.0下fck_docprops.html和phpcompat.php 报错
- BOOST的LU分解求解线性方程以及求逆
- AudioFlinger的SRC造成的谐波失真测试
- TFIDF算法java实现
- 电力负荷控制终端(Q/GDW 376.1)整体解决方案
- Arrays.sort()对自己定义的类进行排序
- Trie树及其应用
- JavaScript细微差别(函数声明与函数表达式及浏览器差异)
- mxml中嵌入代码时出错
- Spring事件驱动模型详解