用Powell法優(yōu)化設(shè)計程序與一維搜索黃金分割法組合

1、用Powell法優(yōu)化設(shè)計程序與一維搜索黃金分割法組合編程求解函數(shù)的極小點x，初始點x0=1,1T，迭代精度=0.001。解：已知，初始點x0=1,1T，迭代精度=0.001。在該優(yōu)化設(shè)計過程中，黃金分割搜索法作為POWELL算法主程序中的一部分。在POWELL算法運行過程中會多次調(diào)用黃金分割搜索算法程序。這樣可以縮短優(yōu)化設(shè)計計算時間。1. MATLAB源程序代碼1.1 關(guān)于的目標函數(shù)源代碼function m = y(x1,x2,d1,d2,alpha)%建立關(guān)于的一元函數(shù)y()m = (x1+alpha*d1)2+2*(x2+alpha*d2)2-4*(x1+alpha*d1)-2*(x1+

2、alpha*d1)*(x2+alpha*d2);1.2 一元函數(shù)最小值區(qū)間函數(shù)源代碼function a,b = section(x1,x2,d1,d2)%采用外推法求解一元函數(shù)的最小值區(qū)間x11 = x1;x22 = x2;d11 = d1;d22 = d2;h0 = 1;h = h0;alpha1 = 0;y1 = y(x11,x22,d11,d22,alpha1);alpha2 = h;y2 = y (x11,x22,d11,d22,alpha2);t = 0;if y2>y1 h = -h;alpha3 = alpha1;y3 = y1;t = 1;endwhile(1) if

3、t = 1 alpha1 = alpha2;y1 = y2; alpha2 = alpha3;y2 = y3; else t = 1; end alpha3 = alpha2+h;y3 = y(x11,x22,d11,d22,alpha3); if y3<y2 h = 2*h; else break; endendif alpha1>alpha3 tem = alpha1;alpha1 = alpha3;alpha3 = tem; a = alpha1;b = alpha3;else a=alpha1;b = alpha3;end1.3黃金分割搜索法函數(shù)源代碼function al

4、pha = ALPHA(x1,x2,d1,d2,A,B)%利用黃金分割法求解關(guān)于的函數(shù)y()的極小點*x11 = x1;x22 = x2;d11 = d1;d22 = d2;a = A;b = B;ep = 0.001;r = 0.618;alpha1 = b-r*(b-a);y1 = y(x11,x22,d11,d22,alpha1);alpha2 = a+r*(b-a);y2 = y(x11,x22,d11,d22,alpha2);while(1) if y1>=y2 a = alpha1;alpha1 = alpha2; y1 = y2; alpha2 = a+r*(b-a); y

5、2 = y(x11,x22,d11,d22,alpha2); else b = alpha2;alpha2 = alpha1; y2 = y1; alpha1 = b-r*(b-a); y1 = y(x11,x22,d11,d22,alpha1); end if abs(b-a)<ep && abs(y2-y1)<ep break; endendalpha = 0.5*(a+b);1.4 POWELL算法程序源代碼function z,fmin = powell(f) %根據(jù)無約束優(yōu)化方法POWELL算法，在二維空間中求解f(x1,x2)的最小值點，求解結(jié)果返回變量

6、坐標(x1,x2)和極小值fmink = 0;n = 2;x = 0;0;ff(1) = f(x(1),x(2);ep = 0.001;d = 1;0;0;1;while(1) x00 = x(1);x(2); for i = 1:n a(i),b(i)=section(x(2*i-1),x(2*i),d(2*i-1),d(2*i); alpha(i)=ALPHA(x(2*i-1),x(2*i),d(2*i-1),d(2*i),a(i),b(i); x(2*i+1)=x(2*i-1)+alpha(i)*d(2*i-1); x(2*i+2)=x(2*i)+alpha(i)*d(2*i); ff(

7、i+1)=f(x(2*i+1),x(2*i+2); end for i=1:n Delta(i)=ff(i)-ff(i+1); end delta=max(Delta); for i = 1:n if delta = Delta(i) m = i; break; end end d(2*n+1) = x(2*n+1)-x(1); d(2*n+2)=x(2*n+2)-x(2); x(2*n+3)=2*x(2*n+1)-x(1); x(2*n+4)=2*x(2*n+2)-x(2); ff(n+2)=f(x(2*n+3),x(2*n+4); f0=ff(1);f2=ff(n+1);f3=ff(n+2

8、); k=k+1; R(k,:)=k,x',d',ff; if f3<f0 && (f0-2*f2+f3)*(f0-f2-delta)2<0.5*delta*(f0-f3)2 a(n+1),b(n+1)=section(x(2*n+1),x(2*n+2),d(2*n+1),d(2*n+2); alpha(n+1)=ALPHA(x(2*n+1),x(2*n+2),d(2*n+1),d(2*n+2),a(n+1),b(n+1); x(1)=x(2*n+1)+alpha(n+1)*d(2*n+1); x(2)=x(2*n+2)+alpha(n+1)*d(2

9、*n+2); for i=m:n d(2*i-1)=d(2*i+1); d(2*i)=d(2*i+2); end else if f2<f3 x(1)=x(2*n+1); x(2)=x(2*n+2); else x(1)=x(2*n+3); x(2)=x(2*n+4); end end RR(k,:)=alpha; ff(1)=f(x(1),x(2); if (x(2*n+1)-x00(1)2+(x(2*n+2)-x00(2)2)(1/2)<ep break; endendz = x(1);x(2);fmin = f(x(1),x(2);2. 運行程序計算>> f =

10、(x1,x2) x12+2*x22-4*x1-2*x1*x2f = (x1,x2)x12+2*x22-4*x1-2*x1*x2>> z,fmin = powell(f)z = 3.9998 1.9998fmin = -8.0000>>經(jīng)過計算可知，極小值點為(3.9998,1.9998)，極小值為-8。3. 驗算通過MATLAB軟件內(nèi)置fminsearch函數(shù)進行驗算。3.1 繪圖>> x = linspace(-10,10,50);y = linspace(-10,10,50);x,y = meshgrid(x,y);z = x.2+2*y.2-4*x-2

11、*x.*y;subplot(1,2,1);%在一幅框中顯示兩張圖，顯示第一張圖cs = contour(x,y,z);clabel(cs); %繪制等高線圖xlabel('x1');ylabel('x2');title('(a)Contour plot');grid;>> subplot(1,2,2); %在一幅框中顯示兩張圖，顯示第二張圖cs = surfc(x,y,z); %繪制空間曲面圖zmin = floor(min(z);zmax = ceil(max(z);xlabel('x1');ylabel('x2');zlabel('f(x1,x2)');title('(b) Mesh plot');>>通過觀察可以發(fā)現(xiàn)極小值點在(4,2)附近。3.2 fminsearch函數(shù)驗算>> f = (x) x(1)2+2*x(2)2-4*x(1)-2*x(1)*x(2)f = (x)x(1)2+2*x(2)2-4*x(1)-2*x(1)*x(2)>> x,fval = fm