0

此代码工作正常，但绘图不正确，因为优化函数fmincon将取决于初始条件x0和迭代次数。对于 alpha ( a) 和 beta ( b) 的每个值，我应该使用不同的初始条件多次运行优化，x0以验证我得到了正确的答案。可能需要更多的迭代才能获得准确的答案。

我希望能够在 和 的不同初始条件下运行x0优化。ab

函数文件

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 回答