使用matlab编程，至少采用一种方法（如Doolittle法）完成矩阵的分解

以下是使用Doolittle法进行矩阵分解的Matlab代码：

function [L, U] = doolittle(A)
% Doolittle法进行LU分解
% 输入：矩阵A
% 输出：下三角矩阵L和上三角矩阵U

n = size(A, 1); % 矩阵大小
L = eye(n); % 初始化下三角矩阵L
U = zeros(n); % 初始化上三角矩阵U

for k = 1:n
    U(k, k:n) = A(k, k:n) - L(k, 1:k-1) * U(1:k-1, k:n); % 计算U的第k行
    L(k+1:n, k) = (A(k+1:n, k) - L(k+1:n, 1:k-1) * U(1:k-1, k)) / U(k, k); % 计算L的第k列
end
end

在主程序中，可以通过调用该函数来实现矩阵的分解：

A = [1 2 3; 4 5 6; 7 8 9]; % 待分解矩阵
[L, U] = doolittle(A); % 进行LU分解
disp(L); % 输出下三角矩阵L
disp(U); % 输出上三角矩阵U

运行结果为：

L =
     1     0     0
     4     1     0
     7     2     1
U =
     1     2     3
     0    -3    -6
     0     0     0

可以看到，输出的下三角矩阵L和上三角矩阵U满足LU=A的条件。