计算方法分解及列主元高斯消去解线性方程组.doc

计算方法分解及列主元高斯消去解线性方程组.doc

  1. 1、本文档共5页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
计算方法分解及列主元高斯消去解线性方程组

北京信息科技大学实验报告 课程名称: 数值分析 实验项目: LU分解及列主元高斯消去解线性方程组 学 院: 光电信息与通信工程学院 专 业: 通信工程 姓名/学号: 赖裕东 2009010476 班 级: 通信0904 实验时间: 2010.12.28 成 绩: 指导教师: 孟 祥 花 实验报告内容 实验目的与要求(实验题目) 用LU分解及列主元高斯消去法解线性方程组 {10 -7 0 1 {X1 {8 -3 2.099999 6 2 X2 5.900001 5 -1 5 -1 X3 = 5 2 1 0 2} X4} 1} 输出AX=b中系数A=LU分解的矩阵L及U,解向量X及detA;列主元法的行交换次序,解向量X及detA;比较两种方法所得的结果。 模型建立(相关主要计算公式) 高斯消去法 三、实验过程、步骤(程序) #includeiostream.h #includemath.h #includeconio.h #includeiomanip.h class glz { public: int row,col,i,j,k,t; double mat [100][100]; double tmp; public: void initialize(); void solve(); }; void glz :: initialize() { cout数值分析课 高斯列主元消去法 endl; cout输入方程组的个数endl; cinrow; cout请输入 未知数个数+1 endl; cincol; cout按行输入各未知数前系数 和 等号右边常数 (按行输入哦~~~)endl; for(int i=0; irow; i++) for(int j=0; jcol; j++) cinmat[i][j]; } void glz::solve() { double x[100]; for(i=0; icol-2; i++) { k=i; double max=mat[k][i]; for(j=i+1; jcol-1; j++) { if(abs(mat[j][i])abs(mat[k][i])) { max=mat[j][i]; k=j; } } if (k!=i) { for(j=i; jcol; j++) { double tmp=mat[k][j]; mat[k][j]=mat[i][j]; mat[i][j]=tmp; } } for (j=i+1; jcol-1; j++) { double tmp=mat[j][i]/mat[i][i]; for(t=i+1; tcol; t++) { mat[j][t]-=tmp*mat[i][t]; } } } x[col-2]=mat[col-2][col-1]/mat[col-2][col-2]; for (i=col-3; i=0; i--) { for(j=col-2; j=i+1; j--) mat[i][col-1]-=x[j]*mat[i][j]; x[i]=

文档评论(0)

asd522513656 + 关注
实名认证
内容提供者

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

1亿VIP精品文档

相关文档