2019年实验课程名称数值分析.docVIP

实验课程名称：数值分析 实验项目名称 运用二分法和迭代法 实验成绩 实验者 张伟 专业班级 软件0703 组别 26 同组者 实验日期 08年 月 日 第一部分：实验分析与设计 一、实验内容描述 用二分法和简单迭代法计算X3-X-1=0在[1，2]内的近似根，要求精度达到10-4和10-5，并给出二分和迭代的次数，分析和比较算法的效率差异，最后与理论估计的二分和迭代次数作比较。 用Aitken算法通过编程计算上述算是的近似根，要求精度达到10-4，并分析算法精度与二分法比较。 用Newton迭代法和弦截法通过编程计算上述算式，精度达到10-5，并分析算法精度及不同方法的计算量。分析和比较算法的效率差异：收敛速度和总计算量等之间的关系。 二、实验基本原理与设计 eq \o\ac(○,1) 二分法： 　　1 确定区间[a,b],验证f(a)·f(b)0,给定精确度pr. 　　2 求区间(a,b)的中点c. 　　3 计算f(c). 　　(1) 若f(c)=0,则c就是函数的零点; 　　(2) 若f(a)·f(c)0,则令b=c; 　　(3) 若f(c)·f(b)0,则令a=c. 4 判断是否达到精确度pr:即若┃a-b┃pr,则得到零点近似值a(或b),否则重复2-4. eq \o\ac(○,2)迭代法： 1 确定初始值a，给定精度pr，定义迭代函数x=f（x） 2 fa=a;a=f（a） 3判断是否达到精确度pr:即若┃a-fa┃pr,则得到零点近似值a,否则重复2-3. eq \o\ac(○,3)Aitken迭代法 1 确定初始值a，给定精度pr，定义迭代函数x=f（x） 2 a1=a;fa=f（a）; faa=f（fa）；a=faa-(faa-fa)*(faa-fa)/(faa-fa*2+a); 3判断是否达到精确度pr:即若┃a-a1┃pr,则得到零点近似值a,否则重复2-3. eq \o\ac(○,4)Newton迭代 1 确定初始值a，给定精度pr，定义迭代函数xk+1= xk-(xk*xk*xk-xk-1)/(3*xk*xk-1) 2 fa=a;a=f（a） 3判断是否达到精确度pr:即若┃a-fa┃pr,则得到零点近似值a,否则重复2-3. eq \o\ac(○,5)弦截法 1 确定初始值a，b给定精度pr，定义函数f()和弦截法函数g（a,b）{return b-f(b)*(b-a)/(f(b)-f(a));} 2 b1=b;b=g(a,b); a=b1; 3判断是否达到精确度pr:即若┃b-a┃pr,则得到零点近似值b,否则重复2-3. 三、主要仪器设备及耗材 Ubuntu8.10 g++ 第二部分：实验调试与结果分析 一、调试过程 eq \o\ac(○,1)二分法 #includeiostream using namespace std; double f(double); //申明函数 int main() { double a,b,c,i=0,fa,fb,pr; coutplease input the start number:; cina; coutplease input the last number:; cinb; //输入边界值 coutplease input precision:; cinpr; //输入精度 if(f(a)*f(b)0) {coutBetween a and b does not have a answer!endl;}//判断是否有根 else { do { c=(a+b)/2; fa=f(a); fb=f(c); if(fa*fb0) b=c; else a=c; i++; } while(b-apr); coutIterative Numbers is iendl; coutThe answer is cendl; //输出迭代次数和结果 } return 0; } double f(double x) //定义函数 { return (x*x*x-x-1); } 结果 eq \o\ac(○,2)简单迭代法 #includeiostream #includemath.h using namespace std; double f(double);