matlab中矩阵lu分解复习进程.docVIP

MATLAB中矩阵LU分解 __________________________________________________ __________________________________________________ __________________________________________________ __________________________________________________ __________________________________________________ __________________________________________________ 题目 编写实现对N阶非奇矩阵A进行LU分解的程序。 算法组织 若n阶方阵的各阶顺序主子行列式不为零则存在唯一的单位上三角矩阵L 和上三角矩阵L式的A=LU。其基本思想是GAUSS消去法。参照《计算方法》第38页L、U各项计算公式编写公式。 输入带分解矩阵A For i=1，2，……n 2.1 将L对角线元素赋值L（i，i）=1； For j=1，2，……n 将U 第一行元素赋值U（1，j）=A（1，j）； For k=2，……n 将L第一列元素赋值L(k,1)=A(k,1)/U(1,1); For i=2，……n For j=i，……n 5.1.1 5.2 For k=i+1，……n 5.2.1 程序实现 clear all clc A=input(请输入一个方阵 );%输入一个n阶方阵 [n,n]=size(A); L=zeros(n,n); U=zeros(n,n); for i=1:n %将L的主对角线元素赋值1 L(i,i)=1; end for j=1:n %求矩阵U的第一行元素 U(1,j)=A(1,j); end for k=2:n %求矩阵L的第一列元素 L(k,1)=A(k,1)/U(1,1); end for i=2:n %求L、U矩阵元素 for j=i:n s=0; for t=1:i-1 s=s+L(i,t)*U(t,j); end U(i,j)=A(i,j)-s; end for k=i+1:n r=0; for t=1:i-1 r=r+L(k,t)*U(t,i); end L(k,i)=(A(k,i)-r)/U(i,i); end end %输出矩阵L、U L U 计算与结果分析 本题中根据GAUSS消去算法对任意一各阶 顺序主子式不为零的n阶方阵A进行L分解。 并对《计算方法》第二章42页例2.2进行了 计算。计算结果如右图。 该计算结果与实际分解结果相同，该算法 可靠。但该算法没有考虑A顺序主子行列式 值为0的情况，因此还有欠缺。