用Powell优化设计程序与一维搜索黄金分割法组合.docVIP

用Powell法优化设计程序与一维搜索黄金分割法组合 编程求解函数 ε=0.001。 解：已知，初始点x0=[1,1]T，迭代精度ε=0.001。 在该优化设计过程中，黄金分割搜索法作为POWELL算法主程序中的一部分。在POWELL算法运行过程中会多次调用黄金分割搜索算法程序。这样可以缩短优化设计计算时间。 1. MATLAB源程序代码 1.1 关于α的目标函数源代码 function m = y(x1,x2,d1,d2,alpha)%建立关于α的一元函数y(α) m = (x1+alpha*d1)^2+2*(x2+alpha*d2)^2-4*(x1+alpha*d1)-2*(x1+alpha*d1)*(x2+alpha*d2); 1.2 一元函数最小值区间函数源代码 function [a,b] = section(x1,x2,d1,d2)%采用外推法求解一元函数的最小值区间 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 y2y1 h = -h;alpha3 = alpha1;y3 = y1;t = 1; end