moead.zip_DTLZ_DTLZ函数_MOEAD_alikesfl_daily4b7

function [pa]=OVcom(v,NO) x=v; Pa=[]; % OVcom: define the test problems,compute the objective values to v % usage: % pa=OVcom(v,NO); % where % NO: the serial number of test problem % v: the solutions in decision space % pa: the solutions in objective space % Authors: Maoguo Gong and Licheng Jiao % April 7, 2006 % Copyright (C) 2005-2006 by Maoguo Gong (e-mail: gongmg@126.com) %-------------------------------------------------------------------------- switch NO %%1-3 are MISA's problems case 1000 r=[0.006 0.0062 0.0104 0.0231 0.0067 0.0089 0.0164 0.0261 0.0078 0.0156 0.0223 0.012]; k=[0.001 0.002 0.003 0.005 0.003 0.002 0.004 0.005 0.002 0.003 0.004 0.002]; x0=[0.02 0.03 0.04 0.15 0.04 0.04 0.1 0.2 0.05 0.1 0.15 0.08]; R=[0.0064 0.00655 0.0109 0.02345 0.00725 0.00935 0.01685 0.02645 0.00825 0.01615 0.0226 0.0124]; Q=[0.026 0.005 0.007;0.028 0.006 0.008;0.034 0.003 0.005;0.042 0.007 0.008;0.029 0.005 0.007;0.032 0.004 0.008;0.039 0.006 0.005;0.044 0.007 0.009;0.03 0.004 0.006;0.038 0.006 0.009;0.04 0.003 0.007;0.036 0.005 0.006]; h=[0.022 0.034;0.032 0.044;0.014 0.025;0.016 0.034;0.012 0.026;0.032 0.096;0.033 0.054; 0.025 0.037;0.031 0.057;0.023 0.046;0.016 0.045;0.023 0.036]; H=[0.04 0.06]; Uj=0.4; x=v; [n,m]=size(x); f1=sum(ones(n,1)*r.*x,2)-sum(ones(n,1)*k.*abs(x-ones(n,1)*x0),2)-abs(sum(x,2)-1)*0; f2=sum(ones(n,1)*R.*x,2)-sum(ones(n,1)*k.*abs(x-ones(n,1)*x0),2)-abs(sum(x,2)-1)*0; f3=abs(sum(x,2)-1)*0; for i=1:12 for j=1:12 f3=f3+x(:,i).*x(:,j).*(Q(i,3)+Q(i,2)).*(Q(j,3)+Q(j,2)); end end pa=[1-f1,1-f2,f3*1000]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% case 81 vg=v(:,2:10); t=(vg-sin(pi*v(:,1)/2)*ones(1,size(vg,2)))/1; % gx=2*sin(pi*v(:,1)).*sum(-0.9*t.^2+(abs(t)).^0.6,2); % gx=20*sum(-0.9*t.^2+(abs(t)).^0.6,2); gx=2*sin(pi*v(:,1)).*(9+sum(t.^2-cos(2*pi*t),2)); f1=v(:,1).*(1+gx); f2=(1+gx).*(1-sqrt(v(:,1))); pa=[f1,f2]; case 82 vg=v(:,2:10); t=vg-sin(pi*v(:,1)/2)*ones(1,size(vg,2)); % gx=2*sin(pi*v(:,1)).*sum(abs(t)./(1+exp(5*abs(t))),2); % gx=2*sum(-0.9*t.^2+(abs(t)).^0.6,2); % gx=10*sum(abs(t)./(1+exp(5*abs(t))),2); gx=2*sin(pi*v(:,1)).*(9+sum(t.^2-cos(2*pi*t),2)); f1=v(:,1).*(1+gx); f2=(1+gx).*(1-v(:,1).^2); pa=[f1,f2]; case 83 vg=v(:,2:8); t=vg-sin(pi*v(:,1)/2)*ones(1,size(vg,2)); % gx=2*sin(pi*v(:,1)).*sum(abs(t)./(1+exp(5*abs(t))),2); % gx=20*sum(-0.9*t.^2+(abs(t)).^0.6,2); gx=2*sin(pi*v(:,1)).*(7+sum(t.^2-cos(2*pi*t),2)); f1=cos(pi*v(:,1)/2).*(1+gx); f2=(1+gx).*sin(pi*v(:,1)/2); pa=[f1,f2]; case 84 vg=v(:,2:8); t=vg-sin(pi*v(:,1)/2)*ones(1,size(vg,2)); % gx=2*sin(pi*v(:,1)).*sum(abs(t)./(1+exp(5*abs(t))),2); % gx=20*sum(-0.9*t.^2+(abs(t)).^0.6,2); gx=2*sin(pi*v(:,1)).*(7+sum(t.^2-cos(2*pi*t),2)); f1=v(:,1).*(1+gx); f2=(1+gx).*(1-v(:,1).^0.5.*(cos(pi*v(:,1)*2)).^2); pa=[f1,f2]; case 85 vg=v(:,2:10); t=vg-sin(pi*v(:,1))*ones(1,size(vg,2)); gx=2*abs(cos(pi*v(:,1))).*sum(-0.9*t.^2+(abs(t)).^0.6,2); % gx=20*sum(-0.9*t.^2+(abs(t)).^0.6,2); gx=2*sin(pi*v(:,1)).*(sum(t.^2,2)); f1=v(:,1).*(1+gx); f2=(1+gx).*(1-sqrt(v(:,1))); pa=[f1,f2]; case 86 vg=v(:,3:10); % t=vg-0.5; t=(vg-v(:,1).*v(:,2)*ones(1,size(vg,2)))/1; % gx=2*(sin(pi*v(:,1))).*sum(-0.9*t.^2+(abs(t)).^0.6,2); % gx=20*sum(-0.9*t.^2+(abs(t)).^0.6,2); gx=2*sin(pi*v(:,1)).*(8+sum(t.^2-cos(2*pi*t),2)); f1=v(:,1).*v(:,2).*(1+gx); f2=(1+gx).*v(:,1).*(1-v(:,2)); f3=(1+gx).*(1-v(:,1)); pa=[f1,f2,f3]; case 87 vg=v(:,3:8); t=(vg-v(:,1).*v(:,2)*ones(1,size(vg,2)))/1; % gx=2*sin(pi*v(:,1)).*sum(-0.9*t.^2+(abs(t)).^0.6,2); % gx=20*sum(-0.9*t.^2+(abs(t)).^0.6,2); gx=2*sin(pi*v(:,1)).*(6+sum(t.^2-cos(2*pi*t),2)); f1=cos(v(:,1)*pi/2).*cos(v(:,2)*pi/2).*(1+gx); f2=(1+gx).*cos(v(:,1)*pi/2).*sin(v(:,2)*pi/2); f3=(1+gx).*sin(v(:,1)*pi/2); pa=[f1,f2,f3]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% case 711 s=0;S=0; x=v;y=v(:,2); for i=1:14 s=s+(x(:,2*i+1)-x(:,1).^(0.5+1.5*(2*i-1)/28)).^2; S=S+(x(:,2*i)-x(:,1).^(0.5+1.5*(2*i-2)/28)).^2; end f1=x(:,1)+2*s/14; f2=(1+10*y).*(2-(x(:,1)./(1+10*y)).^2-(x(:,1)./(1+10*y)).*sin(2*pi*4*x(:,1)))+2*(S+(x(:,30)-x(:,1).^2).^2)/15; pa=[f1 f2]; case 7127 %ZDT3 s=0;S=0; x=v;y=v(:,2); for i=1:5 s=s+(x(:,2*i+1)-x(:,1).^(0.5+1.5*(2*i-1)/10)).^2; S=S+(x(:,2*i)-x(:,1).^(0.5+1.5*(2*i-2)/10)).^2; end f1=v(:,1)+2*s/5; S=S+(x(:,12)-x(:,1).^2).^2; f2=2*S/6+(3-sqrt(v(:,1)./(1+S))-(v(:,1)./(1+S)).*sin(10*pi*v(:,1))); pa=[f1,f2]; case 712 s=0;S=0; x=v;y=v(:,1); for i=1:14 s=s+(x(:,2*i+1)-x(:,1).^(0.5+1.5*(2*i-1)/28)).^2; S=S+(x(:,2*i)-x(:,1).^(0.5+1.5*(2*i-2)/28)).^2; end f1=zeros(size(v,1),1); f2=(y-1).^2; sx1=find(y<=1);x1=y(sx1);f1(sx1)=1-x1; sx2=find(y>1&y<=3);x2=y(sx2);f1(sx2)=1-2+x2; sx3=find(y>3&y<=4);x3=y(sx3);f1(sx3)=5-x3; sx4=find(y>4);x4=y(sx4);f1(sx4)=x4-3; f1=f1+2*s/14; f2=f2+2*(S+(x(:,30)-x(:,1).^2).^2)/15; pa=[f1 f2]; case 2727 %ZDT3 vg=v(:,2:6); gx=1+100*(5+sum((vg-0.5).^2-cos(20*pi*(vg-0.5)),2)); f1=v(:,1); f2=(3-sqrt(v(:,1)./gx)-(v(:,1)./gx).*sin(10*pi*v(:,1))); pa=[f1,f2]; case 101 %% MISA's Example 1 x=v(:,1);y=v(:,2); vg=v(:,3:12); gx=100*(10+sum((vg-0.5).^2-cos(20*pi*(vg-0.5)),2)); f1=(1+gx).*(1+(0.5*sin(1)-2*cos(1)+sin(2)-1.5*cos(2)-(0.5*sin(x)-2*cos(x)+sin(y)-1.5*cos(y))).^2+(1.5*sin(1)-cos(1)+2*sin(2)-0.5*cos(2)-(1.5*sin(x)-cos(x)+2*sin(y)-0.5*cos(y))).^2); f2=(1+gx).*((x+3).^2+(y+1).^2); pa=[f1,f2]; case 1 %% MISA's Example 1 x=v(:,1);y=v(:,2); vg=v(:,3:7); gx=100*(5+sum((vg-0.5).^2-cos(20*pi*(vg-0.5)),2)); f1=(1+gx).*x; f2=(1+gx).*(1+10*y).*(2-(x./(1+10*y)).^2-(x./(1+10*y)).*sin(2*pi*4*x)); pa=[f1,f2]; case 2%% MISA's Example 2 x=v(:,1); vg=v(:,2:6); gx=100*(5+sum((vg-0.5).^2-cos(20*pi*(vg-0.5)),2)); f1=zeros(size(v,1),1); f2=(1+gx).*(x-5).^2; sx1=find(x<=1);x1=x(sx1);f1(sx1)=1-x1; sx2=find(x>1&x<=3);x2=x(sx2);f1(sx2)=1-2+x2; sx3=find(x>3&x<=4);x3=x(sx3);f1(sx3)=5-x3; sx4=find(x>4);x4=x(sx4);f1(sx4)=x4-3; pa=[f1,(1+gx).*f2]; case 212%%DTLZ3 vg=v(:,3:7);x=v(:,1);y=v(:,2); f1=zeros(size(v,1),1); f2=f1; f3=f1; gx=100*(5+sum((vg-0.5).^2-cos(20*pi*(vg-0.5)),2)); sx1=find(x<=1); x1=x(sx1);y1=y(sx1);f1(sx1)=cos(v(sx1,1)*0.5*pi).*cos(v(sx1,2)*0.5*pi); f2(sx1)=cos(v(sx1,1)*0.5*pi).*sin(v(sx1,2)*0.5*pi); f3(sx1)=4+sin(v(sx1,1)*0.5*pi); sx2=find(x>1&x<=2);f1(sx2)=2+cos((v(sx2,1)-1)*0.5*pi).*cos(v(sx2,2)*0.5*pi); f2(sx2)=2+cos((v(sx2,1)-1)*pi*0.5).*sin(v(sx2,2)*0.5*pi); f3(sx2)=2+sin((v(sx2,1)-1)*pi*0.5); sx3=find(x>2);f1(sx3)=4+cos((v(sx3,1)-2)*0.25*pi).*cos(v(sx3,2)*0.5*pi); f2(sx3)=4+cos((v(sx3,1)-2)*pi*0.25).*sin(v(sx3,2)*0.5*pi); f3(sx3)=sin((v(sx3,1)-2)*pi*0.25); % pa=[f1,f2,f3]; % f1=(1+gx).*cos(v(:,1)*0.5*pi).*cos(v(:,2)*0.5*pi); % f2=(1+gx).*cos(v(:,1)*0.5*pi).*sin(v(:,2)*0.5*pi); % f3=(1+gx).*sin(v(:,1)*0.5*pi); pa=[(1+gx).*f1,(1+gx).*f2,(1+gx).*f3]; case 213 %%DTLZ1 y=v(:,1); f1=zeros(size(v,1),1); s=0;S=0; t=0;x=v; for i=1:9