多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
$ l" Z4 }5 s, Q0 P8 U9 p9 ]3 a/ zglobal u;
* a0 ~( ~  e- S%贝塞尔函数求根
* L$ c7 T9 B( h7 t%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
4 I6 X% h7 x  Gu(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
for s=2:maxs
    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
8 w! @! A7 x5 I8 o' l7 yend
NIND=40;      %个体数目40
9 S) z9 n+ X: rMAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
" `3 ?# \% x" O) r# QGGAP=0.9;      %代沟
trace1=[];
' |4 Y; w0 v3 K: I. F6 ftrace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
- e0 B) f" K# b! q) r  x%建立区域描述器
, U9 V2 G# N( hFieldD=[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])];
2 J7 @# K" a* m" ?& K4 OChrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
) q& W2 [( n2 a5 ?5 h( X6 hv=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
* Z( _2 e0 {; R- f) N8 Y% V, }& [5 Zgen=1;                               %代计数器
%遗传开始
while gen<MAXGEN
    [NIND,N]=size(Chrom);
6 S7 n, R$ a& P    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
6 @" o+ N' @3 U; o    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
" k+ a; n# f5 c    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
% _( j3 m% @- E( B, b6 X4 f% K    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
    trace1(gen,1)=max(Highpass(v));
* S, ?  ^# N. M) B- s6 V5 }: f    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
    trace2(gen,1)=min(Lowpass(v));
0 ^( g3 M+ m% F) p    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
0 D# X. z9 A- q    gen=gen+1;
end
- p# [! A2 U: v8 n, Z* Rfigure(1);clf;
; C2 S  I# h& B" bplot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
# j3 o5 t' J) x0 l: ulegend('解的变化','种群均值的变化')
, e+ z" Y/ a2 n, n3 d4 R6 r3 Zxlabel('迭代次数');ylabel('高通目标函数值');
figure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
% z: h" k& g$ ?* W# D0 qplot(trace2(:,1),'.');
4 G  N, J; T9 Bplot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
. s" l! ~5 I6 q( ixlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
plot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;
, G0 F. \' W/ \9 j0 V1 L