An extension of Nelder-Mead method to nonlinear mixed-integer optimization problems

==Resumen==

Este artículo presenta un algoritmo nuevo basado en una extensión del método algorítmico simplex de Nelder Mead para la identificación de al menos un óptimo local cuando se usa en problemas no lineales enteros mixtos irrestrictos. Este método algorítmico, denominado Algoritmo Simplex Entero Mixto (ASEM) por el autor, se basa en una doble estructura de símpleces que está compuesta por una estructura simplex real de dimensión ''n''  (simplex real) y otra estructura simplex entera de igual dimensión ''n''  (simplex entero). Las operaciones propuestas originalmente por Nelder Mead se ejecutan sobre el simplex real. Por otra parte, un grupo de nuevas operaciones presentadas en este artículo se emplean en el simplex entero. Este conjunto de nuevas operaciones junto con las operaciones originales del método de Nelder Mead garantizan que en cada iteración del ASEM se arroje un nuevo punto de prueba definido en el campo numérico de dimensión 2''n  ''  enteros mixtos <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math> , con el objetivo de garantizar la identificación del óptimo local entero mixto, sin convertir las variables enteras en variables reales.

==Abstract==

This article presents a new algorithm based on the Nelder-Mead simplex algorithmic method for identifying a local optimum, at least on unconstrained nonlinear mixed-integer problems. The algorithmic method, Integer Mixed Simplex Algorithm (IMSA), so called by the author, is based on a double simplex structure, which is composed of a real ''n'' -dimensional simplex structure (real simplex) and an integer ''n'' -dimensional simplex structure (integer simplex). The original Nelder-Mead operations are applied on the real simplex. Meanwhile, a novel group of operations are applied on the integer simplex. This new set of operations, together with the original Nelder-Mead operations, guarantee a new trail point at each IMSA iteration in the search of the local optimum in the integer real mixed 2''n  '' -dimensional numerical field <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math>  without the need of integer to real conversions.  

==Palabras clave==

Método simplex ; Método de Nelder Mead ; optimización entera mixta ; búsqueda directa

==Keywords==

Simplex method ; Nelder-Mead method ; Mixed-integer optimization ; Direct search

==1. Introducción==

El desarrollo de métodos de optimización para su empleo en la búsqueda de soluciones a problemas enteros mixtos ocupa un lugar importante en diversos campos de la ingeniería, debido a las amplias posibilidades de aplicación en problemas de diseño óptimo, cuando en este tipo de problemas las variables enteras no pueden ser relajadas a variables reales, que suelen ser abordados con el propósito de emplear algún método algorítmico de optimización desarrollado para variables reales. Además, esta situación de inconvertibilidad de las variables enteras a variables reales se presenta con frecuencia cuando la función objetivo del problema de optimización es representada mediante la función de desempeño del sistema bajo estudio, la cual es medida mediante algún modelo de simulación que depende tanto de variables reales como de variables enteras. Bajo estas circunstancias es inconveniente transformar el problema de optimización entero mixto (OEM) a un problema de optimización de variables reales, debido al hecho de que existen variables de decisión que no pueden ser tratadas como variables reales. Como ejemplo de sistemas que son definidos por variables enteras mixtas se tienen: el número de cajeros en un supermercado y su respectiva tasa de atención; la capacidad de clientes que pueden aguardar a ser atendidos en cada cola correspondiente a su taquilla de una agencia bancaria y la tasa de servicio de cada cajero; el número de dientes de cada engranaje que conforman una caja reductora y la dimensión de cada engranaje; el número de etapas de un sistema amplificador de señales y la ganancia de cada etapa de amplificación entre otros.

Es en esta situación donde los métodos de búsqueda directa, llamados así por primera vez por Hooke y Jeeves [[#bib0005|[1]]] , pueden brindar una solución a este tipo de problemas de optimización, bien sea por no disponer de una expresión matemática de la función objetivo o por no contar con el gradiente de esta, principalmente cuando se requiere optimizar un proceso que depende tanto de variables reales como de variables enteras, lo que exige la aplicación de estrategias de exploración que permitan avanzar a un mejor punto de prueba en cada iteración, con el fin de identificar la solución óptima del problema.      

En este artículo se presenta una extensión del método simplex de Nelder Mead (MSNM) [[#bib0010|[2]]]  a problemas OEM cuando la función objetivo es no lineal y no está sujeta a restricciones. El Algoritmo Simplex Entero Mixto (ASEM), denominado así por el autor debido a que se basa en el patrón de búsqueda conocido como simplex, se fundamenta en las operaciones definidas por el MSNM y en un conjunto de nuevas operaciones propuestas por Brea [[#bib0015|[3]]] , que son ejecutadas únicamente sobre las variables enteras.      

De acuerdo con Brandts et al. [[#bib0020|[4]]] , puede decirse que un simplex (''d  '' -simplex) definido en el espacio euclidiano <math display="inline">{\mathbb{R}}^d</math> , siendo ''d''  un entero positivo, es un casco o envoltorio convexo de ''d''  + 1 puntos, los cuales no todos pertenecen al mismo hiperplano definido en el espacio <math display="inline">{\mathbb{R}}^d</math> , y cuyos puntos son denominados vértices del simplex.      

En consecuencia, un 1-simplex es un segmento de recta en el espacio <math display="inline">{\mathbb{R}}^1</math> ; un 2-simplex es un triángulo en el espacio <math display="inline">{\mathbb{R}}^2</math> ; un 3-simplex es un tetraedro en el espacio <math display="inline">{\mathbb{R}}^3</math> , y así sucesivamente.      

Debido a que los problemas de optimización enteros mixtos son abordados en esta investigación mediante estructuras de simplex, se hace necesario definir 2 símpleces: uno para las variables reales, denominada aquí simplex real y otro que es empleado en el campo entero, denominado en este artículo simplex entero. Sobre estos símpleces es definido un conjunto de operaciones a fin de explorar en el espacio de búsqueda entero mixto y así poder identificar la solución óptima local del problema.

El ASEM ha mostrado ser lo suficientemente eficiente a problemas OEM irrestrictos [[#bib0015|[3]]] . No obstante, también ha sido probado con problemas no lineales enteros mixtos bajo restricciones a través de funciones de penalización, lo cual ha permitido comprobar que el ASEM es lo suficientemente eficaz en la identificación de óptimos locales en problemas de OEM bajo restricciones.      

Mediante la ejecución de un número importante de ejemplos numéricos, Brea [[#bib0025|[5]]]  ha comprobado el buen desempeño del MSNM cuando la función objetivo está sujeta a ruido, y además cuando el problema está sujeto a restricciones lineales. Sin embargo, existen condiciones particulares en donde el método algorítmico que él desarrolló [[#bib0025|[5]]]  converge muy lentamente al óptimo.      

Por otra parte, Barton e Ivey [[#bib0030|[6]]]  y Humphrey y Wilson [[#bib0035|[7]]]  han estudiado los parámetros que ofrecen el mejor desempeño del MSNM, cuando la función objetivo de un problema de optimización irrestricto está influenciada por ruido.      

Estos hechos han motivado al autor a desarrollar un método algorítmico basado en el MSNM en el campo de los números enteros mixtos.

El resto de este artículo está estructurado como sigue: en la sección [[#sec0010|2]]  se describe matemáticamente el tipo de problema de OEM que hay que resolver mediante el ASEM; en la sección [[#sec0015|3]]  se enuncian un conjunto de definiciones preliminares que permitirán exponer claramente las proposiciones que definen las nuevas operaciones que se deben emplear sobre el simplex entero, además de facilitar la explicación del algoritmo; en la sección [[#sec0020|4]]  se enuncian las proposiciones que definen las nuevas operaciones que son aplicadas sobre el simplex entero; en la sección [[#sec0025|5]]  son establecidos los procedimientos empleados en el ASEM; en la sección [[#sec0030|6]]  se ofrece una explicación del ASEM mediante una inédita representación basada en grafo; en la sección [[#sec0035|7]]  se presenta un modo de flexibilizar el ASEM con el propósito de emplearlo en la optimización de problemas, en donde el número de componentes reales es diferente al número de componentes enteras de las variables de decisión; la sintonización o los ajustes de los parámetros del ASEM se muestran en la sección [[#sec0040|8]] ; en la sección [[#sec0045|9]]  se muestra un conjunto de experimentos numéricos con el objeto de mostrar la eficacia del ASEM; una comparación del ASEM con otros algoritmos de optimización en el problema de diseño óptimo de un tanque presurizado se presenta en la sección [[#sec0050|10]] ; finalmente, se presentan conclusiones y propuestas de futuras investigaciones en la sección [[#sec0055|11]] .      

<span id='sec0010'></span>
==2. El problema==

<span id='enun0005'></span>

==Problema 1.                     ==

Sean <math display="inline">\boldsymbol{\mbox{x}}\in {\mathbb{R}}^n</math>  el conjunto de ''n  ''  variables reales independientes e <math display="inline">\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^m</math>  el conjunto de ''m''  variables enteras independientes. Entonces, el problema de minimización irrestricto entero mixto puede ser expresado por:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\underset{\boldsymbol{\mbox{x}}\in {\mathbb{R}}^n,\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^m}{minimice\quad }f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}}),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 1)
|}

donde <math display="inline">f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}}):{\mathbb{R}}^n\times {\mathbb{Z}}^m\rightarrow \mathbb{R}</math>  es una función objetivo no lineal.

<span id='sec0015'></span>
==3. Definiciones preliminares==

En esta sección se establecerá un conjunto de definiciones que facilitarán la explicación del método algorítmico desarrollado.

==Definición 1 Espacio euclidiano entero mixto.                     ==

Se entenderá por espacio euclidiano entero mixto de dimensión ''n''  + ''m  ''  el conjunto de puntos definidos por el producto cartesiano <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^m</math> , donde <math display="inline">\mathbb{R}</math>  es el conjunto de los números reales y <math display="inline">\mathbb{Z}</math>  es el conjunto de los números enteros.

Es importante señalar que, debido a la naturaleza del problema, el espacio de exploración donde deben ejecutarse las operaciones matemáticas es el espacio euclidiano entero mixto <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^m</math> , cuya dimensión corresponde a ''n''  + ''m'' , indicando respectivamente así el número de variables reales y enteras que conforman el problema.

<span id='enun0015'></span>

==Definición 2 Vector precedente.                     ==

Sea una función <math display="inline">f(\boldsymbol{\mbox{w}}):W\rightarrow \mathbb{R}</math>  y sean '''a'''  y '''b'''  2 vectores cualesquiera definidos en un dominio <math display="inline">W</math> . Se dice que '''a'''  precede o es igualmente precedente a '''b''' , y se denota mediante '''a'''  ⪯ '''b''' , si ''f'' ('''a''' ) ≤ ''f'' ('''b''' ).

La [[#enun0015|definición 2]]  también se puede entender en el sentido estricto, esto es: un vector precede estrictamente a otro si la relación entre los valores de la función es estrictamente menor, es decir, '''a'''  ≺ '''b'''  si ''f'' ('''a''' ) < ''f'' ('''b''' ).

<span id='enun0020'></span>

==Definición 3 Vector subsecuente.                     ==

Sea una función <math display="inline">f(\boldsymbol{\mbox{w}}):W\rightarrow \mathbb{R}</math>  y sean '''a'''  y '''b'''  2 vectores cualesquiera definidos en un dominio <math display="inline">W</math> . Se dice que '''a'''  es subsecuente o es igualmente subsecuente a '''b''' , y se denota como '''a'''  ⪰ '''b''' , si ''f'' ('''a''' ) ≥ ''f'' ('''b''' ).

También la [[#enun0020|definición 3]]  se puede entender en su sentido estrictamente subsecuente. En este caso se dice que el vector '''a'''  es estrictamente subsecuente al vector '''b''' , y se denota como '''a'''  ≻ '''b''' , si ''f'' ('''a''' ) > ''f'' ('''b''' ).

==Definición 4 Vector entero mixto.                     ==

Sean <math display="inline">\boldsymbol{\mbox{x}}\in {\mathbb{R}}^n</math>  e <math display="inline">\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^m</math>  2 subvectores que conforman un vector <math display="inline">\boldsymbol{\mbox{v}}={\left[\boldsymbol{\mbox{x}}^t\quad \boldsymbol{\mbox{y}}^t\right]}^t</math> . Se dice que el vector '''v'''  es un vector entero mixto en el espacio euclidiano entero mixto de dimensión ''n''  + ''m  '' , para indicar respectivamente la dimensión de las componentes reales y enteras, si el vector <math display="inline">\boldsymbol{\mbox{v}}\in {\mathbb{R}}^n\times {\mathbb{Z}}^m</math> .

==Definición 5 Subvector real.                     ==

Sea <math display="inline">\boldsymbol{\mbox{v}}\in {\mathbb{R}}^n\times {\mathbb{Z}}^m</math>  un vector entero mixto. Se dice entonces que <math display="inline">\boldsymbol{\mbox{x}}\in {\mathbb{R}}^n</math>  es el subvector real de dimensión ''n'' , si sus componentes representan las cantidades reales del vector entero mixto '''v''' .

==Definición 6 Subvector entero.                     ==

Sea <math display="inline">\boldsymbol{\mbox{v}}\in {\mathbb{R}}^n\times {\mathbb{Z}}^m</math>  un vector entero mixto. Se dice que <math display="inline">\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^m</math>  es el subvector entero de dimensión ''m'' , si sus componentes representan cantidades enteras del vector entero mixto '''v''' .

Debido a que el método algorítmico desarrollado para resolver problemas enteros mixtos requiere que las operaciones sobre cada subvector se ejecuten de forma independiente, es decir, que los nuevos puntos sean definidos por operaciones sobre los subvectores reales y sobre los subvectores enteros de forma separada, y que se permita así la coexistencia de ambos conjuntos de subvectores sin que ellos cambien sus condiciones en cuanto al espacio numérico que lo definen, se hacen necesarias las definiciones de operaciones sobre los subvectores pertenecientes a un mismo campo numérico.

Como consecuencia de las exigencias impuestas, se estudiará el caso particular donde ''n''  = ''m  '' , quedando entonces las siguientes definiciones en el espacio euclidiano <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math> . No obstante, se propondrá una adecuación a los problemas en el espacio euclidiano entero mixto <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^m</math>  para cuando ''n''  ≠ ''m'' , y de esta forma poder aplicar el algoritmo desarrollado a situaciones en donde el número de variables reales y enteras son distintos.

==Definición 7 Centroide.                     ==

Sea un conjunto de ''k  ''  puntos <math display="inline">\boldsymbol{\mbox{p}}_i={\left(\boldsymbol{\mbox{x}}_i^t\quad \boldsymbol{\mbox{y}}_i^t\right)}^t</math>  definidos en el espacio euclidiano entero mixto <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math>  para cada ''i''  = 1, … ''k  '' . Se dice que el centroide <math display="inline">\overline{\boldsymbol{\mbox{p}}}</math>  es el centro de masa de los ''k''  puntos, si a cada uno de los puntos se le ha considerado una masa ''m''<sub>''i''</sub> . En el caso de que la masa de cada uno de los puntos '''p'''<sub>''i''</sub>  sea igual, el centroide es determinado por

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\overline{\boldsymbol{\mbox{p}}}=\frac{1}{k}\sum_{i=1}^k\boldsymbol{\mbox{p}}_i.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 2)
|}

==Observación 1.                     ==

Nótese que el punto <math display="inline">\overline{\boldsymbol{\mbox{p}}}</math>  no necesariamente pertenece al espacio euclidiano <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math> . Más aún, <math display="inline">\overline{\boldsymbol{\mbox{p}}}=</math><math>{\left({\overline{\boldsymbol{\mbox{x}}}}^t,{\overline{\boldsymbol{\mbox{y}}}}^t\right)}^t</math> , donde

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\overline{x}}=\frac{1}{k}\sum_{i=1}^k\boldsymbol{x}_i,\mbox{ }\overline{\boldsymbol{\mbox{y}}}=</math><math>\frac{1}{k}\sum_{i=1}^k\boldsymbol{y}_i</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 3)
|}

<span id='enun0050'></span>

==Definición 8 Simplex entero mixto.                     ==

Se dice que un ''q  '' -ésimo simplex entero mixto en el espacio euclidiano entero mixto <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math>  es un conjunto de diferentes puntos <math display="inline">\boldsymbol{\mbox{v}}_i={\left[\boldsymbol{\mbox{x}}_i^t\quad \boldsymbol{\mbox{y}}_i^t\right]}^t</math>  para todo ''i''  = 1, …, ''ν  '' , donde <math display="inline">\boldsymbol{\mbox{x}}_i\in {\mathbb{R}}^n</math>  corresponde al ''i'' -ésimo subvector real de dimensión ''n  '' , <math display="inline">\boldsymbol{\mbox{y}}_i\in {\mathbb{Z}}^n</math>  denota cada ''i'' -ésimo subvector entero de dimensión ''n'' , ''ν''  = ''n''  + 1 representa el número de vértices del simplex entero mixto, los cuales no todos pertenecen a la misma hipercara definida en el espacio euclidiano <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math> , y cada uno es representado por un vector de dimensión ''n''  + ''n'' . El ''q  '' -ésimo simplex entero mixto puede entonces ser representado en notación matricial como <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}=</math><math>[\boldsymbol{\mbox{v}}_1^{\left[q\right]}\quad ;\boldsymbol{\mbox{v}}_2^{\left[q\right]}\quad ;\cdots \quad ;\boldsymbol{\mbox{v}}_{\nu }^{\left[q\right]}]</math> .

De acuerdo con la [[#enun0050|definición 8]] , se establece que

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}=\left[\begin{array}{c}
\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}\\
\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}\\

\end{array}\right]=\left[\begin{array}{cccc}
\boldsymbol{\mbox{x}}_1 & \boldsymbol{\mbox{x}}_2 & \cdots  & \boldsymbol{\mbox{x}}_{\nu }\\
\boldsymbol{\mbox{y}}_1 & \boldsymbol{\mbox{y}}_2 & \cdots  & \boldsymbol{\mbox{y}}_{\nu }\\
 &  &  & 
\end{array}\right],</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 4)
|}

donde <math display="inline">\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}\in {\mathbb{R}}^{n\times (n+1)}</math>  e <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}\in {\mathbb{Z}}^{n\times (n+1)}</math>  representan respectivamente el simplex real y el simplex entero.      

Un concepto importante que debe conocerse es el referente al significado de simplex colapsado o degenerado, el cual se identifica cuando todos los vértices del simplex pertenecen a un mismo hiperplano.

==Definición 9 Menor simplex entero.                     ==

Sea la matriz <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}=</math><math>[\boldsymbol{\mbox{y}}_1^{\left[q\right]}\quad ;\boldsymbol{\mbox{y}}_2^{\left[q\right]}\quad ;\cdots \quad ;\boldsymbol{\mbox{y}}_{\nu }^{\left[q\right]}]</math>  un ''q  '' -ésimo simplex entero. Se dice que <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  es el más pequeño simplex entero no colapsado contenido en <math display="inline">{\mathbb{Z}}^{n\times (n+1)}</math> , si existe un vértice que equidista de sus otros vértices del simplex entero <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  en 1.

<span id='enun0060'></span>

'''Definición 10 Aristas de <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math> .                                                         '''

Una ''p'' -ésima matriz de aristas de un ''q  '' -ésimo simplex entero mixto <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  en el espacio euclidiano entero mixto <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math>  se define como una matriz de dimensión 2''n''  × ''n'' , cuya ''j  '' -ésima columna representa la arista del simplex entero mixto <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  entre un vértice referencial '''v'''<sub>''p''</sub>  y '''v'''<sub>''j''</sub>  para todo ''j''  = 1, …, ''ν''  y ''j''  ≠ ''p'' , es decir,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{A}}_p^{\left[q\right]}=[\boldsymbol{\mbox{v}}_1-</math><math>\boldsymbol{\mbox{v}}_p;\cdots ;\boldsymbol{\mbox{v}}_{j\not =p}-</math><math>\boldsymbol{\mbox{v}}_p;\cdots ;\boldsymbol{\mbox{v}}_{\nu }-</math><math>\boldsymbol{\mbox{v}}_p],\forall p=1,\ldots ,\nu .</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 5)
|}

Nótese que de acuerdo con la [[#enun0060|definición 10]] , una matriz <math display="inline">\boldsymbol{\mbox{A}}_p^{\left[q\right]}</math>  representa el conjunto de aristas del ''q  '' -ésimo simplex entero mixto <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  que comparten un mismo ''p'' -ésimo vértice, es decir, son todas las aristas que se intersectan en el ''p'' -ésimo vértice del ''q  '' -ésimo simplex <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math> .

'''Definición 11 Aristas de <math display="inline">\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}</math> .                                                         '''

Una ''p'' -ésima submatriz de aristas de un ''q  '' -ésimo simplex real <math display="inline">\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}</math>  en el espacio euclidiano <math display="inline">{\mathbb{R}}^n</math>  se define como una matriz de dimensión ''n''  × ''n'' , cuya ''j  '' -ésima columna representa la arista del simplex real <math display="inline">\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}</math>  entre un vértice referencial '''x'''<sub>''p''</sub>  y '''x'''<sub>''j''</sub>  para todo ''j''  = 1, …, ''ν''  y ''j''  ≠ ''p'' , es decir,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{A}}_{x_p}^{\left[q\right]}=[\boldsymbol{\mbox{x}}_1-</math><math>\boldsymbol{\mbox{x}}_p\quad ;\cdots :\boldsymbol{\mbox{x}}_{j\not =p}-</math><math>\boldsymbol{\mbox{x}}_p\quad ;\cdots \quad ;\boldsymbol{\mbox{x}}_{\nu }-</math><math>\boldsymbol{\mbox{x}}_p],</math>
|-
|<math>\forall p=1,\ldots ,\nu .</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 6)
|}

De igual manera puede ser definida la matriz de aristas del simplex entero, la cual se puede enunciar como:

'''Definición 12 Aristas de <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math> .                                                         '''

Una ''p'' -ésima submatriz de aristas de un ''q  '' -ésimo simplex entera <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  en el espacio euclidiano <math display="inline">{\mathbb{Z}}^n</math> , se define como una matriz de dimensión ''n''  × ''n'' , cuya ''j  '' -ésima columna representa la arista del simplex entero <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  entre un vértice referencial '''y'''<sub>''p''</sub>  e '''y'''<sub>''j''</sub>  para todo ''j''  = 1, …, ''ν''  y ''j''  ≠ ''p'' , es decir,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{A}}_{y_p}^{\left[q\right]}=[\boldsymbol{\mbox{y}}_1-</math><math>\boldsymbol{\mbox{y}}_p\quad ;\cdots \quad ;\boldsymbol{\mbox{y}}_{j\not =p}-</math><math>\boldsymbol{\mbox{y}}_p\quad ;\cdots \quad ;\boldsymbol{\mbox{y}}_{\nu }-</math><math>\boldsymbol{\mbox{y}}_p],</math>
|-
|<math>\forall p=1,\ldots ,\nu .</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 7)
|}

==Observación 2.                     ==

Nótese que un simplex real <math display="inline">\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}</math>  es colapsado si la matriz de arista <math display="inline">\boldsymbol{\mbox{A}}_{x_p}^{\left[q\right]}</math>  es singular para todo ''j''  = 1, …, ''ν  '' . De igual modo se puede afirmar que un simplex entero <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  es colapsado o degenerado si la matriz <math display="inline">\boldsymbol{\mbox{A}}_{y_p}^{\left[q\right]}</math>  es singular para todo ''j''  = 1, …, ''ν'' .

'''Definición 13 Simplex <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  jerarquizado.                                                         '''

Se dice que un ''q  '' -ésimo simplex entero mixto <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  es jerarquizado de acuerdo con ''f'' ('''v''' ) si '''v'''<sub>1</sub>  ⪯ '''v'''<sub>2</sub>  ⪯ ⋯ ⪯ '''v'''<sub>''ν''</sub> , lo que significa que ''f'' ('''v'''<sub>1</sub> ) ≤ ''f'' ('''v'''<sub>2</sub> ) ≤ ⋯ ≤ ''f'' ('''v'''<sub>''ν''</sub> ), en el espacio euclidiano entero mixto <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^n</math> .

==Observación 3.                     ==

Cuando el simplex no esté jerarquizado según algún criterio, esto se indica con el símbolo tilde, es decir, como <math display="inline">{\tilde{\boldsymbol{\mbox{Y}}}}_{\nu }^{\left[q\right]}</math> . Por otra parte, nótese que cuando un ''q  '' -ésimo simplex mixto está jerarquizado, tanto su correspondiente simplex real <math display="inline">\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}</math>  como su simplex entero <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  están jerarquizados de acuerdo con ''f'' ('''v''' ).

==Definición 14 Función signo de un entero.                     ==

Sea un número entero <math display="inline">k\in \mathbb{Z}</math> . Entonces la función signo entero de ''k''  y que es denotada por sgnd (''k'' ) adquiere el valor de 1, si ''k''  > 0; 0, si ''k''  = 0; o −1, si ''k''  < 0.

<span id='sec0020'></span>
==4. Proposiciones de nuevas operaciones==

Debido a que el ASEM requiere contar con un conjunto de operaciones que han de ser aplicadas sobre la submatriz de elementos enteros <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math> , se requiere establecer un conjunto de proposiciones que formarán parte de los operadores que actuarán sobre el simplex entero mixto jerarquizado de acuerdo con ''f'' ('''v''' ), junto con las ya conocidas operaciones de Nelder Mead [[#bib0010|[2]]] , las cuales son aplicadas sobre el simplex <math display="inline">\boldsymbol{\mbox{X}}_{\nu }^{\left[q\right]}</math> .

==Proposición 1 Punto reflejado entero.                     ==

Sea un ''q  '' -ésimo simplex entero jerarquizado <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  de acuerdo a ''f'' ('''v''' ), en el espacio euclidiano entero <math display="inline">{\mathbb{Z}}^n</math> . Si el punto reflejado entero '''v'''<sub>''r''</sub>  a través del centroide de su hipercara <math display="inline">H_{\nu }^{\left[q\right]}=\lbrace \boldsymbol{\mbox{v}}_i\in \boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}\vert \quad i=</math><math>1,\ldots ,\nu -1\rbrace </math>  opuesta al ''ν'' -ésimo vértice '''v'''<sub>''ν''</sub>  es calculado por

<span id='eq0040'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{y}}_r=\boldsymbol{\mbox{y}}_{\nu }+</math><math>{\alpha }_e\quad \mu \quad sgnd\quad (\overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu }),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 8)
|}

donde para un <math display="inline">{\alpha }_e\in \mathbb{N}(\geq 2)</math> , el cual es el denominado coeficiente de reflexión entero, <math display="inline">\mu =\lceil \parallel \overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu }\parallel \rceil </math>  es el valor entero próximo superior de la distancia entre '''y'''<sub>''ν''</sub>  e <math display="inline">\overline{\boldsymbol{\mbox{y}}}</math> , dado que el operador ⌈· ⌉ aproxima el resultado a su próximo superior entero, e <math display="inline">\overline{\boldsymbol{\mbox{y}}}</math>  es determinado por <math display="inline">\overline{\boldsymbol{\mbox{y}}}=</math><math>\frac{1}{\nu -1}{\sum }_{i=1}^{\nu -1}\boldsymbol{\mbox{y}}_i</math> . Entonces, se cumple que:
* '''y'''<sub>''r''</sub>  es vector definido en el espacio entero <math display="inline">{\mathbb{Z}}^n</math> ;                                  
* '''y'''<sub>''r''</sub>  e '''y'''<sub>''ν''</sub>  están a sendos lados de la hipercara <math display="inline">H_{\nu }^{\left[q\right]}</math> .          

''Demostración''  Parte I) Denote (''y''<sub>1,''ν''</sub> , …, ''y''<sub>''n'' ,''ν''</sub> )<sup>''t''</sup>  el vértice '''y'''<sub>''ν''</sub> , es decir, '''y'''<sub>''ν''</sub>  = (''y''<sub>1,''ν''</sub> , …, ''y''<sub>''n'' ,''ν''</sub> )<sup>''t''</sup> . Entonces, de acuerdo con la ecuación [[#eq0040|(8)]] , se tiene que

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{y}}_r=\boldsymbol{\mbox{y}}_{\nu }+</math><math>{\alpha }_e\quad \mu \quad {\left(k_1,\ldots ,k_n\right)}^t=</math><math>{\left(y_{1,\nu }+{\alpha }_e\quad \mu \quad k_1,\ldots ,y_{n,\nu }+{\alpha }_e\quad \mu \quad k_n\right)}^t,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}

donde ''k''<sub>''i''</sub>  ∈ { −1, 0, 1} para todo ''i''  = 1, …, ''n''  debido a la función signo discreta que es aplicada sobre cada elemento del vector. En consecuencia, cada una de las componentes de '''y'''<sub>''r''</sub>  es entera, en virtud de que ''y''<sub>''i'' ,''ν''</sub> , ''α''<sub>''e''</sub> , ''μ''  y ''k''<sub>''i''</sub>  son todos números enteros, y para cada ''i''  = 1, …, ''n''  su operación arroja un número entero.      

Parte II) Si tanto el vértice '''y'''<sub>''ν''</sub>  como el punto '''y'''<sub>''r''</sub>  están ubicados en sendos lados de la hipercara <math display="inline">H_{\nu }^{\left[q\right]}</math> , entonces el ángulo entre los vectores formados por <math display="inline">\boldsymbol{\mbox{y}}_{\nu }-</math><math>\overline{\boldsymbol{\mbox{y}}}</math>  e <math display="inline">\boldsymbol{\mbox{y}}_r-\overline{\boldsymbol{\mbox{y}}}</math>  debe ser obtuso.      

Este hecho se demuestra cuando, al realizar el producto escalar <math display="inline">{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t(\boldsymbol{\mbox{y}}_r-</math><math>\overline{\boldsymbol{\mbox{y}}})</math> , este es negativo. De esta manera se tiene que de la ecuación [[#eq0040|(8)]]  se puede escribir

<span id='eq0050'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{y}}_r-\overline{\boldsymbol{\mbox{y}}}=</math><math>(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}})-</math><math>{\alpha }_e\quad \mu \quad sgnd\quad (\boldsymbol{\mbox{y}}_{\nu }-</math><math>\overline{\boldsymbol{\mbox{y}}}).</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 9)
|}

Entonces, al premultiplicar la ecuación [[#eq0050|(9)]]  por <math display="inline">{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t</math> , se obtiene que

<span id='eq0055'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t(\boldsymbol{\mbox{y}}_r-</math><math>\overline{\boldsymbol{\mbox{y}}})=</math>
|-
|<math>{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t(\boldsymbol{\mbox{y}}_{\nu }-</math><math>\overline{\boldsymbol{\mbox{y}}})-{\alpha }_e\quad \mu \quad {\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^tsgnd\quad (\boldsymbol{\mbox{y}}_{\nu }-</math><math>\overline{\boldsymbol{\mbox{y}}}).</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 10)
|}

Simplificando la ecuación [[#eq0055|(10)]] , se tiene que

<span id='eq0060'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t(\boldsymbol{\mbox{y}}_r-</math><math>\overline{\boldsymbol{\mbox{y}}})=</math>
|-
|<math>\parallel \boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}{\parallel }^2-</math><math>{\alpha }_e\quad \mu \quad {\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^tsgnd\quad (\boldsymbol{\mbox{y}}_{\nu }-</math><math>\overline{\boldsymbol{\mbox{y}}}).</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 11)
|}

Si se denota a <math display="inline">{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t=</math><math>{\left(d_1,d_2,\ldots ,d_n\right)}^t</math>  se tiene que

<span id='eq0065'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^tsgnd\quad (\boldsymbol{\mbox{y}}_{\nu }-</math><math>\overline{\boldsymbol{\mbox{y}}})=\sum_{i=1}^n\vert d_i\vert .</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 12)
|}

Sustituyendo la ecuación [[#eq0065|(12)]]  en la ecuación [[#eq0060|(11)]]  se puede afirmar que

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t(\boldsymbol{\mbox{y}}_r-</math><math>\overline{\boldsymbol{\mbox{y}}})=\sum_{i=1}^n\left(d_i^2-\right. </math><math>\left. {\alpha }_e\quad \mu \quad \vert d_i\vert \right).</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 13)
|}

Dado que <math display="inline">\mu =\lceil \parallel \overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu }\parallel \rceil \geq \vert d_i\vert </math>  para todo ''i''  = 1, …, ''n'' , y que ''α''<sub>''e''</sub>  ≥ 2 se puede asegurar que <math display="inline">d_i^2-{\alpha }_e\quad \mu \quad \vert d_i\vert <0</math>  para todo ''i''  = 1, …, ''n  '' , y en consecuencia <math display="inline">{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t(\boldsymbol{\mbox{y}}_r-</math><math>\overline{\boldsymbol{\mbox{y}}})<0</math> . Solo cuando ''α''<sub>''e''</sub>  = 1 el producto escalar <math display="inline">{\left(\boldsymbol{\mbox{y}}_{\nu }-\overline{\boldsymbol{\mbox{y}}}\right)}^t(\boldsymbol{\mbox{y}}_r-</math><math>\overline{\boldsymbol{\mbox{y}}})</math>  puede ser igual a cero, lo cual provocaría el colapso del simplex entero. □

<span id='enun0100'></span>

==Lema 1.                     ==

Sean: un vector '''y'''  = (''y''<sub>1</sub> , …, ''y''<sub>''n''</sub> )<sup>''t''</sup>  de componentes enteros definido en el espacio euclidiano <math display="inline">{\mathbb{Z}}^n</math> , ''ξ  ''  un número entero positivo, y <math display="inline">{\tilde{\boldsymbol{\mbox{k}}}}_n=</math><math>{\left(k_1,\ldots ,k_n\right)}^t</math>  un vector cuyas componentes ''k''<sub>''i''</sub>  ∈ { −1, 0, 1}, para todo ''i''  = 1, …, ''n  '' . Entonces, <math display="inline">\boldsymbol{\mbox{y}}+\xi {\tilde{\boldsymbol{\mbox{k}}}}_n</math>  también pertenece al espacio euclidiano entero <math display="inline">{\mathbb{Z}}^n</math> , y su distancia a '''y'''  es de <math display="inline">\xi \sqrt{n-l}</math> , donde ''l''  ≤ ''n  ''  es el número de componentes nulas en el vector <math display="inline">{\tilde{\boldsymbol{\mbox{k}}}}_n</math> .

''Demostración.  ''  En virtud de que <math display="inline">\boldsymbol{\mbox{y}}={\left(y_1,\ldots ,y_n\right)}^t\in {\mathbb{Z}}^n</math> , entonces,

<span id='eq0075'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{y}}_k=\boldsymbol{\mbox{y}}+\xi {\tilde{\boldsymbol{\mbox{k}}}}_n=</math><math>{\left(y_1+\xi k_1,\ldots ,y_n+\xi k_n\right)}^t.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 14)
|}

Debido al hecho de que cada ''i'' -ésima componente de '''y'''<sub>''k''</sub>  indicada en la ecuación [[#eq0075|(14)]] , es decir, ''y''<sub>''i''</sub>  + ''ξk''<sub>''i''</sub>  para todo ''i''  = 1, …, ''n  ''  es un número entero, puede asegurarse que <math display="inline">\boldsymbol{\mbox{y}}_k\in {\mathbb{Z}}^n</math> . Por otra parte, la norma euclidiana del vector definido por la diferencia entre '''y'''<sub>''k''</sub>  e '''y'''  es

<span id='eq0080'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\parallel \boldsymbol{\mbox{y}}_k-\boldsymbol{\mbox{y}}\parallel =</math><math>\sqrt{\underset{n-l}{\underbrace{{\left(\pm \quad \xi \right)}^2+\cdots +{\left(\pm \quad \xi \right)}^2}}+\underset{l}{\underbrace{0^2+\cdots +0^2}}},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 15)
|}

en virtud de que para todo ''i''  = 1, …, ''n''  los números escalares ''k''<sub>''i''</sub>  ∈ { −1, 0, 1}. La ecuación [[#eq0080|(15)]]  determina la distancia de '''y'''  a un punto vecino '''y'''<sub>''k''</sub> . □

<span id='enun0105'></span>

==Proposición 2 Punto expandido entero.                     ==

Sea <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  un ''q'' -ésimo simplex entero jerarquizado de acuerdo con ''f'' ('''v''' ), definido exclusivamente en el espacio euclidiano entero <math display="inline">{\mathbb{Z}}^n</math> . Si el punto expandido entero '''y'''<sub>''e''</sub>  es definido por

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{y}}_e=\boldsymbol{\mbox{y}}_r+{\beta }_e\quad \mu \quad sgnd\quad (\overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu }),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 16)
|}

donde <math display="inline">{\beta }_e\in \mathbb{N}(>0)</math>  representa el coeficiente de expansión entero e '''y'''<sub>''r''</sub>  es determinado por la ecuación [[#eq0040|(8)]] . Entonces, se cumple que:
* '''y'''<sub>''e''</sub>  es un punto en el espacio de los números enteros <math display="inline">{\mathbb{Z}}^n</math> ;                                  
* con respecto al vértice '''y'''<sub>''ν''</sub> , '''y'''<sub>''e''</sub>  nunca está más cerca de '''y'''<sub>''ν''</sub> , de lo que está '''y'''<sub>''r''</sub>  de '''y'''<sub>''ν''</sub> .          

''Demostración.''  Parte I) Represente '''y'''<sub>''r''</sub>  = (''y''<sub>1,''r''</sub> , …, ''y''<sub>''n'' ,''r''</sub> )<sup>''t''</sup>  un punto reflejado producto de la operación definida por la ecuación [[#eq0040|(8)]] . Entonces

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{y}}_e=\boldsymbol{\mbox{y}}_r+{\beta }_e\quad \mu \quad {\left(k_1,\ldots ,k_n\right)}^t,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 17)
|}

donde ''k''<sub>''i''</sub>  ∈ { −1, 0, 1} para todo ''i''  = 1, …, ''n''  debido a la función signo discreta. En consecuencia, cada una de las componentes de '''y'''<sub>''e''</sub>  es entera, debido al hecho de que ''y''<sub>''i'' ,''r''</sub> , ''β''<sub>''e''</sub>  y ''k''<sub>''i''</sub>  son todos números enteros, para cada ''i''  = 1, …, ''n''  su operación arroja un número entero.      

Parte II) Claramente el vector <math display="inline">sgnd\quad (\overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu })</math>  tiene la misma orientación del vector <math display="inline">\overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu }</math> , y este último va desde el punto '''y'''<sub>''ν''</sub>  hasta el punto <math display="inline">\overline{\boldsymbol{\mbox{y}}}</math> . Como resultado de esto, el vector <math display="inline">sgnd\quad (\overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu })</math>  tiene como orientación un vector que va desde el punto <math display="inline">\overline{\boldsymbol{\mbox{y}}}</math>  hasta un punto que se encuentra en la región discreta entre <math display="inline">\overline{\boldsymbol{\mbox{y}}}</math>  e '''y'''<sub>''r''</sub> , y la norma euclidiana de <math display="inline">{\beta }_e\quad sgnd\quad (\overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu })</math> , según el [[#enun0100|lema 1]] , es <math display="inline">{\beta }_e\quad \sqrt{n-l}\not<0</math> , lo cual, al sumarse con '''y'''<sub>''r''</sub> , se ubica en un punto que no está más cerca de lo que está '''y'''<sub>''r''</sub>  de '''y'''<sub>''ν''</sub> .

==Proposición 3 Punto contraído entero.                     ==

Sea un ''q  '' -ésimo simplex entero jerarquizado <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  de acuerdo a ''f'' ('''v''' ), definido únicamente en el espacio euclidiano entero <math display="inline">{\mathbb{Z}}^n</math> . Si el punto contraído entero '''y'''<sub>''c''</sub>  es determinado por

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{y}}_c=\boldsymbol{\mbox{y}}_r-{\gamma }_e\quad \mu \quad sgnd\quad (\overline{\boldsymbol{\mbox{y}}}-</math><math>\boldsymbol{\mbox{y}}_{\nu }),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 18)
|}

donde <math display="inline">{\gamma }_e\in \mathbb{N}(>0)</math>  es el llamado coeficiente de contracción entero, '''y'''<sub>''r''</sub>  es determinado por [[#eq0040|(8)]] . Entonces, se tiene que:
* '''y'''<sub>''c''</sub>  es un punto en el espacio de los números enteros <math display="inline">{\mathbb{Z}}^n</math> ;                                  
* con respecto al vértice '''y'''<sub>''ν''</sub>  está más cerca, de lo que está '''y'''<sub>''r''</sub>  de '''y'''<sub>''ν''</sub> .          

''Demostración.''  Aplicando el mismo método empleado en la demostración de la  [[#enun0105|proposición 2]] , se demuestra la parte I) y parte II) de la proposición.

<span id='enun0115'></span>

==Proposición 4 Simplex entero encogido.                     ==

Sea un ''q  '' -ésimo simplex entero jerarquizado <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  de acuerdo a ''f'' ('''v''' ), definido solamente en el espacio euclidiano entero <math display="inline">{\mathbb{Z}}^n</math> . Si 0 < ''δ''<sub>''e''</sub>  < 1 es el coeficiente de encogimiento entero, y además,

<span id='eq0100'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{D}}^{\left[q\right]}=[\boldsymbol{\mbox{d}}_1^{\left[q\right]}\quad ;\boldsymbol{\mbox{d}}_2^{\left[q\right]}\quad ;\cdots \quad ;\boldsymbol{\mbox{d}}_{\nu -1}^{\left[q\right]}]=</math><math>\lceil {\delta }_eA_1^{\left[q\right]}\rceil ,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 19)
|}

donde <math display="inline">\boldsymbol{\mbox{A}}_1^{\left[q\right]}</math>  es la matriz de arista de <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math> , y el operador ⌈· ⌉ aproxima cada elemento de la matriz <math display="inline">\boldsymbol{\mbox{A}}_1^{\left[q\right]}</math>  al próximo entero por encima. Entonces, se cumple que:
* la matriz formada por

<span id='eq0105'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\tilde{\boldsymbol{\mbox{S}}}}_{\nu }^{\left[q+1\right]}=</math><math>[\boldsymbol{\mbox{y}}_1^{\left[q\right]};\quad \boldsymbol{\mbox{y}}_1^{\left[q\right]}+</math><math>\boldsymbol{\mbox{d}}_1^{\left[q\right]}\quad ;\cdots \quad ;\boldsymbol{\mbox{y}}_1^{\left[q\right]}+</math><math>\boldsymbol{\mbox{d}}_{\nu -1}^{\left[q\right]}];</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 20)
|}

es un simplex compuesto de elementos enteros.                                  
* las aristas del simplex <math display="inline">{\tilde{\boldsymbol{\mbox{Y}}}}_{\nu }^{\left[q+1\right]}</math>  con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math>  nunca son mayores a las aristas del simplex <math display="inline">\boldsymbol{\mbox{Y}}_{\nu }^{\left[q\right]}</math>  con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math> , en consecuencia el nuevo simplex es igual o menor en tamaño.          

''Demostración.''  Parte I) Del análisis de la operación definida por la ecuación  [[#eq0100|(19)]]  se puede afirmar que esta arroja como resultado una matriz de dimensión ''n''  × ''n''  compuesta de elementos enteros, en virtud de que el operador ⌈· ⌉ aproxima cada uno de los elementos de la matrix al próximo entero por encima.      

Parte II) Para demostrar esta afirmación, se tiene primero que, de acuerdo con la ecuación [[#eq0105|(20)]] , sendas primeras columnas de <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  y <math display="inline">{\tilde{\boldsymbol{\mbox{S}}}}_{\nu }^{\left[q+1\right]}</math>  son iguales al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math> . Ahora se deberán estudiar las normas de las aristas de sendos símpleces con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math> .      

Se tiene que la matriz de aristas del simplex <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math>  está definida por

<span id='eq0110'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{A}}_1^{\left[q\right]}=[\boldsymbol{\mbox{a}}_1^{\left[q\right]}\quad ;\boldsymbol{\mbox{a}}_2^{\left[q\right]}\quad ;\cdots \quad ;\boldsymbol{\mbox{a}}_n^{\left[q\right]}],</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 21)
|}

donde <math display="inline">\boldsymbol{\mbox{a}}_i^{\left[q\right]}=</math><math>{\left(a_{1i},a_{2i},\ldots ,a_{ni}\right)}^t</math>  para todo ''i''  = 1, …, ''n'' .      

Por otra parte, la matriz de arista del simplex <math display="inline">{\tilde{\boldsymbol{\mbox{S}}}}_{\nu }^{\left[q+1\right]}</math>  con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q+1\right]}</math> , de acuerdo a la ecuación [[#eq0105|(20)]]  viene dada por

<span id='eq0115'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{A}}_1^{\left[q+1\right]}=[\lceil {\delta }_e\boldsymbol{\mbox{a}}_1^{\left[q\right]}\rceil \quad ;\lceil {\delta }_e\boldsymbol{\mbox{a}}_2^{\left[q\right]}\rceil ;\cdots \quad ;\lceil {\delta }_e\boldsymbol{\mbox{a}}_n^{\left[q\right]}\rceil ].</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 22)
|}

Ahora, al comparar la norma euclidiana de cada arista de ambas matrices, es decir, las matrices dadas por las ecuaciones [[#eq0110|(21)]]  y [[#eq0115|(22)]] , se tiene que para todo ''j''  = 1, …, ''n''

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>+\sqrt{\sum_{i=1}^n\lceil {\delta }_ea_{ij}{\rceil }^2}\leq +</math><math>\sqrt{\sum_{i=1}^na_{ij}^2}.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 23)
|}

Debido al hecho de que 0 < ''δ''<sub>''e''</sub>  < 1, para cada ''i''  = 1, …, ''n''  y ''j''  = 1, …, ''n  '' , siempre se cumple que <math display="inline">\lceil {\delta }_ea_{ij}{\rceil }^2\leq a_{ij}^2</math> , lo cual permite asegurar que para todo ''i''  = 1, …, ''n'' , la norma euclidiana de cada ''j  '' -ésima arista del simplex <math display="inline">{\tilde{\boldsymbol{\mbox{S}}}}_{\nu }^{\left[q+1\right]}</math>  con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math>  nunca es mayor a la norma euclidiana para esa misma ''j  '' -ésima arista del simplex <math display="inline">\boldsymbol{\mbox{S}}_{\nu }^{\left[q\right]}</math>  con respecto al vértice <math display="inline">\boldsymbol{\mbox{y}}_1^{\left[q\right]}</math> . □

<span id='enun0120'></span>

==Observación 4.                     ==

La construcción de un nuevo simplex entero encogido mediante la [[#enun0115|proposición 4]]  puede producir un simplex entero colapsado <math display="inline">{\tilde{\boldsymbol{\mbox{S}}}}_{\nu }^{\left[q+1\right]}</math> . Como ejemplo, nótese que la [[#enun0115|proposición 4]]  aplicada al simplex entero <math display="inline">\boldsymbol{\mbox{S}}_3^{\left[q\right]}=</math><math>[\boldsymbol{\mbox{y}}_1\quad ;\boldsymbol{\mbox{y}}_2\quad ;\boldsymbol{\mbox{y}}_3]=</math><math>\left[\begin{array}{ccc}
1 & 2 & 3\\
3 & 1 & 1\\
 &  & 
\end{array}\right]</math> , con una valor de ''δ''<sub>''e''</sub>  = 1/2, arroja como resultado el simplex entero colapsado <math display="inline">{\tilde{\boldsymbol{\mbox{S}}}}_3^{\left[q+1\right]}=</math><math>\left[\begin{array}{ccc}
1 & 2 & 2\\
3 & 2 & 2\\
 &  & 
\end{array}\right]</math> .

La [[#fig0005|figura 1]]  muestra las operaciones de reflexión, expansión, contracción y construcción de un simplex encogido que pueden ser ejecutadas sobre un ''q  '' -ésimo simplex entero jerarquizado <math display="inline">\boldsymbol{\mbox{S}}_3^{\left[q\right]}=</math><math>[\boldsymbol{\mbox{y}}_1\quad ;\boldsymbol{\mbox{y}}_2\quad ;\boldsymbol{\mbox{y}}_3]=</math><math>\left[\begin{array}{ccc}
1 & 3 & 1\\
5 & 2 & 1\\
 &  & 
\end{array}\right]</math> , definido en el espacio <math display="inline">{\mathbb{Z}}^2</math> , cuando los coeficientes tomaron los valores ''α''<sub>''e''</sub>  = 2, y ''β''<sub>''e''</sub>  = ''γ''<sub>''e''</sub>  = 1, y ''δ''<sub>''e''</sub>  = 1/2.

<span id='fig0005'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_404207060-1-s2.0-S0213131513000345-gr1.jpg|center|376px|Operaciones sobre un simplex entero jerarquizado.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 1.

Operaciones sobre un simplex entero jerarquizado.

</span>
|}

Obsérvese que el punto <math display="inline">\overline{\boldsymbol{\mbox{y}}}=</math><math>{\left(2,\frac{7}{2}\right)}^t</math>  claramente no pertenece a los números <math display="inline">{\mathbb{Z}}^2</math> , sin embargo, el resto de los puntos arrojados por las distintas operaciones sobre el simplex entero <math display="inline">\boldsymbol{\mbox{S}}_3^{\left[q\right]}</math>  sí pertenecen al espacio de los números <math display="inline">{\mathbb{Z}}^2</math> .      

Por último, nótese además que en la [[#fig0005|figura 1]]  se muestra el simplex entero encogido, el cual es denotado por <math display="inline">{\tilde{\boldsymbol{\mbox{S}}}}_3^{\left[q+1\right]}=</math><math>[\boldsymbol{\mbox{y}}_1\quad ;{\acute{\boldsymbol{\mbox{y}}}}_2\quad ;{\acute{\boldsymbol{\mbox{y}}}}_3]=</math><math>\left[\begin{array}{ccc}
1 & 2 & 1\\
5 & 4 & 3\\
 &  & 
\end{array}\right]</math> , para lo cual se empleó la [[#enun0115|proposición 4]] , la cual podría arrojar un simplex entero colapsado según la [[#enun0120|observación 4]] . No obstante, el algoritmo propuesto cuenta con etapas, las cuales deben iniciar con un nuevo simplex entero mixto.      

<span id='sec0025'></span>

==5. Procedimientos==

El ASEM se basa, entre otros pasos, en 2 procedimientos principales, que se presentan en esta sección. Es importante señalar que el superíndice [''q'' ] es en algunas oportunidades suprimido con el propósito de simplificar la notación, sin que pierda su significado.      

El procedimiento mostrado en la [[#fig0010|figura 2]] , llamado Procedimiento de Construcción de un Simplex Entero Mixto (CSEM), permite construir un simplex entero mixto basándose, entre otros parámetros, en los valores de Δ<sub>''r''</sub>  y Δ<sub>''e''</sub>  que son definidos inicialmente por el usuario del ASEM a los efectos de tener control de los parámetros del algoritmo. Nótese que en la figura se emplea la notación '''0'''<sub>2''n''</sub> , la cual representa un vector nulo de dimensión 2''n'' , es decir, un vector de elementos todos nulos.

<span id='fig0010'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_404207060-1-s2.0-S0213131513000345-gr2.jpg|center|320px|Construcción de un Simplex Entero Mixto.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 2.

Construcción de un Simplex Entero Mixto.

</span>
|}

Otro procedimiento empleado en el ASEM es el Mejoramiento de un Simplex Entero Mixto (MSEM), el cual tiene como propósito redefinir el simplex entero mixto inicial a través de la orientación de sus aristas. Este procedimiento se presenta en la [[#fig0015|figura 3]]  con el objetivo de buscar la mejor orientación del simplex construido por el procedimiento CSEM. Los experimentos numéricos realizados por Brea [[#bib0015|[3]]]  probaron la eficacia del MSEM.

<span id='fig0015'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_404207060-1-s2.0-S0213131513000345-gr3.jpg|center|331px|Mejoramiento de un Simplex Entero Mixto.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 3.

Mejoramiento de un Simplex Entero Mixto.

</span>
|}

<span id='sec0030'></span>
==6. El algoritmo==

Para iniciar la búsqueda de un mínimo local del [[#enun0005|problema 1]] , el ASEM, mostrado en la [[#fig0020|figura 4]] , arranca con sus parámetros definidos por el usuario, el número de la etapa ''s''  = 1 y un punto de inicio <math display="inline">\boldsymbol{\mbox{v}}_0\in {\mathbb{R}}^w\times {\mathbb{Z}}^w</math> , donde <math display="inline">w=max(n,m)</math> . Con el punto de inicio '''v'''<sub>0</sub>  y los parámetros, el algoritmo construye un simplex entero mixto mediante la ejecución del procedimiento CSEM. Si el ASEM es ejecutado bajo la modalidad de selección del mejor simplex entero mixto construido, el cual se indica con ''m''  = 1, el algoritmo ejecuta el procedimiento MSEM e inicia el contador de iteraciones ''q''  en cero.

<span id='fig0020'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_404207060-1-s2.0-S0213131513000345-gr4.jpg|center|565px|Diagrama del Algoritmo Simplex Entero Mixto.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 4.

Diagrama del Algoritmo Simplex Entero Mixto.

</span>
|}

Una vez construido un simplex entero mixto inicial, el ASEM evalúa la función objetivo para cada uno de los vértices del simplex actual de los que no se conoce su valor, mediante el paso (3) Evaluación . Seguidamente, se determina el orden de precedencia de sus vértices a través del paso (4) Ordenamiento , lo cual permite determinar el punto de prueba reflejado mediante la ejecución del paso (5) Reflexión . Al final de este último paso, el ASEM selecciona uno de los 4 posibles casos, en función de la precedencia del punto de prueba reflejado, denotado por '''v'''<sub>''r''</sub>  con relación a: '''v'''<sub>1</sub> , '''v'''<sub>''ν'' −1                            </sub>  y '''v'''<sub>''ν''</sub> . Si '''v'''<sub>''r''</sub>  ≺ '''v'''<sub>1</sub> , es decir, si ''f'' ('''v'''<sub>''r''</sub> ) < ''f'' ('''v'''<sub>1</sub> ), el ASEM ejecuta una operación de expansión, para luego decidir si el vértice '''v'''<sub>''ν''</sub>  será sustituido por el punto de prueba '''v'''<sub>''r''</sub>  o '''v'''<sub>''e''</sub> . Por otra parte, si '''v'''<sub>1</sub>  ⪯ '''v'''<sub>''r''</sub>  ⪯ '''v'''<sub>''ν'' −1                            </sub> , el ASEM reemplaza el vértice '''v'''<sub>''ν''</sub>  por el punto de prueba '''v'''<sub>''r''</sub> . No obstante, si '''v'''<sub>''ν'' −1                            </sub>  ≺ '''v'''<sub>''r''</sub>  ≺ '''v'''<sub>''ν''</sub>  el algoritmo sustituye el vértice '''v'''<sub>''ν''</sub>  del simplex actual por el punto de prueba '''v'''<sub>''r''</sub>  para luego ejecutar el paso (7) Contracción . Finalmente, en el caso de que '''v'''<sub>''ν''</sub>  ⪯ '''v'''<sub>''r''</sub> , el ASEM ejecuta una contracción del simplex actual por medio del paso (7) Contracción .      

Después de que el ASEM tome una de las opciones de acuerdo con la precedencia de '''v'''<sub>''r''</sub>  con relación a '''v'''<sub>1</sub> , '''v'''<sub>''ν'' −1                            </sub>  y '''v'''<sub>''ν''</sub> , el algoritmo determina la cualidad del nuevo simplex de acuerdo al diámetro del simplex '''X''' , el cual es denotado por ℓ<sub>''x''</sub> . Este valor de ℓ<sub>''x''</sub>  es comparado con un criterio de parada <math display="inline">{\kappa }_x^{\left[s\right]}</math>  de la ''s  '' -ésima etapa, el cual adquiere distintos valores según la etapa que se esté ejecutando. Si <math display="inline">{l}_x\geq {\kappa }_x^{\left[s\right]}</math> , el ASEM ejecuta una nueva iteración. En caso contrario, el algoritmo detiene las iteraciones dentro de la etapa actual para decidir si continúa la búsqueda o no, de acuerdo con el valor de <math display="inline">\Delta =\vert \vert {\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}}}^{\left[s\right]}-</math><math>{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}}}^{\left[s-1\right]}\vert \vert </math> . Si Δ ≥ ''ɛ  '' , el ASEM asigna a <math display="inline">{\kappa }_x^{\left[s\right]}=</math><math>\varphi {\kappa }_x^{\left[s-1\right]}</math> ; <math display="inline">{\Delta }_r^{\left[s\right]}=</math><math>\rho {\Delta }_r^{\left[s-1\right]}</math>  y con '''v'''<sub>1</sub>  va al paso (2) Construcción de un Simplex Entero Mixto  a fin de rearrancar una nueva etapa. En caso contrario, el ASEM reporta a <math display="inline">{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}}}^{\left[s\right]}=</math><math>{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}}}_1</math>  como el mínimo.      

<span id='sec0035'></span>
==7. Flexibilidad del Algoritmo Ximplex Entero Mixto==

Una característica particular que tiene el ASEM es su flexibilidad con respecto a problemas en donde el número de componentes reales y enteras son diferentes, es decir, problemas del tipo <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^m</math> , donde <math display="inline">n\in {\mathbb{N}}_+</math>  y <math display="inline">m\in {\mathbb{N}}_+</math>  pueden ser diferentes.      

Para tratar este tipo de problemas se debe convertir el problema de <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^m</math>  a <math display="inline">{\mathbb{R}}^{max(n,m)}\times {\mathbb{Z}}^{max(n,m)}</math> , para así permitir definir un simplex entero mixto con igual número de componentes reales y enteras.      

Es importante destacar que la función objetivo definida en el espacio <math display="inline">{\mathbb{R}}^n\times {\mathbb{Z}}^m</math>  es tratada con su número original de variables. No obstante, el simplex es definido con max(''n'' , ''m'' ) + 1 vértices enteros mixtos todos definidos en espacio euclidiano entero mixto <math display="inline">{\mathbb{R}}^{max(n,m)}\times {\mathbb{Z}}^{max(n,m)}</math> .      

Esta generalización de los problemas enteros mixtos se muestra a través de las funciones de pruebas presentadas en el apéndice [[#sec0060|A]] .      

<span id='sec0040'></span>
==8. Sintonización del Algoritmo Ximplex Entero Mixto==

Para el proceso de sintonización del ASEM se consideraron los parámetros obtenidos en la sintonización del Algoritmo Simplex Entero Relajado (ASER) desarrollado por Brea [[#bib0015|[3]]] , además de los valores recomendados y apoyados en las investigaciones de Brea sobre el MSNM cuando el problema de optimización no lineal está bajo restricciones lineales [[#bib0025|[5]]]  ([[#tbl0005|tabla 1]] ).      

En consecuencia, los parámetros ''α''<sub>''r''</sub> , ''β''<sub>''r''</sub> , ''γ''<sub>''r''</sub>  y ''δ''<sub>''r''</sub>  indicados en la [[#tbl0005|tabla 1]]  fueron ajustados de acuerdo a los valores recomendados en Nelder y Mead [[#bib0010|[2]]] , mientras los parámetros Δ<sub>''e''</sub> , ''α''<sub>''e''</sub> , ''β''<sub>''e''</sub> , ''γ''<sub>''e''</sub>  y ''δ''<sub>''e''</sub>  se fijaron de acuerdo con los valores obtenidos en el proceso de sintonización del ASER (véase Brea [[#bib0015|[3, cap. 5]]] ). Además, los valores de los parámetros ''φ  ''  y <math display="inline">k_x^{\left[1\right]}</math>  mostrados en la [[#tbl0005|tabla 1]]  se obtuvieron mediante pruebas preliminares, las cuales permitieron mostrar un desempeño satisfactorio del ASEM. Finalmente, el parámetro ''ɛ''  representa el criterio de comparación de los resultados obtenidos en 2 etapas sucesivas y, en consecuencia, debe ser fijado por el usuario dependiendo del máximo error que espera obtener en la comparación de 2 etapas sucesivas.

<span id='tbl0005'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabla 1.

Resultado de la sintonización del ASEM

|-

| Parámetro
| 
| 
| ''ɛ''
| ''φ''
| <math display="inline">k_x^{\left[1\right]}</math>
|-

| Valores
| 
| 
| 0,1
| 0,3
| 1
|-

| Parámetro
| 
| ''α''<sub>''r''</sub>
| ''β''<sub>''r''</sub>
| ''γ''<sub>''r''</sub>
| ''δ''<sub>''r''</sub>
|-

| Valores
| 
| 1
| 2
| 0,5
| 0,5
|-

| Parámetro
| Δ<sub>''e''</sub>
| ''α''<sub>''e''</sub>
| ''β''<sub>''e''</sub>
| ''γ''<sub>''e''</sub>
| ''δ''<sub>''e''</sub>
|-

| Valores
| 1
| 2
| 2
| 1
| 0,4
|}

Estas consideraciones permitieron simplificar significativamente la experimentación para la obtención de los valores apropiados de ''ρ  ''  y <math display="inline">{\Delta }_r^{\left[1\right]}</math> . Es importante señalar que durante la realización de los experimentos, el número máximo de iteraciones que se le permitió ejecutar al ASEM por etapa fue de 15.000, debido a que el ASEM podía quedar en un proceso iterativo, sin encontrar un avance significativo en la minimización de la función objetivo.      

La [[#tbl0010|tabla 2]]  muestra el número de evaluaciones (NE) ejecutadas sobre la función objetivo FC<math display="inline">\mathbb{R}\mathbb{Z}</math> [''d''  + ''d'' ] con ''a''<sub>''i''</sub>  = ''b''<sub>''i''</sub>  = 1 para todo ''i''  = 1, …, ''d'' , de acuerdo con los valores de ''ρ  ''  y <math display="inline">{\Delta }_r^{\left[1\right]}</math>  indicados en la tabla.

<span id='tbl0010'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabla 2.

Resultados obtenidos del experimento a 3 niveles sobre la función FC<math display="inline">\mathbb{R}\mathbb{Z}</math> [''d''  + ''d'' ]                  

|-

! ''d''  = 2                                                    
! <math display="inline">\rho \setminus {\Delta }_r^{\left[1\right]}</math>
! 0,5
! 1,0
! 1,5
|-

| 
| 0,7
| 175 (288e-6)
| 127 (163e-6)
| 133 (169e-6)
|-

| 
| 0,8
| 192 (10e-6)
| 115 (14e-6)
| 146 (19e-6)
|-

| 
| 0,9
| 107 (214e-6)
| 165 (71e-6)
| 169 (2)
|-

| ''d''  = 4                                                    
| <math display="inline">\rho \setminus {\Delta }_r^{\left[1\right]}</math>
| 0,5
| 1,0
| 1,5
|-

| 
| 0,7
| 529 (5e-6)
| 457 (1)
| 425 (100e-6)
|-

| 
| 0,8
| 727 (2e-6)
| 435 (174e-6)
| 544(107e-6)
|-

| 
| 0,9
| 883 (15e-6)
| 379 (174e-6)
| 544 (3)
|-

| ''d''  = 8                                                    
| <math display="inline">\rho \setminus {\Delta }_r^{\left[1\right]}</math>
| 0,5
| 1,0
| 1,5
|-

| 
| 0,7
| 12498 (0)
| 3123 (1e-6)
| 5452 (3e-6)
|-

| 
| 0,8
| 14192 (0)
| 7651 (0)
| 4835 (2e-6)
|-

| 
| 0,9
| 5686 (0)
| 2761 (2)
| 3642 (1)
|-

| ''d''  = 14                                                    
| <math display="inline">\rho \setminus {\Delta }_r^{\left[1\right]}</math>
| 0,5
| 1,0
| 1,5
|-

| 
| 0,7
| 57251 (90)
| 96558 (0)
| 45253 (1)
|-

| 
| 0,8
| 17107 (1e-6)
| 80936 (0)
| 38234 (0)
|-

| 
| 0,9
| 5686 (0)
| 2761 (2)
| 3642 (1)
|-

| ''d''  = 18                                                    
| <math display="inline">\rho \setminus {\Delta }_r^{\left[1\right]}</math>
| 0,5
| 1,0
| 1,5
|-

| 
| 0,7
| 35378 (0)
| 270236 (596e-6)
| 40550 (3)
|-

| 
| 0,8
| 188136 (0)
| 72880 (0)
| 58245 (2)
|-

| 
| 0,9
| 83971 (0)
| 146017 (3e-6)
| 33041 (5)
|-

| ''d''  = 20                                                    
| <math display="inline">\rho \setminus {\Delta }_r^{\left[1\right]}</math>
| 0,5
| 1,0
| 1,5
|-

| 
| 0,7
| 405221 (3952e-6)
| 253326 (1009)
| 71239 (2)
|-

| 
| 0,8
| 281525 (0)
| 161640 (0)
| 57712 (1)
|-

| 
| 0,9
| 95165 (0)
| 160039 (0)
| 68075 (4)
|}

Además, en la [[#tbl0010|tabla 2]]  se muestra entre paréntesis el mejor valor de la función objetivo para el NE indicado, siendo en todos los casos el punto óptimo en <math display="inline">\overset{\mbox{ˆ}}{\boldsymbol{\mbox{x}}}=</math><math>\boldsymbol{0}_d</math>  e <math display="inline">\overset{\mbox{ˆ}}{\boldsymbol{\mbox{y}}}=</math><math>\boldsymbol{0}_d</math> , con un valor de la función objetivo de <math display="inline">f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})=</math><math>0</math> .      

Como ejemplo de lectura de la [[#tbl0010|tabla 2]]  se tiene que para ''d''  = 8, ''ρ''  = 0, 8 y <math display="inline">{\Delta }_r^{\left[1\right]}=</math><math>1</math>  el NE fue igual a 7.651 y el valor mínimo obtenido de la función objetivo fue entre 0 y un valor menor a 10<sup>−6</sup> , en particular fue igual a 7.67546044785799 × 10<sup>−8</sup> .      

De la [[#tbl0010|tabla 2]]  se puede apreciar que unos valores adecuados para ''ρ  ''  y <math display="inline">{\Delta }_r^{\left[1\right]}</math>  son 0,8 y 1 respectivamente, debido al hecho de que para esos valores el ASEM identificó el punto mínimo de la función objetivo para todos los casos de dimensión entera mixta considerados.      

Es importante destacar que estos valores de sintonización del ASEM son referenciales y que no deben tomarse como valores que deben satisfacer la totalidad de los casos de estudio, puesto que fueron obtenidos de casos particulares.

<span id='sec0045'></span>
==9. Experimentos numéricos==

La [[#tbl0015|tabla 3]]  muestra los valores mínimos obtenidos por el ASEM, el mínimo teórico y el orden <math display="inline">O</math>  de magnitud del error ''e'' , definido como

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>e=\left|\tilde{f}(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})-\right. </math><math>\left. f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})\right|,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 24)
|}

donde <math display="inline">\tilde{f}(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})</math>  representa el valor mínimo determinado por el ASEM <math display="inline">f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})</math>  representa el valor mínimo teórico para cada una de las funciones de prueba definidas en el Apéndice [[#sec0060|A]] .

<span id='tbl0015'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabla 3.

Experimentos numéricos del ASEM sobre diversas funciones objetivo

|-

! Problema
! colspan="3" | Condiciones
! colspan="3" | Mínimo
! colspan="2" | Total
|-

! 
! ''ρ''
! '''x'''<sub>0</sub>
! '''y'''<sub>0</sub>
! ASEM
! Teórico
! <math display="inline">O(e)</math>
! NE
! NI
|-

| FC<math display="inline">\mathbb{R}\mathbb{Z}[5+10]</math>
| 0, 85
| 10 · '''1'''<sub>5</sub>
| 10 · '''1'''<sub>10</sub>
| 3, 85399595970641 × 10<sup>−17</sup>
| 0
| 10<sup>−16</sup>
| 2.687
| 1.525
|-

| FC<math display="inline">\mathbb{R}\mathbb{Z}[10+5]</math>
| 0, 85
| 10 · '''1'''<sub>10</sub>
| 10 · '''1'''<sub>5</sub>
| 3, 19278384117684 × 10<sup>−7</sup>
| 0
| 10<sup>−6</sup>
| 22.718
| 14.854
|-

| FC<math display="inline">\mathbb{R}\mathbb{Z}[10+10]</math>
| 0, 85
| 10 · '''1'''<sub>10</sub>
| 10 · '''1'''<sub>10</sub>
| 3, 2938056148257 × 10<sup>−11</sup>
| 0
| 10<sup>−10</sup>
| 27.550
| 18.956
|-

| FC<math display="inline">\mathbb{R}\mathbb{Z}[20+20]</math>
| 0, 85
| 10 · '''1'''<sub>20</sub>
| 10 · '''1'''<sub>20</sub>
| 1, 42023676961727 × 10<sup>−9</sup>
| 0
| 10<sup>−8</sup>
| 156.785
| 121.154
|-

| FCD<math display="inline">\mathbb{R}\mathbb{Z}[5+10]</math>
| 0, 85
| 10 · '''1'''<sub>5</sub>
| 10 · '''1'''<sub>10</sub>
| 0,777777777837533
| <math display="inline">\frac{7}{9}</math>
| 10<sup>−9</sup>
| 15.428
| 10.553
|-

| FCD<math display="inline">\mathbb{R}\mathbb{Z}[10+5]</math>
| 0, 85
| 10 · '''1'''<sub>10</sub>
| 10 · '''1'''<sub>5</sub>
| 0,444444465663501
| <math display="inline">\frac{4}{9}</math>
| 10<sup>−7</sup>
| 5.511
| 3.551
|-

| FCD<math display="inline">\mathbb{R}\mathbb{Z}[10+10]</math>
| 0, 85
| 10 · '''1'''<sub>10</sub>
| 10 · '''1'''<sub>10</sub>
| 0,777777777837533
| <math display="inline">\frac{7}{9}</math>
| 10<sup>−9</sup>
| 15.428
| 10.553
|-

| FCD<math display="inline">\mathbb{R}\mathbb{Z}[20+20]</math>
| 0, 85
| 10 · '''1'''<sub>20</sub>
| 10 · '''1'''<sub>20</sub>
| 1, 55555555557834
| <math display="inline">\frac{14}{9}</math>
| 10<sup>−10</sup>
| 68.884
| 53.399
|-

| FGM<math display="inline">\mathbb{R}\mathbb{Z}[5+10]</math>
| 0, 85
| 10 · '''1'''<sub>5</sub>
| 10 · '''1'''<sub>10</sub>
| 3, 1137516964225
| -
| -
| 1.639
| 863
|-

| FGM<math display="inline">\mathbb{R}\mathbb{Z}[10+5]</math>
| 0, 85
| 10 · '''1'''<sub>10</sub>
| 10 · '''1'''<sub>5</sub>
| −0, 9998056161237
| -
| -
| 3.036
| 1.942
|-

| FGM<math display="inline">\mathbb{R}\mathbb{Z}[10+10]</math>
| 0, 85
| 10 · '''1'''<sub>10</sub>
| 10 · '''1'''<sub>10</sub>
| 3, 91853917817428
| -
| -
| 2.982
| 1.703
|-

| FGM<math display="inline">\mathbb{R}\mathbb{Z}[20+20]</math>
| 0, 85
| 10 · '''1'''<sub>20</sub>
| 10 · '''1'''<sub>20</sub>
| 0, 995491554751617
| -
| -
| 55.667
| 43.607
|-

| FPM<math display="inline">\mathbb{R}\mathbb{Z}[5+5]</math>
| 0, 85
| 10 · '''1'''<sub>5</sub>
| 10 · '''1'''<sub>5</sub>
| 1, 91104307081247 × 10<sup>−5</sup>
| 0
| 10<sup>−4</sup>
| 1.276
| 371
|-

| FPM<math display="inline">\mathbb{R}\mathbb{Z}[10+10]</math>
| 0, 85
| 10 · '''1'''<sub>10</sub>
| 10 · '''1'''<sub>10</sub>
| 0, 00596898501698284
| 0
| 10<sup>−2</sup>
| 2.172
| 927
|-

| FPM<math display="inline">\mathbb{R}\mathbb{Z}[20+20]</math>
| 0, 85
| 10 · '''1'''<sub>20</sub>
| 10 · '''1'''<sub>20</sub>
| 3, 07 × 10<sup>−5</sup>
| 0
| 10<sup>−4</sup>
| 27.897
| 17.457
|-

| FR<math display="inline">\mathbb{R}\mathbb{Z}[10+10]</math>
| 0, 99
| 5 · '''1'''<sub>10</sub>
| 5 · '''1'''<sub>10</sub>
| 6, 1531 × 10<sup>−13</sup>
| 0
| 10<sup>−12</sup>
| 10.213
| 6.677
|-

| FR<math display="inline">\mathbb{R}\mathbb{Z}[20+20]</math>
| 0, 99
| 5 · '''1'''<sub>20</sub>
| 5 · '''1'''<sub>20</sub>
| 4, 75 × 10<sup>−6</sup>
| 0
| 10<sup>−5</sup>
| 48.777
| 36.038
|-

| FAD<math display="inline">\mathbb{R}\mathbb{Z}[2+1]</math>
| 0, 85
| 10 · '''1'''<sub>2</sub>
| 10 · '''1'''<sub>1</sub>
| −13, 9942
| -14
| 10<sup>−2</sup>
| 272
| 67
|}

Adicionalmente, la [[#tbl0015|tabla 3]]  presenta el número de evaluaciones (NE) de la función objetivo y el número de iteraciones (NI) totales ejecutadas por el ASEM, para las funciones de prueba. La tabla también suministra las condiciones con que fue ejecutado el ASEM, tales como el factor ''ρ''  que significa la tasa de reducción del paso de cada simplex de arranque por etapa, así como el punto inicial '''v'''<sub>0</sub> , donde '''1'''<sub>''d''</sub>  es el vector unitario de dimensión ''d'' . Además, el resto de los parámetros fueron ajustados de acuerdo con los valores indicados en la  [[#tbl0005|tabla 1]] .      

Nótese que para los tipos de problemas denotados como FGM<math display="inline">\mathbb{R}\mathbb{Z}[n+m]</math>  no se especificó el valor teórico debido a que esos tipos de problemas presentan múltiples mínimos locales y se ha supuesto que el mínimo identificado por el ASEM es un mínimo local al problema, mas no el mínimo global.      

Para los casos en donde fue considerado un punto de inicio '''v'''<sub>0</sub>  con componentes todas iguales a 10, se puede inferir que, de haberse identificado sus respectivos mínimos por búsqueda exhaustiva, solo por las variables enteras, se hubiesen requerido 10<sup>20</sup>  evaluaciones, es decir, 100 × 10<sup>6</sup>  T evaluaciones de la función objetivo, lo cual permite inferir la alta eficiencia que posee el ASEM desde el punto de vista del NE, en el caso de funciones objetivo convexas.      

La [[#fig0025|figura 5]]  muestra la gráfica de <math display="inline">f(\boldsymbol{\mbox{v}}_1^{\left[q\right]})</math>  del problema FC<math display="inline">\mathbb{R}\mathbb{Z}[5+10]</math>  en función de la iteración ''q''  del ASEM, para un valor de ''ρ''  = 0, 8, dejando ajustado el resto de los parámetros según lo indicado en la [[#tbl0005|tabla 1]] , y punto <math display="inline">\boldsymbol{\mbox{v}}_0^t=(10\cdot \boldsymbol{1}_5^t,10\cdot \boldsymbol{1}_{10}^t)</math> . La figura muestra al final de cada ''s'' -ésima etapa, sobre la escala denotada por ''q'' , el valor alcanzado del número de iteraciones ''q'' . Por otro lado, en la parte inferior de la figura se muestra la escala correspondiente al número de la etapa ''s''  del ASEM, y es denotada por ''s  '' . Nótese además que en la figura se representa el valor de <math display="inline">f(\boldsymbol{\mbox{v}}_1^{\left[1\right]})</math> , es decir, al final de la primera iteración, cuyo valor fue de 1462.

<span id='fig0025'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_404207060-1-s2.0-S0213131513000345-gr5.jpg|center|376px|Evolución de la función objetivo para el Problema FCℝℤ[5+10], para el caso ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 5.

Evolución de la función objetivo para el Problema FC<math display="inline">\mathbb{R}\mathbb{Z}[5+10]</math> , para el caso cuando ''ρ''  = 0, 8.                  

</span>
|}

Por otra parte, la [[#tbl0020|tabla 4]]  presenta un resumen por etapa del ASEM, cuando este fue aplicado al problema FC<math display="inline">\mathbb{R}\mathbb{Z}[5+10]</math>  bajo las mismas condiciones con que fue reportada la [[#fig0025|figura 5]] . Nótese que los datos numéricos que se muestran en la columna ''q''  corresponden al número máximo de iteraciones empleadas por el ASEM en cada etapa, así como el NE de la función objetivo.

<span id='tbl0020'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabla 4.

Reporte por etapa del Problema FC<math display="inline">\mathbb{R}\mathbb{Z}[5+10]</math> , para el caso cuando ''ρ''  = 0, 8                  

|-

! ''s''
! ℓ<sub>''x''</sub>
! <math display="inline">{\kappa }_x^{\left[s\right]}</math>
! <math display="inline">{\Delta }_r^{\left[s\right]}</math>
! <math display="inline">f({\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}}}^{\left[s\right]})</math>
! ''q''
! NE
|-

| 1
| 0,55972
| 1
| 1
| 397,8225
| 44
| 109
|-

| 2
| 0,25699
| 0,3
| 0,8
| 27,1793
| 385
| 671
|-

| 3
| 0,07269
| 0,09
| 0,064
| 9,0016
| 124
| 250
|-

| 4
| 0,01749
| 0,027
| 0,512
| 2,488×10<sup>−6</sup>
| 170
| 334
|-

| 5
| 0,00766
| 0,0081
| 0,4096
| 7,379×10<sup>−13</sup>
| 267
| 491
|}

<span id='sec0050'></span>
==10. Comparación con otros métodos==

El problema del diseño óptimo de un tanque presurizado ha servido de ejemplo de aplicación y comparación de diversos métodos algorítmicos de optimización, entre los cuales se pueden señalar los métodos de optimización estudiados por Zahara y Kao [[#bib0040|[8]]]  y por Tahera et al. [[#bib0045|[9]]] .      

La [[#fig0030|figura 6]]  muestra la concha semiesférica izquierda y el cilindro de un tanque presurizado, mas no la concha semiesférica derecha que es también parte de la estructura del tanque. La gráfica presenta la geometría del tanque que hay que diseñar junto a las variables de decisión: ''R  ''  (<math display="inline">x_1\in \mathbb{R}</math> , radio interno del cilindro y de las conchas semiesféricas medido en pulgadas); ''L  ''  (<math display="inline">x_2\in \mathbb{R}</math> , longitud del cilindro en pulgadas); ''E''<sub>''p''</sub>  (<math display="inline">y_1\in \mathbb{Z}</math> , espesor de la pared del cilindro en cantidades enteras de la fracción <math display="inline">\frac{1}{16}</math>  pulgadas, debido a los distintos calibres con que es fabricado el material) y ''E''<sub>''c''</sub>  (<math display="inline">y_2\in \mathbb{Z}</math> , espesor de la pared de las conchas semiesféricas en cantidades enteras de la fracción <math display="inline">\frac{1}{16}</math>  pulgadas).

<span id='fig0030'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_404207060-1-s2.0-S0213131513000345-gr6.jpg|center|376px|Tanque cilíndrico presurizado.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 6.

Tanque cilíndrico presurizado.

</span>
|}

El problema en términos matemáticos presentado por Zahara y Kao [[#bib0040|[8]]]  and [[#bib0045|[9]]] , que expresa en dólares americanos el costo de producir un tanque cilíndrico presurizado en el que se adecuan las variables a la notación empleada en esta investigación, viene dada por

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\underset{\boldsymbol{\mbox{x}}\in {\mathbb{R}}^2,\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^2}{minimice\quad }f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}}),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}

donde para todo <math display="inline">\boldsymbol{\mbox{x}}\in {\mathbb{R}}^2</math>  e <math display="inline">\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^2</math> ,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})=0,6224y_1x_1x_2+</math><math>1,7781y_2x_1^2+3,1661y_1^2x_2+19,84y_2^2x_1,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 25)
|}

sujeto a:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>g_1:0,0193x_1-y_1\leq 0;</math>
|-
|<math>g_2:0,00954x_1-y_2\leq 0;</math>
|-
|<math>g_3:1296000-\pi x_1^2x_2-\frac{4}{3}\pi x_1^3\leq 0;</math>
|-
|<math>g_4:x_2-240\leq 0;</math>
|-
|<math>b_1:0\leq x_1\leq 200;</math>
|-
|<math>b_2:0\leq x_2\leq 200;</math>
|-
|<math>b_3:0\leq y_1\leq 100;</math>
|-
|<math>b_4:0\leq y_2\leq 100.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 26)
|}

Con el propósito de redefinir el problema como un problema de minimización irrestricto, se definieron las siguientes funciones de penalización:

<span id='eq0145'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>p_1(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^2max(0,0193x_1-</math><math>y_1,0);</math>
|-
|<math>p_2(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^2max(0,00954x_1-</math><math>y_2,0);</math>
|-
|<math>p_3(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^2max(1296000-</math><math>\pi x_1^2x_2-\frac{4}{3}\pi x_1^3,0);</math>
|-
|<math>p_4(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^2max(x_2-</math><math>240,0);</math>
|-
|<math>p_5(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^6(max(-</math><math>x_1,0)+max(x_1-200,0));</math>
|-
|<math>p_6(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^6(max(-</math><math>x_2,0)+max(x_2-200,0));</math>
|-
|<math>p_7(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^6(max(-</math><math>y_1,0)+max(y_1-100,0));</math>
|-
|<math>p_8(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})={10}^6(max(-</math><math>y_2,0)+max(y_2-100,0)),</math>
|-
|<math></math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 27)
|}

lo que permitió redefinir el problema a través de las funciones de penalización dadas por la ecuación [[#eq0145|(27)]]  como:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\underset{\boldsymbol{\mbox{x}}\in {\mathbb{R}}^2,\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^2}{minimice\quad }\left[f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})+\right. </math><math>\left. \sum_{i=1}^8\quad p_i(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})\right].</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 28)
|}

Para la identificación del mínimo se empleó el ASEM con sus parámetros ajustados de acuerdo con lo especificado por la [[#tbl0005|tabla 1]]  y con valores de ''ρ''  = 0, 8 y <math display="inline">{\Delta }_r^{\left[1\right]}=</math><math>1</math> , y como punto de inicio se fijaron todas las componentes en 10.      

El resultado obtenido a través del ASEM se comparó con los resultados de otros métodos tales como: el GADYM, acrónimo del inglés ''Gender-Age structure DYnamic parameter tuning and Mandatory self perfection scheme'' , propuesto por Tahera et al.  [[#bib0045|[9]]] ; el GRTA, del inglés ''Generalized Random Tunneling Algorithm'' , desarrollado por Kitayama y Yamazaki  [[#bib0050|[10]]] ; el NM-PSO, acrónimo del inglés ''Nelder-Mead simplex search method and Particle Swarm Optimization'' , de Zahara y Kao  [[#bib0040|[8]]]  y el GRG, conocido del inglés ''Generalised Reduced Gradient'' , reportado también en Tahera et al.  [[#bib0045|[9]]] , los cuales son indicados en la [[#tbl0025|tabla 5]] .

<span id='tbl0025'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabla 5.

Resultados comparativos en el problema del diseño óptimo de un envase presurizado

|-

! 
! 
! 
! colspan="4" | Variables medidas en pulgadas [in]
|-

! Algoritmo
! ''f'' ('''x''' , '''y''' )                                                    
! ''S'' (''L'' , ''R'' ) [in<sup>2</sup> ]                                                    
! ''R'' (''x''<sub>1</sub> )                                                    
! ''L'' (''x''<sub>2</sub> )                                                    
! ''E''<sub>''p''</sub> (''y''<sub>1</sub> )                                                    
! ''E''<sub>''c''</sub> (''y''<sub>2</sub> )                                                    
|-

| GADYM
| 6.062,65
| 69.026,731
| 42,098
| 176,769
| 0,4357
| 0,8125
|-

| GRTA
| 6.059,58
| 68.992,022
| 42,098
| 176,634
| 0,4375
| 0,8125
|-

| NM-PSO
| 5.930,31
| 69.511,951
| 41,639
| 182,412
| 0,3972
| 0,8036
|-

| GRG
| 9.664,09
| 74.691,075
| 37,692
| 240,000
| 1 (<math display="inline">\frac{1}{16}</math> )                                                    
| 1 (<math display="inline">\frac{1}{16}</math> )                                                    
|-

| ASEM
| 8.796,92
| 60.623,882
| 51,807
| 84,627
| 1 (<math display="inline">\frac{1}{16}</math> )                                                    
| 1 (<math display="inline">\frac{1}{16}</math> )                                                    
|}

La [[#tbl0025|tabla 5]]  muestra los valores obtenidos de la función objetivo y sus correspondientes valores de las variables de decisión para el mejor punto identificado por cada algoritmo. Además, incluye el valor de la superficie, ''S'' (''L'' , ''R'' ) = 2''πRL''  + 4''πR''<sup>2</sup> , de material requerido para la construcción del envase presurizado, puesto que la comparación debe hacerse según la superficie de material requerido, en virtud de que únicamente el GRG y el ASEM tomaron en cuenta la naturaleza entera de las variables ''E''<sub>''c''</sub>  y ''E''<sub>''p''</sub> , las cuales están indicadas en la tabla mediante el número de fracciones de <math display="inline">\frac{1}{16}</math>  pulgadas. El hecho de no poder asignar valores reales a las variables ''E''<sub>''c''</sub>  y ''E''<sub>''p''</sub>  produce un aumento en la función objetivo con relación a los valores que hubieran tomado ''E''<sub>''c''</sub>  y ''E''<sub>''p''</sub> , en caso de haberse considerado variables reales. Nótese que de acuerdo con los resultados mostrados en la [[#tbl0025|tabla 5]] , el ASEM alcanzó a identificar la menor superficie de material para la fabricación del envase presurizado, con sólo 440 NE, mientras que el resto de los algoritmos emplearon decenas de miles de evaluaciones de la función objetivo.      

Más aún, el ASEM identificó un mejor resultado en comparación al algoritmo GRG, el cual es un algoritmo para problemas de OEM.

<span id='sec0055'></span>
==11. Conclusiones==

Como conclusiones de este trabajo de investigación se tiene que:
* El empleo de la estructura del simplex dentro del campo de los números enteros, junto a las redefinidas operaciones de reflexión, extensión, contracción y encogimiento, result&rsquo;0 muy apropiado para definir una estrategia de identificación de óptimos de problemas enteros mixtos irrestrictos, sin el relajamiento de las variables enteras de decisión.
* La principal característica que tienen los problemas de variables enteras mixtas en cuanto a la derivabilidad de sus funciones objetivo es que estas no cuentan con expresiones de derivadas con respecto a las variables enteras, desde el punto de vista estricto de la definición de derivada. Este hecho ha llevado a desarrollar un método libre de derivadas, además de haber desarrollado métodos algorítmicos que requieren un relativo bajo número de evaluaciones de la función objetivo.

Por otra parte, el desarrollo de otros enfoques para la identificación de óptimos así como de estructuras teóricas son temas de investigación que podrían constituir interesantes investigaciones, debido al reto que imponen. En específico, se pueden sugerir como temas de investigación:
* El estudio teórico de las condiciones de optimalidad dentro de la programación no lineal entera mixta.
* El desarrollo de funciones de penalización para el empleo del ASEM con problemas bajo restricciones. Este último puede constituir un importante trabajo de investigación, debido a las innumerables aplicaciones que tiene este tipo de problemas dentro de la ingeniería.

==Agradecimientos==

El autor desea expresar su gratitud al Consejo de Desarrollo Científico y Humanístico de la Universidad Central de Venezuela por la financiación de esta investigación a través del proyecto PI-08-7506-2009/1. Él desea realmente agradecer también a los árbitros por sus valiosos comentarios y sugerencias, los cuales permitieron mejorar el artículo.

<span id='sec0060'></span>

==Appendix A. Funciones de prueba==

En este apartado son mostrados los problemas enteros mixtos empleados en los experimentos numéricos, para lo cual se empleó una notación que permite definir la dimensión del problema de optimización. Además, para cada uno de los problemas mostrados se establece su respectiva solución óptima.

==Función cuadrática entera mixta: FC<math display="inline">\mathbb{R}\mathbb{Z}</math> [''n''  + ''m'' ]                            ==

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})=\sum_{i=1}^na_ix_i^2+</math><math>\sum_{i=1}^mb_iy_i^2,\quad \forall \boldsymbol{\mbox{x}}\in {\mathbb{R}}^n,\quad \boldsymbol{\mbox{y}}\in {\mathbb{Z}}^m,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 29)
|}

donde ''d''  es la dimensión del problema, y ''a''<sub>''i''</sub>  y ''b''<sub>''i''</sub>  para ''i''  = 1, …, ''d''  representan coeficientes reales positivos.      

'''Solución.'''  El mínimo global está ubicado en <math display="inline">{\boldsymbol{\overset{\mbox{ˆ}}{\mbox{x}}}}^t=</math><math>(\underset{n}{\underbrace{0,0,\ldots ,0}})</math>  e <math display="inline">{\boldsymbol{\overset{\mbox{ˆ}}{\mbox{y}}}}^t=</math><math>(\underset{m}{\underbrace{0,0,\ldots ,0}})</math>  con un valor de <math display="inline">f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})=</math><math>0</math> .

==Función cuadrática entera mixta desplazada del origen: FCD<math display="inline">\mathbb{R}\mathbb{Z}</math> [''n''  + ''m'' ]                            ==

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})=\sum_{i=1}^na_i{\left(x_i-x_{0i}\right)}^2+</math><math>\sum_{i=1}^mb_i{\left(y_i-y_{0i}\right)}^2,</math>
|-
|<math>\forall \boldsymbol{\mbox{x}}\in {\mathbb{R}}^n,\boldsymbol{\mbox{y}}\in {\mathbb{Z}}^m,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 30)
|}

donde ''n''  es la dimensión del subvector real, ''m  ''  es la dimensión del subvector entero, y <math display="inline">a_i\in \mathbb{R}</math>  para ''i''  = 1, …, ''n  '' , y <math display="inline">b_i\in \mathbb{R}</math>  para ''i''  = 1, …, ''m  ''  representan coeficientes reales positivos, y <math display="inline">x_{0i}\in \mathbb{R}</math>  para ''i''  = 1, …, ''n  '' , y <math display="inline">y_{0i}\in \mathbb{R}</math>  para ''i''  = 1, …, ''m''  representan las coordenadas del vértice de la función cuadrática.      

Para el caso particular de los experimentos numéricos ejecutados, se consideró ''a''<sub>''i''</sub>  = 1 y ''x''<sub>0''i''</sub>  = ''i'' /3 para todo ''i''  = 1, …, ''n'' , ''b''<sub>''i''</sub>  = 1 y ''y''<sub>0''i''</sub>  = ''i'' /3 para todo ''i''  = 1, …, ''m'' .      

'''Solución.'''  El mínimo global está ubicado en <math display="inline">{\boldsymbol{\overset{\mbox{ˆ}}{\mbox{x}}}}^t=</math><math>(\underset{n}{\underbrace{x_{01},x_{02},\ldots ,x_{0n}}})</math>  e <math display="inline">{\boldsymbol{\overset{\mbox{ˆ}}{\mbox{y}}}}^t=</math><math>(\underset{m}{\underbrace{{\overset{\mbox{ˆ}}{y}}_1,{\overset{\mbox{ˆ}}{y}}_2,\ldots ,{\overset{\mbox{ˆ}}{y}}_m}})</math> , donde

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\overset{\mbox{ˆ}}{y}}_i=\begin{array}{ll}
\lfloor y_{0i}\rfloor , & \forall y_{0i}-\lfloor y_{0i}\rfloor <\lceil y_{0i}\rceil -y_{0i};\\
\lfloor y_{0i}\rfloor \mbox{ }\mbox{ó}\mbox{ }\lceil y_{0i}\rceil , & \forall y_{0i}-\lfloor y_{0i}\rfloor =\lceil y_{0i}\rceil -y_{0i};\\
\lceil y_{0i}\rceil , & \forall y_{0i}-\lfloor y_{0i}\rfloor >\lceil y_{0i}\rceil -y_{0i},
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 31)
|}

con un valor de <math display="inline">f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})=</math><math>{\sum }_{i=1}^mb_i{\left({\overset{\mbox{ˆ}}{y}}_i-y_{0i}\right)}^2</math> .

==Función producto de módulos entero mixto: FPM<math display="inline">\mathbb{R}\mathbb{Z}[d+d]</math>==

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})=-d+</math><math>\sum_{i=1}^da_i(1+\vert x_i\vert )(1+\vert y_i\vert ),\forall \boldsymbol{\mbox{x}}\in {\mathbb{R}}^d,\quad \boldsymbol{\mbox{y}}\in {\mathbb{Z}}^d,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 32)
|}

donde ''d''  es la dimensión del problema y ''a''<sub>''i''</sub>  para todo ''i''  = 1, …, ''d''  representan coeficientes reales positivos.      

'''Solución.'''  El mínimo global está ubicado en <math display="inline">{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{x}}}}^t=</math><math>(\underset{d}{\underbrace{0,0,\ldots ,0}})</math>  e <math display="inline">{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{y}}}}^t=</math><math>(\underset{d}{\underbrace{0,0,\ldots ,0}})</math>  con un valor de <math display="inline">f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})=</math><math>0</math> .

==A.1. Función Rosenbrock: FR<math display="inline">\mathbb{R}\mathbb{Z}[d+d]</math>==

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{y}})=\sum_{i=0}^{d-1}{\left(x_{i+1}-y_{i+1}^2\right)}^2+</math><math>{\left(1-y_{i+1}\right)}^2,\quad \forall \boldsymbol{\mbox{x}}\in {\mathbb{R}}^d,\quad \boldsymbol{\mbox{y}}\in {\mathbb{Z}}^d,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 33)
|}

donde ''d''  es la dimensión del problema y es un número par positivo.      

'''Solución.'''  El mínimo global está ubicado en <math display="inline">{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{x}}}}^t=</math><math>\left(\underset{d}{\underbrace{1,1,\ldots ,1}}\right)</math>  e <math display="inline">{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{y}}}}^t=</math><math>\left(\underset{d}{\underbrace{1,1,\ldots ,1}}\right)</math>  con un valor de <math display="inline">f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{y}}})=</math><math>0</math> .

==Función Audet & Dennis: FAD<math display="inline">\mathbb{R}\mathbb{Z}[2+1]</math>==

El problema original de [[#bib0055|[11]]]  es:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\underset{\boldsymbol{\mbox{x}}\in {\mathbb{R}}^2,\boldsymbol{\mbox{y}}\in \mathbb{Z}}{minimice\quad }f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 34)
|}

sujeto a:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>-2\leq x_1\leq 2;\quad -2\leq x_2\leq 2;\quad y_1\in \lbrace 0,1\rbrace ,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 35)
|}

donde

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})=g(x_1,x_2)(1-</math><math>y_1)+h(x_1,x_2)y_1,\quad \forall \boldsymbol{\mbox{x}}\in {\mathbb{R}}^2,\quad \boldsymbol{\mbox{y}}\in \mathbb{Z};</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 36)
|}

y las funciones ''g'' (''x''<sub>1</sub> , ''x''<sub>2</sub> ) y ''h'' (''x''<sub>1</sub> , ''x''<sub>2</sub> ) para todo <math display="inline">\boldsymbol{\mbox{x}}\in {\mathbb{R}}^2</math>  vienen dadas por

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>g(x_1,x_2)=x_1^2+x_2^2;</math>
|-
|<math>h(x_1,x_2)=x_1^2x_2+x_1(1-x_2).</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 37)
|}

El planteamiento del problema, empleando el concepto de función de penalización, queda entonces de la siguiente forma:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})=g(x_1,x_2)(1-</math><math>y_1)+(x_1,x_2)y_1+p_1(x_1)+p_2(x_2)+p_3(y_1),\quad \forall \boldsymbol{\mbox{x}}\in {\mathbb{R}}^2,\quad \boldsymbol{\mbox{y}}\in \mathbb{Z},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 38)
|}

donde para todo <math display="inline">x_1\in \mathbb{R}</math> ; <math display="inline">x_2\in \mathbb{R}</math> ; <math display="inline">y_1\in \mathbb{Z}</math>  y ''k''  = 10<sup>3</sup>

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>p_1(x_1)=kmax(-2-x_1,0)+kmax(x_1-2,0);</math>
|-
|<math>p_2(x_2)=kmax(-2-x_2,0)+kmax(x_2-2,0);</math>
|-
|<math>p_3(y_1)=kmax(0-y_1,0)+kmax(y_1-1,0).</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 39)
|}

'''Solución.'''  El mínimo global está ubicado en el punto <math display="inline">{\overset{\mbox{ˆ}}{\boldsymbol{\mbox{x}}}}^t=</math><math>(-2,-2)</math>  e <math display="inline">\overset{\mbox{ˆ}}{\boldsymbol{\mbox{y}}}=</math><math>(0)</math>  con un valor de <math display="inline">f(\overset{\mbox{ˆ}}{\boldsymbol{\mbox{v}}})=</math><math>-14</math> .

==Función Griewank Modificada: FGM<math display="inline">\mathbb{R}\mathbb{Z}[n+m]</math>==

Una versión modificada en el campo de los números enteros mixtos al problema de optimización de la función de Griewank presentada por Fu et al. [[#bib0060|[12]]]  se muestra a continuación:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(\boldsymbol{\mbox{x}},\boldsymbol{\mbox{y}})=1+\frac{1}{20}\sum_{i=1}^nx_i^2+</math><math>\frac{1}{20}\sum_{i=1}^my_i^2-\prod_{i=1}^ncos(\frac{2\pi }{5}x_i)-</math><math>\prod_{i=1}^mcos(\frac{2\pi }{5}y_i),</math>
|-
|<math>\forall \boldsymbol{\mbox{x}}\in {\mathbb{R}}^n,\quad \boldsymbol{\mbox{y}}\in {\mathbb{Z}}^m.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 40)
|}

==References==

<ol style='list-style-type: none;margin-left: 0px;'><li><span id='bib0005'></span>
[[#bib0005|[1]]] R. Hooke, T.A. Jeeves; “Direct Search” Solution of Numerical and Statistical Problems; Journal of the Association for Computing Machinery, 8 (2) (1961), pp. 212–229</li>
<li><span id='bib0010'></span>
[[#bib0010|[2]]] J.A. Nelder, R. Mead; A simplex method for function minimization; The Computer Journal, 7 (1965), pp. 308–313</li>
<li><span id='bib0015'></span>
[[#bib0015|[3]]] E. Brea, (2009) Extensiones del método de Nelder Mead a problemas de variables enteras y enteras mixtas. Trabajo de Ascenso a la categoría de Profesor Titular, Universidad Central de Venezuela. Caracas, Venezuela.</li>
<li><span id='bib0020'></span>
[[#bib0020|[4]]] J. Brandts, S. Korotov, M. Křížek, J. Šolc; On Nonobtuse Simplicial Partitions; SIAM Review, 51 (2) (2009), pp. 317–335</li>
<li><span id='bib0025'></span>
[[#bib0025|[5]]] E. Brea, (2004) Nelder-Mead Optimization Under Linear Constraints. Ph.D. Thesis, University of Southampton, Southampton, England.</li>
<li><span id='bib0030'></span>
[[#bib0030|[6]]] R.R. Barton, J.S. Ivey Jr.; Nelder-Mead simplex modifications simulation optimization; Management Science, 42 (7) (1996), pp. 954–973</li>
<li><span id='bib0035'></span>
[[#bib0035|[7]]] D.G. Humphrey, J.R. Wilson; A revised simplex procedure for stochastic simulation response surface optimization; INFORMS Journal on Computing, 12 (4) (2000), pp. 272–283</li>
<li><span id='bib0040'></span>
[[#bib0040|[8]]] E. Zahara, Y.T. Kao; Hybrid Nelder-Mead simplex search and particle swarm optimization for constrained engineering design problems; Expert Systems with Applications, 36 (2, Part 2) (2009), pp. 3880–3886</li>
<li><span id='bib0045'></span>
[[#bib0045|[9]]] K. Tahera, R.N. Ibrahim, P.B. Lochert; CADYM - A novel genetic algorithm in mechanical design problems; Journal of Universal Computer Science, 14 (15) (2008), pp. 2566–2581</li>
<li><span id='bib0050'></span>
[[#bib0050|[10]]] S. Kitayama, K. Yamazaki; Global Optimization by Generalized Random Tunneling Algorithm (4th Report Application to the Nonlinear Optimum Design Problem of the Mixed Design Variables); Journal of Computational Science and Technology, 2 (1) (2008), pp. 258–267</li>
<li><span id='bib0055'></span>
[[#bib0055|[11]]] C. Audet, J.E. Dennis Jr.; Pattern Search Algorithms for Mixed Variable Programming; SIAM Journal on Optimization, 11 (3) (2001), pp. 573–594</li>
<li><span id='bib0060'></span>
[[#bib0060|[12]]] M. Fu, S. Marcus, J. Hu; Model-Based Randomized Methods for Global Optimization; 17th International Symposium on Mathematical Theory of Networks and Systems, vol. 1, Kyoto, Japan (2006), pp. 355–363</li>
</ol>