Gauss消元法解线形方程组的程序_第1页
Gauss消元法解线形方程组的程序_第2页
Gauss消元法解线形方程组的程序_第3页
全文预览已结束

下载本文档

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

文档简介

1、数值分析实验报告实验序号： 题目名称：gauss消元法解方程组学号：姓名:任课教师: 专业班级:1、实验目的：编写用gauss消元法解线形方程组的程序2、实验分析：计算方法分析：gauss消元法的基本做法就是把方程组转化成为一个如下图的等价的三角方程组，这个 过程叫做消元。得到三角方程组后，就可以逐个求出xn, xnl,，xi,这个过程叫回代。程序代码分析：建立两个数组a和b,通过循环语句将n阶增广矩阵输入进去，通过对列的循环对每一 列进行消去未知数，通过n小步n大步把矩阵化简成上三角形矩阵，最后通过迭代法解得方 程组得解。3、函数分析：程序中只用到一个主函数，求解线形方程组得算法都放在主函数

2、中， 利用以下函数进行求解：aiu=aiu-aik/akk*akj;bi=bi-aik/akk*bk;迭代:for(i=n-l;i>0;i)double sum=0; for(j=i+l;j<=n;j+) sum=sum+ai u *x j ; xi=(bi-sum)/aii;4、实验数据:int main(int argc.char* argv)int n,i,j,k,l;double a100100,b100j;while(cout«"请输入矩阵的阶数:"«endl&&cin»n) cout«h请输入第

3、增广矩阵：m«endl;fbr(i=l;i<=n;i+)for(j=l;j<=n+l;j+)for(k= 1 ;k<=n-1 ;k+)fbr(i=k+l;i<=n;i+)for(j=k+l ;j<=n+l ;j+)aiu=aiu-(aik/akk)*akju；1aik=o;bn=ann+/ann;for(i=n-l;i>=l;i)double sum=0.0;for(j=n;j>=i+l;j-)sum=sum+aij*bj;bi=(afin + l-sum)/aii;cout«n方程的根为：m«endl;for(int loop=l;loop<=n;loop+) cout«"xm«loop«h = "« bloop&

人人文库> 全部分类> 生活休闲 > 科普知识

温馨提示

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

最新文档

评论

0/150

提交评论