%%%程序编写者 西北工业大学自动化学院 Email: yincwxa2013@mail.nwpu.edu.cn
%% All rights reserved
clear
a=[4,5,6;3,6,7;8,6,9]
n=length(a(1,:));
for j=1:n
U(1,j)=a(1,j);
L(j,1)=a(j,1)/U(1,1);
end
for k=1:n
for j=1:n
t=0;
for m=1:k-1
t=t+L(k,m)*U(m,j);
end
U(k,j)=a(k,j)-t;
end
for i=k+1:n
t=0;
for m=1:k-1
t=t+L(i,m)*U(m,k);
end
L(i,k)=a(i,k)-t;
end
end
for i=1:n
L(i,i)=1;
end
L
U
%a =
%4 5 6
%3 6 7
%8 6 9
%L =
%1 0 0
%3 1 0
%8 -34 1
%U =
%4 5 6
%-9 -9 -11
%-330 -340 -413
转载本文请联系原作者获取授权,同时请注明本文来自殷春武科学网博客。
链接地址:http://blog.sciencenet.cn/blog-748007-751155.html
上一篇:数值分析----三对角方程组的追赶法matlab程序
下一篇:数值分析---雅克比求特征值matlab程序