多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
global u;
%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
; v/ p& n5 R2 a, v1 n( g/ @/ _8 nmaxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
6 t4 b% u3 Z9 R  r+ ku(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附近求零阶贝塞尔函数的第二个及以上的解
end
NIND=40;      %个体数目40
MAXGEN=100;   %最大遗传代数100
/ {/ Y0 L2 J: a  m. V8 X3 MNVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
4 k( Z# t# r1 z9 J$ tPRECI=20;      %变量的二进制位数
. N5 m" O; W/ MGGAP=0.9;      %代沟
3 J0 ~0 a& d( S4 P2 a& Ztrace1=[];
trace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
$ K$ N. ?+ ^& c+ U' T! }4 m( C%建立区域描述器
, w. L2 |) T5 C4 f+ Y% U  sFieldD=[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])];
  e1 k* e2 ?8 B' z' BChrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
v=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
gen=1;                               %代计数器
%遗传开始
& K  [% l. s+ Lwhile gen<MAXGEN
    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
; k2 K" C" u# U& K7 {0 c) s3 ]    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
& W. @  [! P7 f; o3 J4 W# h    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
9 `# L$ f" ^! M  N, n    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
4 B2 ]6 R- w: r( K* F; C: J+ y    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
: `$ w# G- ], K; }% O0 H- L  c    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
, |$ ]- n9 Z8 Y# J    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
    trace1(gen,1)=max(Highpass(v));
    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
: t  `  D/ I  o9 x5 o: R0 N    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
( k6 ~, h, |6 d9 a8 H6 D$ E    gen=gen+1;
2 s, T- O* @# ~3 Y# K- K0 r2 p8 Gend
figure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
" F3 ~0 `/ k" K$ {, ilegend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('高通目标函数值');
figure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
plot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
: p' D7 A$ ~3 a2 F0 f9 ulegend('解的变化','种群均值的变化')
4 `" V5 P+ K  ~5 W3 F! fxlabel('迭代次数');ylabel('低通目标函数值');
: a. k0 o" }- r/ Xfigure(3);clf;
3 c% B' r1 ^  F3 f; t) @plot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;
7 l) a; e$ R( t2 G2 r) Y
# \. ]7 K; ~6 w' h$ w1 O3 B