<h1>Four heuristic optimization algorithms applied to wind energy: determination of Weibull curve parameters for three Brazilian sites</h1>
			<ul class="item-list">
	        	        <li>
	            Carla Freitas de Andrade	            	            	            	            <sup aria-label="Affiliated with Mechanical Engineering Department, Federal University of Ceará, Fortaleza, Ceará, BR">
	                1	            </sup>
	            	        </li>
	        	        <li>
	            Lindemberg Ferreira dos Santos	            	            	            	            <sup aria-label="Affiliated with Mechanical Engineering Department, Federal University of Ceará, Fortaleza, Ceará, BR">
	                1	            </sup>
	            	        </li>
	        	        <li>
	            Marcus V. Silveira Macedo	            	            <abbr title="This is the corresponding author for this article">*</abbr>
	            	            	                <a href="mailto:marcusmacedo87@gmail.com" class="tiny-icon email-link mx-1" title="Email Marcus V. Silveira Macedo">
	                    Email
	                </a>
	            	            	            <sup aria-label="Affiliated with Mechanical Engineering Department, Federal University of Ceará, Fortaleza, Ceará, BR">
	                1	            </sup>
	            	        </li>
	        	        <li>
	            Paulo A. Costa Rocha	            	            	            	            <sup aria-label="Affiliated with Mechanical Engineering Department, Federal University of Ceará, Fortaleza, Ceará, BR">
	                1	            </sup>
	            	        </li>
	        	        <li>
	            Felipe Ferreira Gomes	            	            	            	            <sup aria-label="Affiliated with Mechanical Engineering Department, Federal University of Ceará, Fortaleza, Ceará, BR">
	                1	            </sup>
	            	        </li>
	        	    </ul>
	    	    <ul class="affiliations" aria-hidden="true">
	        <li>Mechanical Engineering Department, Federal University of Ceará, Fortaleza, Ceará, BR</li>	    </ul>
	    
<h2>Abstract</h2>
<p>Minimizing errors in wind resource analysis brings significant reliability gains for any wind power generation project. The characterization of the wind regime is one of fundamental importance, and the two parameters Weibull distribution is the most applied function for it. This study aims to determine the scale and shape factor in an attempt to establish acceptable criteria to a better utilization of wind power in the states of Pernambuco and Rio Grande do Sul, which is a national prominence in the use of renewable sources for electricity generation in Brazil. The following heuristic optimization algorithms were applied: Harmony Search, Cuckoo Search Optimization, Particle Swarm Optimization and Ant Colony Optimization. The fit tests were performed with data from the Brazilian Federal Government’s SONDA (National System of Environmental Data Organization) project, referring to Triunfo, Petrolina and São Martinho da Serra, states of Pernambuco and Rio Grande do Sul, cities in the northeast and south regions of Brazil, during the period of 1 year. The tests were made in 2006 and 2010, all at 50 m from ground level. The results were analyzed and compared with those obtained by the maximum likelihood method, moment method, empirical method and equivalent energy method, methods that presented significant results in regions with characteristics similar to the regions studied in this study. The performance of each method was evaluated by the RMSE (root mean square error), MAE (mean absolute error), 
<em>R</em>
<span id="IEq1"><mml:math id="IEq1_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq1_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq1.gif"/></span>
 (coefficient of determination) and WPD (wind production deviation) tests . The statistical tests showed that ACO is the most efficient method for determining the parameters of the Weibull distribution for Triunfo and São Martinho da Serra and CSO is the most efficient for Petrolina.</p><hr/><section><h2>Introduction</h2>
<p>Search for energy forms that reduce or eliminate the carbon dioxide emission to the atmosphere has encouraged the renewable energy sector development, with the wind energy being highlighted. According to the World Wind Energy Association, the installed capacity of wind power in the world reached 486.661 MW at the end of 2016, 54.846 MW more than in 2015, representing a growth rate of 11.8%.</p>
<p>Wind resource analysis is a key step in the wind power generation projects development. Reducing errors in this step brings significant reliability gains for the project. One of the most important information in the wind resource analysis is the characterization of the wind regime according to a probability distribution, which aims to transform the discrete data collected in a measurement campaign into continuous data. In this procedure, velocities are grouped in intervals and a probability distribution function is fitted to this histogram. Depending on the wind conditions, the curve to be adjusted may follow the Gauss, Rayleigh or, more commonly, two parameters 
<em>k</em>
 and 
<em>c</em>
 Weibull distributions [
<a href="#CR35"><sup>35</sup></a>
].</p>
<p>One of the challenges in applying the Weibull distribution to represent the region wind regime is the estimation of the parameters 
<em>k</em>
 and 
<em>c</em>
, and an adjustment must be obtained with the smallest possible error. Dorvlo [
<a href="#CR12"><sup>12</sup></a>
] used the Chi-square method to determine the Weibull parameters in four locations in Oman and Saudi Arabia. Silva [
<a href="#CR32"><sup>32</sup></a>
] presented the equivalent energy method, where the parameters are found from the square error minimization power. Akdag and Dinler [
<a href="#CR1"><sup>1</sup></a>
] proposed the energy pattern factor method, with which it would be possible to determine the 
<em>k</em>
 and 
<em>c</em>
 parameters from the power density and average velocity. Rocha et al. [
<a href="#CR28"><sup>28</sup></a>
] dealt with the analysis and comparison of seven numerical methods for the assessment of effectiveness in determining the parameters for the Weibull distribution, using wind data collected for Camocim and Paracuru cities in the northeast region of Brazil. Also in the Brazilian northeastern region, Andrade et al [
<a href="#CR3"><sup>3</sup></a>
] compared the graphical method, moment, pattern energy, maximum likelihood, empirical and equivalent energy and evaluated the efficiency through the predicted and measured power available.</p>
<p>However, in some cases, these methods cannot represent satisfactorily the wind speed distribution. Therefore, a favorable condition for the study of the heuristic method applications has been applied in more recent studies in the field of wind energy. Rahmani et al. [
<a href="#CR27"><sup>27</sup></a>
] estimated, applying the Particle Swarm Optimization, the wind speed and the power produced in the Binaloud Wind Farm. Barbosa [
<a href="#CR6"><sup>6</sup></a>
] estimated the Weibull curve parameters through the Harmony Search for two Brazilian regions. Wang et al. [
<a href="#CR36"><sup>36</sup></a>
] used the Cuckoo Search Optimization and Ant Colony Optimization methods to evaluate wind potential and predict wind speed in four locations in China. Gonzlez et al. [
<a href="#CR17"><sup>17</sup></a>
] presented a new approach for optimizing the layout of offshore wind farms comparing the behavior of two metaheuristic optimization algorithms, the genetic algorithm and Particle Swarm Optimization. Hajibandeh et al. [
<a href="#CR18"><sup>18</sup></a>
] used the multicriteria multi-objective heuristic method to propose a new model for wind energy and DR integration, optimizing supply and demand side operations by the time to use (TOU) or incentive with the emergency DR program (EDRP), as well as combining TOU and EDRP together. Salcedo-Sanz et al. [
<a href="#CR29"><sup>29</sup></a>
] addressed a problem of representative selection of measurement points for long-term wind energy analysis, as the objective of selecting the best set of 
<em>N</em>
 measurement points, such that a measure of wind energy error reconstruction is minimized considering a monthly average wind energy field, for which the metaheuristic algorithm, Coral Reef Optimization with Substrate Layer, was used, which is an evolutionary type method capable of combining different search procedures within a single population. Faced with the inconsistent relationship between China’s economy and the distribution of wind power potential that caused unavoidable difficulties in wind power transport and even network integration, Jiang et al. [
<a href="#CR19"><sup>19</sup></a>
] studied, by optimization methods, among them the Cuckoo Search and the Particle Swarm, the establishment of an integrated electric energy system with low-speed wind energy. Marzband et al. [
<a href="#CR23"><sup>23</sup></a>
] used four heuristically optimized optimization algorithms to implement a market structure based on transactional energy, to ensure that market participants can obtain a higher return.</p>
<p>Considering the presented works, this study aims to analyze four heuristic optimization methods and compare them with four other deterministic numerical methods, to suggest which is the most efficient to determine the parameters of the Weibull probability distribution curve for Petrolina, Triunfo and São Martinho da Serra regions.</p></section>
<section><h2>Weibull distribution</h2>
<p>Wind speed is a random variable, and it is useful to use statistical analysis to determine the wind potential of a region [
<a href="#CR2"><sup>2</sup></a>
, 
<a href="#CR9"><sup>9</sup></a>
, 
<a href="#CR35"><sup>35</sup></a>
]. Commonly, the two parameters Weibull distribution is the one that presents the best fit and is therefore the most used to estimate this potential. [
<a href="#CR8"><sup>8</sup></a>
, 
<a href="#CR22"><sup>22</sup></a>
].</p>
<p>The Weibull distribution for the velocity 
<em>v</em>
 is expressed by the probability density function, wind velocity frequency curve, shown in Eq. 
<a href="#Equ1"><sup>1</sup></a>
. Equation 
<a href="#Equ2"><sup>2</sup></a>
 expresses its cumulative probability function [
<a href="#CR10"><sup>10</sup></a>
, 
<a href="#CR24"><sup>24</sup></a>
].
<section id="Equ1"><mml:math display="block" id="Equ1_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>v</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:mfenced close=")" open="(" separators=""><mml:mfrac><mml:mi>k</mml:mi><mml:mi>c</mml:mi></mml:mfrac></mml:mfenced><mml:mo>·</mml:mo><mml:msup><mml:mfenced close=")" open="(" separators=""><mml:mfrac><mml:mi>v</mml:mi><mml:mi>c</mml:mi></mml:mfrac></mml:mfenced><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msup><mml:mo>·</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:msup><mml:mfenced close=")" open="(" separators=""><mml:mfrac><mml:mi>v</mml:mi><mml:mi>c</mml:mi></mml:mfrac></mml:mfenced><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ1_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} f(v)= & {} \left( \frac{k}{c}\right) \cdot \left( \frac{v}{c}\right) ^{(k-1)}\cdot e^{-\left( \frac{v}{c}\right) ^{k}}, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ1.gif"/></section>
<section id="Equ2"><mml:math display="block" id="Equ2_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>F</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>v</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:msubsup><mml:mo>∫</mml:mo><mml:mn>0</mml:mn><mml:mi>v</mml:mi></mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>v</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>v</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mfrac><mml:mi>v</mml:mi><mml:mi>c</mml:mi></mml:mfrac><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:msup><mml:mspace width="1em"/><mml:mi>v</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.166667em"/><mml:mi>k</mml:mi><mml:mspace width="0.166667em"/><mml:mi mathvariant="normal">and</mml:mi><mml:mspace width="0.166667em"/><mml:mi>c</mml:mi><mml:mo>&gt;</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ2_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} F(v)= & {} \int _0^v f(v) \mathrm {d}v = 1 - e^{-(\frac{v}{c})^{k}} \quad v, \, k \, \mathrm{and} \, c > 0, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ2.gif"/></section>
where 
<em>c</em>
 is the scaling factor with unit 
<span id="IEq2"><mml:math id="IEq2_Math"><mml:mrow><mml:mi>m</mml:mi><mml:mo>·</mml:mo><mml:msup><mml:mi>s</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math><tex-math id="IEq2_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$m \cdot s^{-1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq2.gif"/></span>
, 
<em>k</em>
 is the shape factor (dimensionless) and 
<em>F</em>
(
<em>v</em>
) denotes the probability of velocities smaller than or equal to 
<em>v</em>
.</p>
<p>Among the methods already studied with the purpose of Weibull curves estimating parameters for the regions studied in this paper, or with similar characteristics, the maximum likelihood method, moment method, empirical method and the equivalent energy method were shown to be the most effective [
<a href="#CR3"><sup>3</sup></a>
, 
<a href="#CR5"><sup>5</sup></a>
, 
<a href="#CR28"><sup>28</sup></a>
].</p>
<section><h2>Maximum likelihood method (MLM)</h2>
<p>In the maximum likelihood method, numerical iterations are required to determine the Weibull distribution parameters [
<a href="#CR15"><sup>15</sup></a>
]. In this method [
<a href="#CR28"><sup>28</sup></a>
], the parameters 
<em>k</em>
 and 
<em>c</em>
 are determined according to the Eqs. 
<a href="#Equ3"><sup>3</sup></a>
 and 
<a href="#Equ4"><sup>4</sup></a>
.
<section id="Equ3"><mml:math display="block" id="Equ3_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:msup><mml:mfenced close="]" open="[" separators=""><mml:mfrac><mml:mrow><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msubsup><mml:mi>v</mml:mi><mml:mi>i</mml:mi><mml:mi>k</mml:mi></mml:msubsup><mml:mi mathvariant="normal">ln</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msubsup><mml:mi>v</mml:mi><mml:mi>i</mml:mi><mml:mi>k</mml:mi></mml:msubsup></mml:mrow></mml:mfrac><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mi mathvariant="normal">ln</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mi>n</mml:mi></mml:mfrac></mml:mfenced><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ3_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} k= & {} \left[ \frac{\sum _{i=1}^{n}{v_i^k}\mathrm{ln}{(v_i)}}{\sum _{i=1}^{n}{v_i^k}}- \frac{\sum _{i=1}^{n}\mathrm{ln}{(v_i)}}{n}\right] ^{-1}, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ3.gif"/></section>
<section id="Equ4"><mml:math display="block" id="Equ4_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>c</mml:mi><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:msup><mml:mfenced close=")" open="(" separators=""><mml:mfrac><mml:mn>1</mml:mn><mml:mi>n</mml:mi></mml:mfrac><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:msubsup><mml:mi>v</mml:mi><mml:mi>i</mml:mi><mml:mi>k</mml:mi></mml:msubsup></mml:mfenced><mml:mfrac><mml:mn>1</mml:mn><mml:mi>k</mml:mi></mml:mfrac></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ4_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} c= & {} \left( \frac{1}{n}\sum _{i=1}^{n}{v_i^k}\right) ^{\frac{1}{k}}, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ4.gif"/></section>
where 
<em>n</em>
 is the number of observed data and 
<span id="IEq3"><mml:math id="IEq3_Math"><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:math><tex-math id="IEq3_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$v_{i}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq3.gif"/></span>
 is the wind speed measured in the interval 
<em>i</em>
.</p></section>
<section><h2>Moment method (MM)</h2>
<p>The moment method may be used as an alternative to the maximum likelihood method and it is recommended when the mean and standard deviation of the elements are known and are initially on an appropriate scale [
<a href="#CR21"><sup>21</sup></a>
]. In this case [
<a href="#CR28"><sup>28</sup></a>
], the 
<em>k</em>
 and 
<em>c</em>
 parameters are determined by Eqs. 
<a href="#Equ5"><sup>5</sup></a>
 and 
<a href="#Equ6"><sup>6</sup></a>
.
<section id="Equ5"><mml:math display="block" id="Equ5_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>σ</mml:mi><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:mi>c</mml:mi><mml:mo>·</mml:mo><mml:msqrt><mml:mrow><mml:mi>Γ</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>2</mml:mn><mml:mi>k</mml:mi></mml:mfrac></mml:mfenced><mml:mo>-</mml:mo><mml:msup><mml:mi>Γ</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mfenced close=")" open="(" separators=""><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>2</mml:mn><mml:mi>k</mml:mi></mml:mfrac></mml:mfenced></mml:mrow></mml:msqrt><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ5_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \sigma= & {} c\cdot \sqrt{\varGamma \left( 1+\frac{2}{k}\right) -\varGamma ^2\left( 1+\frac{2}{k}\right) }, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ5.gif"/></section>
<section id="Equ6"><mml:math display="block" id="Equ6_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mover><mml:mi>v</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:mi>c</mml:mi><mml:mo>·</mml:mo><mml:mi>Γ</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>k</mml:mi></mml:mfrac></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ6_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} {\overline{v}}= & {} c\cdot \varGamma \left( 1+\frac{1}{k}\right) , \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ6.gif"/></section>
where 
<span id="IEq4"><mml:math id="IEq4_Math"><mml:mover accent="true"><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover></mml:math><tex-math id="IEq4_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\bar{v}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq4.gif"/></span>
, 
<span id="IEq5"><mml:math id="IEq5_Math"><mml:mi>σ</mml:mi></mml:math><tex-math id="IEq5_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\sigma$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq5.gif"/></span>
, 
<span id="IEq6"><mml:math id="IEq6_Math"><mml:mi>Γ</mml:mi></mml:math><tex-math id="IEq6_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\varGamma$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq6.gif"/></span>
 are, respectively, the average wind speed, the standard deviation of the observed wind speed data, and the gamma function.</p></section>
<section><h2>Empirical method (EM)</h2>
<p>The empirical method [
<a href="#CR10"><sup>10</sup></a>
, 
<a href="#CR28"><sup>28</sup></a>
] is considered a simplified form of the moment method, in which the determination of the 
<em>k</em>
 parameter follows Eq. 
<a href="#Equ7"><sup>7</sup></a>
 and the 
<em>c</em>
 parameter Eq. 
<a href="#Equ8"><sup>8</sup></a>
.
<section id="Equ7"><mml:math display="block" id="Equ7_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:msup><mml:mfenced close=")" open="(" separators=""><mml:mfrac><mml:mi>σ</mml:mi><mml:mover><mml:mi>v</mml:mi><mml:mo>¯</mml:mo></mml:mover></mml:mfrac></mml:mfenced><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>086</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ7_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} k= & {} \left( \frac{\sigma }{{\overline{v}}}\right) ^{-1,086} \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ7.gif"/></section>
<section id="Equ8"><mml:math display="block" id="Equ8_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mover><mml:mi>v</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mrow/><mml:mi>c</mml:mi><mml:mo>·</mml:mo><mml:mi>Γ</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>k</mml:mi></mml:mfrac><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ8_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} {\overline{v}}= {} c\cdot \varGamma (1+\frac{1}{k}), \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ8.gif"/></section>
where 
<span id="IEq7"><mml:math id="IEq7_Math"><mml:mover accent="true"><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover></mml:math><tex-math id="IEq7_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\bar{v}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq7.gif"/></span>
 and 
<span id="IEq8"><mml:math id="IEq8_Math"><mml:mi>σ</mml:mi></mml:math><tex-math id="IEq8_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\sigma$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq8.gif"/></span>
 are, respectively, the mean wind speed and the standard deviation of the observed wind speed data.</p></section>
<section><h2>Equivalent energy method (EEM)</h2>
<p>The equivalent energy method seeks the equivalence between the energy density of the observations and the theoretical Weibull curve. For this, the 
<em>k</em>
 parameter is estimated from the third moment of the velocity, by minimizing the square error related to the adjustment, represented by Eq. 
<a href="#Equ9"><sup>9</sup></a>
 and the 
<em>c</em>
 parameter is adjusted by using Eq. 
<a href="#Equ10"><sup>10</sup></a>
 [
<a href="#CR3"><sup>3</sup></a>
, 
<a href="#CR32"><sup>32</sup></a>
].
<section id="Equ9"><mml:math display="block" id="Equ9_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msup><mml:mi>ϵ</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:msup><mml:mfenced close="}" open="{" separators=""><mml:msub><mml:mi>W</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:msup><mml:mfenced close="]" open="[" separators=""><mml:mfrac><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>Γ</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>3</mml:mn><mml:mi>k</mml:mi></mml:mfrac><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo stretchy="false">/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mover><mml:msup><mml:mi>v</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo stretchy="false">/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:mfrac></mml:mfenced><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:msup><mml:mfenced close="]" open="[" separators=""><mml:mfrac><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>Γ</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>3</mml:mn><mml:mi>k</mml:mi></mml:mfrac><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo stretchy="false">/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mover><mml:msup><mml:mi>v</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo stretchy="false">/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:mfrac></mml:mfenced><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:msup></mml:mfenced><mml:mn>2</mml:mn></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ9_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \epsilon ^2= & {} \sum \limits _{i=1}^{n}\left\{ {W_i}-e^{- \left[ \frac{(v_i- 1)(\varGamma (1+\frac{3}{k}))^{1/3}}{(\overline{v^3})^{1/3}}\right] ^{k}}+e^{-\left[ \frac{(v_i)(\varGamma (1+\frac{3}{k}))^{1/3}}{(\overline{v^3})^{1/3}}\right] ^{k}}\right\} ^2, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ9.gif"/></section>
<section id="Equ10"><mml:math display="block" id="Equ10_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>c</mml:mi><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:msup><mml:mfenced close="]" open="[" separators=""><mml:mfrac><mml:mover><mml:msup><mml:mi>v</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>Γ</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>3</mml:mn><mml:mi>k</mml:mi></mml:mfrac></mml:mfenced></mml:mrow></mml:mfrac></mml:mfenced><mml:mrow><mml:mn>1</mml:mn><mml:mo stretchy="false">/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ10_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} c= & {} \left[ \frac{\overline{v^3}}{\varGamma \left( 1+\frac{3}{k}\right) }\right] ^{1/3}. \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ10.gif"/></section></p></section></section>
<section><h2>Heuristic methods</h2>
<p>Heuristics encompasses a set of methods where, to solve a problem, the variables in question use the experience gained over the iterations. Heuristic methods combine different concepts intelligently to explore the search space, so that learning strategies are used to structure information and find efficient and almost optimal solutions [
<a href="#CR25"><sup>25</sup></a>
]. Many of the heuristic approaches depend on probabilistic decisions made during the algorithm run. The main difference against pure random search is that in heuristic algorithms, randomness is not used blindly but intelligently and biased [
<a href="#CR34"><sup>34</sup></a>
]. It is valid to emphasize that every optimization procedure searches for the best result of a function for the desired scenario. This function is called the objective function. In this paper, the objective function is the one presented in Eq. 
<a href="#Equ11"><sup>11</sup></a>
, which represents the minimization of the square error sum applied to the frequency of occurrence values found by the curve adjusted by the method and the observed frequency of occurrence in the histogram of the data.
<section id="Equ11"><mml:math display="block" id="Equ11_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msup><mml:mi>ϵ</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi mathvariant="normal">adjustment</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi mathvariant="normal">observed</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ11_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \epsilon ^2 = \sum \limits _{i=1}^{n} (f_{\mathrm{adjustment}} - f_{\mathrm{observed}})^2, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ11.gif"/></section>
where 
<em>n</em>
 is the number of histogram velocity intervals and 
<span id="IEq9"><mml:math id="IEq9_Math"><mml:msub><mml:mi>f</mml:mi><mml:mi mathvariant="normal">adjustment</mml:mi></mml:msub></mml:math><tex-math id="IEq9_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$f_{\mathrm{adjustment}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq9.gif"/></span>
 and 
<span id="IEq10"><mml:math id="IEq10_Math"><mml:msub><mml:mi>f</mml:mi><mml:mi mathvariant="normal">observed</mml:mi></mml:msub></mml:math><tex-math id="IEq10_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$f_{\mathrm{observed}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq10.gif"/></span>
 are the occurrence frequencies from the adjusted curve and that observed in the histogram, respectively.</p>
<section><h2>Harmony search (HS)</h2>
<p>The Harmony Search is a heuristic algorithm based on the analogy of the artificial phenomenon of a musical group in search of the best harmony. This search occurs by the combination of existing elements and the generation of new elements that are combined to form possible solutions [
<a href="#CR16"><sup>16</sup></a>
]. The search process begins with the formation of a harmony memory (HM), by the memorization of a series of possible solutions, denominated harmonies. At each iteration, a new harmony is formed and compared to the harmonies stored in the HM. The algorithm is presented by the following steps [
<a href="#CR4"><sup>4</sup></a>
]:
<ul><li><p>initialize the HARMONY MEMORY;</p></li>
<li><p>improvise a new harmony from HM;</p></li>
<li><p>if the new harmony is better than the minimum harmony in HM, include the new harmony in HM and exclude the minimum harmony from HM. If not, the new harmony is excluded;</p></li>
<li><p>if the stopping criteria are not satisfied, go to step 2.</p></li></ul>
Figure 
<a href="#Fig1"><sup>1</sup></a>
 shows the flowchart of the Harmony Search algorithm, summarizing all steps described previously.
<figure id="Fig1"><h3>Fig. 1</h3>
<figcaption><p>Flowchart of the HS algorithm</p></figcaption>
<img src="40095_2018_285_Fig1_HTML.png" /></figure></p></section>
<section><h2>Cuckoo search optimization (CSO)</h2>
<p>Cuckoos are birds with an aggressive breeding strategy. Some species such as Ani and Guira cuckoos place their eggs in communal nests, and sometimes remove other species’ eggs to increase the hatching probability of their own eggs. Other species lay their eggs in nesting host birds (often of other species). New World brood-parasitic Tapera species have evolved in such a way that the female parasitic cuckoos are often very specialized in the mimicry of the color and pattern of the eggs of a few chosen host species. This ability reduces the probability of their eggs being abandoned and thus increasing their reproductivity [
<a href="#CR26"><sup>26</sup></a>
].</p>
<p>The CSO has its origin inspired by the behavior of the cuckoo in the process of finding nests, in which a nest is a possible solution. First, an initial population of nests is randomly generated. Later, new solutions are generated via Lévy flights. and from these the best solutions are stored in comparison to the current solutions. There are several ways to implement the distribution of Lévy distribution, the simplest is the Mantegna
<span id="IEq11"><mml:math id="IEq11_Math"><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:math><tex-math id="IEq11_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$'$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq11.gif"/></span>
s algorithm [
<a href="#CR37"><sup>37</sup></a>
], and its distribution takes the form presented by Eq. 
<a href="#Equ12"><sup>12</sup></a>
 [
<a href="#CR20"><sup>20</sup></a>
].
<section id="Equ12"><mml:math display="block" id="Equ12_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:msub><mml:mi>α</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>·</mml:mo><mml:mfrac><mml:mrow><mml:mi>ϕ</mml:mi><mml:mo>·</mml:mo><mml:mi>u</mml:mi></mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">|</mml:mo><mml:mi>v</mml:mi><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mi>β</mml:mi></mml:mfrac></mml:msup></mml:mfrac><mml:mo>·</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi mathvariant="normal">best</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ12_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} x_i^{t+1}=x_i^{(t)} + \alpha _0 \cdot \frac{\phi \cdot u}{|v|^\frac{1}{\beta }} \cdot (x_i^{(t)}-x_{\mathrm{best}}^{(t)} ), \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ12.gif"/></section>
where 
<span id="IEq12"><mml:math id="IEq12_Math"><mml:msubsup><mml:mi>X</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup></mml:math><tex-math id="IEq12_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$X_i^{(t)}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq12.gif"/></span>
 is the previous solution from which the new solution 
<span id="IEq13"><mml:math id="IEq13_Math"><mml:msubsup><mml:mi>X</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup></mml:math><tex-math id="IEq13_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$X_i^{(t+1)}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq13.gif"/></span>
 has been generated, 
<span id="IEq14"><mml:math id="IEq14_Math"><mml:msub><mml:mi>α</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:math><tex-math id="IEq14_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\alpha _0$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq14.gif"/></span>
 is a constant, usually 0.01, 
<span id="IEq15"><mml:math id="IEq15_Math"><mml:msubsup><mml:mi>X</mml:mi><mml:mrow><mml:mi mathvariant="normal">best</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup></mml:math><tex-math id="IEq15_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$X_{\mathrm{best}}^{(t)}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq15.gif"/></span>
 represents the best actual solution, 
<em>u</em>
 and 
<em>v</em>
 are drawn from normal distributions, 
<span id="IEq16"><mml:math id="IEq16_Math"><mml:mi>β</mml:mi></mml:math><tex-math id="IEq16_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\beta$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq16.gif"/></span>
 is the scale factor which has an assigned value of 1.5 and 
<span id="IEq17"><mml:math id="IEq17_Math"><mml:mi>ϕ</mml:mi></mml:math><tex-math id="IEq17_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\phi$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq17.gif"/></span>
 is calculated according to Eq. 
<a href="#Equ13"><sup>13</sup></a>
, where 
<span id="IEq18"><mml:math id="IEq18_Math"><mml:mi>Γ</mml:mi></mml:math><tex-math id="IEq18_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\varGamma$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq18.gif"/></span>
 is the gamma function.
<section id="Equ13"><mml:math display="block" id="Equ13_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>ϕ</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mfenced close="}" open="{" separators=""><mml:mfrac><mml:mrow><mml:mi>Γ</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>β</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>·</mml:mo><mml:mo>sin</mml:mo><mml:mfenced close=")" open="(" separators=""><mml:mfrac><mml:mrow><mml:mi>π</mml:mi><mml:mo>·</mml:mo><mml:mi>β</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:mfrac></mml:mfenced></mml:mrow><mml:mrow><mml:mi>Γ</mml:mi><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">[</mml:mo></mml:mrow><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>β</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:mfrac><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">]</mml:mo></mml:mrow><mml:mo>·</mml:mo><mml:mi>β</mml:mi><mml:mo>·</mml:mo><mml:msup><mml:mn>2</mml:mn><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>β</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:mfrac></mml:msup></mml:mrow></mml:mfrac></mml:mfenced><mml:mfrac><mml:mn>1</mml:mn><mml:mi>β</mml:mi></mml:mfrac></mml:msup><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ13_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \phi = \left\{ \frac{\varGamma (1+\beta ) \cdot \sin \left( \frac{\pi \cdot \beta }{2}\right) }{\varGamma \bigg [\frac{(1+ \beta )}{2}\bigg ]\cdot \beta \cdot 2^\frac{(\beta -1)}{2}} \right\} ^\frac{1}{\beta } . \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ13.gif"/></section>
Then, the solution subset is discarded according to the probability of detection 
<span id="IEq19"><mml:math id="IEq19_Math"><mml:msub><mml:mi>P</mml:mi><mml:mi>a</mml:mi></mml:msub></mml:math><tex-math id="IEq19_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$P_a$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq19.gif"/></span>
<span id="IEq20"><mml:math id="IEq20_Math"><mml:mo>∈</mml:mo></mml:math><tex-math id="IEq20_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\in$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq20.gif"/></span>
 [0,1] and new solutions are obtained, according to Eq. 
<a href="#Equ14"><sup>14</sup></a>
, with the same quantity of solutions abandoned [
<a href="#CR20"><sup>20</sup></a>
, 
<a href="#CR30"><sup>30</sup></a>
].
<section id="Equ14"><mml:math display="block" id="Equ14_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:mi>r</mml:mi><mml:mo>·</mml:mo><mml:mfenced close=")" open="(" separators=""><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ14_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} x_i^{t+1}=x_i^{(t)} + r \cdot \left( x_{i,c}^{(t)}-x_{i,k}^{(t)} \right) . \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ14.gif"/></section>
In Eq. 
<a href="#Equ14"><sup>14</sup></a>
, 
<em>r</em>
 is a uniformly distributed random number from 0 to 1, and 
<span id="IEq21"><mml:math id="IEq21_Math"><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msub></mml:math><tex-math id="IEq21_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$X_{(i,c)}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq21.gif"/></span>
 and 
<span id="IEq22"><mml:math id="IEq22_Math"><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msub></mml:math><tex-math id="IEq22_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$X_{(i,k)}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq22.gif"/></span>
 denote the two random solutions of the 
<em>i</em>
th generation.</p>
<p>Figure 
<a href="#Fig2"><sup>2</sup></a>
 shows the flowchart of the Cuckoo Search Optimization algorithm, summarizing all steps described previously.
<figure id="Fig2"><h3>Fig. 2</h3>
<figcaption><p>Flowchart of the CSO algorithm</p></figcaption>
<img src="40095_2018_285_Fig2_HTML.png" /></figure></p></section>
<section><h2>Particle swarm optimization (PSO)</h2>
<p>In a PSO system, each particle “flies” through the multidimensional search space, adjusting its position in space according to its own experience, however, also considering the experience of the neighboring particle. A particle uses the best position found by itself and the best position of its neighbors to position itself toward an ideal solution. The effect is that the particles “fly” toward a global optimum, while still investigating an area around the best current solution [
<a href="#CR14"><sup>14</sup></a>
]. For each particle 
<em>k</em>
 positioned in a two-dimensional plane and for each iteration i, the positions and the best individual results 
<span id="IEq23"><mml:math id="IEq23_Math"><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>k</mml:mi><mml:mi mathvariant="normal">best</mml:mi></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>k</mml:mi><mml:mi mathvariant="normal">best</mml:mi></mml:msubsup><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math><tex-math id="IEq23_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$(x_k^{\mathrm{best}},y_k^{\mathrm{best}})$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq23.gif"/></span>
 are recorded. Then, the best result among the k particles is recorded 
<span id="IEq24"><mml:math id="IEq24_Math"><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi mathvariant="normal">global</mml:mi></mml:mrow><mml:mi mathvariant="normal">best</mml:mi></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mrow><mml:mi mathvariant="normal">global</mml:mi></mml:mrow><mml:mi mathvariant="normal">best</mml:mi></mml:msubsup><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math><tex-math id="IEq24_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$(x_{\mathrm{global}}^{\mathrm{best}},y_{\mathrm{global}}^{\mathrm{best}})$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq24.gif"/></span>
. Each particle’s movement will be proportional to the distance between the current position of the particle and the resulting point of the weighted average between the best individual position of the particle and the best position of the swarm, according to Eqs. 
<a href="#Equ15"><sup>15</sup></a>
 and 
<a href="#Equ16"><sup>16</sup></a>
 [
<a href="#CR13"><sup>13</sup></a>
].
<section id="Equ15"><mml:math display="block" id="Equ15_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msubsup><mml:mi>x</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mrow/><mml:msubsup><mml:mi>x</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ15_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} x_k^{(i+1)}= {} x_k^{(i)}+V_{x,k}^{(i)} , \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ15.gif"/></section>
<section id="Equ16"><mml:math display="block" id="Equ16_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msubsup><mml:mi>y</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mrow/><mml:msubsup><mml:mi>y</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ16_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} y_k^{(i+1)}= {} y_k^{(i)}+V_{y,k}^{(i)} , \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ16.gif"/></section>
where
<section id="Equ17"><mml:math display="block" id="Equ17_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mrow/><mml:msup><mml:mi>ω</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msup><mml:mo>·</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>·</mml:mo><mml:mi>λ</mml:mi><mml:mo>·</mml:mo><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="normal">best</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>·</mml:mo><mml:mi>μ</mml:mi><mml:mo>·</mml:mo><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="normal">global</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="normal">best</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ17_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} V_{x,k}^{(i+1)}= {} \omega ^{(i)}\cdot V_{x,k}^{(i)} +c_1\cdot \lambda \cdot [x_k^{\mathrm{(best)}} - x_k^{(i)}] + c_2\cdot \mu \cdot [x_{\mathrm{(global)}}^{\mathrm{(best)}} - x_k^{(i)}], \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ17.gif"/></section>
<section id="Equ18"><mml:math display="block" id="Equ18_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mrow/><mml:msup><mml:mi>ω</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msup><mml:mo>·</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>·</mml:mo><mml:mi>η</mml:mi><mml:mo>·</mml:mo><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="normal">best</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>·</mml:mo><mml:mi>ϵ</mml:mi><mml:mo>·</mml:mo><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="normal">global</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="normal">best</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msubsup><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ18_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} V_{y,k}^{(i+1)}= {} \omega ^{(i)}\cdot V_{y,k}^{(i)} +c_1\cdot \eta \cdot [y_k^{\mathrm{(best)}} - y_k^{(i)}] + c_2\cdot \epsilon \cdot [x_{\mathrm{(global)}}^{\mathrm{(best)}} - y_k^{(i)}]. \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ18.gif"/></section>
In these, 
<span id="IEq25"><mml:math id="IEq25_Math"><mml:mi>λ</mml:mi></mml:math><tex-math id="IEq25_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\lambda$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq25.gif"/></span>
, 
<span id="IEq26"><mml:math id="IEq26_Math"><mml:mi>μ</mml:mi></mml:math><tex-math id="IEq26_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\mu$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq26.gif"/></span>
, 
<span id="IEq27"><mml:math id="IEq27_Math"><mml:mi>η</mml:mi></mml:math><tex-math id="IEq27_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\eta$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq27.gif"/></span>
 and 
<span id="IEq28"><mml:math id="IEq28_Math"><mml:mi>ϵ</mml:mi></mml:math><tex-math id="IEq28_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\epsilon$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq28.gif"/></span>
 are random numbers belonging to the set [0, 1] and 
<span id="IEq29"><mml:math id="IEq29_Math"><mml:mi>ω</mml:mi></mml:math><tex-math id="IEq29_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\omega$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq29.gif"/></span>
 is the particle inertia term, defined by Eq. 
<a href="#Equ19"><sup>19</sup></a>
 [
<a href="#CR31"><sup>31</sup></a>
].
<section id="Equ19"><mml:math display="block" id="Equ19_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msup><mml:mi>ω</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msub><mml:mi>ω</mml:mi><mml:mi mathvariant="normal">initial</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>ω</mml:mi><mml:mi mathvariant="normal">final</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>ω</mml:mi><mml:mi mathvariant="normal">initial</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>·</mml:mo><mml:mrow><mml:mo maxsize="1.623em" minsize="1.623em" stretchy="true">(</mml:mo></mml:mrow><mml:mfrac><mml:mi>i</mml:mi><mml:mi>m</mml:mi></mml:mfrac><mml:mrow><mml:mo maxsize="1.623em" minsize="1.623em" stretchy="true">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ19_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \omega ^{(i)}=\omega _{\mathrm{initial}}+(\omega _{\mathrm{final}}-\omega _{\mathrm{initial}})\cdot \Big (\frac{i}{m}\Big ), \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ19.gif"/></section>
where 
<em>m</em>
 is the maximum number of iterations.</p>
<p>Figure 
<a href="#Fig3"><sup>3</sup></a>
 shows the flowchart of the Particle Search Optimization algorithm, summarizing all steps described previously.
<figure id="Fig3"><h3>Fig. 3</h3>
<figcaption><p>Flowchart of the PSO algorithm</p></figcaption>
<img src="40095_2018_285_Fig3_HTML.png" /></figure></p></section>
<section><h2>Ant colony optimization (ACO)</h2>
<p>In an ant colony, the communication between individuals, or between the individuals and the environment, is based on the pheromone produced by them. The trail pheromone is a specific type of pheromone that some ant species use to mark paths on the ground. When detecting pheromone trails, forage ants may follow the path trodden by other ants to the food source. The first ants when sniffing the pheromone tend to choose, probabilistically, the trails marked with stronger concentrations of pheromone. The second group of ants will notice more intense the shortest path due to the shorter evaporation time. With the continuation of this procedure by all the ants, at one point in this process, one of the paths stands out for being the most frequented, being indicated by the intensity of ants’ pheromone and density superior to the others. At this point, the best path found by the ants is defined. This behavior inspired the optimization method by ant colonies [
<a href="#CR11"><sup>11</sup></a>
].</p>
<p>In the ACO method, the parameters 
<em>k</em>
 and 
<em>c</em>
 of the Weibull curve form a Cartesian plane that is divided into 
<em>N</em>
 equal parts. The center point of each new area will be an ordered pair (
<em>k</em>
, 
<em>c</em>
) where the curve fit will be evaluated [
<a href="#CR5"><sup>5</sup></a>
]. The probability of occurrence of each reticulum is defined by Eq. 
<a href="#Equ20"><sup>20</sup></a>
.
<section id="Equ20"><mml:math display="block" id="Equ20_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>r</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:msub><mml:mi>τ</mml:mi><mml:mi>r</mml:mi></mml:msub><mml:mrow><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>R</mml:mi></mml:msubsup><mml:msub><mml:mi>τ</mml:mi><mml:mi>r</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ20_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} P_{r}= \frac{\tau _r}{\sum _{r=1}^{R} \tau _r}, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ20.gif"/></section>
where 
<span id="IEq30"><mml:math id="IEq30_Math"><mml:msub><mml:mi>τ</mml:mi><mml:mi>r</mml:mi></mml:msub></mml:math><tex-math id="IEq30_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tau _r$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq30.gif"/></span>
 is the pheromone intensity for the reticle 
<em>r</em>
.</p>
<p>Each ant is then randomly positioned in the plane through a roulette draw, where each slice of the roulette represents a reticle and is defined by the probability of occurrence. The visited quadrants are indicated by the pheromone deposit according to Eq. 
<a href="#Equ21"><sup>21</sup></a>
. At each iteration, the amount of the hormone is also reduced at a constant rate to simulate the hormone volatility, according to Eq. 
<a href="#Equ22"><sup>22</sup></a>
.
<section id="Equ21"><mml:math display="block" id="Equ21_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msub><mml:mi>τ</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>r</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:msub><mml:mi>τ</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi>r</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">(</mml:mo></mml:mrow><mml:mfrac><mml:mi>μ</mml:mi><mml:msub><mml:mi mathvariant="normal">err</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:mfrac><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">)</mml:mo></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ21_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \tau _{i, r}= & {} \tau _{i-1, r} + \bigg (\frac{\mu }{{\mathrm{err}}_f} \bigg ) \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ21.gif"/></section>
<section id="Equ22"><mml:math display="block" id="Equ22_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msub><mml:mi>τ</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>r</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:mi>ρ</mml:mi><mml:mo>·</mml:mo><mml:msub><mml:mi>τ</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>r</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ22_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \tau _{i, r}= & {} \rho \cdot \tau _{i,r}, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ22.gif"/></section>
where 
<span id="IEq31"><mml:math id="IEq31_Math"><mml:msub><mml:mi>τ</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>r</mml:mi></mml:mrow></mml:msub></mml:math><tex-math id="IEq31_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tau _{i, r}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq31.gif"/></span>
 is the pheromone intensity for the reticle 
<em>r</em>
, at iteration 
<em>i</em>
, 
<span id="IEq32"><mml:math id="IEq32_Math"><mml:msub><mml:mi mathvariant="normal">err</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:math><tex-math id="IEq32_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\mathrm{err}}_f$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq32.gif"/></span>
 is the error evaluated by the ant 
<em>f</em>
, 
<span id="IEq33"><mml:math id="IEq33_Math"><mml:mi>μ</mml:mi></mml:math><tex-math id="IEq33_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\mu$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq33.gif"/></span>
 is the deposition constant and 
<span id="IEq34"><mml:math id="IEq34_Math"><mml:mi>ρ</mml:mi></mml:math><tex-math id="IEq34_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\rho$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq34.gif"/></span>
 is the evaporation constant [
<a href="#CR33"><sup>33</sup></a>
].</p>
<p>While the iterations follow up, some reticles will be more attractive to ants because they have a large amount of pheromone, this attraction being symbolized by the larger slices of the roulette, until most of the ants follow the same path.</p>
<p>Figure 
<a href="#Fig4"><sup>4</sup></a>
 shows the flowchart of the Ant Colony Optimization algorithm, summarizing all steps described previously.
<figure id="Fig4"><h3>Fig. 4</h3>
<figcaption><p>Flowchart of the ACO algorithm</p></figcaption>
<img src="40095_2018_285_Fig4_HTML.png" /></figure></p></section>
<section><h2>Parameters applied to the heuristic methods</h2>
<p>Each heuristic method depends on a certain number of parameters, with its adjustment being necessary to reduce the computational time response that leads to convergence to the optimal values. The parameters applied here, presented in Table 
<a href="#Tab1"><sup>1</sup></a>
, were extracted from the works whose authors used the proposed methods in wind energy applications [
<a href="#CR6"><sup>6</sup></a>
, 
<a href="#CR7"><sup>7</sup></a>
, 
<a href="#CR27"><sup>27</sup></a>
, 
<a href="#CR36"><sup>36</sup></a>
].
<figure id="Tab1"><h3>Table 1</h3>
<figcaption><p>Parameters applied to the heuristic methods</p></figcaption>
<table frame="hsides" rules="groups"><thead><tr><th align="left" colspan="2"><p>HS</p></th><th align="left" colspan="2"><p>CSO</p></th><th align="left" colspan="2"><p>PSO</p></th><th align="left" colspan="2"><p>ACO</p></th></tr></thead><tbody><tr><td align="left"><p><inline-formula id="IEq35"><alternatives><mml:math id="IEq35_Math"><mml:msub><mml:mi>N</mml:mi><mml:mi>h</mml:mi></mml:msub></mml:math><tex-math id="IEq35_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$N_h$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq35.gif"/></alternatives></inline-formula><inline-formula id="IEq36"><alternatives><mml:math id="IEq36_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">a</mml:mi></mml:msup></mml:math><tex-math id="IEq36_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{a}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq36.gif"/></alternatives></inline-formula></p></td><td align="left"><p>6</p></td><td align="left"><p><inline-formula id="IEq37"><alternatives><mml:math id="IEq37_Math"><mml:msub><mml:mi>N</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:math><tex-math id="IEq37_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$N_n$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq37.gif"/></alternatives></inline-formula><inline-formula id="IEq38"><alternatives><mml:math id="IEq38_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">b</mml:mi></mml:msup></mml:math><tex-math id="IEq38_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{b}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq38.gif"/></alternatives></inline-formula></p></td><td align="left"><p>50</p></td><td align="left"><p><inline-formula id="IEq39"><alternatives><mml:math id="IEq39_Math"><mml:msub><mml:mi>N</mml:mi><mml:mi>p</mml:mi></mml:msub></mml:math><tex-math id="IEq39_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$N_p$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq39.gif"/></alternatives></inline-formula><inline-formula id="IEq40"><alternatives><mml:math id="IEq40_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">d</mml:mi></mml:msup></mml:math><tex-math id="IEq40_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{d}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq40.gif"/></alternatives></inline-formula></p></td><td align="left"><p>30</p></td><td align="left"><p><inline-formula id="IEq41"><alternatives><mml:math id="IEq41_Math"><mml:msub><mml:mi>N</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:math><tex-math id="IEq41_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$N_f$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq41.gif"/></alternatives></inline-formula><inline-formula id="IEq42"><alternatives><mml:math id="IEq42_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">i</mml:mi></mml:msup></mml:math><tex-math id="IEq42_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{i}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq42.gif"/></alternatives></inline-formula></p></td><td align="left"><p>100</p></td></tr><tr><td align="left"/><td align="left"/><td align="left"><p><inline-formula id="IEq43"><alternatives><mml:math id="IEq43_Math"><mml:msub><mml:mi>P</mml:mi><mml:mi>a</mml:mi></mml:msub></mml:math><tex-math id="IEq43_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$P_a$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq43.gif"/></alternatives></inline-formula><inline-formula id="IEq44"><alternatives><mml:math id="IEq44_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">c</mml:mi></mml:msup></mml:math><tex-math id="IEq44_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{c}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq44.gif"/></alternatives></inline-formula></p></td><td align="left"><p>0.25</p></td><td align="left"><p><inline-formula id="IEq45"><alternatives><mml:math id="IEq45_Math"><mml:msub><mml:mi>w</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:math><tex-math id="IEq45_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$w_i$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq45.gif"/></alternatives></inline-formula><inline-formula id="IEq46"><alternatives><mml:math id="IEq46_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">e</mml:mi></mml:msup></mml:math><tex-math id="IEq46_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{e}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq46.gif"/></alternatives></inline-formula></p></td><td align="left"><p>1.8</p></td><td align="left"><p><inline-formula id="IEq47"><alternatives><mml:math id="IEq47_Math"><mml:mi>μ</mml:mi></mml:math><tex-math id="IEq47_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\mu$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq47.gif"/></alternatives></inline-formula><inline-formula id="IEq48"><alternatives><mml:math id="IEq48_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">j</mml:mi></mml:msup></mml:math><tex-math id="IEq48_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{j}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq48.gif"/></alternatives></inline-formula></p></td><td align="left"><p>0.2</p></td></tr><tr><td align="left"/><td align="left"/><td align="left"/><td align="left"/><td align="left"><p><inline-formula id="IEq49"><alternatives><mml:math id="IEq49_Math"><mml:msub><mml:mi>w</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:math><tex-math id="IEq49_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$w_f$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq49.gif"/></alternatives></inline-formula><inline-formula id="IEq50"><alternatives><mml:math id="IEq50_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">f</mml:mi></mml:msup></mml:math><tex-math id="IEq50_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{f}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq50.gif"/></alternatives></inline-formula></p></td><td align="left"><p>0.2</p></td><td align="left"><p><inline-formula id="IEq51"><alternatives><mml:math id="IEq51_Math"><mml:mi>ρ</mml:mi></mml:math><tex-math id="IEq51_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\rho$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq51.gif"/></alternatives></inline-formula><inline-formula id="IEq52"><alternatives><mml:math id="IEq52_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">k</mml:mi></mml:msup></mml:math><tex-math id="IEq52_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{k}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq52.gif"/></alternatives></inline-formula></p></td><td align="left"><p>0.1</p></td></tr><tr><td align="left"/><td align="left"/><td align="left"/><td align="left"/><td align="left"><p><inline-formula id="IEq53"><alternatives><mml:math id="IEq53_Math"><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:math><tex-math id="IEq53_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$c_1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq53.gif"/></alternatives></inline-formula><inline-formula id="IEq54"><alternatives><mml:math id="IEq54_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">g</mml:mi></mml:msup></mml:math><tex-math id="IEq54_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{g}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq54.gif"/></alternatives></inline-formula></p></td><td align="left"><p>1.0</p></td><td align="left"/><td align="left"/></tr><tr><td align="left"/><td align="left"/><td align="left"/><td align="left"/><td align="left"><p><inline-formula id="IEq55"><alternatives><mml:math id="IEq55_Math"><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:math><tex-math id="IEq55_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$c_2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq55.gif"/></alternatives></inline-formula><inline-formula id="IEq56"><alternatives><mml:math id="IEq56_Math"><mml:msup><mml:mrow/><mml:mi mathvariant="normal">h</mml:mi></mml:msup></mml:math><tex-math id="IEq56_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{\mathrm{h}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq56.gif"/></alternatives></inline-formula></p></td><td align="left"><p>1.0</p></td><td align="left"/><td align="left"/></tr></tbody></table>
<footer><p><span id="IEq57"><mml:math id="IEq57_Math"><mml:msup><mml:mrow/><mml:mi>a</mml:mi></mml:msup></mml:math><tex-math id="IEq57_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{a}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq57.gif"/></span>
 Harmony number</p>
<p><span id="IEq58"><mml:math id="IEq58_Math"><mml:msup><mml:mrow/><mml:mi>b</mml:mi></mml:msup></mml:math><tex-math id="IEq58_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{b}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq58.gif"/></span>
 Nest number</p>
<p><span id="IEq59"><mml:math id="IEq59_Math"><mml:msup><mml:mrow/><mml:mi>c</mml:mi></mml:msup></mml:math><tex-math id="IEq59_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{c}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq59.gif"/></span>
 Probability of detection</p>
<p><span id="IEq60"><mml:math id="IEq60_Math"><mml:msup><mml:mrow/><mml:mi>d</mml:mi></mml:msup></mml:math><tex-math id="IEq60_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{d}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq60.gif"/></span>
 Particle number</p>
<p><span id="IEq61"><mml:math id="IEq61_Math"><mml:msup><mml:mrow/><mml:mi>e</mml:mi></mml:msup></mml:math><tex-math id="IEq61_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{e}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq61.gif"/></span>
 Inicial particle inertia</p>
<p><span id="IEq62"><mml:math id="IEq62_Math"><mml:msup><mml:mrow/><mml:mi>f</mml:mi></mml:msup></mml:math><tex-math id="IEq62_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{f}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq62.gif"/></span>
 Final particle inertia</p>
<p><span id="IEq63"><mml:math id="IEq63_Math"><mml:msup><mml:mrow/><mml:mi>g</mml:mi></mml:msup></mml:math><tex-math id="IEq63_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{g}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq63.gif"/></span>
 Local learning factor</p>
<p><span id="IEq64"><mml:math id="IEq64_Math"><mml:msup><mml:mrow/><mml:mi>h</mml:mi></mml:msup></mml:math><tex-math id="IEq64_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{h}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq64.gif"/></span>
 Global learning factor</p>
<p><span id="IEq65"><mml:math id="IEq65_Math"><mml:msup><mml:mrow/><mml:mi>i</mml:mi></mml:msup></mml:math><tex-math id="IEq65_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{i}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq65.gif"/></span>
 Ant number</p>
<p><span id="IEq66"><mml:math id="IEq66_Math"><mml:msup><mml:mrow/><mml:mi>j</mml:mi></mml:msup></mml:math><tex-math id="IEq66_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{j}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq66.gif"/></span>
 Pheromone deposit</p>
<p><span id="IEq67"><mml:math id="IEq67_Math"><mml:msup><mml:mrow/><mml:mi>k</mml:mi></mml:msup></mml:math><tex-math id="IEq67_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^{k}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq67.gif"/></span>
 Evaporation constant</p></footer></figure></p></section></section>
<section><h2>Statistical tests</h2>
<p>The performance evaluation of the applied methods was realized by the following tests:</p>
<p>Root mean square error (RMSE) (Eq, 
<a href="#Equ23"><sup>23</sup></a>
):
<section id="Equ23"><mml:math display="block" id="Equ23_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi mathvariant="normal">RMSE</mml:mi><mml:mspace width="0.166667em"/><mml:mo>=</mml:mo><mml:msqrt><mml:mfrac><mml:mrow><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mtext>calculated</mml:mtext></mml:msubsup><mml:mo>-</mml:mo><mml:msup><mml:mi>y</mml:mi><mml:mtext>measured</mml:mtext></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mi>n</mml:mi></mml:mfrac></mml:msqrt><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ23_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} \mathrm{RMSE}\, = \sqrt{\frac{\sum _{i=1}^{n}(y_i^{\text {calculated}} - y^{\text {measured}} )^2}{n}}. \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ23.gif"/></section>
Mean absolute error (MAE) (Eq. 
<a href="#Equ24"><sup>24</sup></a>
):
<section id="Equ24"><mml:math display="block" id="Equ24_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>M</mml:mi><mml:mi>A</mml:mi><mml:mi>E</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>n</mml:mi></mml:mfrac><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:mfenced close="|" open="|" separators=""><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mtext>calculated</mml:mtext></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mtext>measured</mml:mtext></mml:msubsup></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ24_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} MAE = \frac{1}{n}\displaystyle \sum _{i=1}^{n} {\left| y_i^{\text {calculated}} - y_i^{\text {measured}}\right| }. \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ24.gif"/></section>
Determination coefficient 
<em>R</em>
<span id="IEq68"><mml:math id="IEq68_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq68_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq68.gif"/></span>
 (Eq. 
<a href="#Equ25"><sup>25</sup></a>
):
<section id="Equ25"><mml:math display="block" id="Equ25_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mtext>measured</mml:mtext></mml:msubsup><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover></mml:mrow><mml:mtext>measured</mml:mtext></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mo>-</mml:mo><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mtext>measured</mml:mtext></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mtext>calculated</mml:mtext></mml:msubsup><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mtext>measured</mml:mtext></mml:msubsup><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover></mml:mrow><mml:mtext>measured</mml:mtext></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ25_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} R^2 = \frac{\sum _{i=1}^{n}(y_i^{\text {measured}} - {\bar{y}}^{\text {measured}} )^2 - \sum _{i=1}^{n}(y_i^{\text {measured}} - y_i^{\text {calculated}})^2 }{\sum _{i=1}^{n}(y_i^{\text {measured}} - {\bar{y}}^{\text {measured}})^2}, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ25.gif"/></section>
where 
<em>n</em>
 is the number of observations, 
<span id="IEq69"><mml:math id="IEq69_Math"><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mi mathvariant="normal">calculated</mml:mi></mml:msubsup></mml:math><tex-math id="IEq69_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$y_i^{\mathrm{calculated}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq69.gif"/></span>
 is the frequency of Weibull, 
<span id="IEq70"><mml:math id="IEq70_Math"><mml:msubsup><mml:mover accent="true"><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover><mml:mi>i</mml:mi><mml:mi mathvariant="normal">measured</mml:mi></mml:msubsup></mml:math><tex-math id="IEq70_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\bar{y}}_i^{\mathrm{measured}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq70.gif"/></span>
 is the mean wind speed and 
<span id="IEq71"><mml:math id="IEq71_Math"><mml:msubsup><mml:mi>y</mml:mi><mml:mi>i</mml:mi><mml:mi mathvariant="normal">measured</mml:mi></mml:msubsup></mml:math><tex-math id="IEq71_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$y_i^{\mathrm{measured}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq71.gif"/></span>
 is the frequency of observations.</p>
<p>The percentage value of the wind production deviation (WPD) between the obtained Weibull probability distribution curve and the data histogram was evaluated as in Eq. 
<a href="#Equ26"><sup>26</sup></a>
.
<section id="Equ26"><mml:math display="block" id="Equ26_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:mi>D</mml:mi><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">(</mml:mo></mml:mrow><mml:mfrac><mml:mrow><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">estimated</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">measured</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">measured</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">)</mml:mo></mml:mrow><mml:mo>·</mml:mo><mml:mn>100</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ26_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} WPD= & {} \bigg (\frac{WPD_{\mathrm{estimated}} -WPD_{\mathrm{measured}}}{WPD_{\mathrm{measured}}}\bigg )\cdot 100, \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ26.gif"/></section>
<section id="Equ27"><mml:math display="block" id="Equ27_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">measured</mml:mi></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mo>·</mml:mo><mml:mi>ρ</mml:mi><mml:mo>·</mml:mo><mml:msup><mml:mi>v</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ27_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} WPD_{\mathrm{measured}}= & {} \frac{1}{2}\cdot \rho \cdot v^3 \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ27.gif"/></section>
<section id="Equ28"><mml:math display="block" id="Equ28_Math"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">estimated</mml:mi></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd columnalign="left"><mml:mrow><mml:mrow/><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mo>·</mml:mo><mml:mi>ρ</mml:mi><mml:mo>·</mml:mo><mml:msup><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>·</mml:mo><mml:mi>Γ</mml:mi><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">(</mml:mo></mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mfrac><mml:mn>3</mml:mn><mml:mi>k</mml:mi></mml:mfrac><mml:mrow><mml:mo maxsize="2.047em" minsize="2.047em" stretchy="true">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math><tex-math id="Equ28_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\begin{aligned} WPD_{\mathrm{estimated}}= & {} \frac{1}{2}\cdot \rho \cdot c^3 \cdot \varGamma \bigg (1 + \frac{3}{k} \bigg ), \end{aligned}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_Equ28.gif"/></section>
where 
<span id="IEq72"><mml:math id="IEq72_Math"><mml:mi>ρ</mml:mi></mml:math><tex-math id="IEq72_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\rho$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq72.gif"/></span>
 is the specific mass of the air, 
<em>v</em>
 is the wind speed, 
<span id="IEq73"><mml:math id="IEq73_Math"><mml:mi>Γ</mml:mi></mml:math><tex-math id="IEq73_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\varGamma$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq73.gif"/></span>
 is the gamma function and 
<em>k</em>
 and 
<em>c</em>
 are the estimated Weibull curve parameters.</p></section>
<section><h2>Wind site data processing</h2>
<p>The data of each location were separated into intervals with a variation of 1 m/s, and to fit the interval, the velocity should be higher than the lower value of the interval and less than or equal to the upper value, except the first interval where: 0 m/s 
<span id="IEq74"><mml:math id="IEq74_Math"><mml:mo>≤</mml:mo></mml:math><tex-math id="IEq74_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\le$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq74.gif"/></span>
 
<em>V</em>
  
<span id="IEq75"><mml:math id="IEq75_Math"><mml:mo>≤</mml:mo></mml:math><tex-math id="IEq75_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\le$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq75.gif"/></span>
 1 m/s. Once separated, the data size within each interval was evaluated, and this amount of each interval was divided by the data size, thus generating a relative frequency value for each interval. The data is validated by a SONDA project methodology, which does not change the databases, eliminating data considered invalid by the process. However, this only indicates the data considered as suspicious for the user to decide whether or not to use them. Data collected by the SONDA project for the Triunfo, Petrolina and São Martinho da Serra accounted for a total of 52,560 for the three stations, although, after the processing, a total of 52,560, 52,514 and 52,366 data were considered, representing a use of 100%, 99.91% and 99.63%. (Figs. 
<a href="#Fig5"><sup>5</sup></a>
, 
<a href="#Fig6"><sup>6</sup></a>
, 
<a href="#Fig7"><sup>7</sup></a>
, 
<a href="#Fig8"><sup>8</sup></a>
, 
<a href="#Fig9"><sup>9</sup></a>
, 
<a href="#Fig10"><sup>10</sup></a>
)</p></section>
<section><h2>Results and discusssion</h2>
<p>The results of the statistical tests for the TRI23 station located in Triunfo, PTR11 station located in Petrolina and SMS08 located in São Martinho da Serra are presented in Tables  
<a href="#Tab2"><sup>2</sup></a>
, 
<a href="#Tab3"><sup>3</sup></a>
 and 
<a href="#Tab4"><sup>4</sup></a>
. Figures 
<a href="#Fig5"><sup>5</sup></a>
, 
<a href="#Fig8"><sup>8</sup></a>
, 
<a href="#Fig11"><sup>11</sup></a>
, 
<a href="#Fig6"><sup>6</sup></a>
, 
<a href="#Fig9"><sup>9</sup></a>
 and 
<a href="#Fig12"><sup>12</sup></a>
 present the Weilbull distribution curves for deterministic numerical methods and heuristic methods. Figures  
<a href="#Fig7"><sup>7</sup></a>
 and 
<a href="#Fig13"><sup>13</sup></a>
compare the results obtained by ACO and the EM for Triunfo and São Martinho da Serra and Figure 
<a href="#Fig10"><sup>10</sup></a>
 compares the results obtained by CSO and the EM for Petrolina
<figure id="Tab2"><h3>Table 2</h3>
<figcaption><p>Statistical analysis: Triunfo, year 2010</p></figcaption>
<table frame="hsides" rules="groups"><thead><tr><th align="left"><p>Method</p></th><th align="left"><p><italic>k</italic></p></th><th align="left"><p><italic>c</italic></p></th><th align="left"><p>RMSE</p></th><th align="left"><p>MAE</p></th><th align="left"><p><italic>R</italic><inline-formula id="IEq76"><alternatives><mml:math id="IEq76_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq76_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq76.gif"/></alternatives></inline-formula></p></th><th align="left"><p>WPD</p></th></tr></thead><tbody><tr><td align="left"><p>MLM</p></td><td align="left"><p>3.3337</p></td><td align="left"><p>15.2254</p></td><td align="left"><p>0.000775</p></td><td align="left"><p>0.002964</p></td><td align="left"><p>0.980183</p></td><td align="left"><p>0.080741</p></td></tr><tr><td align="left"><p>MM</p></td><td align="left"><p>3.3361</p></td><td align="left"><p>15.2103</p></td><td align="left"><p>0.000775</p></td><td align="left"><p>0.002973</p></td><td align="left"><p>0.980174</p></td><td align="left"><p><inline-formula id="IEq77"><alternatives><mml:math id="IEq77_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.241481</mml:mn></mml:mrow></mml:math><tex-math id="IEq77_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.241481$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq77.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>EM</p></td><td align="left"><p>3.3290</p></td><td align="left"><p>15.2119</p></td><td align="left"><p>0.000769</p></td><td align="left"><p>0.002956</p></td><td align="left"><p>0.980511</p></td><td align="left"><p><inline-formula id="IEq78"><alternatives><mml:math id="IEq78_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.140927</mml:mn></mml:mrow></mml:math><tex-math id="IEq78_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.140927$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq78.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>EEM</p></td><td align="left"><p>3.0004</p></td><td align="left"><p>15.0249</p></td><td align="left"><p>0.000843</p></td><td align="left"><p>0.003614</p></td><td align="left"><p>0.976576</p></td><td align="left"><p><inline-formula id="IEq79"><alternatives><mml:math id="IEq79_Math"><mml:mrow><mml:mn>2.22</mml:mn><mml:mo>·</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>14</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math><tex-math id="IEq79_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$2.22\cdot 10^{-14}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq79.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>ACO</p></td><td align="left"><p>3.1936</p></td><td align="left"><p>15.1312</p></td><td align="left"><p>0.000694</p></td><td align="left"><p>0.002821</p></td><td align="left"><p>0.984081</p></td><td align="left"><p><inline-formula id="IEq80"><alternatives><mml:math id="IEq80_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.322802</mml:mn></mml:mrow></mml:math><tex-math id="IEq80_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.322802$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq80.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>CSO</p></td><td align="left"><p>3.1930</p></td><td align="left"><p>15.1307</p></td><td align="left"><p>0.000694</p></td><td align="left"><p>0.002822</p></td><td align="left"><p>0.984081</p></td><td align="left"><p><inline-formula id="IEq81"><alternatives><mml:math id="IEq81_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.325351</mml:mn></mml:mrow></mml:math><tex-math id="IEq81_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.325351$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq81.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>HS</p></td><td align="left"><p>3.1827</p></td><td align="left"><p>15.1247</p></td><td align="left"><p>0.000695</p></td><td align="left"><p>0.002831</p></td><td align="left"><p>0.984060</p></td><td align="left"><p><inline-formula id="IEq82"><alternatives><mml:math id="IEq82_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.327501</mml:mn></mml:mrow></mml:math><tex-math id="IEq82_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.327501$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq82.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>PSO</p></td><td align="left"><p>3.3012</p></td><td align="left"><p>15.7376</p></td><td align="left"><p>0.001034</p></td><td align="left"><p>0.004310</p></td><td align="left"><p>0.964724</p></td><td align="left"><p>10.875590</p></td></tr></tbody></table></figure></p>
<p><figure id="Fig5"><h3>Fig. 5</h3>
<figcaption><p>Weibull curve adjustment for TRI23-deterministic numerical methods</p></figcaption>
<img src="40095_2018_285_Fig5_HTML.png" /></figure>
<figure id="Fig6"><h3>Fig. 6</h3>
<figcaption><p>Weibull curve adjustment for TRI23-heuristic methods</p></figcaption>
<img src="40095_2018_285_Fig6_HTML.png" /></figure>
<figure id="Fig7"><h3>Fig. 7</h3>
<figcaption><p>Weibull curve adjustment for TR23 station in Triunfo</p></figcaption>
<img src="40095_2018_285_Fig7_HTML.png" /></figure>
<figure id="Tab3"><h3>Table 3</h3>
<figcaption><p>Statistical analysis: PTR11, year 2006</p></figcaption>
<table frame="hsides" rules="groups"><thead><tr><th align="left"><p>Method</p></th><th align="left"><p><italic>k</italic></p></th><th align="left"><p><italic>c</italic></p></th><th align="left"><p>RMSE</p></th><th align="left"><p>MAE</p></th><th align="left"><p><italic>R</italic><inline-formula id="IEq83"><alternatives><mml:math id="IEq83_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq83_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq83.gif"/></alternatives></inline-formula></p></th><th align="left"><p>WPD</p></th></tr></thead><tbody><tr><td align="left"><p>MLM</p></td><td align="left"><p>3.0258</p></td><td align="left"><p>5.4536</p></td><td align="left"><p>0.002440</p></td><td align="left"><p>0.006081</p></td><td align="left"><p>0.984086</p></td><td align="left"><p><inline-formula id="IEq84"><alternatives><mml:math id="IEq84_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.009484</mml:mn></mml:mrow></mml:math><tex-math id="IEq84_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.009484$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq84.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>MM</p></td><td align="left"><p>3.0593</p></td><td align="left"><p>5.4665</p></td><td align="left"><p>0.002283</p></td><td align="left"><p>0.005948</p></td><td align="left"><p>0.986071</p></td><td align="left"><p>0.252696</p></td></tr><tr><td align="left"><p>EM</p></td><td align="left"><p>3.0595</p></td><td align="left"><p>5.4665</p></td><td align="left"><p>0.002282</p></td><td align="left"><p>0.005946</p></td><td align="left"><p>0.986082</p></td><td align="left"><p>0.249015</p></td></tr><tr><td align="left"><p>EEM</p></td><td align="left"><p>2.8169</p></td><td align="left"><p>5.3952</p></td><td align="left"><p>0.003640</p></td><td align="left"><p>0.008796</p></td><td align="left"><p>0.964585</p></td><td align="left"><p><inline-formula id="IEq85"><alternatives><mml:math id="IEq85_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>5.55</mml:mn><mml:mo>·</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>14</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math><tex-math id="IEq85_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-5.55\cdot 10^{-14}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq85.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>ACO</p></td><td align="left"><p>3.2928</p></td><td align="left"><p>5.4768</p></td><td align="left"><p>0.001727</p></td><td align="left"><p>0.004985</p></td><td align="left"><p>0.992025</p></td><td align="left"><p>−1.861296</p></td></tr><tr><td align="left"><p>CSO</p></td><td align="left"><p>3.2924</p></td><td align="left"><p>5.4774</p></td><td align="left"><p>0.001727</p></td><td align="left"><p>0.004988</p></td><td align="left"><p>0.992025</p></td><td align="left"><p><inline-formula id="IEq86"><alternatives><mml:math id="IEq86_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>1.823873</mml:mn></mml:mrow></mml:math><tex-math id="IEq86_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-1.823873$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq86.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>HS</p></td><td align="left"><p>3.2989</p></td><td align="left"><p>5.4643</p></td><td align="left"><p>0.001732</p></td><td align="left"><p>0.004923</p></td><td align="left"><p>0.991972</p></td><td align="left"><p><inline-formula id="IEq87"><alternatives><mml:math id="IEq87_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>2.588528</mml:mn></mml:mrow></mml:math><tex-math id="IEq87_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-2.588528$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq87.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>PSO</p></td><td align="left"><p>3.2921</p></td><td align="left"><p>5.4846</p></td><td align="left"><p>0.001728</p></td><td align="left"><p>0.005007</p></td><td align="left"><p>0.992011</p></td><td align="left"><p><inline-formula id="IEq88"><alternatives><mml:math id="IEq88_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>1.433228</mml:mn></mml:mrow></mml:math><tex-math id="IEq88_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-1.433228$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq88.gif"/></alternatives></inline-formula></p></td></tr></tbody></table></figure></p>
<p><figure id="Fig8"><h3>Fig. 8</h3>
<figcaption><p>Weibull curve adjustment for PTR11-deterministic numerical methods</p></figcaption>
<img src="40095_2018_285_Fig8_HTML.png" /></figure>
<figure id="Fig9"><h3>Fig. 9</h3>
<figcaption><p>Weibull curve adjustment for PTR11-heuristic methods</p></figcaption>
<img src="40095_2018_285_Fig9_HTML.png" /></figure>
<figure id="Fig10"><h3>Fig. 10</h3>
<figcaption><p>Weibull curve adjustment for PTR11 station in Petrolina</p></figcaption>
<img src="40095_2018_285_Fig10_HTML.png" /></figure></p>
<p><figure id="Tab4"><h3>Table 4</h3>
<figcaption><p>Statistical analysis: SMS08, year 2006</p></figcaption>
<table frame="hsides" rules="groups"><thead><tr><th align="left"><p>Method</p></th><th align="left"><p><italic>k</italic></p></th><th align="left"><p><italic>c</italic></p></th><th align="left"><p>RMSE</p></th><th align="left"><p>MAE</p></th><th align="left"><p><italic>R</italic><inline-formula id="IEq89"><alternatives><mml:math id="IEq89_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq89_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq89.gif"/></alternatives></inline-formula></p></th><th align="left"><p>WPD</p></th></tr></thead><tbody><tr><td align="left"><p>MLM</p></td><td align="left"><p>2.7188</p></td><td align="left"><p>3.6600</p></td><td align="left"><p>0.002666</p></td><td align="left"><p>0.006216</p></td><td align="left"><p>0.990161</p></td><td align="left"><p>0.059538</p></td></tr><tr><td align="left"><p>MM</p></td><td align="left"><p>2.7760</p></td><td align="left"><p>3.6717</p></td><td align="left"><p>0.002251</p></td><td align="left"><p>0.005308</p></td><td align="left"><p>0.992983</p></td><td align="left"><p><inline-formula id="IEq90"><alternatives><mml:math id="IEq90_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.080201</mml:mn></mml:mrow></mml:math><tex-math id="IEq90_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.080201$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq90.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>EM</p></td><td align="left"><p>2.7833</p></td><td align="left"><p>3.6713</p></td><td align="left"><p>0.002215</p></td><td align="left"><p>0.005208</p></td><td align="left"><p>0.993209</p></td><td align="left"><p><inline-formula id="IEq91"><alternatives><mml:math id="IEq91_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.244181</mml:mn></mml:mrow></mml:math><tex-math id="IEq91_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-0.244181$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq91.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>EEM</p></td><td align="left"><p>2.0573</p></td><td align="left"><p>3.4135</p></td><td align="left"><p>0.010163</p></td><td align="left"><p>0.021280</p></td><td align="left"><p>0.856988</p></td><td align="left"><p><inline-formula id="IEq92"><alternatives><mml:math id="IEq92_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mn>1.11</mml:mn><mml:mo>·</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>14</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math><tex-math id="IEq92_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$-1.11\cdot 10^{-14}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq92.gif"/></alternatives></inline-formula></p></td></tr><tr><td align="left"><p>ACO</p></td><td align="left"><p>2.9104</p></td><td align="left"><p>3.7168</p></td><td align="left"><p>0.001692</p></td><td align="left"><p>0.003699</p></td><td align="left"><p>0.996033</p></td><td align="left"><p>1.302106</p></td></tr><tr><td align="left"><p>CSO</p></td><td align="left"><p>2.9102</p></td><td align="left"><p>3.7173</p></td><td align="left"><p>0.001693</p></td><td align="left"><p>0.003700</p></td><td align="left"><p>0.996033</p></td><td align="left"><p>1.347609</p></td></tr><tr><td align="left"><p>HS</p></td><td align="left"><p>2.9100</p></td><td align="left"><p>3.7185</p></td><td align="left"><p>0.001692</p></td><td align="left"><p>0.003701</p></td><td align="left"><p>0.996032</p></td><td align="left"><p>1.454803</p></td></tr><tr><td align="left"><p>PSO</p></td><td align="left"><p>2.8539</p></td><td align="left"><p>3.7561</p></td><td align="left"><p>0.001904</p></td><td align="left"><p>0.004346</p></td><td align="left"><p>0.994979</p></td><td align="left"><p>5.516611</p></td></tr></tbody></table></figure></p>
<p><figure id="Fig11"><h3>Fig. 11</h3>
<figcaption><p>Weibull curve adjustment for SMS08—deterministic numerical methods</p></figcaption>
<img src="40095_2018_285_Fig11_HTML.png" /></figure>
<figure id="Fig12"><h3>Fig. 12</h3>
<figcaption><p>Weibull curve adjustment for SMS08—heuristics methods</p></figcaption>
<img src="40095_2018_285_Fig12_HTML.png" /></figure>
<figure id="Fig13"><h3>Fig. 13</h3>
<figcaption><p>Weibull curve adjustment for SMS08 station in São Martinho da Serra</p></figcaption>
<img src="40095_2018_285_Fig13_HTML.png" /></figure></p>
<p>Graphically, it was observed that the methods EM, MLM, MM, ACO, CSO, PSO and HS, to determine the shape parameter 
<em>k</em>
 and the scale parameter 
<em>c</em>
 of the Weibull distribution, presented a better curve fit with the histogram of the wind speed for the cities Triunfo, Petrolina and São Martinho da Serra. Moreover, it was further observed that the heuristic methods Ant Colony Optimization and Cuckoo Search Optimization were completely adequate to estimate the Weibull parameters. This fact was clearly validated by means of the statistical tests, i.e., RMSE, MAE and R
<span id="IEq93"><mml:math id="IEq93_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq93_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq93.gif"/></span>
, and by the WPD test. Tables 
<a href="#Tab2"><sup>2</sup></a>
, 
<a href="#Tab3"><sup>3</sup></a>
 and 
<a href="#Tab4"><sup>4</sup></a>
 show the statistical tests results for all deterministic and heuristic methods and WPD test considered in the analysis. It was also observed from the statistical and wind production deviation analysis that the values of RMSE, MAE, 
<em>R</em>
<span id="IEq94"><mml:math id="IEq94_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq94_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq94.gif"/></span>
 and WPD have low variation magnitudes to each other for all the methods.</p>
<p>It can be concluded that the ACO method for Triunfo and São Martinho da Serra and the CSO method for Petrolina have a good performance, since the results among all the used methods obtained the lowest values of RMSE, MAE and WPD, highlighting the WPD test values less than 2%, which was below the acceptable limit for the wind production deviation. It can also be concluded that the EEM for Triunfo, Petrolina and São Martinho da Serra has the worst performance, since it obtained the highest values of RMSE and MAE, and the lowest value of 
<em>R</em>
<span id="IEq95"><mml:math id="IEq95_Math"><mml:msup><mml:mrow/><mml:mn>2</mml:mn></mml:msup></mml:math><tex-math id="IEq95_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$^2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40095_2018_285_Article_IEq95.gif"/></span>
 among all methods, although this method presented great performance of the WPD test, since it obtained negligible values of wind production deviation. Among the heuristic methods, PSO for Triunfo and São Martinho da Serra had the worst performance, since it obtained WPD value higher than 2%.</p></section>
<section><h2>Conclusion</h2>
<p>The following conclusions can be drawn from the preceding analysis:
<ul><li><p>Graphically, the EEM method was the least effective to fit Weibull distribution curves for wind speed data from the region of Pernambuco and Rio Grande do Sul, respectively, using the data analyzed for the cities of Triunfo, Petrolina and São Martinho da Serra.</p></li>
<li><p>Regarding the parameter 
<em>k</em>
, it was observed that its values range from 2 to 3 for the cities of Triunfo, Petrolina and São Martinho da Serra, showing less constancy of the wind speed for that location. The values of 
<em>c</em>
 for Petrolina and São Martinho da Serra cities range from 3 to 6 and for Triunfo range 13 to 16 for the mean wind speed occurring in those aforementioned places.</p></li>
<li><p>Ant Colony Optimization was an efficient method to determine the Weibull distribution parameters, 
<em>k</em>
 and 
<em>c</em>
, for Triunfo and São Martinho da serra, and Cuckoo Search Optimization was an efficient method for Petrolina.</p></li>
<li><p>A suggestion for future work is to evaluate more periods of time and use the predicted values for 
<em>k</em>
 and 
<em>c</em>
 to calculate the average wind speed and its standard deviation to achieve a rank for each method.</p></li></ul></p></section><hr/><h2>Acknowledgements</h2>
<p>This research was jointly supported by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Brazilian governmental agencies.</p><hr/><h2>Publisher's Note</h2>
<p>Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
<hr/><h2>References</h2>
<ol><li>Akdag (2009) <em>Dinler: a new method to estimate weibull parameters for wind energy applications</em> 50(7) (pp. 1761-1766) <a href="https://doi.org/10.1016/j.enconman.2009.03.020" target="_blank">10.1016/j.enconman.2009.03.020</a></li><li>Akpinar and Akpinar (2004) <em>Determination of the wind energy potential for maden-elazig</em> (pp. 2901-2914) <a href="https://doi.org/10.1016/j.enconman.2003.12.016" target="_blank">10.1016/j.enconman.2003.12.016</a></li><li>Andrade et al. (2014) <em>An efficiency comparison of numerical methods for determining weibull parameters for wind energy applications:a new approach applied to the northeast region of brazil</em> (pp. 801-808) <a href="https://doi.org/10.1016/j.enconman.2014.06.046" target="_blank">10.1016/j.enconman.2014.06.046</a></li><li>Askarzadeh and Zebarjadi (2014) <em>Wind power modeling using harmony search with a novel parameter setting approach</em> (pp. 70-75) <a href="https://doi.org/10.1016/j.jweia.2014.10.012" target="_blank">10.1016/j.jweia.2014.10.012</a></li><li>Unknown () <em></em></li><li>Unknown () <em></em></li><li>Benhala and Bouattane (2014) <em>Ga and aco techniques for the analog circuits design optimization</em> (pp. 413-419)</li><li>Burton et al. (2001) <em></em> Wiley <a href="https://doi.org/10.1002/0470846062" target="_blank">10.1002/0470846062</a></li><li>Celik (2003) <em>A statistical analysis of wind power density based on the weibull and rayleigh models at the southern region of turkey</em> (pp. 593-604) <a href="https://doi.org/10.1016/j.renene.2003.07.002" target="_blank">10.1016/j.renene.2003.07.002</a></li><li>Chang (2011) <em>Estimation of wind energy potential using different probability density functions</em> (pp. 1848-1856) <a href="https://doi.org/10.1016/j.apenergy.2010.11.010" target="_blank">10.1016/j.apenergy.2010.11.010</a></li><li>Dorigo and stützle (2004) <em></em> MIT Press <a href="https://doi.org/10.7551/mitpress/1290.001.0001" target="_blank">10.7551/mitpress/1290.001.0001</a></li><li>Dorvlo (2002) <em>Estimating wind speed distribution</em> 43(17) (pp. 2311-2318) <a href="https://doi.org/10.1016/S0196-8904(01)00182-0" target="_blank">10.1016/S0196-8904(01)00182-0</a></li><li>Unknown () <em></em></li><li>Engelbrecht (2007) <em></em> Wiley <a href="https://doi.org/10.1002/9780470512517" target="_blank">10.1002/9780470512517</a></li><li>Fisher (1915) <em>Frequency distribution of the values of the correlation coefficient in samples from an indefinitely large population</em> (pp. 507-521)</li><li>Geem et al. (2001) <em>A new heuristic optimization algorithm: harmony search</em> 76(2) (pp. 60-68) <a href="https://doi.org/10.1177/003754970107600201" target="_blank">10.1177/003754970107600201</a></li><li>González et al. (2017) <em>Optimal wind-turbine micro-siting of offshore wind farms: A grid-like layout approach</em> (pp. 28-38) <a href="https://doi.org/10.1016/j.apenergy.2017.05.071" target="_blank">10.1016/j.apenergy.2017.05.071</a></li><li>Hajibandeh et al. (2018) <em>A heuristic multi-objective multi-criteria demand response planning in a system with high penetration of wind power generators</em> (pp. 721-732) <a href="https://doi.org/10.1016/j.apenergy.2017.12.076" target="_blank">10.1016/j.apenergy.2017.12.076</a></li><li>Jiang et al. (2017) <em>Comparison of numerical methods and metaheuristic optimization algorithms for estimating parameters for wind energy potential assessment in low wind regions</em> (pp. 1199-1217) <a href="https://doi.org/10.1016/j.rser.2016.11.241" target="_blank">10.1016/j.rser.2016.11.241</a></li><li>Jiang et al. (2017) <em>Short-term wind speed forecasting using a hybrid model</em> (pp. 561-577) <a href="https://doi.org/10.1016/j.energy.2016.10.040" target="_blank">10.1016/j.energy.2016.10.040</a></li><li>Justus and Mikhail (1976) <em>Height variation of wind speed and wind distribution statistics</em> (pp. 261-264) <a href="https://doi.org/10.1029/GL003i005p00261" target="_blank">10.1029/GL003i005p00261</a></li><li>Manwell et al. (2009) <em></em> Wiley <a href="https://doi.org/10.1002/9781119994367" target="_blank">10.1002/9781119994367</a></li><li>Marzband et al. (2018) <em>Framework for smart transactive energy in home-microgrids considering coalition formation and demand side management</em> <a href="https://doi.org/10.1016/j.scs.2018.04.010" target="_blank">10.1016/j.scs.2018.04.010</a></li><li>Ohunakin et al. (2011) <em>Wind energy evaluation for electricity generation using wecs in seven selected locations in nigeria</em> (pp. 3197-3206) <a href="https://doi.org/10.1016/j.apenergy.2011.03.022" target="_blank">10.1016/j.apenergy.2011.03.022</a></li><li>Osman and Laporte (1996) <em>Metaheuristics: a bibliography</em> (pp. 513-628) <a href="https://doi.org/10.1007/BF02125421" target="_blank">10.1007/BF02125421</a></li><li>Payne et al. (2005) <em></em> Oxford University Press</li><li>Rahmani et al. (2013) <em>Hybrid technique of ant colony and particle swarm optimization for short term wind energy forecasting</em> (pp. 163-170) <a href="https://doi.org/10.1016/j.jweia.2013.10.004" target="_blank">10.1016/j.jweia.2013.10.004</a></li><li>Rocha et al. (2012) <em>Comparison of seven numerical methods for determining weibull parameters for wind energy generation in the northeast region of Brazil</em> 89(1) (pp. 395-400) <a href="https://doi.org/10.1016/j.apenergy.2011.08.003" target="_blank">10.1016/j.apenergy.2011.08.003</a></li><li>Salcedo-Sanz et al. (2018) <em>Wind power field reconstruction from a reduced set of representative measuring points</em> (pp. 1111-1121) <a href="https://doi.org/10.1016/j.apenergy.2018.07.003" target="_blank">10.1016/j.apenergy.2018.07.003</a></li><li>Sanajaoba and Fernandez (2016) <em>Maiden application of cuckoo search algorithm for optimal sizing of a remote hybrid renewable energy system</em> (pp. 1-10) <a href="https://doi.org/10.1016/j.renene.2016.04.069" target="_blank">10.1016/j.renene.2016.04.069</a></li><li>Unknown () <em></em></li><li>Unknown () <em></em></li><li>Unknown () <em></em></li><li>Unknown () <em></em></li><li>Wais (2017) <em>Two and three-parameter weibull distribution in available windpower analysis</em> (pp. 15-29) <a href="https://doi.org/10.1016/j.renene.2016.10.041" target="_blank">10.1016/j.renene.2016.10.041</a></li><li>Wang et al. (2016) <em>Wind energy potential assessment and forecasting research based on the data pre-processing technique and swarm intelligent optimization algorithms</em> 8(11) <a href="https://doi.org/10.3390/su8111191" target="_blank">10.3390/su8111191</a></li><li>Yang (2010) <em></em> Luniver Press</li></ol>