cholesky分解java代码_cholesky分解的实现
Cholesky 分解是把一个对称正定的矩阵表示成一个下三角矩阵L和其转置的乘积的分解。
它要求矩阵的所有特征值必须大于零,故分解的下三角的对角元也是大于零的。
Cholesky分解法又称平方根法,是当A为实对称正定矩阵时,LU三角分解法的变形。
通过直接比较A=L*L^T两边的对应元素来计算L,其中L^T为L的转置。
思路如下:
L为一实下三角矩阵,求L的步骤如下:
1、Amn = Lm1*Ln1 + Lm2*Ln2 + ... + Lmx*Lnx其中x = min(m,n)
2、Umn = Amn - sum(Lmk*Lnk) 其中k ~ (0, min(m,n)),Umn包含Lmn*Lnn
3、当m < n时,由于是下三角矩阵,所以Umn为0,仅求m >= n的情况
4、如果 m == n 时 Lmn = sqrt(Umn),否则,Lmn = Umn / Lnn
5、即求出Lmn的值
根据此思路的代码实现如下:
public class MyCholeskyDecomposition {
/**
* 2.0000000000 0.0000000000 0.0000000000
* 0.5000000000 1.3228756555 0.0000000000
* 0.5000000000 2.0788046016 1.1952286093
* @param args
*/
public static void main(String[] args) {
double[][] A = {{4.,1.,1.},{1.,2.,3.},{1.,3.,6.}};
double[][] L = new double[3][3];
for(int m = 0; m < A.length; m++){
for(int n = 0; n <= m; n++){
L[m][n] = A[m][n];
for(int k = 0; k < n; k++){
L[m][n] -= L[m][k] * L[n][k];
}
if(m == n){
L[m][n] = Math.sqrt(L[m][n]);
}else{
L[m][n] = L[m][n] / L[n][n];
}
}
for(int x = m + 1; x < A.length; x++){
L[m][x] = 0.0;
}
}
for(int i = 0; i < L.length; i++){
for(int j = 0; j < L.length; j++){
System.out.print(L[i][j] + " ");
}
System.out.println();
}
}
}
————————————————
版权声明:本文为CSDN博主「expleeve」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/expleeve/article/details/51933762
cholesky分解java代码_cholesky分解的实现相关推荐
- cholesky分解java代码_Cholesky 分解(转)
Cholesky 分解是把一个对称正定的矩阵表示成一个下三角矩阵L和其转置的乘积的分解. 它要求矩阵的所有特征值必须大于零,故分解的下三角的对角元也是大于零的. Cholesky分解法又称平方根法,是 ...
- c# lu分解的代码_LU分解(1)
1/6 LU 分解 LU 分解可以写成A = LU,这里的L代表下三角矩阵,U代表上三角矩阵.对应的matlab代码如下: function[L, U] =zlu(A) % ZLU - LU deco ...
- qr分解java代码_[转载]【原创】基于Givens变换QR分解Matlab代码
如果有任何问题.建议,或者更多资源.代码.视频,欢迎您访问专业的Matlab技术交流平台--Matlab技术论坛http://www.matlabsky.com Wish my dear girl,h ...
- 递归回溯--数字分解java代码
/* * 5.数字分解 * * 6 * 6 * 5+1 * 4+2 * 4+1+1 * 3+3 * 3+2+1 * 3+1+1+1 * 2+2+2 ...
- 小波分解与小波包分解代码_分解的功能参数和代码可维护性
小波分解与小波包分解代码 Code keeps changing, there's no doubt about that. We always do our best to set some roc ...
- 面对困难的代码,分解困难各个击破
怎么说呢,最近遇到一大堆程序的时候真的有点手足无措,感觉很被动阿. 迎难而上才能有所得.这次的经验是面对一大堆要写的代码该怎么办呢? 答案:先搭框架! 其实上面这四个字着实有点太忽悠人,写代码少的人很 ...
- c语言素数筛法与分解素因数,质因数分解及代码:
计算方法 短除法 求一个数分解质因数,要从最小的质数除起,一直除到结果为质数为止.分解质因数的算式的叫短除法,和除法的性质差不多,还可以用来求多个个数的公因式: 求最大公因数的一种方法,也可用来求最小 ...
- LMD改进的局部均值分解matlab代码模版
LMD改进的局部均值分解matlab代码模版 ID:9520675130845877比邻星机灵的薯片
- 颜色分类识别代码matlab——分解RGB通道通过阈值来判断
颜色分类识别代码matlab--分解RGB通道通过阈值来判断 代码下载链接 代码下载链接 代码下载链接 首先有一张包含多种颜色的图片: 然后可以编写代码,来实现分别提取出不同颜色的操作:
- 【二】Java编程之分解质因素
质因数百度百科解释 质因数(素因数或质因子)在数论里是指能整除给定正整数的质数.除了1以外,两个没有其他共同质因子的正整数称为互质.因为1没有质因子,1与任何正整数(包括1本身)都是互质.正整数的因数 ...
最新文章
- ORB-SLAM2从理论到代码实现(八):Tracking.cc程序详解(下)
- 注意事项:二、sqlserver2008服务配置
- Python subprocess
- nginx-rtmp message
- ppt上的倒计时小工具_办公小技巧:轻松玩转PPT秒针倒计时
- 802.11 n wlan linux驱动下载,802.11n无线网卡驱动
- python编辑器中文字体倒立的_matplotlib的安装和允许中文及几种字体
- 【Qt串口调试助手】1.5 - 发送/接收字节数量统计
- python实践项目(五)
- Facebook开源项目:我们为什么要用Fresco框架?
- java中reg=new,JavaScript replace new RegExp使用介绍
- HTML右下角音乐播放器,利用HTML5实现网页底部音乐播放器代码
- linux dstat,dstat 用法详解
- Python坐标系转换
- 【Arcgis】球面坐标系转投影坐标
- notion函数_最适合设计师的笔记软件 Notion
- 记录一次面试经历(深刻)
- java编程题身高排队_编程初学者入门4_从键盘输入5个人的身高(米),求他们的平均身...
- JVM之jstat命令
- 再生龙盘对盘拷贝Linux
热门文章
- ubuntu14.04安装QQ
- WIN10_用户获取最高的管理员权限(关闭UAC控制)
- 我的账号 小米云服务器地址,小米云服务登录
- 远程登陆cisco防火墙
- cisco 防火墙安全策略
- 仿《91创业网》网站源码 招商加盟致富商机网站 帝国cms模版+采集
- 用Python回忆QQ空间里的青春
- 诺基亚x6 android one,诺基亚X6手机推送固件更新:提升安卓8.1系统流畅性,新增后台锁定...
- 创联网络宽带上网助手的脱壳
- JDBC————PreparedStatement批处理