1、function y=ff(x1,x2)%y=10*(x1+x2-5)2+(x1-x2)2; y=x12+x22-x1*x2-10*x1-4*x2+60;function alpha=goldmethod2(x01,x02,d,h0)% 輸入分割比gamagama=0.618; % 確定搜索區(qū)間a,ba,b=search2(x01,x02,d,h0); %計算a1,a2和函數(shù)值y1,y2, 確定最優(yōu)步長 alphaa1=b-gama*(b-a);y1=ff(x01+d(1)*a1,x02+d(2)*a1);a2=a+gama*(b-a);y2=ff(x01+d(1)*a2,x02+d(2)*a

2、2);for n=1:100 if y1=y2 x1(n)=a1;x2(n)=a2;yp1(n)=y1;yp2(n)=y2; a=a1;a1=a2;y1=y2; a2=a+gama*(b-a);y2=ff(x01+d(1)*a2,x02+d(2)*a2); else x1(n)=a1;x2(n)=a2;yp1(n)=y1;yp2(n)=y2; b=a2;a2=a1;y2=y1; a1=b-gama*(b-a);y1=ff(x01+d(1)*a1,x02+d(2)*a1); end aa(n)=(a+b)/2; %輸出極值點 y(n)=ff(x01+d(1)*aa(n),x02+d(2)*aa(

3、n); %輸出極值點函數(shù)值 e(n)=abs(b-a); alpha=(a+b)/2; if abs(b-a)1e-5 break; endend clear;h0=0.1;x0=0,0;e=1,0;0,1;xp=0,0;for k=1:20 k for i=1:2 i x01=x0(1);x02=x0(2);d=e(i,:); alpha=goldmethod2(x01,x02,d,h0); % 確定最優(yōu)步長 alpha alpha x1=x0+alpha*d; x0=x1; xp end xn=x1; if abs(xn-xp)1e-1 break; end xp=xn;endclear;

4、% 給定初始值h0=0.1;x0=0,0;e=1,0;0,1;n=2;for k=1:20 k % k 代表第輪計算 % 第 1 到 n 次搜索計算 for i=1:1 % 第1次搜索 i x01=x0(1);x02=x0(2);d=e(i,:); F0=ff(x0(1),x0(2); alpha=goldmethod2(x01,x02,d,h0); % 確定最優(yōu)步長 alpha alpha x(i,:)=x0+alpha*d; F(i)=ff(x(i,1),x(i,2); Delta(i)=F(i)-F0; end for i=2:n % 第2到n次搜索 i x01=x(i-1,1);x02

5、=x(i-1,2);d=e(i,:); alpha=goldmethod2(x01,x02,d,h0); % 確定最優(yōu)步長 alpha alpha x(i,:)=x(i-1,:)+alpha*d; F(i)=ff(x(i,1),x(i,2); Delta(i)=F(i-1)-F(i); end % 111 % 精度是否滿足要求 2 if abs(x(n,:)-x0)1e-3 break; end % 222 % 第n次結(jié)束參數(shù)計算 d=x(n,:)-x0; x(n+1,:)=2*x(n,:)-x0; Delta_m,j=max(Delta); F2=F(n);F3=ff(x(n+1,1),x(n+1,2); % 333 % 搜索方向是否進行替換 if F3=Delta_m/2*(0-F3)2 x01=x(n,1);x02=x(n,2); alpha=goldmethod2(x01,x02,d,h0); % 確定最優(yōu)步長 alpha alpha x(n+1,:)=x(n,:)+alpha*d; e(j,:)=; e(n,:)=d; x0=x(n+1,:); else if F2y1 %反向搜索 h=-h; a3=a1;y3=y1; a1=a2;y1=y2; a2=a3;y2=y3;en