数值实验报告侯瑞祥.docVIP

本科实验报告 课程名称： 计算机数值方法 实验项目： 方程求根 实验地点： 致远楼B301 专业班级： 软件1216 学号：2012005403 学生姓名： 侯瑞祥 指导教师： 田华 2014年 5 月 21 日 学生姓名 侯瑞祥 实验成绩 实验名称 实验一 方程求根 实验目的和要求 熟悉使用、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：二分法f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度|x*-xn|0.5×10-5 实验内容和原理 二分法：函数f(x)在区间（x，y）上连续先（x，y），(a)，(b)异号，说明在区间(ab)内零点，然后求[(a+b)/2]。 　　假设(a)0,F(b)0,ab， 如果[(a+b)/2]=0，该点零点 如果[(a+b)/2]0，则区间（(a+b)/2，b)内零点，(a+b)/2a；　　如果[(a+b)/2]0，则区间(a,(a+b)/2)内零点，(a+b)/2b； 返回①重新循环，不断接近零点每次把(x)的零点所在区间收缩一半的方法，使区间两个端点逐步函数零点，求得零点近似值’(x0)，再把f(x)在x1处展开为泰勒级数，取其线性部分f(x)=0的近似方程，如此进行下去，得到牛顿法的迭代公式：Xn+1=Xn-f(Xn)/f’(Xn) 主要仪器设备 Vc++、WIN-TC、HP6470b 实验记录 牛顿法： #includestdio.h void main() { float x[10]={1.5}; float y,w; int n; y=x[n]*x[n]*x[n]+4*x[n]*x[n]-10; w=3*x[n]*x[n]+8*x[n];//对原函数进行求导 x[1]=1.5-(1.5*1.5*1.5+4*1.5*1.5-10)/(3*1.5*1.5+8*1.5); while(n=0,n=0,n++){ while(x[n+1]-x[n]0.000005) { x[n+1]=x[n]-y/w;//牛顿法的迭代式 } } printf(%f\n,x[n]); getch(); } 二分法： #includestdio.h void main() { float a=1.0,b=2.0; int n; float mid; for(n=0;n18;n++)//根据1/2^n=0.5*10^(-3)可得迭代18次后停止迭代 { mid=(a+b)/2; if ((a*a*a+4*a*a-10)*(mid*mid*mid+4*mid*mid-10)0)//二分法的中心思想 { b=mid; } else { a=mid; } } printf(%f\n,mid); } 实验结果和分析 牛顿法： 二分法： 分析： 使用二分法和牛顿法均能计算出方程的根，但利用牛顿法要比二分法计算的次数少，并且能够较早的达到精度要求。 心得体会 本次数值计算方法程序设计实验是在不断的习题练习中跳脱出来，直接面对实用性较强的程序代码编写。效果很好，不仅加深对二分法、牛顿法的理解，还加强了实际用运能力。将理论成功地转化成实践结果，虽然在实验过程中有很多地方并不是一步到位，但是经过不断的改进，最终达到了目的，这也是对自己能力的一次锻炼。 本科实验报告 课程名称： 计算机数值方法 实验项目： 线性方程组的直接求法 实验地点： 致远楼B301 专业班级： 软件1216 学号：2012005403 学生姓名： 侯瑞祥 指导教师： 田华 2014年 5 月 21 日 实验名称