多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
8 ]: l* j0 W: i' kglobal u;
%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
4 o! m: m2 U/ C1 G! s  ?maxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
u(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
for s=2:maxs
" V( N6 ]3 E3 W) S    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
end
NIND=40;      %个体数目40
) S. Z! m$ [, a. E0 W0 {/ h- L3 lMAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
; v4 P% t, F% B% @1 w. m% APRECI=20;      %变量的二进制位数
! r" ]$ c+ h4 v7 E; bGGAP=0.9;      %代沟
1 T, U" C( r& e1 U; ]trace1=[];
trace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
" @6 I5 [; h2 j+ a%建立区域描述器
& Y# T' J4 u# v. ^$ R% LFieldD=[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);        %初始化种群
v=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
8 J) S0 X, O+ ^gen=1;                               %代计数器
+ L! S; ^+ f& |1 D%遗传开始
+ |2 i  b, A1 qwhile gen<MAXGEN
- D/ A2 P% i1 P, U; T  }    [NIND,N]=size(Chrom);
; e- s. Q; n7 B) a8 R    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
9 a; `# l( }4 }7 P    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
  b4 }$ ^  B3 k" A& b9 B& }    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
6 i$ ]& r; Q  M& ~# f& S* C% W    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的交叉概率进行单点交叉重组
0 z6 ^+ S" }1 C( b% r& ]& w    Chrom=mut(SelCh);                 %变异
, X+ O+ L7 |7 l  E8 }    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
" r  w. a& i( D" l6 B- ~  X    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
    trace1(gen,1)=max(Highpass(v));
0 J) A! z5 a( J9 W% B; F- C; ?    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
& W# c1 I6 m2 }, \    trace2(gen,1)=min(Lowpass(v));
6 b6 E$ Z! P) h) i7 X. H+ O    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
& M2 V9 A$ I  i( `7 O    gen=gen+1;
end
figure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
legend('解的变化','种群均值的变化')
  J* t4 C: n: W; J9 `5 x7 f; fxlabel('迭代次数');ylabel('高通目标函数值');
8 Q6 n5 k2 @' X. S! O6 ofigure(2);clf;
/ g# A' j# Q9 w  o( C! ^+ M  vplot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
5 n" o8 Z+ I; N- d9 Kplot(trace2(:,1),'.');
' y- D& y% I) |! v# R; |* e, Rplot(trace2(:,2),'.');grid;
& p. x2 v- p* U# k( u! S; zlegend('解的变化','种群均值的变化')
7 y* V/ N0 J" r0 U$ zxlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
plot(Highpass(v));hold on;
' d" ~( A# I0 x5 Dplot(Lowpass(v),'r-.');grid;
1 a5 V9 g' ]6 |; ~