计算机实习课作业文件_第1页
计算机实习课作业文件_第2页
计算机实习课作业文件_第3页
计算机实习课作业文件_第4页
免费预览已结束，剩余1页可下载查看

付费下载

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、实习一：用高斯消元法及高斯列主元消去法求解线性方程组 序并讨论结果。其中59.1746.78 b=1解：一、高斯消元法：首先，编写.m文件如下：fun cti on X = GaussXiaoyua n(n, A,b)B=A b;for k= 1:nif abs(B(k,k)v=0.1*10八(-4)Bretur n ;elsefor i=k+1:1: nL(i,k)= B(i,k)./ B(k,k);B(i, n+1)= B(i, n+1)-L(i,k).* B(k, n+1); for j=1: nif k+1aa=B(i,k);r=i;endj=1:n+1C(k,j)=B(r,j); B

2、(r,j)=B(k,j); B(k,j)=C(k,j);end for i=k+1:1:n L(i,k)= B(i,k)./ B(k,k) ;B(i,n+1)= B(i,n+1)-L(i,k).* B(k,n+1); for j=1:n ifk+1=j=nB(i,j)=B(i,j)-L(i,k).*B(k,j); else B(i,j)=0; endend end b=B(1:n,n+1) A=B(1:n,1:n)X(n)=b(n)./A(n,n); for q=n-1:-1:1X(q)=(b(q)-sum(A(q,q+1:4).*X(q+1:4)./A(q,q);end end之后在Comm

3、a nd Win dOW输入：A二0.3*10八(-15) 59.14 3 1;5.291 -6.13 -1 2;11.3 9 5 2;1 2 1 1;b = 59.1746.7812;Liezhuyuan(4,A,b)输出结果为：b =1.000046.311859.17001.911511.3000 9.0000 5.0000 2.0000 0-10.3441-3.3412 0 59.14003.0000 0 1.2035 0.55751.06351.00000.82301.000059.170056.66110.707411.3000 9.0000 5.0000 2.000000059.1400 3.0000 1.00000 -2.8164 1.23840 0.4965 0.80271.000059.170056.661110.695411.3000 9.0000 5.0000 2.000000059.1400 3.0000 1.00000 -2.8164 1.23840 0 1.02101.000059.170056.661110.695411.3000 9.0000 5.0000 2.000059.1400 3.0000 1.0000 0 -2.8164 1.2384 0 0 1.0210ans =列主元消去法避

人人文库> 全部分类> 应用文书 > 作业报告

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论