数值实验报告..docVIP

（实验报告的首页） 本科实验报告 课程名称： 计算机数值方法 实验项目： 方程求根、线性方程组的直接求解和迭代求解、 代数差值和最小二乘法拟合 实验地点： 专业班级： 学号： 学生姓名： 指导教师： 年 月 日  学生姓名 实验成绩 实验名称 实验一 方程求根 实验目的和要求（必填） 熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度|x*-xn|0.5×10-5 实验内容和原理（必填） 函数f(x)在区间（x，y）上连续，先在区间（x，y）确定a与b，若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)的零点所在区间收缩一半的方法，使区间内的两个端点逐步逼近函数零点，最终求得零点近似值。 主要仪器设备 台式或笔记本计算机 实验记录(写出实验内容中的程序代码和运行结果)(可分栏或加页) 割线法 程序代码 #include stdio.h #include math.h int main() { float c,a=1.0,b=2.0; 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.00001)break; b=c; printf(%f\n,b); } printf(%f\n,c); getch(); } 2.二分法 程序代码 #include stdio.h #include stdlib.h #include math.h int main() { double a=1.0, b=2.0; double c,d; while(1) { c=(a+b)/2; d=pow(c,3)+4*c*c-10; if (-0.000005 d d 0.000005)break; else if(d 0){ a=c;} else if(d 0) { b=c;} printf(%f\t%f\n,a,b);} printf(%f\n,c); printf(%f\n,d); getch(); return 0;} 截图 实验结果和分析 用二分法和割线法均能计算出方程的根，但割线计算次数少于二分法，且能够较早达到精度要求，割线法的代码也较为简单明了，实用方便。 心得体会（遇到的问题和解决方法） 通过本次试验，编写代码，运行程序，使我复习了c语言的同时也加深了对二分法和割线法的理解，并认识了一些新的函数。  实验名称 实验二 线性方程组的直接求解 实验目的和要求 合理选择利用Gauss消元法、LU分解法、追赶法求解下列方程组： ① ② ③ ④（n=5,10,100,…） 实验内容 LU分解法：将系数矩阵A转化为A=L*U，L为单位下三角矩阵，U为普通上三角矩阵，然后通过解方程组l*y=b，u*x=y,来求解x。 高斯消元法：将原方程组化为三角形方阵的方程组：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 ） xn= ann+1/ ann xk=( akn+1-∑akj xj)/ akk 完全主元素消元法和列主元素消元法 主要仪器设备 台式或笔记本计算机 实验记录(写出实验内容中的程序代码和运行结果)(可分栏或加页) 1.列主元素消元法 程序代码 #include stdio.h int main() { float a[3][4]={{1,2,3,14},{0,1,2,8},{2,4,1,13}};