高斯分解LU矩阵之MATLAB实现

高斯分解LU矩阵之MATLAB实现

数值代数是从这个星期开起来的，但是由于老师的影响感觉貌似已经上了两个星期一样，其实才上三大节课，这个星期把绪论和第一章的第一节：三角形方程组和三角分解讲完了。老师说：“加强编程能力”。20221111

1. 高斯变换

函数

function [L,U]=jiaLU(A)

n=length(A);

for k=1:n-1

    A(k+1:n,k)=A(k+1:n,k)/A(k,k);

    A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);

end

L=eye(n)+tril(A,-1);  %tril是对下三角矩阵的抽取，抽取到主对角线下一条对角线

U=triu(A);  %triu是对上三角矩阵的抽取

1. 运行窗口
   就是自己随便找特殊的矩阵

 A=[1 2 3;4 5 6;7 8 9]
jiaLU(A)

运行结果

感悟：20221111对于编程来说还是需要把一般性的过程写出来才能够实现。