2025年数值分析线性方程组求解高效直接法实践探索报告_第1页
2025年数值分析线性方程组求解高效直接法实践探索报告_第2页
2025年数值分析线性方程组求解高效直接法实践探索报告_第3页
2025年数值分析线性方程组求解高效直接法实践探索报告_第4页
2025年数值分析线性方程组求解高效直接法实践探索报告_第5页
已阅读5页，还剩12页未读， 继续免费阅读

下载本文档

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

文档简介

计算措施试验二试验汇报专业班级：姓名：学号：试验成绩：1.【试验题目】2.【试验目的】●掌握三角分解法与追赶法3.【试验内容】。4.【试验规定】(2)观测小主元并分析对计算成果的影响。赶法公式。A=G₁/b₁β=c/(b-aβ)y=(f-ay_)/(b-aβ_₁)a6.【源程序(带注释)】(1)一：列主元素消去法(cout<<"下面请输入未知数的个数m=";cout<<"请按次序输入增广矩阵a:"<<endl;c=(fabs(a[j][i])>fabs(a[i]/*找列最大元素*/}/*将列最大数防在对角线上*/a[k][r]=a[k][r]-1[k][i]}x[m-1]=a[m-1][m]/a[m-1][x[i]=(a[i][m]-d)/a[icout<<"x["<<i<<"]="<<}j<m+1;j++)}一般消去法voidsolve(floatI[][100],floatu[][100],flfor(i=1;i<=n;i++)/*第一次回代过程开始*/{{y[i]=(b[i]+s1)/L[i]}for(i=n;i>=1;i--)/*第二次回代过程开始*/}}}{floata[100][100],I[100][100],u[100][100],x[100],b[100];for(i=1;i<=99;i++)/*将所有的数组置零，同步将L矩阵的对角值设为1*/{}printf("inputn:\n");/*输入方程组的个数*/printf("inputarrayA:\n");/*读取原矩阵A*/for(r=1;r<=n;r++)/*求解矩阵L和U*/{{}{}printf("arrayL:\n");/*输出矩阵L*/for(i=1;i<=n;i++)printf("%7.3f",I[i][j]);}printf("arrayU:\n");/*输出矩阵U*/{printf("%7.3f",u[i]lj]);printf("\n");}printf("}2)消元过程中，即使,用其作除数会导致计算中间结果数量级严重增长和舍入误差的累积、扩大，最后使得计算结果不可靠。3)应避免采用绝对值很小的主元素ak);对一般的系数矩阵，最好保持乘数|≤1,因此，在高斯消去法中应引进选主元技巧，以便减少计算过程中舍入误差对求解的影响。doublea[N],b[N],c[N-1],f[N],voidLUDecomposevoidbackSubs();//回代{{printf("由于该维数过于犀利，导致程序退出!");}printf("\n输入下三角元素\n");printf("\n输入主对角线元素\n");printf("\n输入上三角元素\n");printf("x%

人人文库> 全部分类> 教育资料 > 英语等级

温馨提示

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

最新文档

评论

0/150

提交评论