多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
global u;
%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
4 Q6 b" J9 I4 E# h- l- kincr=4.0;
u(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
& d+ S! W& l: E' ]% o. xfor s=2:maxs
    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
end
NIND=40;      %个体数目40
2 B. {" D6 D8 l& m! |MAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
trace1=[];
trace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
+ @4 _. N# V1 \%建立区域描述器
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])];
Chrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
8 ^* j0 i: \$ Y$ }* Av=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
9 W5 \- H5 I5 P3 i; Mgen=1;                               %代计数器
) b0 ?/ q0 E3 B/ E7 J) w- J%遗传开始
% B7 M, R6 ?, d( h3 o- F+ f& Qwhile gen<MAXGEN
    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
' y: J3 L1 y# e! A( w    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
- Y* J9 m/ B1 B8 g: Y1 h1 \) U    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
9 W# _+ p. D* J# w' c- |- m    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
$ ]+ A9 Y8 F7 K    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
4 v3 n/ ^+ _. Z% {6 t    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
* p# j& N- ~* `- _2 T" z    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
$ W4 _, [1 O' k% d) m    Chrom=mut(SelCh);                 %变异
& k& u3 {1 ]! {    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
1 k: [2 P) |0 G: m- K    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
# V+ a# }8 n4 |( R3 l    trace1(gen,1)=max(Highpass(v));
* u1 |8 H) \0 e4 V, f7 w7 W8 @! u    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
! I* ~; D* k, c# y, K! p# K. K4 L    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
    gen=gen+1;
) A" a0 ^6 K% U7 V' Z/ cend
/ W$ q5 d$ }9 |( ffigure(1);clf;
- A, A3 K- R/ Dplot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
! m4 {) Q4 F. q- a+ y  g- f. P" Rplot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
legend('解的变化','种群均值的变化')
6 H) r5 j' m7 X5 U% U0 q" @' F5 p; |7 Oxlabel('迭代次数');ylabel('高通目标函数值');
figure(2);clf;
) M( R* ]$ u2 H" hplot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
& _5 U* V# e1 }: _0 tplot(trace2(:,1),'.');
: K; i7 K8 x2 R5 a/ `( ?plot(trace2(:,2),'.');grid;
" ?3 z; C, y! Q4 y( t" R/ h$ a7 l/ ~% Glegend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('低通目标函数值');
4 e' G9 A9 Z  Bfigure(3);clf;
! ~0 B0 B$ v* ^( W  {0 o2 Uplot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;

( t& v! T) G" s: h3 a