多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
2 o3 e( b+ D8 @( O0 r0 s7 @global u;
( ^  \7 V/ w4 P& q+ Z+ t% L0 |%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
3 H% a. K6 V4 b; X$ P9 Omaxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
; \/ Y9 E+ E) d/ E& Su(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
for s=2:maxs
# j' W6 W+ H$ b% R( N6 r    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
end
$ {# n) H" i' i: ]9 X, uNIND=40;      %个体数目40
MAXGEN=100;   %最大遗传代数100
& z" k" e* W9 u* QNVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
: `2 C1 w# e, o9 {6 Ctrace1=[];
3 ]. w9 U. S$ D% Z% m" Z4 n; ]trace2=[];   %性能跟踪器，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])];
! z4 E) K; l& |& u$ B" K2 ^0 RChrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
# J& q1 S" d" T6 V3 h# r" N: [8 |v=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
gen=1;                               %代计数器
%遗传开始
/ B/ h9 l. T8 J% o( }while gen<MAXGEN
8 g& S% \( Z: [1 [2 x7 \+ Q5 L9 b    [NIND,N]=size(Chrom);
2 k) G* k7 j! f- U: c! a    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
9 f& D, j% W) t    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
* Q( @; k' ?1 F! E' z( F1 l4 J. G    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对第二组种群基于适应度值进行选择
    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
. ]0 i1 ^. U, b- ?7 X8 s$ x& _    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
# w3 V; B5 e' W8 }8 U; l7 I, [5 `9 ?    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
    trace1(gen,1)=max(Highpass(v));
+ p# w0 N: ~! X  i    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
    trace2(gen,1)=min(Lowpass(v));
3 U. S; |) Q% J5 h    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
: G; S/ m/ e" J: H3 B    gen=gen+1;
end
% A1 G# D: `2 c! zfigure(1);clf;
* S+ |. ]4 i9 ]: U4 H" Kplot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
8 H  k+ m1 Y" z" [( j2 U- U& L* ]( ^plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
3 b; K; N% I+ }  r8 h! Z) Tlegend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('高通目标函数值');
! {7 l$ l* _1 X6 K) C. h  ]+ G; qfigure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
4 e4 S$ z' O: ]$ a- cplot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
( f2 L2 Z) `) m2 X, |! E2 Uplot(Highpass(v));hold on;
) U$ C% v" |' l1 G. @1 `1 \plot(Lowpass(v),'r-.');grid;