太原理工大学数值计算实验报告材料

《太原理工大学数值计算实验报告材料》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、实用文案本科实验报告课程名称：计算机数值方法实验项目：方程求根,线性方程组的直接求解,线性方程组的迭代求解,代数插值和最小二乘法拟合实验地点：软件楼208专业班级：学号：学生姓名：指导教师：李志崔冬华2015年6月5日标准文档实用文案学生姓名实验成绩实验名称实验一方程求根实验内容和要求（必填）目的：《数值计算方法》是一门实践性和实用性都很强的课程，学习的目的在于应用已有的理论基础和技术方法，实现对实际计算问题的计算机求解。本实验配合《计算数值方法》课堂教学，通过上机实践实验环节，巩固课堂和书本上所学的知识、加强实践

2、能力、提高解决实际计算问题的水平、启发创新思想。其基本目的是：(1)培养理论联系实际的思路，训练综合运用已经学过的理论和实际知识去分析和解决实际问题的能力。(2)帮助学生全面消化已学的相关课程内容，深刻理解计算数值方法课程的内涵，培养使用电子计算机进行科学计算和解决问题的能力。(3)进行基本技能训练和巩固。使学生得到选择算法、编写程序、分析数值结果、写数值试验报告、课堂讨论等环节的综合训练。要求：(1)应用结构化程序设计编出通用程序，源程序要有详细的注释和说明；(2)比较计算结果，分析数值解误差的原因；(3)实验完

3、成，要求提交实验结果并写出报告，分析计算结果是否符合问题的要求，找出计算成功的原因或计算失败的教训。实验原理（必填）选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度

4、x*-xn

5、<0.5×10-5主要仪器设备笔记本计算机实验记录(写出实验内容中的程序代码和运行结果)(可分栏或加页)标准文档实用文案迭代法：#include"stdafx.h"#include"stdio.h"#include"math.h"#include"iostream"usingnames

6、pacestd;floatmain(){floata;cin>>a;floatt,x;x=a;do{x=sqrt((10-x*x*x)/4);t=a;a=x;}while(fabs(a-t)>0.5*1e-5);printf("x=%f",a);system("pause");}割线法:#include"stdafx.h"#include"stdio.h"#include"math.h"#include"iostream"usingnamespacestd;floatmain(){floatc,a=1.0,b=2.

7、0;//cin>>a>>b;while(1){c=b-(b*b*b+4*b*b-10)*(b-a)/(b*b*b+4*b*b-(a*a*a+4*a*a));if(fabs(b-c)<0.5*0.000001)break;b=c;}cout<

8、LU分解法、追赶法。（2）加深对线性方程组求解方法的认识，掌握算法。（3）会进行误差分析，并能对不同方法进行比较。实验原理合理利用Gauss消元法、LU分解法或追赶法求解下列方程组：1、2、3、4、（n=5,10,100,…）主要仪器设备台式或笔记本计算机实验记录(写出实验内容中的程序代码和运行结果)(可分栏或加页)①Gauss消元法：#include"stdafx.h"#include"stdio.h"#include"iostream"usingnamespacestd;标准文档实用文案floatmain(){

9、floata[3][4]={{1,2,3,14},{0,1,2,8},{2,4,1,13}};floatx[3];floatsum=0;intk,i,j;for(k=0;k<2;k++)for(i=k+1;i<3;i++)for(j=k+1;j<4;j++)a[i][j]=a[i][j]-a[i][k]/a[k][k]*a[k][j];for(i=0;i<3;i++)for(j=0;j<4;j++)printf("a[%d][%d]=%f,",i,j,a[i][j]);cout<

10、/a[2][2];for(k=1;k>=0;k--){sum=0;for(j=k+1;j<3;j++){sum+=a[k][j]*x[j];}x[k]=(a[k][3]-sum)/a[k][k];}for(i=0;i<3;i++)printf("x[%d]=%f,",i+1,x[i]);}③LU分解法：#include"stdafx.h"#include