矩阵LU分解

来源：要发发知识网

%矩阵的LU分解计算

function x=jzlufj(A,b)

n=length(A(:,1));

L=eye(n);

L1=eye(n);

for k=1:n-1

for t=k+1:n

m=A(t,k)/A(k,k);

L1(t,k)=-m;

L(t,k)=m;

A=L1*A;

L1=eye(n);

end

U=A;

y=zeros(n,1);

y(1)=b(1);

for v=2:n

y(v)=b(v)-L(v,1:v-1)*y(1:v-1);

end

x=zeros(n,1);

x(n)=y(n)/U(n,n);

for s=n-1:-1:1

x(s)=(y(s)-U(s,s+1:n)*x(s+1:n))/U(s,s);

end

验证结果

>> A=[1 1 1;0 4 -1;2 -2 1]

A =

1 1 1

0 4 -1

2 -2 1

>> b=[6;5;1]

b =

>> x=jzlufj(A,b)

x =

因篇幅问题不能全部显示，请点此查看更多更全内容

查看全文

违法及侵权请联系：TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务