多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
global u;
%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
8 b7 Z+ c! F0 l5 l8 _& {$ I1 @incr=4.0;
u(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
for s=2:maxs
% d# J! J( W' F/ `) k8 K' o8 Z    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
end
NIND=40;      %个体数目40
MAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
: Q. s% I) |' r8 D" J. M! xtrace1=[];
% {# f& Z3 o7 _! V6 Strace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
%建立区域描述器
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])];
  b, f$ h1 \/ ]% M2 q  U! RChrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
( m* i4 K& |! U. O( @" ^  |$ K0 R) sv=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
gen=1;                               %代计数器
%遗传开始
  O, |. E+ ^6 J. S' I/ mwhile gen<MAXGEN
. p9 f. D* ?: D    [NIND,N]=size(Chrom);
: q7 Y- B4 i% c( }# p) `& K9 p  l    M=fix(NIND/2);
# V) z6 n, M6 d( ~8 T    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
/ ~4 `  H: A. `5 q" _+ T$ B    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
1 u( q7 R1 U7 O/ \7 {    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
) `7 o9 d, Z2 i$ T2 q    trace1(gen,1)=max(Highpass(v));
* I* M: o& b: n9 H/ z, |    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
; y+ ~5 ?9 Z; G- r/ x# _# L    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
7 y1 ]: r0 }# |    gen=gen+1;
& G1 U* X$ p& i( Y# a" a) Lend
: J  w1 g  \' m" O2 D1 [figure(1);clf;
) q5 c2 s' V* m, n8 ?, ?) q) n- Xplot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
legend('解的变化','种群均值的变化')
% H! V* k5 }5 }5 y/ q. Oxlabel('迭代次数');ylabel('高通目标函数值');
. t  x, D0 K7 z6 V" Kfigure(2);clf;
plot(trace2(:,1));hold on;
$ o% V, q- x, f! rplot(trace2(:,2),'-.');
: B4 H! G1 T  Iplot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
7 b' O  N: S+ o3 s4 @xlabel('迭代次数');ylabel('低通目标函数值');
8 V. e# s  u1 I# Rfigure(3);clf;
: D2 A( s( T% V* ?0 T, @- U7 R3 C. Splot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;
, q  Q* U$ K9 c- B; O7 {) _# O
# @4 e; b: {6 Y- Q