多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
2 H: ^* ~' J3 A' p2 o- nglobal u;
9 l  q2 t# ]  w5 F2 O6 Q%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
( M( r' ~4 {3 W+ |& K  smaxs=20;%要求的贝塞尔函数的非零根个数
# ?8 `1 v, S- S- K: X* k$ l0 M# Qincr=4.0;
8 e( Z9 f5 N3 R9 Mu(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
4 z# ?& }' \' s7 j  c( s0 q2 qfor s=2:maxs
    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
end
# v( C7 @: ~' O3 cNIND=40;      %个体数目40
5 z2 X1 G! ~! ^3 `( gMAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
- C! O: ~7 d9 ~: DGGAP=0.9;      %代沟
trace1=[];
0 |( e5 ]6 }9 B' v! w# Q* F' f( `trace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
%建立区域描述器
/ u8 [! f7 q/ PFieldD=[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)
gen=1;                               %代计数器
%遗传开始
' L/ K) H. y& k) d$ A$ F- m, gwhile gen<MAXGEN
    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
' q; j' z# _; ^" \: G% ^    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
: _, I/ T9 S; t! ~" L    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
1 P5 w( W$ x( z' a0 }; A    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
! i2 d: f- J' v; }: q    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
7 L% S6 x6 e/ T+ h7 W    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
# H" u8 {2 {- j5 ^% |/ T0 y$ m; M    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
9 A6 w: r' w6 Q& Z/ `: ?" K    trace1(gen,1)=max(Highpass(v));
- \. v" O  s7 w! z    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
! Y9 l, e! ~4 Z* H! Y, ?    gen=gen+1;
end
- k1 b, d! {7 h, Ffigure(1);clf;
& y* {0 A- ]4 `7 ]) Splot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
legend('解的变化','种群均值的变化')
2 b4 A" j! C0 R# ~  v1 u( m8 cxlabel('迭代次数');ylabel('高通目标函数值');
figure(2);clf;
plot(trace2(:,1));hold on;
' N* d- I# W6 F6 ^9 M% F7 Aplot(trace2(:,2),'-.');
; r4 k4 x, b5 L$ a. bplot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
3 o" U9 K9 ]/ Y' t. `. l) [7 ]xlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
plot(Highpass(v));hold on;
$ F5 {6 c# W. `3 E7 Splot(Lowpass(v),'r-.');grid;
& V' E% _2 a7 I: q