多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
global u;
%贝塞尔函数求根
& {. ^, W9 B' U1 q: T2 U, U: ^, S%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
6 p: m0 B- Z  }' {" ^8 k% x1 Ou(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
! A# j) N* ]( P# lfor s=2:maxs
    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
1 d" j! q; X6 W& ~0 tend
0 n2 Y, l* ~$ a% x1 WNIND=40;      %个体数目40
MAXGEN=100;   %最大遗传代数100
; e; V1 z1 G3 [' _: Q4 V% y9 yNVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
# w( D( T8 q' TPRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
+ `5 h( I& x7 g5 W  p, q) y, a& u8 Utrace1=[];
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])];
Chrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
0 M# w$ I  j8 G( @0 Cv=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
' S( s% q4 e8 F" t' S, B' o1 C' Ggen=1;                               %代计数器
%遗传开始
while gen<MAXGEN
, V( ?% V% N6 @( Q. p1 ?4 J    [NIND,N]=size(Chrom);
- Q- [+ E/ T1 M0 |! j    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
% {1 e# D3 S! M+ K: r, s9 ~    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
8 n7 o' S6 ]& [6 G2 g0 T) K8 }    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
4 g" U( M# v/ n/ F    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
( }, `4 u$ l# C) r4 X0 y    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
1 ?3 Y7 H) L* f    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
8 N+ a8 T! C0 A+ `- m! }    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
' H( S) e; ~- |; f7 }    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
( {9 _4 }( y4 L8 k, x, `    trace1(gen,1)=max(Highpass(v));
    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
; _& |0 _( U# S/ k    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
& r- N4 T5 ]& K: J9 m    gen=gen+1;
end
figure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
: H( B, b, u) \% i9 uplot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
" W# {* j" U# F8 M. vlegend('解的变化','种群均值的变化')
' i/ g2 T6 c; r3 J/ [2 B: x/ Xxlabel('迭代次数');ylabel('高通目标函数值');
0 h' _% z0 q& g1 U0 K; K  g2 ]figure(2);clf;
; n: Z# Y" \9 _9 g/ d) ?5 G$ Tplot(trace2(:,1));hold on;
& @- l5 s6 C; l  k# uplot(trace2(:,2),'-.');
! L; n! s2 V+ U% E9 v6 C/ }6 dplot(trace2(:,1),'.');
- H% n: p, {  x1 Mplot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
) o, w) s, V. K# H* rxlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
plot(Highpass(v));hold on;
$ G2 c' Y( b  \2 ]plot(Lowpass(v),'r-.');grid;
" \: b; e# q( p* j; E6 a
  |( c, ]2 c! N" D1 d2 d