GA Vectorized do not return the best FF value

1 次查看（过去 30 天）

显示更早的评论

Lucas Carneiro 2018-2-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/384967-ga-vectorized-do-not-return-the-best-ff-value

编辑： Lucas Carneiro 2018-3-18

在 MATLAB Online 中打开

Hello, I am having a similar problem to this:

https://www.mathworks.com/matlabcentral/answers/367573-ga-didn-t-return-the-x-for-the-best-fitness-value

In my case, I am running GA Vectorized both with ga and gamultiobj. All is fine when Vectorization is off but the solver does not return the best solution with Vectorization on. I am asking myself how to solve this because I need to speed up my optimization.

Here is how my FF is implemented:

function FF=multiobj_mics(pop,n,r,max_d,max_o)
popSize=size(pop,1);
FF=zeros(popSize,2);
k=nchoosek(n,2);
%%Variables
%rad
%meters
%phi=2*pi*rand(1,n);      %azimuth
%theta=2*pi*rand(1,n);    %elevation
%phi=x(1)....x(n)         %-pi pi
%theta=x(n+1)....x(2*n)   %0 pi
%%Orientation
vet(1:n-1,1:3,1:n-1,1:popSize)=NaN;
for i=1:1:n-1
    x1=r*sin(pop(1:popSize,n+i)).*cos(pop(1:popSize,i));
    y1=r*sin(pop(1:popSize,n+i)).*sin(pop(1:popSize,i));
    z1=-r*cos(pop(1:popSize,n+i));
      p=1;
      for j=i+1:1:n
          x2=r*sin(pop(1:popSize,n+j)).*cos(pop(1:popSize,j));
          y2=r*sin(pop(1:popSize,n+j)).*sin(pop(1:popSize,j));
          z2=-r*cos(pop(1:popSize,n+j));
          %Normalized mic orientation vector
          vet(p,:,i,:)=reshape(([x1 y1 z1]/norm([x1 y1 z1])-[x2 y2 z2]/norm([x2 y2 z2])),[1,3,1,popSize]);
          %Cannot use ./
          for ii=1:popSize 
              vet(p,:,i,ii)=vet(p,:,i,ii)/norm(vet(p,:,i,ii));
          end
          p=p+1;
      end
  end
%Convert 4D matrix to 3D matrix (vet(3 lines) x pop)
Vet=reshape(vet,[],1);
Vet=reshape(Vet,n-1,3,popSize*(n-1));
C=permute(Vet,[1 3 2]);
C = reshape(C,[],size(Vet,2),1);
[rr,c] = size(C);
Vet   = permute(reshape(C',[c,rr/popSize,popSize]),[2,1,3]);
null=any(isnan(Vet),2);
Vet(find(null(:,1,1)==1),:,:) = [];
% Minimum angle btw pairs of mics
for i=1:k
      angle=zeros(popSize,k-1);
      p=1;
      for ii=1:k
          if ii~=i
              for j=1:popSize
                  angle(j,p)=atan2(norm(cross(Vet(i,:,j),Vet(ii,:,j))), dot(Vet(i,:,j),Vet(ii,:,j)));
                  if angle(j,p)>pi/2
                      angle(j,p)=pi-angle(j,p);
                  end
              end
              p=p+1;
          end
      end
      FF(1:popSize,1)=FF(1:popSize,1)+min(angle(:,:),[],2);
end
FF(1:popSize,1)=-max_o*FF(1:popSize,1);
}
    if true
      % code
    end