罚函数法MATLAB程序.doc

《罚函数法MATLAB程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、一、进退法、0.618法、Powell法、罚函数法的Matlab程序设计罚函数法（通用）functiony=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+95

2、17.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，以及允许的误差rfunctionG=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);

3、%调用不等式约束函数，将其值%存入数组gh(1)=hh1(x),h(2)=hh2(x)...h(t)=hht(x);%调用等式约束函数，将其值%存入数组hfori=1:pifg(i)

4、a)d=eye(n);%n个线性无关的初始搜索方向k=1;kk=1;xx(1,1:n)=x0;while(kk)y(1,1:n)=xx(k,1:n);forj=1:ns(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);endd(n+1,1:n)=y(n+1,1:n)-y(1,1:n);if(norm(d(n+1,1:n),2)=ww)

5、m=i;endendcha=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

6、el=y(n+1,1:n)functionw=HJ(x0,h,d,dd,q)%0.618算法[a,b]=JTF(x0,h,d,dd,q);%调用进退法算法，确定范围r=0.618;r1=a+(1-r)*(b-a);r2=a+r*(b-a);y1=ff(x0+r1.*dd,q);y2=ff(x0+r2.*dd,q);k=1;while(abs(r1-r2)>=0.1)ify1

7、d,q);endendw=(r1+r2)/2%进退法function[a,b]=JTF(x0,h,d,dd,q)r0=0;y0=ff(x0+r0.*dd,q);k=0;l=1;while(l)r1=r0+h;y1=ff(x0+r1.*dd,q);ify1