多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
( z! [4 m; v& F. L0 e- e" eglobal u;
/ [' f* N! S- \- x/ E%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
: t& K4 @- S/ f& ^6 I3 bmaxs=20;%要求的贝塞尔函数的非零根个数
1 U' s0 v5 o% K/ f3 Uincr=4.0;
! \1 V* ~: P* h3 r) au(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
( f, c0 d: b5 w  |; G& O7 |1 Efor s=2:maxs
, F3 }- \+ p3 E    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
end
NIND=40;      %个体数目40
MAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
. y4 K5 g* G; r' A4 A1 n# @1 bPRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
trace1=[];
: _2 l- e' H- d6 r8 Z, f; G: Ftrace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
; h1 V. p+ q# ?  F& @" k%建立区域描述器
- N- a! P5 z7 _4 ~$ m) Y8 b  ^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);        %初始化种群
1 ?7 n5 r* I3 yv=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
  H- L2 _$ o* \gen=1;                               %代计数器
& |( f1 B) O# ?5 k% G%遗传开始
* w7 E9 M$ o! e4 \/ Zwhile gen<MAXGEN
    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
/ Y' |" s) e! P& v( {    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
2 s  l+ A- E6 u1 f    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
! F! e+ A' j& N/ D    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
6 G6 G* U2 g4 i9 v6 U* ]' {    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
    trace1(gen,1)=max(Highpass(v));
    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
    trace2(gen,1)=min(Lowpass(v));
- i) D/ }# R! K+ d8 D    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
1 u) @4 R) Y/ m9 w    gen=gen+1;
end
figure(1);clf;
. ]: l( O/ E! d8 l0 }7 Q/ yplot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
legend('解的变化','种群均值的变化')
% I9 C5 f/ ^0 ?4 z% y, s5 @2 bxlabel('迭代次数');ylabel('高通目标函数值');
; L6 j5 x) l# z* u2 r. ]: o. z, E# jfigure(2);clf;
" y: X- M. b. O3 M6 J3 f, Pplot(trace2(:,1));hold on;
# P4 n6 H, t4 ]' m+ `plot(trace2(:,2),'-.');
plot(trace2(:,1),'.');
2 r. R4 I+ X% k0 A- s! F0 I( Aplot(trace2(:,2),'.');grid;
* W: a( C. X4 q% Plegend('解的变化','种群均值的变化')
# a, I4 D% n& |! N% |3 Q, y+ f( P0 Axlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
4 D. V! J  O. N! L2 A$ G) {' e; L1 t' Jplot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;
: Q* O8 n/ V0 y$ E+ G$ G2 I8 ~
' _$ N' q5 o4 E! e( v