多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
" h% y& J8 v2 u# hglobal u;
%贝塞尔函数求根
5 O5 h2 v+ d7 y%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
maxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
u(1,1)=fzero(@(x)besselj(1,x),3.9);%一阶贝塞尔函数的第一个解（在3.9附近）
2 i) E1 ?9 P, zfor s=2:maxs
2 {. Z$ j  |$ d; M- d; |5 v    u(1,s)=fzero(@(x)besselj(1,x),u(1,s-1)+incr);%在第一个解加4附近求零阶贝塞尔函数的第二个及以上的解
: ^& s1 u3 q; M  Oend
+ g1 x2 z; s& p. GNIND=40;      %个体数目40
$ @/ F% k" j$ g$ x7 oMAXGEN=100;   %最大遗传代数100
- q, L( L- B$ b& DNVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
: U* Z! X7 Z0 iPRECI=20;      %变量的二进制位数
GGAP=0.9;      %代沟
* g0 i* {+ K  E5 C7 |; Ttrace1=[];
8 ?' l3 o# `( K/ a$ d( Vtrace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
%建立区域描述器
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);        %初始化种群
v=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
gen=1;                               %代计数器
: ]6 k- ^4 [$ m9 v. n2 f9 i. P%遗传开始
while gen<MAXGEN
    [NIND,N]=size(Chrom);
    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
6 t) Y  [% N7 |: y, @    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
: u" Q- I# \" g( q9 f    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
: b: Z7 `6 k9 i& O    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
/ p- m4 {( ~: H! s    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
$ s4 }. [; ^* N7 c    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
    trace1(gen,1)=max(Highpass(v));
# \+ d1 [+ L! m) b2 l. k; C    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
% a: C) G' m+ X    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
    gen=gen+1;
  X, ?# ]- N0 S0 U! n; ?1 q7 ~( Qend
# o6 a/ {/ @; v' I) _7 |; m5 jfigure(1);clf;
! W' \/ h) ^# |; \+ {plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
+ u! ?& L" f3 G2 n' Hplot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
legend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('高通目标函数值');
9 j3 g6 h9 N: Y6 tfigure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
  P* y% \1 B3 Q3 A' \, y6 N0 `plot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
8 M! i8 Q/ {: ]( ?7 a6 zlegend('解的变化','种群均值的变化')
3 o, n! Z; \8 {7 o6 z, |xlabel('迭代次数');ylabel('低通目标函数值');
) [( B$ ?- v" afigure(3);clf;
$ S& f* C  V7 z" s2 Q. A. a5 _plot(Highpass(v));hold on;
plot(Lowpass(v),'r-.');grid;
# n9 Y9 q: I; r