多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
( M4 A6 R/ F/ ]; Uglobal u;
  u& V/ `% X  W0 {%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
u(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
1 z  J  o. x6 ~  O0 |$ @for s=2:maxs
    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
: N/ M( I7 H; A4 }; Y, g2 mend
9 E* F# |" c& K3 `NIND=40;      %个体数目40
# u' S; C0 i3 j! k$ ]MAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
, N- }$ ?/ J6 A6 s. e, M9 l* YGGAP=0.9;      %代沟
) j, x" ~" j. d2 |trace1=[];
- o5 g9 b" Q1 ?- B. u1 qtrace2=[];   %性能跟踪器，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])];
9 C6 }7 {% S9 G2 N9 z9 E3 p" LChrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
1 C2 t, f; b4 Z# h- cv=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
gen=1;                               %代计数器
3 S/ I* Q5 b( f$ y; k. H# D+ s+ S%遗传开始
while gen<MAXGEN
$ Y' j' l5 y2 z  |( L    [NIND,N]=size(Chrom);
/ o8 K6 L( J; ^    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
/ P" s1 b- P/ |! @' k' a3 q7 W    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
# K; K, Q! `1 N    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
& ~$ Y3 {: S7 w0 w* g    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
6 u6 t7 t/ x  [7 ?8 k) P6 P! w    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
9 G" T& A$ V! G# y  F    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
6 I7 P& _2 d8 @" r: t    trace1(gen,1)=max(Highpass(v));
: B: p" R1 Y5 O2 B  Q9 N    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
8 ]9 }, i# ^0 Y    gen=gen+1;
end
figure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
! t) J0 H; T2 l# Y, V" clegend('解的变化','种群均值的变化')
; Z) l9 t! U+ z* kxlabel('迭代次数');ylabel('高通目标函数值');
7 T6 P* N" ?8 b& Gfigure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
. A) u; ?; o) J" F4 O0 d5 v; Iplot(trace2(:,1),'.');
* a2 b" i+ d' F1 Z$ z! }& Q! nplot(trace2(:,2),'.');grid;
. O& m# |/ s/ ]! E/ e$ zlegend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('低通目标函数值');
6 r, w& P6 `. U7 I1 jfigure(3);clf;
plot(Highpass(v));hold on;
4 \6 \9 R' K+ @$ J  r7 K, nplot(Lowpass(v),'r-.');grid;