matlab求doolittle,矩阵计算-MATLAB-doolittle分解

《矩阵计算-MATLAB-doolittle分解》由会员分享，可在线阅读，更多相关《矩阵计算-MATLAB-doolittle分解(1页珍藏版)》请在人人文库网上搜索。

1、姓 名 学 号 专业班级 课程名称 计算矩阵 实验名称 doolittle分解 实验日期 同组人员 指导教师 得分 一、实验目的1了解三角分解与doolittle分解的基本概念及其性质；2掌握doolittle分解算法及matlab实现；3学习、掌握MATLAB软件有关的命令。二、实验准备熟悉doolittle分解，了解MATLAB矩阵的输出。三、实验内容利用doolittle分解下面矩阵四、实验步骤1、在实验时，双击matlab软件，首先建立M文件。具体操作：FilenewM-file粘贴代码保存。代码如下：function f=Doolittle(A)m,n=size(A)L=eye(n)。

2、;U=zeros(n);flag=ok;for i=1:nU(1,i)=A(1,i);endfor r=2:nL(r,1)=A(r,1)/U(1,1);Endfor i=2:nfor j=i:nz=0;for r=1:i-1z=z+L(i,r)*U(r,j);endU(i,j)=A(i,j)-z;endif abs(U(i,i)epsflag=failurereturn;endfor k=i+1:nm=0;for q=1:i-1m=m+L(k,q)*U(q,i);endL(k,i)=(A(k,i)-m)/U(i,i);endendL Uend2、在command window窗口中输入以下程序代码指令，回车一下即可。Doolittle(3 2 2 4;7 3 2 7;1 3 6 9;2 5 9 2)五、实验结果五、实验心得: 通过本次实验我会利用doolittle分解成LU形式是有条件的，首先矩阵必须是非奇异的矩阵，其次矩阵的全部顺序主子式非零的时候才能完全保证矩阵可分解成LU且分解唯一。Doolittle分解法是不选主元的三角分解法，若U(k,k)绝对值很小时，按此分解法可能引起舍入误差的积累。可以通过交换矩阵A的行实现矩阵PA的LU分解。在编写程序的时候也遇到了好多问题，最终在我坚持不懈的努力下，一个个问题被我解决了。本次试验增加了我的MATLAB编程能力。