多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
global u;
; W+ y; S2 j. {# n0 w%贝塞尔函数求根
7 O. L+ Z, i) a0 b( b%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
& d; X( ?3 ?" e# Jincr=4.0;
' c- @9 [& H  _' M, `. ku(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
' z- W& O9 B; d0 H5 J( l; c" pfor s=2:maxs
    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
3 P; }. _/ u& V/ t/ j" aend
/ o0 h+ F8 C& aNIND=40;      %个体数目40
6 I" P. j% W- Z; a5 QMAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
# r; q: C' V8 x( D# Ntrace1=[];
9 H2 r* E* a( ptrace2=[];   %性能跟踪器，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])];
% J) O# ~8 i2 @" t) GChrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
. s; n2 X3 l$ e, P( wv=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
gen=1;                               %代计数器
5 N1 q( r  F, T  F%遗传开始
2 w% {8 ?' G& v) w( l  v; Hwhile gen<MAXGEN
    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
7 U$ g0 d) n8 E    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
5 D, D% {9 Q: k  N. I/ Q    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
- M! B% I% g4 F8 I    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
' ]  Z9 u: a" O8 U3 a6 x9 _# |    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
/ u. K# w* ~' \$ W! ~    trace1(gen,1)=max(Highpass(v));
2 S0 _9 W0 J5 z9 w3 H    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
4 v7 _! X  v8 Z( B1 H0 w- S    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
    gen=gen+1;
end
: L; x1 |  p1 I7 X( pfigure(1);clf;
, L, n* ?2 {, M3 |: `/ G2 p; Xplot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
, F! [  x8 o8 t! H5 o, Vlegend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('高通目标函数值');
, [8 D8 R+ }& K$ Afigure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
plot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
8 P$ ]/ Q- k1 h2 Lxlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
: V3 J1 y- E2 c- ^9 w) eplot(Highpass(v));hold on;
2 F. v: \! S* J, `. Jplot(Lowpass(v),'r-.');grid;

4 F! H& r: M  y0 W/ R( c4 n