矩阵的三角分解法之LU分解之Doolittle分解

function [L,U]=Doolittle(A) % 矩阵的三角分解法之LU分解之Doolittle分解 A=LU % Doolittle分解：LU分解中L为单位下三角阵，U为上三角阵 % 说明：n阶(n>=2)矩阵A有惟一杜里特尔分解的充要条件是：A的前n-1个顺序主子式不为0 % 编程思想：使用待定系数法 % % 作者：野渡无人 % 最后修改日期：2008.4.9 % % >> A=[2 1 2;4 5 4;6 -3 5] % % A = % % 2 1 2 % 4 5 4 % 6 -3 5 % % >> [L,U]=Doolittle(A) % % L = % % 1 0 0 % 2 1 0 % 3 -2 1 % % % U = % % 2 1 2 % 0 3 0 % 0 0 -1 % % >> L*U-A % % ans = % % 0 0 0 % 0 0 0 % 0 0 0

[n,m]=size(A); if n~=m error('请输入方阵'); end L=eye(n); % 初始化L为单位矩阵，即根据Doolittle分解的特性待定L的对角线上的元素为1 U=zeros(n); % 初始化U为全零矩阵 for i=1:n U(1,i)=A(1,i); end for j=2:n L(j,1)=A(j,1)/U(1,1); end for k=2:n for j=k:n U(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j); end for i=k+1:n L(i,k)=(A(i,k)-L(i,1:k-1)*U(1:k-1,k))/U(k,k); end end

矩阵的三角分解法之LU分解之Doolittle分解相关推荐

矩阵三角分解法（LU分解）
矩阵分解法是高斯消元法的变形,它的复杂度和高斯消元法一样都是O(n^3),但是矩阵分解法在处理线性方程组系(具有相同的系数矩阵,但是右端项不同的方程组)时,运算比较方便. 下面是矩阵分解原理的原理: ...
矩阵的三角分解(LU分解)
矩阵的三角分解将矩阵分解成一个下三角矩阵L和一个上三角矩阵U的乘积. 定义:如果n阶矩阵A能够分解成一个下三角矩阵L和一个上三角矩阵U的乘积,则称这种分解为三角分解或LU分解,如果n阶矩阵A能够分解为 ...
【Matlab】矩阵三角分解法求解方程组
%直接三角分解法 function my_LU(a, b)n = length(a);l = zeros(n, n);%初始化u = zeros(n, n);for i=1:nl(i,i) = 1;e ...
从高斯消元到矩阵的三角分解（LU）
本文主要介绍以 Gauss 消去法为根据导出的矩阵的三角(或 LU)分解. 为建立矩阵的三角分解理论,假定化 A\boldsymbol{A} 为
Python实现部分主元法下LU分解
Python实现部分主元法下LU分解 ''' <矩阵分析与应用>小作业1 实现部分主元法下的LU分解 by苗栋程序大体介绍: 引入了numpy便于对数组的操作 ①寻找出一列中绝对值最大的 ...
用MATLAB实现plu分解,编制计算给定矩阵 A 的 LU 分解和 PLU 分解的通用程序
用VB编写一个程序,计算出给定的10*10矩阵(存放在二维数组A中)每行元素的最大值和每列元素的最小值 ModuleModule1SubMain()DimA(,)AsInteger={{1,2,3,4 ...
解方程AX=b与矩阵分解：奇异值分解（SVD分解）特征值分解 QR分解三角分解 LLT分解
文章目录 1. 前言 2. LU三角分解 3. Cholesky分解 - LDLT分解 4. Cholesky分解 - LLT分解 5. QR分解 6. 奇异值分解 7. 特征值分解 1. 前言本博 ...
视觉SLAM中的数学——解方程AX=b与矩阵分解：奇异值分解（SVD分解）特征值分解 QR分解三角分解 LLT分解
前言本博客主要介绍在SLAM问题中常常出现的一些线性代数相关的知识,重点是如何采用矩阵分解的方法,求解线性方程组AX=B.主要参考了<计算机视觉--算法与应用>附录A以及Eigen库的方 ...
机器学习（十一）——机器学习中的矩阵方法（1）LU分解、QR分解
http://antkillerfarm.github.io/ 因子分析的EM估计(续) 去掉和各参数无关的部分后,可得: ∑i=1mE[logp(x(i)|z(i);μ,Λ,Ψ)]=∑i=1mE[1 ...

矩阵的三角分解法之LU分解之Doolittle分解

矩阵的三角分解法之LU分解之Doolittle分解相关推荐

最新文章

热门文章