多目标优化

海贼骑士9235

在多目标优化时，优化结果没有同时让目标满足要求，下面是我的主代码和所得结果，请大神帮忙解答：
global u;
; ]! s, k+ u( u/ G9 b9 L. T5 U%贝塞尔函数求根
1 \, L9 V8 I' ]1 Q  R3 J* [" ~%设定模式数为20，对于TE模式需要求得一阶贝塞尔函数的前20个非零根
8 \" a' Q7 ?/ H, vmaxs=20;%要求的贝塞尔函数的非零根个数
incr=4.0;
- o, U( `( ?- Xu(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附近求零阶贝塞尔函数的第二个及以上的解
+ I% |9 Y, |- r* }! j/ y0 c( dend
NIND=40;      %个体数目40
MAXGEN=100;   %最大遗传代数100
NVAR=10;       %变量个数为10，前五个为半径，后五个位直波导长度
PRECI=20;      %变量的二进制位数
2 s) k  E" A4 OGGAP=0.9;      %代沟
trace1=[];
trace2=[];   %性能跟踪器，1、2分别为第一目标函数、第二目标函数的跟踪
2 k! y2 J. F/ A2 p8 ]+ x%建立区域描述器
3 M! w3 D  S" ^& j; ~) d1 _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])];
3 H8 E! K2 i/ t4 x5 @Chrom=crtbp(NIND,NVAR*PRECI);        %初始化种群
v=bs2rv(Chrom,FieldD)/1000;               %初始种群十进制转换(单位由mm转换为m)
gen=1;                               %代计数器
%遗传开始
( V; j* |8 s# E1 m; o' k/ Uwhile gen<MAXGEN
7 h/ R2 v8 P# z- G    [NIND,N]=size(Chrom);
: y: K) v6 @: u    M=fix(NIND/2);
    ObjV1=Highpass(v(1:M,);         %分组后高通目标函数值（单位m）
# x6 y7 h. h$ }% z; ?0 X: U    FitnV1=ranking(-ObjV1);           %分配高通适应度值（求最大值）
    SelCh1=select('sus', Chrom(1:M,,FitnV1,GGAP);      %调用随机遍历抽样函数sus对第一组种群基于适应度值进行选择
    ObjV2=Lowpass(v(M+1:NIND,);     %分组后低通目标函数值（单位m）
    FitnV2=ranking(ObjV2);            %分配低通适应度值（求最小值）
    SelCh2=select('sus', Chrom((M+1):NIND,,FitnV2,GGAP);      %调用随机遍历抽样函数sus对第二组种群基于适应度值进行选择
    SelCh=[SelCh1;SelCh2];            %将两组选择后种群进行合并
8 a& l7 H' |5 N3 p4 d6 |' z9 E0 }5 s    SelCh=recombin('xovsp',SelCh,0.7);%以0.7的交叉概率进行单点交叉重组
    Chrom=mut(SelCh);                 %变异
    %Chrom=reins(Chrom,SelCh);         %随机均匀选择重插入
    v=bs2rv(Chrom,FieldD)/1000;            %将进化后种群十进制转换(单位由mm转换为m)
# Q2 @8 s# W  S; J3 @    trace1(gen,1)=max(Highpass(v));
, l5 a. f; q, s* ^    trace1(gen,2)=sum(Highpass(v))/length(Highpass(v));
( K% A/ J+ ?0 o    trace2(gen,1)=min(Lowpass(v));
    trace2(gen,2)=sum(Lowpass(v))/length(Lowpass(v));
    gen=gen+1;
3 A6 D! y; L- |8 }6 V. z! n% send
figure(1);clf;
6 C3 X$ p2 `  _1 V# B8 Jplot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
) n" b" P5 q' C' `* y4 Uplot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
) i: S# o. o+ @$ Nlegend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('高通目标函数值');
: I+ h8 `& P, E7 bfigure(2);clf;
2 T2 ]+ L1 e4 i. z' tplot(trace2(:,1));hold on;
+ t- X, @5 c7 {3 E9 `- [plot(trace2(:,2),'-.');
0 J( k! r0 v' t3 A7 X3 m; W! kplot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
/ }# [3 w: o. a8 xlegend('解的变化','种群均值的变化')
xlabel('迭代次数');ylabel('低通目标函数值');
, w, z5 U' `# V" t# m6 x; ^$ q/ nfigure(3);clf;
plot(Highpass(v));hold on;
( Q' q+ |2 G. ^$ H$ Oplot(Lowpass(v),'r-.');grid;
2 s+ ?0 \0 m/ i# o/ H% {4 ^5 e