MATLAB报告用LU分解法求解线性方程组.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《MATLAB》期末大报告 报告内容:用LU分解法求解线性方程组 学院(系): 专 业: 班 级: 学 号: 学生姓名: 2014 年 6 月 23日 报告目的 学会用LU分解法解线性方程组,并且为了保证乘子的绝对值小于1,需要对主元数值很小的时候进行方程变换。 报告内容 1 :测试基本线性方程: 可见程序运行结果正确。 2 :测试主元有接近于0的线性方程: 可见程序运行结果正确。 3:测试矩阵维度大于5情况: 可见程序运行结果正确。 源程序 function X=LU1(A,B) B=B; A=[A;B],n=length(B); X=zeros(n,1); y=zeros(n,1); U=zeros(n); L=eye(n); for k=1:n U(1,k)=A(1,k); L(k,1)=A(k,1)/U(1,1); end for i=2:n for k=i:n lu=0; lu1=0; for j=1:i-1 lu=lu+L(i,j)*U(j,k); lu1=lu1+L(k,j)*U(j,i); end U(i,k)=A(i,k)-lu; L(k,i)=(A(k,i)-lu1)/U(i,i); end end L U for i=1:n ly=0; for j=1:i ly=ly+L(i,j)*y(j); end y(i)=B(i)-ly; end for i=n:-1:1 ly1=0; for j=i+1:n ly1=ly1+U(i,j)*X(j); end X(i)=(y(i)-ly1)/U(i,i); end 四、报告分析与心得 MATLAB编程和C编程一样需要仔细的逻辑和十分的细心,自己还有很多需要学习的地方。

文档评论(0)

153****9595 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档