Cholesky分解算法Matlab代码

n=input('请输入矩阵阶数：\n');
A=input('请输入待求解的对称正定矩阵:\n');
for k =1:n-1A(k,k)=sqrt(A(k,k))A(k+1:n,k)=A(k+1:n,k)/A(k,k);for j =k+1:nA(j:n,j)=A(j:n,j)-A(j:n,k)*A(j,k);end
end
A(n,n)=sqrt(A(n,n));
A

function Cholesky%n=input('请输入矩阵阶数：\n');A=input('请输入待求解的对称正定矩阵:\n');[N,N]=size(A);X=zeros(N,1);Y=zeros(N,1);for i=1:NA(i,i)=sqrt(A(i,i)-A(i,1:i-1)*A(i,1:i-1)');if A(i,i)==0'A is singular, no unique solution';break;endfor j=i+1:NA(j,i)=(A(j,i)-A(j,1:i-1)*A(i,1:i-1)')/A(i,i);endend
for x=1:Nfor y=1:NB(x,y)=A(x,y);      %右上角元素归0if x<yB(x,y)=0;break;endend
end
L=B

Cholesky分解算法Matlab代码相关推荐

【MIMO】两种空间相关信道生成方式的记录（公式+MATLAB代码）
文章目录前言一. Kronecker相关信道模型二.生成方式1 1.公式 2.MATLAB代码三.生成方式2-complex correlation 1.公式 2.MATLAB代码四.生成方 ...
【熵与特征提取】基于“信息熵”的特征指标及其MATLAB代码实现（功率谱熵、奇异谱熵、能量熵）
<三体>中对监听员的日常工作有这样一段描述: ...1379 号监听站已经存在了上千年,像这样的监听站,在三体世界中有⼏千个,它们全神贯注地聆听着宇宙间可能存在的智慧⽂明的信息... .. ...
计及新能源出力不确定性的电气设备综合能源系统协同优化（Matlab代码实现）
目录第一部分文献一<计及新能源出力不确定性的电气设备综合能源系统协同优化> 0 引言 1 新能源出力不确定性处理 1.1 新能源出力预测误差分布 1.2 新能源出力的时间相关性 1.3 ...
【信号分解】基于LMD算法和ELMD算法实现管道泄漏信号处理附matlab代码
1 内容介绍在科技水平相当发达的今天,互联网+.大数据慢慢渗透进人们的生活当中,但科技的进步不仅仅要体现在生活质量水平的提高.经济的快速发展,更应该体现在对社会资源的合理利用.自建国以来,我国管道 ...
最优化方法matlab代码(一) 牛顿类方法
代码:最优化方法牛顿方法matlab代码-从零开始-专业指导文档类资源-CSDN文库https://download.csdn.net/download/benchuspx/58807913目录 1. ...
【时间序列模型】AR模型（原理剖析+MATLAB代码）
目录一.AR模型的原理剖析 1.1 AR模型原理 1.2 模型的参数估计 1.3 模型的定阶方法 (1)AIC (2)BIC 二.MATLAB代码三.AR相关论文参考文献: 前言时间序列分析方 ...
龙格-库塔法(runge-kutta)matlab代码及含义,龙格-库塔法(Runge-Kutta)matlab代码及含义...
龙格-库塔法(Runge-Kutta)matlab代码及含义龙格-库塔法(Runge-Kutta) 数值分析中,龙格-库塔法(Runge-Kutta)是用于模拟常微分方程的解的重要的一类隐式或显式迭 ...
arima模型matlab代码_PSTR面板平滑转换模型简介（附Matlab代码分享）
写论文的时候用到的~相关的资料太少了,做一些简单内容和资料的分享.(PSTR模型的Matlab代码分享在最后)本文主要为简单理论和粗暴实操~ 有用的话可以点个赞哟(知乎小白卑微求赞) 嘻嘻下面进入正题 ...
matlab数值分析拟合实例,数值分析函数拟合matlab代码.doc
数值分析函数拟合matlab代码.doc 第一题MATLAB代码用SPLINE作图XI0204060810YI098092081064038X10012Y1NEWTON3XI,YI,X源代码见M文件Y ...

Cholesky分解算法Matlab代码

Cholesky分解算法Matlab代码相关推荐

最新文章

热门文章