MAT之SA：T1編寫主函數法和T2Matlab自帶的SA工具箱GUI法，兩種方法實作對二進制函數優化求解

輸出結果

%SA：T2法利用Matlab自帶的SA工具箱optimtool通過GUI調用@Jason_niu函數實作對二進制函數優化求解—Jason niu

function fitnessVal = Jason_niu( x )

https://imgconvert.csdnimg.cn/aHR0cHM6Ly9pbWFnZXMyMDE4LmNuYmxvZ3MuY29tL2Jsb2cvMTMxMzQ3NS8yMDE4MDMvMTMxMzQ3NS0yMDE4MDMwMTIxMTQwNjU3MS05NjI4MTkwNTUuZ2lm

實作代碼

%SA：T1法利用Matlab編寫主函數實作對定義域[-5,5]上的二進制函數求最優解—Jason niu

[x,y] = meshgrid(-5:0.1:5,-5:0.1:5);

z = x.^2 + y.^2 - 10*cos(2*pi*x) - 10*cos(2*pi*y) + 20;

figure

mesh(x,y,z)

hold on

xlabel('x')

ylabel('y')

zlabel('z')

title('SA：利用SA最優化，定義域[-5,5]上的二進制函數z =  x^2 + y^2 - 10*cos(2*pi*x) - 10*cos(2*pi*y) + 20的最大值—Jason niu')

maxVal = max(z(:));

[maxIndexX,maxIndexY] = find(z == maxVal);

for i = 1:length(maxIndexX)  

   plot3(x(maxIndexX(i),maxIndexY(i)),y(maxIndexX(i),maxIndexY(i)), maxVal, 'r*','linewidth',2)

    text(x(maxIndexX(i),maxIndexY(i)),y(maxIndexX(i),maxIndexY(i)), maxVal, {['    X: ' num2str(x(maxIndexX(i),maxIndexY(i)))];['    Y: ' num2str(y(maxIndexX(i),maxIndexY(i)))];['    Z: ' num2str(maxVal)]})

   hold on

end