优化设计作业--传动比分配(付廷辉)(精选).docVIP

用matlab建立一个m文件（gnew.m） function f=gnew(f,a,b,epsilon) la=0.618; x1=b-la*(b-a);f1=f(x1); x2=a+la*(b-a);f2=f(x2); for k=1:(log(epsilon/(b-a))/log(0.618)+2) fprintf(1, 迭代次数 k= %3.0f\n,k) if f1=f2 b=x2;x2=x1;f2=f1; x1=b-la*(b-a);f1=f(x1); else a=x1;x1=x2;f1=f2; x2=a+la*(b-a);f2=f(x2); end x=0.5*(b+a); fprintf(1, 迭代区间-左端 a=% 3.4f\n,a) fprintf(1, 试点1坐标值 x1=% 3.4f\n,x1) fprintf(1, 函数值 f1 =% 3.4f\n,f1) fprintf(1, 迭代区间-右端 b=% 3.4f\n,b) fprintf(1, 试点2坐标值 x2=% 3.4f\n,x2) fprintf(1, 函数值 f2=% 3.4f\n,f2) fprintf(1, 区间中点 x=% 3.4f\n,x) disp end 其中f=(0.85-x*((1+sqrt(x*i))/(.97*x*(i+sqrt(30*i))))^(1/3))^2, a=1 , b=2 , epsilon=0.00001, i=（8--40） 运行m文件，并在matlab命令窗口键入： Format long;gnew(inline((0.85-x*((1+sqrt(x*i))/(.97*x*(i+sqrt(30*i))))^(1/3))^2),1,2,0.00001),format short 然后回车 其中i取值为8--40 以i=20为例：(迭代次数25次，x=1.5265) 运行m文件，并在matlab命令窗口键入： Format long;gnew(inline((0.85-x*((1+sqrt(x*20))/(.97*x*(20+sqrt(30*20))))^(1/3))^2),1,2,0.00001),format short 然后回车 Matlab运行结果如下： format long;gnew(inline((0.85-x*((1+sqrt(x*20))/(.97*x*(20+sqrt(20*x))))^(1/3))^2),1,2,0.00001),format short 迭代次数 k= 1 迭代区间-左端 a= 1.3820 试点1坐标值 x1= 1.6180 函数值 f1 = 0.0015 迭代区间-右端 b= 2.0000 试点2坐标值 x2= 1.7639 函数值 f2= 0.0101 区间中点 x= 1.6910 迭代次数 k= 2 迭代区间-左端 a= 1.3820 试点1坐标值 x1= 1.5279 函数值 f1 = 0.0000 迭代区间-右端 b= 1.7639 试点2坐标值 x2= 1.6180 函数值 f2= 0.0015 区间中点 x= 1.5730 迭代次数 k= 3 迭代区间-左端 a= 1.3820 试点1坐标值 x1= 1.4722 函数值 f1 = 0.0006 迭代区间-右端 b= 1.6180 试点2坐标值 x2= 1.5279 函数值 f2= 0.0000 区间中点 x= 1.5000 迭代次数 k= 4 迭代区间-左端 a= 1.4722 试点1坐标值 x1= 1.5279 函数值 f1 = 0.0000 迭代区间-右端 b= 1.6180 试点2坐标值 x2= 1.5623 函数值 f2= 0.0002 区间中点 x= 1.5451 迭代次数 k= 5