0

このコードは問題なく動作しますが、最適化関数は初期条件と反復回数にfmincon依存するため、プロットは正しくありません。alpha ( ) と beta ( )x0の各値に対して、異なる初期条件で何度も最適化を実行して、正しい答えが得られていることを確認する必要があります。正確な答えを得るには、さらに反復が必要になる場合があります。abx0

x0、aおよびの異なる初期条件で最適化を実行できるようにしたいと考えていますb。

関数ファイル

function f = threestate2(x,a,b)
c1 = cos(x(1))*(cos(x(5))*(cos(x(9))+cos(x(11)))+cos(x(7))*(cos(x(9))-cos(x(11))))...
   +cos(x(3))*(cos(x(5))*(cos(x(9))-cos(x(11)))-cos(x(7))*(cos(x(9))+cos(x(11))));
c2=sin(x(1))*(sin(x(5))*(sin(x(9))*cos(x(2)+x(6)+x(10))+sin(x(11))*cos(x(2)+x(6)+x(12)))...
    +sin(x(7))*(sin(x(9))*cos(x(2)+x(8)+x(10))-sin(x(11))*cos(x(2)+x(8)+x(12))))...
    +sin(x(3))*(sin(x(5))*(sin(x(9))*cos(x(4)+x(6)+x(10))-sin(x(11))*cos(x(4)+x(6)+x(12)))...
    -sin(x(7))*(sin(x(9))*cos(x(4)+x(8)+x(10))+sin(x(11))*cos(x(4)+x(8)+x(12))));
f=(a*a-b*b)*c1+2*a*b*c2;

メインファイル

%x=[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8),x(9),x(10),x(11),x(12)]; % angles;

lb=[0,0,0,0,0,0,0,0,0,0,0,0];
ub=[pi,2*pi,pi,2*pi,pi,2*pi,pi,2*pi,pi,2*pi,pi,2*pi];
x0=[pi/8;0;pi/3;0;0.7*pi;.6;0;pi/2;.5;0;pi/4;0];
xout=[];
fout=[];
options = optimoptions(@fmincon,'Algorithm','interior-point','TolX',10^-10,'MaxIter',1500);
a=0:0.01:1;
w=NaN(length(a));

for i=1:length(a)
     bhelp=(1-a(i)*a(i));
if bhelp>0
    b=sqrt(bhelp);
       [x,fval]=fmincon(@(x)threestate2(x,a(i),b),x0,[],[],[],[],lb,ub,[],options);
w(i)=fval;
w(i)=-w(i);
B(i)=b;
else
w(i)=NaN;
B(i)=b;
end 
end
%surface(b,a,w)
%view(3)
%meshc(b,a,w)
x=a.^2;
plot(x,w)
grid on
ylabel('\fontname{Times New Roman} S_{max}(\Psi_{gs})')
xlabel('\fontname{Times New Roman}\alpha^2')
%ylabel('\fontname{Times New Roman}\beta')
title('\fontname{Times New Roman} Maximum of the Svetlichny operator(\alpha|000>+\beta|111>)')
4

1 に答える 1

1

Matlab Parallel Toolbox がある場合は、通常のループに似ているが並列で実行されるparforを使用できます。

それを使用するには、関数ですべての大きな厄介なスクリプトを作成する必要があります。A(i)初期条件を に保存し、結果を に保存すると仮定すると、次のB(i)ようなものを使用できます。

parfor i=1:length(B)
    B(i)=optimise(A(i));
end

ツールボックスがない場合は、他の方法 (MEX ファイルなど) がありますが、基本的に自分でスレッドを管理する必要があるため、お勧めしません。

于 2015-12-18T05:29:01.470 に答える