计算三角分解 matlab,科学网—数值分析----三角分解matlab程序 - 殷春武的博文

%%%程序编写者  西北工业大学自动化学院    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程序