罚函数法MATLAB程序.doc

一、进退法、0.618法、Powell法、罚函数法的Matlab程序设计罚函数法（通用）function y=ff(x,k)y=-17.86*0.42*x(1)/(0.8+0.42*x(1))*(1-exp(-2*(0.8+0.42*x(1))/3))*exp(-1.6)*x(2)-22.99*x(1)/(0.8+x(1))*(1-exp(-2*(0.8+x(1))/3))*x(3)+k*(x(2)-(1.22*10^2*(9517.8*exp(-1.6-2*0.42*x(1)/3)*x(2)+19035.6*exp(-2*x(1)/3)*x(3)))/(1.22*10^2+9517.8*exp(-1.6-2*0.42*x(1)/3)*x(2)+19035.6*exp(-2*x(1)/3)*x(3)))^2+k*(x(3)-exp(-0.8-2*x(1)/3)*x(3)-exp(-2.4-2*0.42*x(1)/3)*x(2))^2;% 主函数，参数包括未知数的个数n，惩罚因子q，惩罚因子增长系数k，初值x0，以及允许的误差rfunction G=FHS(x0,q,k,n,r,h,a) l=1;while (l) x=powell(x0,n,q,r(1),h,a); %调用powell函数 g(1)=ff1(x),g(2)=ff2(x) . . . g(p)=ffp(x); %调用不等式约束函数，将其值 %存入数组g h(1)=hh1(x),h(2)=hh2(x) . . . h(t)=hht(x); %调用等式约束函数，将其值%存入数组h for i=1:p if g(i)r for j=1:t if abs(h(j))r(2) continue; else break; end end else break; end end if (i==p)(j==t) %如果所有约束条件的值都在允许的范围内就结束%搜索，否则加大惩罚力度继续搜索 l=0; else x0=x; q=k*q; endend G=x%powell算法，用于寻找无约束最优值点 function powel=powell(x0,n,q,r,h,a)  d=eye(n); %n个线性无关的初始搜索方向 k=1; kk=1; xx(1,1:n)=x0; while (kk) y(1,1:n)=xx(k,1:n); for j=1:n s(j)=HJ(y(j,1:n),d(j,1:n),q); %调用0.618算法 y(j+1,1:n)=y(j,1:n)+s(j).*d(j,1:n); end d(n+1,1:n)=y(n+1,1:n)-y(1,1:n); if (norm(d(n+1,1:n),2)r) kk=0; break; else ww=0; m=1; for i=1:n gg=ff(y(i,1:n),q)-ff(y(i+1,1:n),q); if (gg=ww) m=i; end end cha=ff(y(1,1:n),q)-2*ff(y(n+1,1:n),q)+ff(2*y(n+1,1:n)-y(1,1:n),q); cha1=2*(ff(y(m,1:n),q)-ff(y(m+1,1:n),q));