矩阵的Cholesky分解的Matlab简洁实现

做完了数值分析的作业，然后上csdn发现大部分人写的都很复杂，于是想着自己也发一个简洁一点的，更易理解

function [outputArray] = chuliesijifenjie(inputArray)
n=size(inputArray,1);%得到输入矩阵的阶数
outputArray=zeros(n);%清零输出矩阵
outputArray(1,1:n)=inputArray(1,1:n);%输出矩阵的第一行先为输入矩阵的第一行
outputArray(2:n,2:n)=inputArray(2:n,2:n);%输出矩阵的n-i阶下主子式先为输入矩阵的n-i阶下主子式
for i=1:1:noutputArray(i,i)=sqrt(outputArray(i,i));%第i行的主元开方u=zeros(1,n-i);%每一次使用u之后都要清空ufor j=i+1:1:noutputArray(i,j)=outputArray(i,j)/outputArray(i,i);%对第i行的主元后面的元素都除于主元u(1,j-i)=outputArray(i,j);%将第i行的主元后面的元素都存在一个行向量里end
U=u'*u;%求外积
outputArray(i+1:n,i+1:n)=outputArray(i+1:n,i+1:n)-U;%在输出矩阵的n-i阶下主子矩阵中减去外积U
end
outputArray=triu(outputArray);%取上三角
end

矩阵的Cholesky分解的Matlab简洁实现相关推荐

矩阵的Cholesky分解的Matlab实现
版权声明:本文为博主原创文章,未经博主允许不得转载.https://blog.csdn.net/weixin_38451800/article/details/88933683 1.Cholesky分 ...
矩阵的Cholesky分解
原文链接:矩阵的Cholesky分解首先来复习线性代数中几个重要的概念. 1)如果一个复矩阵A = A*(共轭转置),则A称为Hermitian矩阵.(注意,矩阵A转置后仍为其本身,显然A一定是方阵 ...
matlab qr分解作用,MATLAB论文_矩阵的QR分解及其MATLAB实现.doc
您所在位置:网站首页 > 海量文档 &nbsp>&nbsp计算机&nbsp>&nbspmatlab MATLAB论文_矩阵的QR分解及其MATLAB实 ...
Matlab与线性代数--矩阵的Cholesky分解
本图文介绍了Matlab对正交矩阵的Cholesky分解操作.
matlab将矩阵分解成lu,10行代码实现矩阵的LU分解（matlab）
最近由于数值分析实验课要求,需要通过matlab实现矩阵的LU分解.但是看了很多网友写的程序,基本上都是通过循环嵌套循环来实现矩阵的LU分解.略感琐碎,因此最近两天便一直在思考能否利用矩阵的乘v法,来 ...
对矩阵进行QR分解的Matlab代码
摘自Introduction to Linear Algebra by Gilbert Strang 结合课后习题进行分析命名初值 n = 3; a = [2;2;-1]; b = [0;-3;3] ...
Cholesky分解（Matlab代码实现）
三十分钟理解：矩阵Cholesky分解，及其在求解线性方程组、矩阵逆的应用
写一篇关于Cholesky分解的文章,作为学习笔记,尽量一文看懂矩阵Cholesky分解,以及用Cholesky分解来求解对称正定线性方程组,以及求对称正定矩阵的逆的应用. 文章目录直接Choles ...
C++实现Cholesky分解
题目: 编制程序求解矩阵 A 的 Cholesky 分解,并用程序求解方程组 Ax=b,其中代码实现: #include <iostream> #include <math.h&g ...
matlab矩阵特征分解,用MATLAB实现矩阵分解
MATLAB求解线性方程的过程基于三种分解法则: (1)Cholesky分解,针对对称正定矩阵: (2)高斯消元法, 针对一般矩阵: (3)正交化, 针对一般矩阵(行数≠列数) 这三种分解 ...

矩阵的Cholesky分解的Matlab简洁实现

矩阵的Cholesky分解的Matlab简洁实现相关推荐

最新文章

热门文章