广东工业大学实验报告广东工业大学实验报告.docVIP

广东工业大学实验报告 成绩评定_______ 学号 3109008538 姓名 刘艳明 教师签名_______ 实验 题目 P150E9 第 9 周 星期 五 第1—2节 编写：用牛顿、平分、成功—失败法、黄金分割法和进退法求函数f(x)的极小点的计算程序，并求解 1.牛顿法： #includeiostream using namespace std; class Newton { double x; double f; double e; public: Newton(double x1,double e1) { x=x1; e=e1; f=x*x*x*x+-16*x*x*x+30*x*x+-24*x+8; } void set(double x1,double e1) { x=x1; e=e1; f=3*x*x*x*x-16*x*x*x+30*x*x-24*x+8; } double getfirstf(double x) { return (12*x*x*x-48*x*x+60*x-24); } double getsecondf(double x) { return (36*x*x-96*x+60); } }; int main() { double x0=3; int i=0; Newton a=Newton(x0,0.01); double f1=a.getfirstf(x0); double f2=a.getsecondf(x0); double x1=x0-f1/f2; while((x1-x0)0.01||(x0-x1)0.01) { i++; x0=x1; a.set(x0,0.01); f1=a.getfirstf(x0); f2=a.getsecondf(x0); coutthe x0 is:x0endl; coutthe x1 is:x1endl; x1=x0-f1/f2; } couteeeeeeeeeeeeeeeeeeeeendl; coutthe last x0 is:x0endl; coutthe last x1 is:x1endl; cout迭代次数为：iendl; return 0; } 2.平分法： #includeiostream using namespace std; double f2(double l,double o,double v,double e); double a,b,eps,sita; int maxstep=30; double x[30]; int main(){ coutPlease enter the initial 区间左端点a: ; cina; coutendl; coutPlease enter the initial 区间左端点b:; cinb; coutendl; coutenter the eps: ; cineps; coutendl; coutenter the min|b-a| : ; cinsita; coutendl; f2(a,b,eps,sita); return 0; } double f2(double l,double o,double v,double e) { for(int i=0;imaxstep;i++){ x[i]=(a+b)/2; cout [a ,b] endl x*=: x[i]endlendl; double df1=5*x[i]*x[i]*x[i]*x[i]+8*x[i]*x[i]*x[i]-12*x[i]*x[i]+2*x[i]+1; if(df10) b=x[i]; else a=x[i]; while(abs(df1)eps||abs(b-a)sita){ coutendl x*= x[i]; return 0; } } } 3.成功—失败法 #includeiostream using namespace std; double funtion1(double b,double u,double g); double funtion2(double x); double chong; double h; double eps; double maxstep=30; int main(){ coutPlease enter the initi