多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
, }% n2 q0 V% Pglobal u;
%贝塞尔函数求根
%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
5 Q* E4 H6 o+ f, E3 L9 d, i7 [incr=4.0;
u(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
& C' U) D$ o; pNIND=40;      %个体数目40
MAXGEN=100;   %最大遗传代数100
- E* N: L4 k. K/ W! wNVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
trace1=[];
/ y* h9 [/ Q; ~: y1 A% otrace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
; a7 q. [* S# k1 P6 {%建立区域描述器
7 Y& O" ^9 ^/ F! Z! lFieldD=[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])];
. A0 g: i5 x- T7 M2 D. Z* Q" U8 TChrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
v=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
8 \) K! e9 {9 e; a% mgen=1;                               %代计数器
%遗传开始
while gen<MAXGEN
# }7 n9 A% P3 h0 J    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
, u$ s$ y" ]: D; M* g/ K! s    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
) i! @9 d9 ~0 c$ I9 ], q    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
; K+ `6 n* y# q7 z- K    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
6 L7 j* }% a3 t2 A' y5 S4 P7 _, q    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
1 l8 y; T% z( t% I5 k3 w    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
6 R9 i6 T1 x5 P( c4 r0 V; S    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
, P$ C! t! J8 L5 `8 `! f7 u; \    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
% {3 H/ @& g, p2 i( S, a% y; c    Chrom=mut(SelCh);                 %变异
! N" B  j9 e! ]( f    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
( L3 t& S0 x5 L* U; a% T    trace1(gen,1)=max(Highpass(v));
    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
    gen=gen+1;
6 G5 D7 D9 q7 \9 ?end
  {* @0 q/ F6 V: K0 ?+ C/ J$ g  ifigure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
& Q* h, h) r2 L1 p+ B7 {legend('解的变化','种群均值的变化')
2 ^( _" N8 [2 `/ V/ A( Gxlabel('迭代次数');ylabel('高通目标函数值');
figure(2);clf;
plot(trace2(:,1));hold on;
6 g( t+ [  D9 @plot(trace2(:,2),'-.');
plot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的变化','种群均值的变化')
; r/ e% \4 R7 ^7 yxlabel('迭代次数');ylabel('低通目标函数值');
figure(3);clf;
plot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;

; |$ k3 x% n4 o