- Parent Category: Unconstrained
- Category: n-Dimensions
- Hits: 14289
Michalewicz's Function
I. Mathematical Expression:
$$f(X)=-\sum_{j=1}^{n}\sin\left(x_j\right)\left[\sin\left(\displaystyle\frac{jx^2_j}{\pi}\right)\right]^{2m} \ \ , \ \ \ m=10$$
where:
\(\bullet\) \(0\leq x_i\leq pi\) , \(i=1,2,\cdots,n\)
\(\bullet\) For \(n=1\):
\( \ \ \ \ \ \ \rhd\) \(f_{min}(x^*)=-0.801303410098552549\)
\( \ \ \ \ \ \ \rhd\) \(x^*=2.20290552017261\)
\(\bullet\) For \(n=2\):
\( \ \ \ \ \ \ \rhd\) \(f_{min}(X^*)=-1.80130341009855321\)
\( \ \ \ \ \ \ \rhd\) \(x^*_i=(2.20290552014618,1.57079632677565)\)
\(\bullet\) These optimal values have been obtained by us using MapleSoft 2015. For ensuring convergence to the global optimal location, the good initial point has been determined by the biogeography-based optimization (BBO) algorithm, which is coded in MATLAB 2011a.
II. Citation Policy:
If you publish material based on databases obtained from this repository, then, in your acknowledgments, please note the assistance you received by using this repository. This will help others to obtain the same data sets and replicate your experiments. We suggest the following pseudo-APA reference format for referring to this repository:
Ali R. Al-Roomi (2015). Unconstrained Single-Objective Benchmark Functions Repository [https://www.al-roomi.org/benchmarks/unconstrained]. Halifax, Nova Scotia, Canada: Dalhousie University, Electrical and Computer Engineering.
Here is a BiBTeX citation as well:
@MISC{Al-Roomi2015,
author = {Ali R. Al-Roomi},
title = {{Unconstrained Single-Objective Benchmark Functions Repository}},
year = {2015},
address = {Halifax, Nova Scotia, Canada},
institution = {Dalhousie University, Electrical and Computer Engineering},
url = {https://www.al-roomi.org/benchmarks/unconstrained}
}
III. 2&3D-Plots:
IV. Controllable 3D Model:
- In case you want to adjust the rendering mode, camera position, background color or/and 3D measurement tool, please check the following link
- In case you face any problem to run this model on your internet browser (it does not work on mobile phones), please check the following link
V. MATLAB M-File:
% Michalewicz's Function
% Range of initial points: 0 < xj < pi , j=1,2,...,n
% For n=1:
% Global minima: x=2.20290552017261
% f(x)=-0.801303410098552549 (determined by Maple 2015)
% For n=2:
% Global minima: (x1,x2)=(2.20290552014618,1.57079632677565)
% f(x1,x2)=-1.80130341009855321 (determined by Maple 2015)
% Coded by: Ali R. Alroomi | Last Update: 22 June 2015 | www.al-roomi.org
clear
clc
warning off
m=10;
x1min=0;
x1max=pi;
x2min=0;
x2max=pi;
R=1500; % steps resolution
x1=x1min:(x1max-x1min)/R:x1max;
x2=x2min:(x2max-x2min)/R:x2max;
for j=1:length(x1)
% For 1-dimensional plotting
f1(j)=-sin(x1(j))*(sin(1*x1(j).^2/pi)).^(2*m);
% For 2-dimensional plotting
for i=1:length(x2)
f(i)=-sin(x1(j))*(sin(1*x1(j).^2/pi)).^(2*m)-sin(x2(i))*(sin(2*x2(i).^2/pi)).^(2*m);
end
fn_tot(j,:)=f;
end
figure(1)
plot(x1,f1);set(gca,'FontSize',12);
xlabel('x','FontName','Times','FontSize',20,'FontAngle','italic');
ylabel('f(x)','FontName','Times','FontSize',20,'FontAngle','italic');
title('2D View','FontName','Times','FontSize',24,'FontWeight','bold');
figure(2)
meshc(x1,x2,fn_tot);colorbar;set(gca,'FontSize',12);
xlabel('x_2','FontName','Times','FontSize',20,'FontAngle','italic');
set(get(gca,'xlabel'),'rotation',25,'VerticalAlignment','bottom');
ylabel('x_1','FontName','Times','FontSize',20,'FontAngle','italic');
set(get(gca,'ylabel'),'rotation',-25,'VerticalAlignment','bottom');
zlabel('f(X)','FontName','Times','FontSize',20,'FontAngle','italic');
title('3D View','FontName','Times','FontSize',24,'FontWeight','bold');
figure(3)
mesh(x1,x2,fn_tot);view(0,90);colorbar;set(gca,'FontSize',12);
xlabel('x_2','FontName','Times','FontSize',20,'FontAngle','italic');
ylabel('x_1','FontName','Times','FontSize',20,'FontAngle','italic');
zlabel('f(X)','FontName','Times','FontSize',20,'FontAngle','italic');
title('X-Y Plane View','FontName','Times','FontSize',24,'FontWeight','bold');
figure(4)
mesh(x1,x2,fn_tot);view(90,0);colorbar;set(gca,'FontSize',12);
xlabel('x_2','FontName','Times','FontSize',20,'FontAngle','italic');
ylabel('x_1','FontName','Times','FontSize',20,'FontAngle','italic');
zlabel('f(X)','FontName','Times','FontSize',20,'FontAngle','italic');
title('X-Z Plane View','FontName','Times','FontSize',24,'FontWeight','bold');
figure(5)
mesh(x1,x2,fn_tot);view(0,0);colorbar;set(gca,'FontSize',12);
xlabel('x_2','FontName','Times','FontSize',20,'FontAngle','italic');
ylabel('x_1','FontName','Times','FontSize',20,'FontAngle','italic');
zlabel('f(X)','FontName','Times','FontSize',20,'FontAngle','italic');
title('Y-Z Plane View','FontName','Times','FontSize',24,'FontWeight','bold');
VI. References:
[1] M. Molga and C. Smutnicki, "Test Functions for Optimization Needs," Apr. 2005, [Accessed March 28, 2013]. [Online]. Available: http://www.zsd.ict.pwr.wroc.pl/files/docs/functions.pdf
[2] S. Rahnamayan, H. R. Tizhoosh, and M. M. A. Salama, "A Novel Population Initialization Method for Accelerating Evolutionary Algorithms," Computers & Mathematics with Applications, vol. 53, no. 10, pp. 1605-1614, May 2007.
[3] Ali R. Alroomi, "The Farm of Unconstrained Benchmark Functions," University of Bahrain, Electrical and Electronics Department, Bahrain, Oct. 2013. [Online]. Available: http://www.al-roomi.org/cv/publications
