罚函数法MATLAB程序工程优化设计中常常用到的惩罚函数法的Matlab实现的源程序代码.zip

一、进退法、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，以及

%存入数组 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;

if (i==p)&(j==t) %如果所有约束条件的值都在允许的范围内就结束

%搜索，否则加大惩罚力度继续搜索

l=0;

else

x0=x;

q=k*q;

end

end

G=x

%powell 算法，用于寻找无约束最优值点

while (kk)

y(1,1:n)=xx(k,1:n);

for j=1:n

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)

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));

if (cha<cha1)

s(n+1)=HJ(y(n+1,1:n),h,a,d(n+1,1:n),q)

xx(k+1,1:n)=y(n+1,1:n);

k=k+1;

end

end

end

powel=y(n+1,1:n)

function w=HJ(x0,h,d,dd,q) %0.618算法