f(1)=W*v*v*1000/(4*pi*F*u*I2*n0*g);
f(2)=x(5);
f(3)=4*3.4265e+005/(4180*pi*0.113*7.8e-6*x(3)*x(3)*x(5));
options=gaoptimset('ParetoFraction',0.5,'PopulationSize',500,'generations',200,'StallGenLimit',200,'TolFun',1e-100,'PlotFcns',@gaplotpareto);
[x,fval]=gamultiobj(fitnessfcn,nvars,A,b,Aeq,beq,lb,ub,@multiconstr,options);
plot3(fval(:,1),fval(:,2),fval(:,3),'*');
title('三維 Pareto 前沿圖');
x1=xlabel('目標(biāo)函數(shù)f(1)');
x2=ylabel('目標(biāo)函數(shù)f(2)');
x3=zlabel('目標(biāo)函數(shù)f(3)'); |