计算机数值计算具体办法.docVIP

本科实验报告 课程名称： 计算机数值方法 实验项目： 二分法 牛顿法 迭代法 实验地点： 学院楼606 专业班级： 学号： 学生姓名： 指导教师： 王丽娟 年 05月 8日 太原理工大学学生实验报告 学院名称 计算机科学与技术 专业班级 学号 学生姓名 实验日期 5月8日 成绩 课程名称 数值计算方法 实验题目 二分法迭代法求解方程 一、实验目的和要求（必填） 熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度|x*-xn|0.5×10-5 实验内容和原理（必填） f(x)在区间（x，y）上连续 　　先找到a、b属于区间（x，y），使f(a)，f(b)异号，说明在区间(a,b)内一定有零点，然后求f[(a+b)/2], 　　现在假设f(a)0,f(b)0,ab ①如果f[(a+b)/2]=0，该点就是零点， 　　如果f[(a+b)/2]0,则在区间（(a+b)/2，b)内有零点，(a+b)/2=a，从①开始继续使用 　　中点函数值判断。 　　如果f[(a+b)/2]0，则在区间(a,(a+b)/2)内有零点，(a+b)/2=b，从①开始继续使用 　　中点函数值判断。 　　这样就可以不断接近零点。 通过每次把f(x)的零点所在小区间收缩一半的方法，使区间的两个端点逐步迫近函数的零点，以求得零点的近似值 三、主要仪器设备（必填） PC,DEV-C 四、操作方法与实验步骤（可选） 1.迭代法 #includestdio.h int main(){ int i=0; double a; double b[100]; printf(请输入任意一个在1与2之间的数\n); scanf(%lf,b[0]); do{ b[i+1]=b[i]-((b[i]*b[i]*b[i]+4*b[i]*b[i]-10)/(3*b[i]*b[i]+8*b[i])); i++; }while(b[i]-b[i+1]0.000005); printf(%.8lf,b[i]); getch(); } 牛顿法 #includestdio.h #includemath.h main() { double x=1.5,a=0; while(pow(pow((x-a),2),0.5)/20.5*0.00001) { a=x; x=0.5*pow(10-pow(x,3),0.5); printf(%f\n,x); } getchar(); } 六、实验结果与分析（必填） 结果与书上的答案一致，编译过程中对do-while循环中条件的存在一定的不合适。。。 实验地点 综合楼六层606室 指导教师 王丽娟 本科实验报告 课程名称： 计算机数值方法 实验项目：高斯消元法，LU分解法，追赶法 实验地点： 学院楼606 专业班级： 学号： 学生姓名： 指导教师： 王丽娟 年 5 月 8日 太原理工大学学生实验报告 学院名称 计算机科学与技术 专业班级 学号 学生姓名 实验日期 5月8日 成绩 课程名称 计算机数值方法 实验题目 高斯消元法，LU分解法，追赶法 实验目的和要求（必填） 2) 线性方程组的直接解法 合理利用Gauss消元法、LU分解法、追赶法求解下列方程组： ① ② ③ ④（n=5,10,100,…） 实验内容和原理（必填） 高斯分解法： = 1 \* GB2 ⑴将原方程组化为三角形方阵的方程组： lik=aik/akk aij= aij- lik* akj k=1,2,…,n-1 i=k+1,k+2, …,n j=k+1,k+2, …,n+1 = 2 \* GB2 ⑵由回代过程求得原方程组的解： xn= ann+1/ an