多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
global u;
%贝塞尔函数求根
% V5 B( r, U7 ]( }& k- A/ x, [6 `%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
4 |! A7 |/ d8 D. O- ]maxs=20;%要求的贝塞尔函数的非零根个数
# [" {2 |3 w9 h6 F) [4 zincr=4.0;
u(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
5 y) M7 T) \3 c* P# Pfor s=2:maxs
" o; J& o7 I9 S7 c- K5 i3 s$ l    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
5 S* T' Z8 X) ?end
& n' m" L( Z5 C" W- V8 ~+ w' |3 L" ZNIND=40;      %个体数目40
9 z; `( h. P8 b6 C8 ]MAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
1 k$ w# {+ _; D  K; c* V1 ^7 E9 FGGAP=0.9;      %代沟
trace1=[];
trace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
+ Z; A8 b* S# \- R" J%建立区域描述器
FieldD=[rep([PRECI],[1,NVAR]);[20,20,20,20,20,10,10,10,10,10;50,50,50,50,50,100,100,100,100,100];rep([1;0;1;1],[1,NVAR])];
Chrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
v=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
: D5 h- l- H8 ?6 C6 D+ R2 ggen=1;                               %代计数器
9 @2 U. H9 G* x7 d  R3 `%遗传开始
while gen<MAXGEN
    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
/ @! c) F! x) I" s3 v* X( S    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
5 m9 |3 S: ^/ h7 @' L2 N7 G    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
4 `% _8 o+ N: \$ s+ N+ ~    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
. M/ }( B: k: x0 e5 B    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
: K) Y% _4 Q. v; c$ S8 p- d: Z6 ^    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
$ F# s  `5 \3 n3 X    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
& s+ e& x, Y- l; h+ R) L0 E    trace1(gen,1)=max(Highpass(v));
: z# S) l( X6 a9 Q    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
2 h% B5 P  a3 N( c4 y5 i    trace2(gen,1)=min(Lowpass(v));
9 H# O% E, f& l" O    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
    gen=gen+1;
end
+ k9 u5 Y$ u2 l! I( e1 Z  Cfigure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
( D' I- {# U, jlegend('解的变化','种群均值的变化')
( J6 H5 I% y* Q& g" M$ [* X: S/ Dxlabel('迭代次数');ylabel('高通目标函数值');
) c# V0 p) |% @! x4 U% sfigure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
% w0 d+ m* B9 Y. u9 ^. ]plot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
plot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;
: e4 E$ K. d" a2 b( g' V9 {5 D& J9 H
  T$ N# y  p4 O