matlab中直接三角分解法,三角分解及牛顿插值matlab算法

一、对非奇异矩阵A, 作A的Doolittle三角分解。

%L为单位下三角阵，U为上三角阵,A为紧凑格式输出矩阵

clear;

A=[9 18 9 -27;18 45 0 -45;9 0 126 9; -27 -45 9 135];%输入非奇异矩阵A=[9 18 9 -27;18 45 0 -45;9 0 126 9; -27 -45 9 135];

n=size(A,1);矩阵A的行数

A(2:n,1)=A(2:n,1)./A(1,1);%求解第一列，L(i,1)

for i=2:n-1

A(i,i)=A(i,i)-sum(A(i,1:i-1)'.*A(1:i-1,i));%求U的对角线上的值

for j=i+1:n

A(i,j)=A(i,j)-sum(A(i,1:i-1)'.*A(1:i-1,j));%求上三角的值

A(j,i)=(A(j,i)-sum(A(j,1:i-1)'.*A(1:i-1,i)))/A(i,i);%求解下三角的值end

end

A(n,n)=A(n,n)-sum(A(n,1:n-1)'.*A(1:n-1,n));%A(4，4)值

A;%输出三角分解后的A矩阵

U=A;L=A;

for i=1:n

L(i,i)=1;

end

for i=1:n-1

for j=i+1:n

L(i,j)=0;

end

End

L %下三角阵

for i=2:n

for j=1:i-1

U(i,j)=0;

end