实验非线性方程求根.docVIP

深 圳 大 学 实 验 报 告 课程名称： 计算方法 实验项目名称： 实验三 非线性方程求根 学院： 计算机与软件学院 专业、班级： 计算机 06 级 3班 指导教师： 张艳 报告人： 卢志敏 学号： 2006131114 实验时间： 09.6.20 实验报告提交时间： 09.6.30 教务处制 要求：提交实验报告（电子版和书面两种形式），内容包含： 实验目的与实验内容 实验目的和实验要求 通过本设计性实验，加深多种非线性方程求根基本原理、相关理论知识的理解，掌握并实现弦截迭代法、Newton迭代法、Newton下山法等算法，培养学生上机操作和实际动手的能力，培养学生分析问题、应用理论知识解决实际问题的能力，使学生逐渐掌握数值算法的设计、分析和应用的综合能力。 要求学生利用极值理论和对分搜索技术，对非线性方程进行分析，设计算法完成根区间隔离，然后实现单根区间上的迭代算法求解方程。 实验题目 求非线性方程x5-3x3+x-1= 0在区间[-10,10〕上的全部实根； 或求cos(x)cosh(x)-1=0在区间[-π,π]上的全部实根. . 试分别用: Newton法;弦截法(割线法); Newton下山法；简化Newton法; 求方程在指定区间上的全部实根.准确到8位有效数字. . 在（1）.的基础上讨论迭代法的加速(选作） 实验内容 对以上非线性方程，分析、设计算法，确定多根区间； 在选取的有根区间内，编程实现用不同方法, 分别求根；（编写通用过程或函数） 对所用算法的收敛阶、局部收敛性、优缺点和精度等作分析及比较.并附运算结果的截图（可用Excel图表或其他绘图软件工具对实验结果分析） 以实验报告的形式提交总结。 学生在完成本实验后，提交实验报告、程序源代码、程序可执行文件以及程序使用说明等文档。 开发环境 Windows xp，vc6.0 模型求解 实验代码 　#includeiostream #includemath.h using namespace std; double a=-10.0,b=10.0,e=0.1,e2=0.000000001; int t1=0,t2=0; int m=0,k=1; //m为解个数,k为对分次数 int* rn=new int[5]; //记录根区间号 double *r=new double[5]; double *r1=new double[5]; // 记录不同方法解得的根 double *r2=new double[5]; double *r3=new double[5]; double *r4=new double[5]; //根区间长度 double h(int k) { int t=1; for(int i=0;ik;i++) t*=2; return((b-a)/t); } //原函数 double f(double x) { return(x*x*x*x*x-3*x*x*x+x-1); } //原函数导数 double fd(double x) { return (5*x*x*x*x-9*x*x+1); } //搜索含根区间 bool search() { while(m5) { if(k4) {break;} k++; double hh=h(k); m=0; for(int i=0;i*hh-10b;i++) { if(f(i*hh-10)*f((i+1)*hh-10)0) {rn[m]=i;m++;} } } return true; } //对分粗略寻根 bool half() {double p,mid,q; for(int i=0;im;i++) { p=(rn[i])*h(k)-10; q=(rn[i]+1)*h(k)-10; mid=(p+q)/2; if(fabs(f(mid))e) r[i]=mid; while(fabs(f(mid))e) { if(f(p)*f