MATLAB Cholesky分解

Cholesky.m文件：

function [L,X] = Cholesky(A,b)
% 先分解
n=size(A,1);
L=zeros(n);
L(1,1:n)=A(1,1:n);
L(2:n,2:n)=A(2:n,2:n);
for i=1:1:nL(i,i)=sqrt(L(i,i));u=zeros(1,n-i);for j=i+1:1:nL(i,j)=L(i,j)/L(i,i);u(1,j-i)=L(i,j);end
U=u'*u;
L(i+1:n,i+1:n)=L(i+1:n,i+1:n)-U;
end
L=triu(L);
% 求X
X = zeros(n,1);
Y = zeros(n,1);
for j=1:nY(j) = (b(j)-L(j,1:j-1)*Y(1:j-1))/L(j,j);
end
L = L';
for k=n:-1:1X(k)=(Y(k)-L(k,k+1:n)*X(k+1:n))/L(k,k);
end

控制台输入：

>> A = [2,1,-5,1;1,-5,2,7;0,2,1,-1;1,7,-1,-4];
>> b = [13;-9;6;0];
>> [L,X] = Cholesky(A,b);

MATLAB Cholesky分解相关推荐

Matlab与线性代数--矩阵的Cholesky分解
本图文介绍了Matlab对正交矩阵的Cholesky分解操作.
matlab 提示秩 1,Cholesky 分解的秩 1 更新
cholupdate Cholesky 分解的秩 1 更新语法 R1 = cholupdate(R,x) R1 = cholupdate(R,x,'+') R1 = cholupdate(R,x,' ...
矩阵的Cholesky分解的Matlab实现
版权声明:本文为博主原创文章,未经博主允许不得转载.https://blog.csdn.net/weixin_38451800/article/details/88933683 1.Cholesky分 ...
matlab矩阵分解
矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积.常见的矩阵分解有LU分解(三角分解).QR分解(正交变换).Cholesky分解,以及Schur分解.Hessenberg分解.奇 ...
矩阵的Cholesky分解
原文链接:矩阵的Cholesky分解首先来复习线性代数中几个重要的概念. 1)如果一个复矩阵A = A*(共轭转置),则A称为Hermitian矩阵.(注意,矩阵A转置后仍为其本身,显然A一定是方阵 ...
【数值分析】Doolittle分解和Cholesky分解的Python实现
Doolittle 分解 import numpy as np # A=[[1.0,2.0,-3.0], # [2.0,-1.0,3.0], # [3.0,-2.0,2.0]]# A=[ [2, 4, ...
matlab 矩阵分解
矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积.常见的矩阵分解有LU分解.QR分解.Cholesky分解,以及Schur分解.Hessenberg分解.奇异分解等. (1) LU ...
3阶Hermitian正定矩阵Cholesky分解通用表达式
pdf文件算法原理将一个nnn阶Hermitian正定矩阵AAA分解为一个下三角矩阵LLL及其转置LTL^TLT的乘积,即: A=LLT(1)A=LL^T\tag{1} A=LLT(1) 推导过程 ...
线性代数笔记： Cholesky分解
1 介绍当一个实矩阵A是对称正定矩阵的时候,它可以分解成一个下三角矩阵L以及它的转置的乘积,即: 1.1 矩阵半正定的情况如果矩阵是正定的话,那么L唯一确定:如果矩阵是半正定的话,那么也可以分解, ...
cholesky分解java代码_Cholesky 分解（转）
Cholesky 分解是把一个对称正定的矩阵表示成一个下三角矩阵L和其转置的乘积的分解. 它要求矩阵的所有特征值必须大于零,故分解的下三角的对角元也是大于零的. Cholesky分解法又称平方根法,是 ...

MATLAB Cholesky分解

MATLAB Cholesky分解相关推荐

最新文章

热门文章