矩阵计算-MATLAB-doolittle分解.doc

《矩阵计算-MATLAB-doolittle分解.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

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

2、A)L=eye(n);U=zeros(n);flag='ok';fori=1:nU(1,i)=A(1,i);endforr=2:nL(r,1)=A(r,1)/U(1,1);Endfori=2:nforj=i:nz=0;forr=1:i-1z=z+L(i,r)*U(r,j);endU(i,j)=A(i,j)-z;endifabs(U(i,i))

3、ommandwindow窗口中输入以下程序代码指令，回车一下即可。Doolittle([3224;7327;1369;2592])五、实验结果五、实验心得:通过本次实验我会利用doolittle分解成LU形式是有条件的，首先矩阵必须是非奇异的矩阵，其次矩阵的全部顺序主子式非零的时候才能完全保证矩阵可分解成LU且分解唯一。Doolittle分解法是不选主元的三角分解法，若U（k,k）绝对值很小时，按此分解法可能引起舍入误差的积累。可以通过交换矩阵A的行实现矩阵PA的LU分解。在编写程序的时候也遇到了好多问题，最终在我坚持不懈的努力下，一个个问题被我解决了。本次试验增

4、加了我的MATLAB编程能力。