矩阵的Doolittle分解(LU分解)

最新推荐文章于 2024-05-02 16:10:04 发布

廣阝

最新推荐文章于 2024-05-02 16:10:04 发布

阅读量4.6k

点赞数

分类专栏：计算机数值方法/数值计算方法文章标签： matlab 线性代数矩阵

本文链接：https://blog.csdn.net/SZU_Kwong/article/details/109563830

版权

计算机数值方法/数值计算方法专栏收录该内容

3 篇文章 0 订阅

订阅专栏

基本三角分解法

设矩阵 $\ A=(a_{ij})_{n \times n}$ 的各阶顺序主子式不为零，A有分解式
$\ \begin{bmatrix} a_{11} \quad a_{12} \quad \cdots \quad a_{1n} \\ a_{21} \quad a_{22} \quad \cdots \quad a_{2n} \\ \vdots \qquad \vdots \qquad \qquad \vdots \\ a_{n1} \quad a_{n2} \quad \cdots \quad a_{nn} \\ \end{bmatrix} = \begin{bmatrix} 1 \qquad \qquad \qquad \\ l_{21} \quad 1 \qquad \qquad \\ \vdots \quad \vdots \quad \ddots \quad \\ l_{n1} \quad l_{n2} \quad \cdots \quad 1 \\ \end{bmatrix} \begin{bmatrix} u_{11} \quad u_{12} \quad \cdots \quad u_{1n} \\ \qquad u_{22} \quad \cdots \quad u_{2n} \\ \qquad \qquad \ddots \qquad \vdots \\ \qquad \qquad \qquad u_{nn} \end{bmatrix} \equiv LU$
可得出关系
$\ a_{ri}= \sum_{k=1}^r l_{rk}u_{ki} (i=r,\cdots,n;r=1,2,\cdots,n)$
$\ a_{ir} = \sum_{k=1}^r l_{ik}u_{kr} (i=r+1,\cdots,n;r=1,\cdots,n-1)$
由这两个关系可以解出L和U：
当r=1时，
$\ u_{1i}=a_{1i} \quad (i=1,2,\cdots,n)$
$\ l_{i1}=\frac{a_{i1}}{u_{11}} \quad (i=2,3,\cdots,n)$
假设求出了U的第1至r-1行，L的1至r-1列，由上面两式得出U，L的第r行，列的计算公式
$\ u_{ri}=a_{ri}-\sum_{k=1}^{r-1}l_{ik}u_{kr} \quad (i=r,\cdots,n;r=1,2,\cdots,n)$
$\ l_{ir}=\frac{a_{ir}-\sum_{k=1}^{r-1}l_{ik}u_{kr}}{u_{rr}} \quad (i=r+1,\cdots,n;r=1,2,\cdots,n-1)$
$\ 规定\sum_{k=1}^0 l_{ik}u_{kr}=0$
称上面所表示的矩阵分解为Doolittle分解，也称LU分解

部分选主元的Doolittle分解

在Doolittle的第r步分解，首先在数组A的第r列主对角元以下（含主对角元）选主元，具体步骤为：
（1） $\ 在r=1,2,\cdots,n-1$ 时
$\ A(i,r)=A(i,r)-\sum_{k=1}^{r-1}A(i,k) \cdot A(k,r) \quad (i=r,\cdots,n)$
（2）挑选绝对值最大的主元，得出对应的行号idx，将第r行的元素与第idx行的元素进行交换
（3）分解计算
$\ A(i,r)=\frac{A(i,r)}{A(r,r)} \quad (i=r+1,\cdots,n;r=1,2,\cdots,n-1)$
当 $\ r \neq n$ 时
$\ A(r,i)=A(r,i)-\sum_{k=1}^{r-1}A(r,k) \cdot A(k,i) \quad (i=r+1,\cdots,n;r=1,2,\cdots,n)$
$\ (\sum_{k=1}^0 A(r,k)\cdot A(k,i)=0)$
当 $\ r=n$ 时
$\ A(r,r)=A(r,r)-\sum_{k=1}^{r-1}A(r,k)\cdot A(k,r)$
按上述方法完成n步分解，便实现了系数矩阵分解 $\ PA=LU$ ，其中， $\ P$ 是排列矩阵
$\ P=I_{i_{n-1}n-1}I_{i_{n-2}n-2} \cdots I_{i_11}$
其中， $I_{i_kk}$ 是单位矩阵在第k步分解中交换了第k行和第 $i_k$ 行(主元中最大绝对值的行数)后的矩阵

Matlab 实现一个矩阵的Doolittle分解

function [L, U, P] = Doolittle(A, pri)
%DOOLITTLE 对矩阵A进行doolittle分解
%   inputs:要分解的矩阵A, pri:是否选主元
%   outputs:单位下三角矩阵L和上三角矩阵U，若选主元，同时输出置换矩阵P

P = eye(size(A));

[row, col] = size(A);

if pri
    for r = 1:row
        % 计算中间量Si，存入A(i,r)
        if r ~= row
            A(r:end, r) = A(r:end, r) - A(r:end, 1:r-1) * A(1:r-1, r);
        end
        % 确定绝对值最大的行号
        [m, idx] = max(abs(A(r:end, r)));
        idx = idx + r - 1;
        % 换行，同时更新排列矩阵P
        I = eye(row, col);
        temp = I(r, :);
        I(r, :) = I(idx, :);
        I(idx, :) = temp;
        A = I * A;
        P = I * P;
        % 分解计算
        if r ~= row
            A(r+1:end, r) = A(r+1:end, r) / A(r, r);
        end
        if r ~= 1
            if r ~= row
                A(r, r+1:end) = A(r, r+1:end) - A(r, 1:r-1) * A(1:r-1, r+1:end);
            else
                A(r, r:end) = A(r, r:end) - A(r, 1:r-1) * A(1:r-1, r:end);
            end
        end
    end
    L = tril(A, -1) + eye(size(A));
    U = triu(A, 0);
else
    U = zeros(size(A));
    L = eye(size(A));
    for r = 1:row
        if r == 1
            U(1, :) = A(1, :);
            L(2:end, 1) = A(2:end, 1) ./ U(1, 1);
        else
            U(r, r:end) = A(r, r:end) - L(r, 1:r-1) * U(1:r-1, r:end);
            L(r+1:end, r) = (A(r+1:end, r) - L(r+1:end, 1:r-1) * U(1:r-1, r)) / U(r, r);
        end
    end
end

end

廣阝

关注

0
点赞
踩
16

收藏

觉得还不错? 一键收藏
0
评论
矩阵的Doolittle分解(LU分解)

基本三角分解法设矩阵 A=(aij)n×n\ A=(a_{ij})_{n \times n} A=(aij)n×n的各阶顺序主子式不为零，A有分解式 [a11a12⋯a1na21a22⋯a2n⋮⋮⋮an1an2⋯ann]=[1l211⋮⋮⋱ln1ln2⋯1][u11u12⋯u1nu22⋯u2n⋱⋮unn]≡LU\ \begin{bmatrix}a_{11} \quad a_{12} \quad \cdots \quad a_{1n} \\a_{21} \quad
复制链接

扫一扫