求下三角矩阵的逆矩阵的详细算法
4星 · 超过85%的资源 需积分: 48 78 浏览量
2012-05-21
16:44:38
上传
评论 7
收藏 56KB DOC 举报 
1. 求下三角矩阵的逆矩阵的详细算法。
解：
function B=nijuzhen(A);
a=input('输入一个整数');
A=unifrnd(0,1000,a,a)
A=tril(A);
n=size(A,1)
B=eye(n);
if rank(A)<n
return
disp('the matrix is singular try another one!')
else
for k=1:n
B(1,k)=B(1,k)/A(1,1);
for i=2:n
B(i,k)=(B(i,k)-A(i,1:i-1)*B(1:i-1,k))/A(i,i);
end
end
end
运行结果：
评论8