对系数矩阵A进行LU分解,没有采用原地存储方式,没有进行选主元
function ex2_7_1
%对系数矩阵A进行LU分解【自己写的版本,没有采用原地存储方式,没有进行选主元】
%将L矩阵和U矩阵输出,
A = [10 -7 0;
-3 2 6;
5 -1 5];
%对系数矩阵A进行LU分解,最终得到L 和 U矩阵
%初始化L和U矩阵.
L = [0 0 0;
0 0 0;
0 0 0]
U = [0 0 0;
0 0 0;
0 0 0]
[n,n] = size(A);
for k = 1:1:n
i = k+1:n;
%计算乘子,将乘子填入到L矩阵的第(2、3)行,第1列
L(i,k) = A(i,k)/A(k,k);
%更新系数矩阵A,就是将A的第一行乘以相应的乘子的相反数,加到第二行和第三行
j = k:n;
A(i,j) = A(i,j) - L(i,k)*A(k,j);
end
U = A;
L
U