?? nalu.m
字號(hào):
function [l,u]=nalu(a)
% 用途:求可逆方陣的Doolittile分解。
% 格式:[l,u]=nalu(a) 其中a為可逆方陣,l返回單位下三角矩陣,u返回上三角矩陣。
n=length(a);
u=zeros(n,n);l=eye(3,3);
u(1,:)=a(1,:);l(2:n,1)=a(2:n,1)/u(1,1);
for k=2:n
u(k,k:n)=a(k,k:n)-l(k,1:k-1)*u(1:k-1,k:n);
l(k+1:n,k)=(a(k+1:n,k)-l(k+1:n,1:k-1)*u(1:k-1,k))/u(k,k);
end
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -