常州大学数值分析第二章作业.docVIP

下载本文档

220
0
约3.94千字
约 5页
2018-10-12 发布于河北
举报
版权申诉

常州大学数值分析第二章作业.doc

1、本文档共5页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

常州大学数值分析第二章作业

专业：材料学姓名：李松彦学号第第 PAGE \* MERGEFORMAT 1 页共 NUMPAGES \* MERGEFORMAT 4 页用高斯消去法求解下列方程组. 解：将方程组写成AX=b 根据下列公式计算：计算结果：解：解：将方程组写成AX=b 根据下列公式计算：计算结果：用LU分解法解方程组。解：将方程组写成AX=b 根据下列公式计算：计算A的LU分解结果：求解Ly=b 得到：求解Ux=y 解得原方程组的解为：设求解：已知计算及。解：由计算得: 根据下列公式计算：经计算得：由方程构造迭代格式取初始值用上述迭代格式进行迭代五次，并讨论迭代法的收敛性。解：根据迭代公式计算得到：次数 0 1 2 3 4 5 X1 0 0.5 0.75 0.815 0.8595 0.86555 X2 0 0.5 0.65 0.735 0.7525 0.76855 X3 0 0.5 0.75 0.815 0.8595 0.86555 用U1*y编写LU分解法、改进平方根法、追赶法的Matlab程序，并进行相关数值实验?。解：将矩阵进行LU分解。解：用Doolittle法分解，matlab程序如下 A=[1 0 2 0;0 1 1 1;2 0 -1 1;0 0 1 1]; n=size(A,1);flag=1; L=eye(n);U=zeros(n,n); if abs(A(1,1))eps U(1,1:n)=A(1,1:n); L(2:n,1)=A(2:n,1)/U(1,1); else flag=0; return end for k=2:n; if abs(A(k,k))eps 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); else flag=0; return end end L U 执行后： L = 1.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.5000 0.0000 1.0000 0.0000 0.0000 0.0000 0.4000 1.0000 U = 2.0000 0.0000 -1.0000 1.0000 0.0000 1.0000 1.0000 1.0000 0.0000 0.0000 2.5000 -0.5000 0.0000 0.0000 0.0000 1.2000 用crout分解法，matlab程序如下： A=[1 0 2 0;0 1 1 1;2 0 -1 1;0 0 1 1]; n=size(A,1);flag=1; U=eye(n);L=zeros(n,n); if abs(A(1,1))eps L(1:n,1)=A(1:n,1); U(1,2:n)=A(1,2:n)/L(1,1); else flag=0; return end for k=2:n; if abs(A(k,k))eps for r=2:k L(k,r)=A(k,r)-L(k,1:r-1)*U(1:r-1,r); end for i=k+1:n U(k,i)=(A(k,i)-L(k,1:k)*U(1:k,i))/L(k,k); end else flag=0; return end end L U 执行后： L = 1.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 2.0000 0.0000 -5.0000 0.0000 0.0000 0.0000 1.0000 1.2000 U = 1.0000 0 2.0000 0 0 1.0000 1.0000 1.0000 0 0 1.0000 -0.2000