优化设计作业及程序.doc

《优化设计作业及程序.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、优化设计程序及作业1进退法求搜索区间程序：functionhu=hunt(x1,h)%定义函数为hunt,x1为初始值，h为步长f=@(x)x^2-7*x+10;f1=f(x1);x2=x1+h;f2=f(x2);iff1

2、1,b=x3end在命令窗口运行hunt（0，0.1），得搜索区间[3.4,3.6]2黄金分割法求极值functiong=gold(a,b,j)%a,b为搜索区间，j为精度f=@(x)x^2-7*x+10;k=0.;x1=a+(1-k)*(b-a);x2=a+k*(b-a);whileabs(a-b)>jf1=f(x1);f2=f(x2);iff1

3、display('所求函数')y=f(x)在命令窗口运行gold（3.4，3.6，0.01）得到所求最优解：x=3.4798所求函数值：y=-2.2496黄金分割总程序functionhj=hjfg(q1,h,j)%定义黄金分割法为hjfg，x1为初始值，h为初始步长，j为初始精度f=@(x)x^2-7*x+10;f1=f(q1);q2=q1+h;f2=f(q2);iff1

4、f3jf1=f(c1);f2=f(q2);iff1

5、1，0.01）得到所求x：x=3.4798所求y：y=-2.24963坐标轮换法求二元函数极值functionzblh(q,j,h)%定义坐标轮换法函数为zblh，q为初始值，j为精度，h为初始步长zb=@(x)60-10*x(1)-4*x(2)+x(1)^2+x(2)^2-x(1)*x(2);p3=q;p1=p3+1;s1=[10];s2=[01];whilenorm(p1-p3)>jfori=1:2ifi==1s=s1;p1=p3;p=p1;elses=s2;p=p2;end%黄金分割法子程序f=@(x)zb(p+x.*

6、s);%进退法子程序x1=0;x2=x1+h;f2=f(x2);f1=f(x1);iff10.1f1=f(c1);f2=f(c2);iff1

7、

8、初始值，j为精度，r为惩罚因子，c为收缩系数，h为初始步长cf=@(x)2*x(1)^2+2*x(2)^2-2*x(1)*x(2)-4*x(1)-6*x(2)+r*(-2+x(1)+x(2)+x(1)-5*x(2)-5+1/x(1)+1/x(2))^2;p3=q;p1=p3+1;s1=[1