% rng(14，'twister') %取消上一行的注释以获得相同的输出ms = MultiStart;opts = optimoptions(@fmincon，'算法'，'内部点');sixmin = @ (x) (4 * x (1) ^ 2 - 2.1 * (1) ^ 4 + x(1) ^ 6/3……+ x(1)*x(2) - 4*x(2)^2 + 4*x(2)^4);问题= createOptimProblem('fmincon'，'x0'，[-1,2]，…sixmin“客观”,“磅”,[3 3],乌兰巴托,3,3,…“选项”,选择);[xminm,fminm,flagm,outptm,manyminsm] = run(ms,problem,50);所有颜色= 'kbgcrm';hold on for i = 1:size(manyminsm,2) %该行颜色cIdx = rem(i-1, length(所有颜色))+ 1; color = possColors(cIdx); % Plot start points u = manyminsm(i).X0; x0ThisMin = reshape([u{:}], 2, length(u)); plot(x0ThisMin(1, :), x0ThisMin(2, :), '.', ... 'Color',color,'MarkerSize',25); % Plot the basin with color i plot(manyminsm(i).X(1), manyminsm(i).X(2), '*', ... 'Color', color, 'MarkerSize',25); end % basin center marked with a *, start points with dots hold off