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

实验课程名称：数值分析 实验项目名称 运用二分法和迭代法 实验成绩实验者张伟专业班级软件0703组别26 同组者实验日期08年月日一、实验内容描述二、实验基本原理与设计 　　1 确定区间[a,b],验证f(a)·f(b)0,给定精确度. 　　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 判断是否达到精确度:即若a-b┃pr,则得到零点近似值a(或b),否则重复2-4.迭代法： 1 确定初始值a，给定精度pr，定义迭代函数x=f（x） 2 fa=a;a=f（a） 3判断是否达到精确度:即若a-fa┃pr,则得到零点近似值a,否则重复2-. 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判断是否达到精确度:即若a-a1┃pr,则得到零点近似值a,否则重复2-. Newton迭代 1 确定初始值a，给定精度pr，定义迭代函数xk+1= xk-(xk*xk*xk-xk-1)/(3*xk*xk-1) 2 fa=a;a=f（a） 3判断是否达到精确度:即若a-fa┃pr,则得到零点近似值a,否则重复2-. 弦截法 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,则得到零点近似值,否则重复2-. 三、主要仪器设备及耗材 Ubuntu8.10 g++ 第二部分：实验调试与结果分析 一、调试过程 #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); } 结果 简单迭代法 #includeiostream #includemath.h using namespace std; double f(double); //申明函数 int main() { double a,i=0,a1,pr; coutplease input the start number:; cina; coutplease input precision:; cinpr; //输入精度 do { a1=a; a=f(a); i++; if(i1000000) {coutIterative is defeated！;return 0;} } while((a-a1)pr||(a-a1)-pr); coutIterative Numbers is iendl; coutThe answer is aendl; //输出迭代次数和结果 } double f(doub