分解
实验日期
同组人员 指导教师
得分
一、实验目的
1.了解三角分解与 doolittle 分解的基本概念及其性质;
2.掌握 doolittle 分解算法及 matlab 实现;
3.学习、掌握 MATLAB 软件有关的命令。
二、实验准备
熟悉 doolittle 分解,了解 MATLAB 矩阵的输出。
三、实验内容
利用 doolittle 分解下面矩阵
for i=2:n
for j=i:n
z=0;
for r=1:i-1
z=z+L(i,r)*U(r,j);
end
U(i,j)=A(i,j)-z;
end
if abs(U(i,i))<eps
&ag='failure'
return;
end
for k=i+1:n
m=0;
for q=1:i-1
m=m+L(k,q)*U(q,i);
end
L(k,i)=(A(k,i)-m)/U(i,i);
end
end
L
U
end
2、在 command window 窗口中输入以下程序代码指令,回车一下即可。
Doolittle([3 2 2 4;7 3 2 7;1 3 6 9;2 5 9 2])
五、实验结果