localización en una red con patrón de elección - revista ...

localización en una red con patrón de elección - revista ... localización en una red con patrón de elección - revista ...

rev.inv.ope.univ.paris1.fr
from rev.inv.ope.univ.paris1.fr More from this publisher
25.07.2013 Views

REVISTA INVESTIGACIÓN OPERACIONAL VOL., 32, NO. 1, 38-52, 2011 LOCALIZACIÓN EN UNA RED CON PATRÓN DE ELECCIÓN DEFINIDO POR UNA DISTANCIA UMBRAL Blas Pelegrín Pelegrín 1 , Saúl Cano Hernández 2 Departamento de Estadística e Investigación Operativa, Universidad de Murcia. 1. INTRODUCCIÓN ABSTRACT We consider the facility location problem on a network for an entering firm, in competition with other already established facilities, with the objective of market share maximization. The consumers choose the facility from which they obtain a maximum utility (binary preference). If a consumer obtains the maximum utility from a preexisting center and a new one, a proportion of his demand is captured by the new facility. The location candidates are the nodes and the points in the edges of the network. If the highest utility is obtained from a facility located to a distance within a certain threshold from the customer, it is proved that the set of candidates to optimal solution is a finite set of points in the network. A procedure to generate the candidates to optimal location is given and a formulation as a mixed integer linear programming problem is presented. A sensitivity analysis related to the proportion, the number of preexisting facilities and the number of new facilities, applied to the Region of Murcia (Spain), is shown. KEYWORDS: Competitive location, binary choice, threshold distance, discrete optimization. MSC: 90C10 RESUMEN Se considera el problema de localización de centros para una firma entrante en una red de transporte, en competencia con otros centros ya establecidos y con el objetivo de maximizar la cuota de mercado. Los consumidores compran en el centro del que obtengan la mayor utilidad (patrón binario). Si un consumidor obtiene la misma utilidad en un centro preexistente y uno nuevo, una proporción de su demanda es capturada por el nuevo. Las posibles localizaciones son los nodos y los puntos en los tramos de la red de transporte. Si la mayor utilidad se obtiene en un centro situado a una distancia dentro de un cierto umbral, se demuestra que el conjunto de candidatos a solución óptima se puede reducir a un conjunto finito de puntos de la red. Se da un procedimiento para generar los posibles candidatos, y se presenta una formulación como problema de programación lineal entera mixta para encontrar las localizaciones óptimas. Se realiza un análisis de sensibilidad respecto de la proporción, el número de centros preexistentes y el número de nuevos centros, aplicado a la Región de Murcia (España). Se considera el problema de localización de una firma entrante que tiene que competir con otras ya establecidas. Los consumidores compran de la firma que les proporciona la máxima utilidad. El objetivo de la firma nueva es encontrar la localización de sus establecimientos que le proporcione la mayor cuota de mercado. En los modelos estudiados en la literatura, la utilidad de un establecimiento depende de su distancia al consumidor y de otros factores, entre los que destacan la calidad y el precio (véase Eiselt y Laporte (1998), y García y otros (2004)). Los modelos que tienen en cuenta la calidad consideran la utilidad como una atracción por el centro, que es definida por alguna función directamente proporcional a la calidad e inversamente proporcional a la distancia (véase Hakimi (1983), Serra y ReVelle (1995), Plastria (2001)). Algunos trabajos consideran un comportamiento binario de los consumidores, que consiste en comprar del centro con mayor utilidad, desarrollando algoritmos para encontrar las localizaciones óptimas (véase Suárez-Vega y otros (2004)). En otros modelos se tiene en cuenta el precio del producto, y la utilidad es mayor cuanto menor es el coste total para adquirir el producto (véase Gabszewicz y Thisse (1992), y Pelegrín y otros (2006)). 1 pelegrin@um.es, 2 saul.cano@um.es 38

REVISTA INVESTIGACIÓN OPERACIONAL VOL., 32, NO. 1, 38-52, 2011<br />

LOCALIZACIÓN EN UNA RED CON PATRÓN DE<br />

ELECCIÓN DEFINIDO POR UNA DISTANCIA<br />

UMBRAL<br />

Blas Pelegrín Pelegrín 1 , Saúl Cano Hernán<strong>de</strong>z 2<br />

Departam<strong>en</strong>to <strong>de</strong> Estadística e Investigación Operativa, Universidad <strong>de</strong> Murcia.<br />

1. INTRODUCCIÓN<br />

ABSTRACT<br />

We <strong>con</strong>si<strong>de</strong>r the facility location problem on a network for an <strong>en</strong>tering firm, in competition with other already<br />

established facilities, with the objective of market share maximization. The <strong>con</strong>sumers choose the facility from which<br />

they obtain a maximum utility (binary prefer<strong>en</strong>ce). If a <strong>con</strong>sumer obtains the maximum utility from a preexisting<br />

c<strong>en</strong>ter and a new one, a proportion of his <strong>de</strong>mand is captu<strong>red</strong> by the new facility. The location candidates are the<br />

no<strong>de</strong>s and the points in the edges of the network. If the highest utility is obtained from a facility located to a distance<br />

within a certain threshold from the customer, it is proved that the set of candidates to optimal solution is a finite set of<br />

points in the network. A procedure to g<strong>en</strong>erate the candidates to optimal location is giv<strong>en</strong> and a formulation as a<br />

mixed integer linear programming problem is pres<strong>en</strong>ted. A s<strong>en</strong>sitivity analysis related to the proportion, the number of<br />

preexisting facilities and the number of new facilities, applied to the Region of Murcia (Spain), is shown.<br />

KEYWORDS: Competitive location, binary choice, threshold distance, discrete optimization.<br />

MSC: 90C10<br />

RESUMEN<br />

Se <strong>con</strong>si<strong>de</strong>ra el problema <strong>de</strong> <strong>localización</strong> <strong>de</strong> c<strong>en</strong>tros para <strong>una</strong> firma <strong>en</strong>trante <strong>en</strong> <strong>una</strong> <strong>red</strong> <strong>de</strong> transporte, <strong>en</strong> compet<strong>en</strong>cia<br />

<strong>con</strong> otros c<strong>en</strong>tros ya establecidos y <strong>con</strong> el objetivo <strong>de</strong> maximizar la cuota <strong>de</strong> mercado. Los <strong>con</strong>sumidores compran <strong>en</strong><br />

el c<strong>en</strong>tro <strong>de</strong>l que obt<strong>en</strong>gan la mayor utilidad (<strong>patrón</strong> binario). Si un <strong>con</strong>sumidor obti<strong>en</strong>e la misma utilidad <strong>en</strong> un c<strong>en</strong>tro<br />

preexist<strong>en</strong>te y uno nuevo, <strong>una</strong> proporción <strong>de</strong> su <strong>de</strong>manda es capturada por el nuevo. Las posibles localizaciones son<br />

los nodos y los puntos <strong>en</strong> los tramos <strong>de</strong> la <strong>red</strong> <strong>de</strong> transporte. Si la mayor utilidad se obti<strong>en</strong>e <strong>en</strong> un c<strong>en</strong>tro situado a <strong>una</strong><br />

distancia d<strong>en</strong>tro <strong>de</strong> un cierto umbral, se <strong>de</strong>muestra que el <strong>con</strong>junto <strong>de</strong> candidatos a solución óptima se pue<strong>de</strong> <strong>red</strong>ucir a<br />

un <strong>con</strong>junto finito <strong>de</strong> puntos <strong>de</strong> la <strong>red</strong>. Se da un procedimi<strong>en</strong>to para g<strong>en</strong>erar los posibles candidatos, y se pres<strong>en</strong>ta <strong>una</strong><br />

formulación como problema <strong>de</strong> programación lineal <strong>en</strong>tera mixta para <strong>en</strong><strong>con</strong>trar las localizaciones óptimas. Se realiza<br />

un análisis <strong>de</strong> s<strong>en</strong>sibilidad respecto <strong>de</strong> la proporción, el número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes y el número <strong>de</strong> nuevos<br />

c<strong>en</strong>tros, aplicado a la Región <strong>de</strong> Murcia (España).<br />

Se <strong>con</strong>si<strong>de</strong>ra el problema <strong>de</strong> <strong>localización</strong> <strong>de</strong> <strong>una</strong> firma <strong>en</strong>trante que ti<strong>en</strong>e que competir <strong>con</strong> otras ya establecidas.<br />

Los <strong>con</strong>sumidores compran <strong>de</strong> la firma que les proporciona la máxima utilidad. El objetivo <strong>de</strong> la firma nueva es<br />

<strong>en</strong><strong>con</strong>trar la <strong>localización</strong> <strong>de</strong> sus establecimi<strong>en</strong>tos que le proporcione la mayor cuota <strong>de</strong> mercado.<br />

En los mo<strong>de</strong>los estudiados <strong>en</strong> la literatura, la utilidad <strong>de</strong> un establecimi<strong>en</strong>to <strong>de</strong>p<strong>en</strong><strong>de</strong> <strong>de</strong> su distancia al<br />

<strong>con</strong>sumidor y <strong>de</strong> otros factores, <strong>en</strong>tre los que <strong>de</strong>stacan la calidad y el precio (véase Eiselt y Laporte (1998), y<br />

García y otros (2004)). Los mo<strong>de</strong>los que ti<strong>en</strong><strong>en</strong> <strong>en</strong> cu<strong>en</strong>ta la calidad <strong>con</strong>si<strong>de</strong>ran la utilidad como <strong>una</strong> atracción<br />

por el c<strong>en</strong>tro, que es <strong>de</strong>finida por alg<strong>una</strong> función directam<strong>en</strong>te proporcional a la calidad e inversam<strong>en</strong>te<br />

proporcional a la distancia (véase Hakimi (1983), Serra y ReVelle (1995), Plastria (2001)). Algunos trabajos<br />

<strong>con</strong>si<strong>de</strong>ran un comportami<strong>en</strong>to binario <strong>de</strong> los <strong>con</strong>sumidores, que <strong>con</strong>siste <strong>en</strong> comprar <strong>de</strong>l c<strong>en</strong>tro <strong>con</strong> mayor<br />

utilidad, <strong>de</strong>sarrollando algoritmos para <strong>en</strong><strong>con</strong>trar las localizaciones óptimas (véase Suárez-Vega y otros (2004)).<br />

En otros mo<strong>de</strong>los se ti<strong>en</strong>e <strong>en</strong> cu<strong>en</strong>ta el precio <strong>de</strong>l producto, y la utilidad es mayor cuanto m<strong>en</strong>or es el coste total<br />

para adquirir el producto (véase Gabszewicz y Thisse (1992), y Pelegrín y otros (2006)).<br />

1 pelegrin@um.es, 2 saul.cano@um.es<br />

38


En este trabajo se g<strong>en</strong>eralizan algunos <strong>de</strong> los mo<strong>de</strong>los anteriores, incluy<strong>en</strong>do todos los casos <strong>en</strong> los que la<br />

máxima utilidad se obti<strong>en</strong>e cuando el <strong>con</strong>sumidor se <strong>en</strong>cu<strong>en</strong>tra a <strong>una</strong> distancia <strong>de</strong> la firma no superior a un<br />

<strong>de</strong>terminado (distancia umbral). Este comportami<strong>en</strong>to es <strong>de</strong>terminista <strong>con</strong> información completa. Sólo <strong>en</strong> el<br />

caso <strong>de</strong> empate <strong>en</strong> máxima utilidad <strong>de</strong> la firma <strong>en</strong>trante <strong>con</strong> otra firma ya establecida, los <strong>con</strong>sumidores se<br />

repartirán <strong>en</strong>tre ellas para adquirir el producto. En los mo<strong>de</strong>los <strong>con</strong> precios, el producto suele ser homogéneo, ya<br />

que <strong>en</strong> esos mo<strong>de</strong>los es el precio el factor que afecta a la <strong>elección</strong> <strong>de</strong> los <strong>con</strong>sumidores. En los mo<strong>de</strong>los que<br />

usan funciones <strong>de</strong> atracción, po<strong>de</strong>mos <strong>con</strong>si<strong>de</strong>rar difer<strong>en</strong>ciación <strong>de</strong> producto o que los <strong>con</strong>sumidores adquier<strong>en</strong><br />

varios productos al <strong>de</strong>splazarse al establecimi<strong>en</strong>to. Esto se ti<strong>en</strong>e <strong>en</strong> cu<strong>en</strong>ta junto <strong>con</strong> otros factores al <strong>de</strong>terminar<br />

la atracción por el c<strong>en</strong>tro (véase Nakanishi y Cooper (1974)).<br />

Se supone que la <strong>de</strong>manda está <strong>con</strong>c<strong>en</strong>trada <strong>en</strong> los nodos <strong>de</strong> <strong>una</strong> <strong>red</strong> <strong>de</strong> transporte (véase Francis y otros (2002)<br />

para la agregación <strong>de</strong> la <strong>de</strong>manda) y es <strong>una</strong> cantidad fija para cada nodo. A difer<strong>en</strong>cia <strong>de</strong> otros mo<strong>de</strong>los, que<br />

<strong>con</strong>si<strong>de</strong>ran sólo los nodos <strong>de</strong> la <strong>red</strong> como posibles candidatos a <strong>localización</strong>, aquí todos los puntos <strong>de</strong> la <strong>red</strong> son<br />

<strong>con</strong>si<strong>de</strong>rados. Se <strong>de</strong>muestra que hay un <strong>con</strong>junto finito <strong>de</strong> puntos <strong>en</strong>tre los cuales se <strong>en</strong>cu<strong>en</strong>tran las<br />

localizaciones óptimas. Se pres<strong>en</strong>ta un algoritmo que g<strong>en</strong>era <strong>de</strong> manera efici<strong>en</strong>te la lista <strong>de</strong> candidatos y <strong>una</strong> vez<br />

discretizado el problema se proce<strong>de</strong> a su resolución mediante su formulación como programación lineal <strong>en</strong>tera<br />

mixta. Se realiza un experim<strong>en</strong>to computacional y unos análisis <strong>de</strong> s<strong>en</strong>sibilidad, aplicados a la región <strong>de</strong><br />

Murcia, España, y finalm<strong>en</strong>te se pres<strong>en</strong>tan alg<strong>una</strong>s <strong>con</strong>clusiones.<br />

2. NOTACIÓN Y DEFINICIONES BÁSICAS<br />

Sea NV ( , El , ) <strong>una</strong> <strong>red</strong> <strong>de</strong> transporte no dirigida, don<strong>de</strong> : 1,2, … es un <strong>con</strong>junto <strong>de</strong> nodos o<br />

vértices, : , ;, es un <strong>con</strong>junto <strong>de</strong> tramos o aristas,: es <strong>una</strong> función que asigna su<br />

longitud a cada tramo e . Sea Γ(i) el <strong>con</strong>junto <strong>de</strong> nodos adyac<strong>en</strong>tes al nodo i . Con dxy ( , ) d<strong>en</strong>otaremos la<br />

distancia (longitud <strong>de</strong>l camino más corto) <strong>en</strong>tre los puntos x y y . Por simplicidad, d<strong>en</strong>otaremos por dix a la<br />

distancia <strong>de</strong>l nodo i al punto . Sea X un <strong>con</strong>junto finito <strong>de</strong> puntos sobre la <strong>red</strong> N , <strong>en</strong>tonces <br />

: d<strong>en</strong>otará la mínima distancia <strong>de</strong>l nodo i a los punto <strong>de</strong> X . Sea D i <strong>una</strong> distancia fija asociada<br />

al nodo i .<br />

Definición 1. Decimos que un punto x <strong>en</strong> N es isodistante relativo a un nodo i si dix = Di.<br />

Como la función distancia d ix es <strong>una</strong> función cóncava lineal <strong>en</strong> dos partes cuando x varía a lo largo <strong>de</strong><br />

cualquier tramo [ uv, , ] resulta que a lo más pue<strong>de</strong> haber dos puntos isodistantes <strong>en</strong> el tramo [ uv , ] relativos al<br />

nodo i . Cada punto isodistante relativo a un nodo i pue<strong>de</strong> ser <strong>de</strong> uno <strong>de</strong> los tres tipos <strong>de</strong> puntos que se<br />

l r lr<br />

muestran <strong>en</strong> Figura 1, d<strong>en</strong>otados por x i , x i y x i , respectivam<strong>en</strong>te. Se verifica que si ,; dix < Disi;<br />

,; dix < Disi;<br />

, ∪ ,<br />

Figura 1<br />

Definición 2. Un punto es llamado isodistante si éste es isodistante relativo a algún nodo.<br />

La <strong>de</strong>finición anterior nos dice que un punto es isodistante si se <strong>en</strong>cu<strong>en</strong>tra a <strong>una</strong> p<strong>red</strong>eterminada distancia<br />

<strong>de</strong> algún nodo <strong>de</strong> la <strong>red</strong>.<br />

39


Un segm<strong>en</strong>to [ x, y ] <strong>en</strong> un tramo [ uv , ] es el sub<strong>con</strong>junto <strong>de</strong> puntos <strong>de</strong> [ uv , ] <strong>en</strong>tre los puntos x y y ,<br />

incluy<strong>en</strong>do a estos puntos. El segm<strong>en</strong>to abierto es , , , y los segm<strong>en</strong>tos semiabiertos<br />

son , , y , , .<br />

Definición 3. Si <strong>con</strong>si<strong>de</strong>ramos todos los puntos isodistantes sobre , , cualquier segm<strong>en</strong>to abierto<br />

<strong>de</strong>finido por dos puntos isodistantes <strong>con</strong>secutivos o por un nodo y su punto isodistante más cercano, es<br />

llamado segm<strong>en</strong>to iso-abierto.<br />

En la Figura 2, pued<strong>en</strong> verse segm<strong>en</strong>tos iso-abiertos <strong>de</strong>terminados por puntos isodistantes.<br />

Figura 2<br />

Un punto isodistante <strong>en</strong> , pue<strong>de</strong> ser isodistante relativo a más <strong>de</strong> un nodo. Sea V x el <strong>con</strong>junto <strong>de</strong><br />

nodos i tales que x es isodistante relativo a i . La sigui<strong>en</strong>te <strong>de</strong>finición <strong>de</strong> puntos isodistantes se usará<br />

para <strong>con</strong>si<strong>de</strong>rar tres clases <strong>de</strong> puntos isodistantes <strong>en</strong> el interior <strong>de</strong> cada tramo , , el cual <strong>de</strong>p<strong>en</strong><strong>de</strong> <strong>de</strong>l<br />

nodo por el cual se pase para obt<strong>en</strong>er D i .<br />

Definición 4. Sea x un punto isodistante <strong>en</strong> el interior <strong>de</strong> algún tramo:<br />

l<br />

i) Si x = xi<br />

para todo , el punto x es llamado isodistante-izquierdo.<br />

r<br />

ii) Si x = xi<br />

para todo , el punto x es llamado isodistante-<strong>de</strong>recho.<br />

iii) En otro caso, el punto x es llamado isodistante-mixto.<br />

Si no hay puntos isodistantes <strong>en</strong> , , este segm<strong>en</strong>to es el único segm<strong>en</strong>to iso-abierto <strong>en</strong> el tramo, .<br />

3. MODELOS DE LOCALIZACIÓN COMPETITIVA CON PATRÓN DE ELECCIÓN BINARIO<br />

Una firma <strong>de</strong>sea <strong>en</strong>trar al mercado t<strong>en</strong>i<strong>en</strong>do que <strong>de</strong>cidir la <strong>localización</strong> <strong>de</strong> r nuevos c<strong>en</strong>tros <strong>en</strong> puntos <strong>de</strong><br />

<strong>una</strong> <strong>red</strong> N , <strong>con</strong> el objetivo <strong>de</strong> maximizar su cuota <strong>de</strong> mercado. La <strong>de</strong>manda se <strong>en</strong>cu<strong>en</strong>tra <strong>en</strong> algunos<br />

nodos <strong>de</strong> la <strong>red</strong>, llamados nodos población, la cual es <strong>con</strong>ocida y fija. Exist<strong>en</strong> otro tipo <strong>de</strong> nodos <strong>en</strong> la<br />

<strong>red</strong>, llamados nodos cruce, los cuales no ti<strong>en</strong><strong>en</strong> <strong>de</strong>manda. Sea w i la cantidad <strong>de</strong> producto requerida por<br />

los <strong>con</strong>sumidores <strong>en</strong> el nodo i . El coste marginal <strong>de</strong> transporte <strong>de</strong> los <strong>con</strong>sumidores que se <strong>en</strong>cu<strong>en</strong>tran <strong>en</strong><br />

el nodo i a un punto localizado a <strong>una</strong> distancia d <strong>de</strong> i , vi<strong>en</strong>e dado por ti( d ) , don<strong>de</strong> ti( d ) es <strong>una</strong><br />

función creci<strong>en</strong>te <strong>de</strong> la distancia y <strong>con</strong> la <strong>con</strong>dición t i(0)<br />

= 0 .<br />

e<br />

En el mercado hay establecidos un <strong>con</strong>junto F = { 1,2,...,q } <strong>de</strong> c<strong>en</strong>tros. Esto es, hay q c<strong>en</strong>tros<br />

preexist<strong>en</strong>tes <strong>en</strong> la <strong>red</strong>, que son propiedad <strong>de</strong> otras firmas <strong>con</strong> las que va a competir la firma <strong>en</strong>trante. El<br />

<strong>con</strong>junto <strong>de</strong> todos los c<strong>en</strong>tros preexist<strong>en</strong>tes y nuevos, es repres<strong>en</strong>tado por 1,<br />

2, … , . D<strong>en</strong>otaremos por X = { xq+ 1, xq+<br />

2,... x q+ r}<br />

al <strong>con</strong>junto <strong>de</strong> r puntos <strong>de</strong> la <strong>red</strong> don<strong>de</strong> los<br />

nuevos c<strong>en</strong>tros son localizados. Para cada , , <strong>de</strong>finimos uij( d ij ) como la utilidad que supone<br />

comprar <strong>en</strong> el c<strong>en</strong>tro j para los <strong>con</strong>sumidores <strong>de</strong>l nodo i . Los <strong>con</strong>sumidores <strong>de</strong> cada nodo i compran<br />

<strong>en</strong> un c<strong>en</strong>tro k , si la utilidad que obti<strong>en</strong><strong>en</strong> es la mayor, es <strong>de</strong>cir, , . Este tipo<br />

<strong>de</strong> comportami<strong>en</strong>to recibe el nombre <strong>de</strong> binario. Si para un nodo i la máxima utilidad se obti<strong>en</strong>e <strong>en</strong> un<br />

c<strong>en</strong>tro <strong>de</strong> la firma <strong>en</strong>trante y otro <strong>de</strong> las firmas preexist<strong>en</strong>tes, <strong>en</strong>tonces supondremos que <strong>una</strong> proporción<br />

q i <strong>de</strong> los <strong>con</strong>sumidores compran <strong>en</strong> la firma <strong>en</strong>trante. La cantidad q i es <strong>de</strong>s<strong>con</strong>ocida y habría que<br />

estimarla para cada situación. Vamos a ver que <strong>en</strong> ocasiones, este comportami<strong>en</strong>to equivale a que los<br />

<strong>con</strong>sumidores <strong>en</strong> el nodo i compran <strong>en</strong> un nuevo c<strong>en</strong>tro si , para algún valor D i .<br />

Ejemplo 1. Un mo<strong>de</strong>lo <strong>con</strong> precios <strong>en</strong> orig<strong>en</strong>.<br />

40


En este tipo <strong>de</strong> mo<strong>de</strong>los los <strong>con</strong>sumidores se <strong>de</strong>splazan al c<strong>en</strong>tro a satisfacer su <strong>de</strong>manda, y por<br />

<strong>con</strong>sigui<strong>en</strong>te pagan los costos <strong>de</strong> transporte. Sea p j el precio <strong>en</strong> orig<strong>en</strong> <strong>de</strong> <strong>una</strong> unidad <strong>de</strong> producto <strong>en</strong> el<br />

c<strong>en</strong>tro . Entonces los <strong>con</strong>sumidores <strong>en</strong> el nodo i compran <strong>en</strong> el c<strong>en</strong>tro k , si <br />

, .<br />

En este caso, la utilidad u ij <strong>de</strong> un <strong>con</strong>sumidor vi<strong>en</strong>e dada <strong>en</strong> términos <strong>de</strong>l valor negativo <strong>de</strong>l costo total<br />

(precio <strong>de</strong>l producto más coste <strong>de</strong> transporte) expresada por:<br />

u d = - ( p + t ( d )) .<br />

( )<br />

ij ij<br />

j i ij<br />

Sean p min el precio mínimo <strong>de</strong> v<strong>en</strong>ta al que se pue<strong>de</strong> ofertar el producto, y p prod el coste <strong>de</strong> producción<br />

<strong>de</strong>l producto, pmin > pprod<br />

. Asumi<strong>en</strong>do que estos precios son comunes a todos los c<strong>en</strong>tros, pj= pmin<br />

, la firma <strong>en</strong>trante ti<strong>en</strong>e máxima utilidad para los <strong>con</strong>sumidores <strong>en</strong> un nodo i si:<br />

: , es <strong>de</strong>cir, si la distancia a algún nuevo c<strong>en</strong>tro no es mayor que la mínima distancia a algún c<strong>en</strong>tro<br />

preexist<strong>en</strong>te.<br />

e<br />

Para este ejemplo, la distancia umbral es D i = di( F ) . Luego, la ganancia <strong>de</strong> la firma <strong>en</strong>trante<br />

<strong>con</strong>seguida <strong>de</strong> los <strong>con</strong>sumidores <strong>en</strong> el nodo i es:<br />

<br />

<br />

<br />

0 <br />

1<br />

Sean M X<br />

2<br />

= { i : di( X) < Di}<br />

y M X = { i : di( X) = Di}<br />

los <strong>con</strong>juntos <strong>de</strong> nodos don<strong>de</strong> los<br />

<strong>con</strong>sumidores se <strong>en</strong>cu<strong>en</strong>tran a <strong>una</strong> distancia m<strong>en</strong>or o igual a la distancia umbral, respectivam<strong>en</strong>te, y pnet<br />

= pmin - p prod el precio neto. Entonces, el problema <strong>de</strong> maximización <strong>de</strong> la ganancia para la firma<br />

<strong>en</strong>trante es:<br />

:⊂, || <br />

<br />

Para este esc<strong>en</strong>ario el problema <strong>de</strong> maximización <strong>de</strong> la ganancia para la firma <strong>en</strong>trante es equival<strong>en</strong>te a<br />

maximizar la cuota <strong>de</strong> mercado:<br />

.<br />

41<br />

<br />

:⊂, || <br />

<br />

A difer<strong>en</strong>cia <strong>de</strong>l mo<strong>de</strong>lo MAXCAP (Serra y Revelle, 1995) que sólo <strong>con</strong>si<strong>de</strong>ra a un <strong>con</strong>junto <strong>de</strong> nodos<br />

como candidatos a <strong>localización</strong>, aquí se <strong>con</strong>si<strong>de</strong>ran nodos y puntos <strong>en</strong> el interior <strong>de</strong> los tramos.<br />

Ejemplo 2. Un mo<strong>de</strong>lo tipo-Huff binario.<br />

Sea a j la calidad <strong>de</strong>l c<strong>en</strong>tro , don<strong>de</strong> la calidad <strong>de</strong> un c<strong>en</strong>tro pue<strong>de</strong> estar dada por su capacidad, su<br />

tamaño, o por alg<strong>una</strong> aglomeración <strong>de</strong> características <strong>de</strong>l c<strong>en</strong>tro (Plastia, 2001). La utilidad u ij que<br />

obti<strong>en</strong>e un <strong>con</strong>sumidor por un c<strong>en</strong>tro vi<strong>en</strong>e <strong>de</strong>terminada <strong>en</strong> términos <strong>de</strong> la atracción que el <strong>con</strong>sumidor<br />

percibe <strong>de</strong> los c<strong>en</strong>tros, dada por el coci<strong>en</strong>te <strong>de</strong>l parámetro calidad y el tiempo <strong>de</strong> <strong>de</strong>splazami<strong>en</strong>to al<br />

c<strong>en</strong>tro:<br />

a j<br />

uij( d ij ) = .<br />

t ( d )<br />

<br />

i ij


Entonces los <strong>con</strong>sumidores <strong>en</strong> el nodo i compran <strong>en</strong> el c<strong>en</strong>tro k , si la atracción que percib<strong>en</strong> por este<br />

c<strong>en</strong>tro es mayor que la atracción que percib<strong>en</strong> por cualquier otro c<strong>en</strong>tro, es <strong>de</strong>cir, si<br />

, La máxima atracción percibida por los <strong>con</strong>sumidores <strong>en</strong> i <strong>de</strong> los c<strong>en</strong>tros preexist<strong>en</strong>tes vi<strong>en</strong>e dada por:<br />

, Si todos los nuevos c<strong>en</strong>tros ti<strong>en</strong>e la misma calidad, esto es, aj= a para j = q + 1,<br />

q + 2 ,...,q + r , la<br />

firma <strong>en</strong>trante sólo pue<strong>de</strong> capturar el mercado <strong>de</strong> los <strong>con</strong>sumidores <strong>en</strong> i si para alguno <strong>de</strong> sus c<strong>en</strong>tros k<br />

se verifica que su atracción es mayor que la máxima atracción <strong>de</strong> los c<strong>en</strong>tros preexist<strong>en</strong>tes:<br />

<br />

Como la función distancia t i es creci<strong>en</strong>te, la <strong>de</strong>sigualdad anterior es equival<strong>en</strong>te a:<br />

<br />

<br />

<br />

<br />

Luego, si X es el <strong>con</strong>junto <strong>de</strong> localizaciones para la firma <strong>en</strong>trante, <strong>en</strong>tonces dicha firma sólo pue<strong>de</strong><br />

capturar mercado <strong>de</strong> los <strong>con</strong>sumidores <strong>en</strong> el nodo i si:<br />

<br />

<br />

<br />

<br />

- 1 a<br />

Para este ejemplo, la distancia umbral es D i = ti<br />

( ) , y la cuota <strong>de</strong> mercado que la firma<br />

e<br />

Ai( F )<br />

<strong>en</strong>trante <strong>con</strong>sigue <strong>en</strong> i es:<br />

<br />

0 <br />

Como <strong>en</strong> el ejemplo anterior, se pued<strong>en</strong> <strong>de</strong>finir los <strong>con</strong>juntos <strong>de</strong> nodos don<strong>de</strong> los <strong>con</strong>sumidores se<br />

1<br />

<strong>en</strong>cu<strong>en</strong>tran a <strong>una</strong> distancia m<strong>en</strong>or o igual a la distancia umbral. Sean M X = { i : di( X) < Di}<br />

y M<br />

= { i : di( X ) = D i}<br />

, <strong>en</strong>tonces el problema <strong>de</strong> maximización <strong>de</strong> la cuota <strong>de</strong> mercado para la firma<br />

<strong>en</strong>trante es:<br />

: ⊂, || <br />

<br />

Observemos que este problema es equival<strong>en</strong>te al problema <strong>de</strong> maximización <strong>de</strong> la cuota <strong>de</strong> mercado <strong>de</strong>l<br />

Ejemplo 1, si bi<strong>en</strong> los niveles umbrales son distintos.<br />

Acabamos <strong>de</strong> ver que la firma <strong>en</strong>trante t<strong>en</strong>drá máxima utilidad para los <strong>con</strong>sumidores <strong>en</strong> un nodo i si se<br />

<strong>en</strong>cu<strong>en</strong>tran a <strong>una</strong> distancia <strong>de</strong> algún nuevo c<strong>en</strong>tro a m<strong>en</strong>or o igual que <strong>una</strong> distancia umbral, que <strong>de</strong>p<strong>en</strong><strong>de</strong><br />

<strong>de</strong> cada mo<strong>de</strong>lo. En lo que sigue vamos a analizar este tipo <strong>de</strong> problemas cuando el objetivo es<br />

maximizar la cuota <strong>de</strong> mercado.<br />

4. LOCALIZACIONES ÓPTIMAS<br />

Vamos a ver que los resultados obt<strong>en</strong>idos para el <strong>de</strong>l problema <strong>de</strong> <strong>localización</strong> <strong>de</strong> un c<strong>en</strong>tro <strong>con</strong> costos <strong>de</strong><br />

transportación lineal <strong>en</strong> Pelegrín y otros (2006), se exti<strong>en</strong>d<strong>en</strong> al problema <strong>de</strong> <strong>localización</strong> <strong>de</strong> varios<br />

c<strong>en</strong>tros <strong>con</strong> costos <strong>de</strong> transportación creci<strong>en</strong>tes. Para ello realizaremos el análisis <strong>de</strong> la función cuota <strong>de</strong><br />

mercado F( X ) fijando todas las localizaciones excepto <strong>una</strong> <strong>de</strong> ellas.<br />

<br />

42<br />

2<br />

X


Propiedad 1. La función F( X ) es <strong>con</strong>stante si X \ {} x i es fijo y x i varía <strong>en</strong> algún segm<strong>en</strong>to isobierto.<br />

Demostración. Supongamos que todas las localizaciones <strong>de</strong> X están fijadas excepto x i , la cual se<br />

<strong>con</strong>si<strong>de</strong>ra que varía <strong>en</strong> algún segm<strong>en</strong>to iso-abierto <strong>con</strong>t<strong>en</strong>ido <strong>en</strong> algún tramo [ uv. , ] De la <strong>de</strong>finición <strong>de</strong><br />

43<br />

1 2<br />

M X y M X cambian cuando x i varía<br />

segm<strong>en</strong>to iso-abierto se <strong>de</strong>spr<strong>en</strong><strong>de</strong> que ninguno <strong>de</strong> los <strong>con</strong>juntos<br />

<strong>en</strong> el segm<strong>en</strong>to iso-abierto. Por lo tanto, F( X ) es <strong>una</strong> función <strong>con</strong>stante si x i varía <strong>en</strong> algún segm<strong>en</strong>to<br />

iso-abierto <strong>de</strong> [ uv, , ] y las otras localida<strong>de</strong>s están fijas. W<br />

La Figura 3, muestra la gráfica <strong>de</strong> F( X ) cuando X \ {} x i es fijo y x i varía <strong>en</strong> un tramo <strong>de</strong> la <strong>red</strong>.<br />

Como <strong>una</strong> <strong>con</strong>secu<strong>en</strong>cia <strong>de</strong>l la propiedad anterior, el problema <strong>de</strong> maximizar F( X ) se <strong>red</strong>uce a un<br />

problema <strong>de</strong> optimización discreta. El sigui<strong>en</strong>te corolario <strong>de</strong>termina qui<strong>en</strong>es son los candidatos a<br />

<strong>localización</strong> óptima.<br />

Figura 3<br />

Corolario 1. Los candidatos a maximizar la función F( X ) son los nodos, los puntos isodistantes <strong>en</strong> el<br />

interior <strong>de</strong> cada tramo y un punto por cada segm<strong>en</strong>to iso-abierto.<br />

El número total <strong>de</strong> candidatos a <strong>localización</strong> es n + 2ISO + m , don<strong>de</strong> n es el número <strong>de</strong> nodos,<br />

ISO es el número puntos isodistantes <strong>en</strong> el interior <strong>de</strong> cada tramo y m es el número <strong>de</strong> tramos <strong>de</strong> la <strong>red</strong><br />

(véase Figura 4).<br />

Figura 4<br />

Veamos que es posible <strong>red</strong>ucir el número total <strong>de</strong> candidatos a <strong>localización</strong>:<br />

a) Los puntos isodistantes-izquierdos e isodistantes-<strong>de</strong>rechos pued<strong>en</strong> ser removidos <strong>de</strong> la lista inicial <strong>de</strong><br />

candidatos, ya que siempre existe un punto la izquierda o la <strong>de</strong>recha, respectivam<strong>en</strong>te, <strong>en</strong> el que la<br />

cuota <strong>de</strong> mercado es mayor o igual al que se obti<strong>en</strong>e <strong>en</strong> el punto isodistante. Por <strong>con</strong>sigui<strong>en</strong>te, el<br />

*<br />

*<br />

número <strong>de</strong> candidatos se <strong>red</strong>uce a n + ISO + m + ISO , don<strong>de</strong> ISO es el número <strong>de</strong> puntos<br />

isodistantes-mixtos <strong>en</strong> el interior <strong>de</strong> los tramos.<br />

b) Si 1 para todo i , los <strong>con</strong>sumidores <strong>en</strong> el nodo i compran <strong>en</strong> un nuevo c<strong>en</strong>tro si este se <strong>en</strong>cu<strong>en</strong>tra<br />

a <strong>una</strong> distancia m<strong>en</strong>or o igual que sus competidores. Los candidatos que se <strong>en</strong>cu<strong>en</strong>tran <strong>en</strong> los segm<strong>en</strong>tos<br />

iso-abiertos pued<strong>en</strong> ser remplazados por uno <strong>de</strong> los extremos <strong>de</strong>l segm<strong>en</strong>to iso-abierto, sin que disminuya<br />

la cuota <strong>de</strong> mercado. Entonces, los puntos <strong>de</strong> los segm<strong>en</strong>tos iso-abiertos pued<strong>en</strong> ser removidos <strong>de</strong> la lista<br />

inicial <strong>de</strong> candidatos, <strong>red</strong>uciéndose el número <strong>de</strong> candidatos a n + ISO , como se ilustra <strong>en</strong> la Figura 5.


Figura 5<br />

c) Si <br />

para todo i , la mayoría <strong>de</strong> los <strong>con</strong>sumidores <strong>en</strong> el nodo i prefier<strong>en</strong> comprar <strong>en</strong> un c<strong>en</strong>tro<br />

<br />

preexist<strong>en</strong>te que <strong>en</strong> un nuevo c<strong>en</strong>tro. Entonces, si x es un punto isodistante mixto, pue<strong>de</strong> comprobarse<br />

que se obti<strong>en</strong>e <strong>una</strong> cuota <strong>de</strong> mercado mayor o igual que <strong>en</strong> x , si sustituimos x por un punto próximo a<br />

*<br />

x situado a la izquierda o a la <strong>de</strong>recha <strong>de</strong> x . Así, los puntos ISO pued<strong>en</strong> ser eliminados y por<br />

<strong>con</strong>sigui<strong>en</strong>te el número <strong>de</strong> candidatos se <strong>red</strong>uce a n + ISO + m , como se muestra <strong>en</strong> la Figura 6.<br />

Figura 6<br />

5. GENERACIÓN DE LA LISTA DE CANDIDATOS<br />

En esta sección se pres<strong>en</strong>ta un algoritmo que g<strong>en</strong>era la lista completa <strong>de</strong> candidatos L , para <strong>una</strong> <strong>red</strong> que<br />

<strong>con</strong>ti<strong>en</strong>e nodos población y nodos cruce. Los elem<strong>en</strong>tos <strong>de</strong> la lista son puntos <strong>de</strong> la <strong>red</strong> <strong>de</strong>finidos <strong>de</strong> la<br />

forma ( iuvx , , , ) , si<strong>en</strong>do x la distancia al nodo u <strong>de</strong>l punto isodistante relativo al nodo i que se<br />

<strong>en</strong>cu<strong>en</strong>tra <strong>en</strong> el tramo [ uv. , ]<br />

Algoritmo<br />

Paso 1. G<strong>en</strong>eración <strong>de</strong> la lista <strong>de</strong> puntos isodistantes:<br />

Determinar la distancia umbral D i para cada nodo población i . Iniciar la lista <br />

Obt<strong>en</strong>er las matrices <strong>de</strong> distancias y p<strong>red</strong>ecesores <strong>en</strong>tre cada par <strong>de</strong> nodos, D = ( d ij ) y P =<br />

( p ij ) .<br />

Para cada nodo población i , hacer:<br />

Para cada elem<strong>en</strong>to u <strong>en</strong> la fila i tal que , hacer:<br />

Si d iu = D i : El nodo u es un punto isodistante, agregar el nodo u a la lista<br />

, , , 0<br />

En otro caso ( iu<br />

d < D i ): Si u es un nodo población, examinar los adyac<strong>en</strong>tes <strong>de</strong> u ,<br />

que sean nodos población y que no estén <strong>en</strong> el camino más corto <strong>de</strong> i a u , es <strong>de</strong>cir,<br />

los nodos población <strong>en</strong> el <strong>con</strong>junto:<br />

Γ\, <br />

Para cada , hacer:<br />

Si ≥2 , el punto , tal que l ux = Di -<br />

d iu es isodistante, agregar dicho punto a la lista, ,,, Paso 2. Agrupar y ord<strong>en</strong>ar los puntos isodistantes <strong>de</strong> la lista L , para cada tramo <strong>de</strong> la<br />

<strong>red</strong>.<br />

Paso 3. Seleccionar el valor <strong>de</strong> la proporción q i .<br />

Paso 4. Si el valor <strong>de</strong> q i = 1, la lista L está completa.<br />

En otro caso, modificamos la lista L <strong>en</strong> cada segm<strong>en</strong>to iso-abierto, <strong>de</strong>p<strong>en</strong>di<strong>en</strong>do <strong>de</strong>l valor <strong>de</strong> : 0 <br />

, ∀, (Figura 7):<br />

<br />

a) Si los extremos <strong>de</strong>l segm<strong>en</strong>to iso-abierto son: un nodo población y un punto isodistante, o<br />

dos puntos isodistantes; si hay un nodo cruce <strong>en</strong> dicho segm<strong>en</strong>to, reemplazar el punto <strong>de</strong>l<br />

extremo <strong>de</strong>recho <strong>de</strong>l segm<strong>en</strong>to por el nodo cruce <strong>en</strong> la lista L .<br />

En otro caso, reemplazar el punto <strong>de</strong>l extremo <strong>de</strong>recho <strong>de</strong>l segm<strong>en</strong>to por el punto medio <strong>de</strong><br />

dicho segm<strong>en</strong>to <strong>en</strong> la lista L .<br />

44


) Si los extremos <strong>de</strong>l segm<strong>en</strong>to iso-abierto son: un punto isodistante y un nodo población; si<br />

hay un nodo cruce <strong>en</strong> dicho segm<strong>en</strong>to, agregar a la lista L el nodo cruce.<br />

En otro caso, agregar a la lista L el punto medio <strong>de</strong> dicho segm<strong>en</strong>to.<br />

Figura 7<br />

Si <br />

1, (Figura 8):<br />

a) Si los extremos <strong>de</strong>l segm<strong>en</strong>to iso-abierto son: un nodo población y un punto isodistante (o<br />

viceversa), o dos puntos isodistantes; si hay un nodo cruce <strong>en</strong> dicho segm<strong>en</strong>to, agregar a la<br />

lista L el nodo cruce.<br />

En otro caso, agregar a la lista L el punto medio <strong>de</strong> dicho segm<strong>en</strong>to.<br />

Figura 8<br />

Una vez obt<strong>en</strong>ida la lista <strong>de</strong> candidatos, se <strong>en</strong><strong>con</strong>traran las soluciones óptimas mediante la resolución <strong>de</strong><br />

un problema <strong>de</strong> programación lineal <strong>en</strong>tera mixta.<br />

6. FORMULACIÓN COMO PROGRAMACIÓN LINEAL ENTERA MIXTA<br />

Consi<strong>de</strong>rando las i<strong>de</strong>as c<strong>en</strong>trales <strong>de</strong> las formulaciones para resolver problemas <strong>de</strong> <strong>localización</strong> discreta<br />

hechas <strong>en</strong> Serra y ReVelle (1995) y Fernán<strong>de</strong>z y otros (2007). Sean L el <strong>con</strong>junto <strong>de</strong> todos los<br />

candidatos a <strong>localización</strong>, : el <strong>con</strong>junto <strong>de</strong> candidatos que se <strong>en</strong>cu<strong>en</strong>tran a <strong>una</strong><br />

distancia m<strong>en</strong>or a la distancia umbral, : el <strong>con</strong>junto <strong>de</strong> candidatos que se<br />

<strong>en</strong>cu<strong>en</strong>tran a <strong>una</strong> distancia igual a la distancia umbral y : el <strong>con</strong>junto <strong>de</strong> nodos<br />

i que se <strong>en</strong>cu<strong>en</strong>tran a <strong>una</strong> distancia <strong>de</strong> un candidato no mayor a su distancia umbral.<br />

Las variables <strong>de</strong> <strong>localización</strong> son:<br />

1 <br />

<br />

0 <br />

Las variables <strong>de</strong> asignación total <strong>de</strong> la <strong>de</strong>manda son:<br />

1 ó <br />

<br />

0 <br />

Las variables <strong>de</strong> asignación parcial <strong>de</strong> la <strong>de</strong>manda son:<br />

1 <br />

<br />

0 <br />

Entonces el problema <strong>de</strong> <strong>localización</strong> se pue<strong>de</strong> formular como:<br />

<br />

Sujeto a:<br />

<br />

<br />

<br />

1, <br />

, <br />

, 2<br />

3<br />

<br />

45<br />

1


0, 1 5<br />

0, 0 6<br />

don<strong>de</strong> el primer <strong>con</strong>junto <strong>de</strong> restricciones impone que la <strong>de</strong>manda <strong>de</strong> los <strong>con</strong>sumidores <strong>en</strong> el nodo i<br />

pue<strong>de</strong> ser capturada total o parcialm<strong>en</strong>te por la nueva firma o bi<strong>en</strong> por sus competidores si la suma <strong>de</strong> las<br />

variables <strong>de</strong> asignación es cero. El segundo <strong>con</strong>junto <strong>de</strong> restricciones obliga que sólo se pue<strong>de</strong> asignar la<br />

<br />

<strong>de</strong>manda totalm<strong>en</strong>te <strong>de</strong> i a , si se localiza un c<strong>en</strong>tro <strong>en</strong> j ( x j = 1) . El tercer <strong>con</strong>junto <strong>de</strong><br />

restricciones impone que sólo se pue<strong>de</strong> asignar parcialm<strong>en</strong>te la <strong>de</strong>manda <strong>de</strong> i a la firma <strong>en</strong>trante si se<br />

localiza algún c<strong>en</strong>tro <strong>en</strong> . La cuarta restricción obliga a que el número <strong>de</strong> c<strong>en</strong>tros <strong>de</strong> la nueva firma<br />

sea r . El quinto <strong>con</strong>junto <strong>de</strong> restricciones impone que las variables <strong>de</strong> <strong>localización</strong> x j sean binarias. El<br />

último <strong>con</strong>junto <strong>de</strong> restricciones obligan a que las variable <strong>de</strong> asignación y ij y z i sean no negativas, <strong>en</strong><br />

vez <strong>de</strong> binarias. Esto se hace para <strong>red</strong>ucir el número <strong>de</strong> variables binarias <strong>en</strong> el mo<strong>de</strong>lo y permitir que el<br />

optimizador resuelva problemas <strong>de</strong> mayor tamaño. Ello no altera el planteami<strong>en</strong>to, ya que si se fijan los<br />

valores <strong>de</strong> las variables x j <strong>en</strong> la formulación pres<strong>en</strong>tada, se pue<strong>de</strong> comprobar que siempre hay<br />

soluciones óptimas <strong>con</strong> valores binarios para las restantes variables.<br />

7. EXPERIMENTO COMPUTACIONAL<br />

Para analizar la s<strong>en</strong>sibilidad respecto <strong>de</strong> los parámetros r y q , utilizaremos datos <strong>de</strong> la <strong>red</strong> <strong>de</strong> transporte<br />

i<br />

<strong>de</strong> la Región <strong>de</strong> Murcia, España. Esta <strong>red</strong> <strong>con</strong>ti<strong>en</strong>e 84 nodos (los primeros 45 son nodos población y 39<br />

son nodos cruce) y 132 tramos, como se muestra <strong>en</strong> Figura 9. La <strong>de</strong>manda w i , i = 1, 2,..., 45 <strong>de</strong> los<br />

nodos población correspon<strong>de</strong> al c<strong>en</strong>so <strong>de</strong> 2007 (Tabla 1) <strong>de</strong> acuerdo al C<strong>en</strong>tro Regional <strong>de</strong> Estadística <strong>de</strong><br />

Murcia.<br />

Se <strong>con</strong>si<strong>de</strong>ran cuatro <strong>con</strong>juntos <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes situados <strong>en</strong> las ciuda<strong>de</strong>s <strong>con</strong> mayor población:<br />

1<br />

F = { Murcia } ,<br />

2<br />

F = { Murcia, Cartag<strong>en</strong>a } ,<br />

3<br />

F = { Murcia, Cartag<strong>en</strong>a, Lorca } ,<br />

4<br />

F = { Murcia, Cartag<strong>en</strong>a, Lorca, Molina <strong>de</strong> Segura}<br />

.<br />

La distancia umbral asociada a los nodos población está dada por: : i = 1,2,...,45.<br />

1 1 3<br />

El rango <strong>de</strong> valores que hemos tomado para la proporción es q i = { 0, , , , 1<br />

4 2 4 } , i = 1,2,.., 45 y el<br />

número <strong>de</strong> nuevos c<strong>en</strong>tros a localizar es r = 1,2,3,4.<br />

Se han <strong>con</strong>si<strong>de</strong>rado 16 esc<strong>en</strong>arios distintos correspondi<strong>en</strong>tes a las difer<strong>en</strong>tes combinaciones <strong>en</strong>tre número<br />

<strong>de</strong> c<strong>en</strong>tros a localizar y número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes. Utilizando el optimizador Xpress-MP (Xpresse<br />

MP, 2007) se resolvieron 80 problemas <strong>con</strong> <strong>patrón</strong> <strong>de</strong> <strong>elección</strong> binario que resultan <strong>de</strong> variar , F y<br />

r . El algoritmo que g<strong>en</strong>era la lista completa <strong>de</strong> candidatos ha sido diseñado <strong>en</strong> Borland C, versión 5.02<br />

y utilizando un equipo Dell Latitu<strong>de</strong> D810, <strong>con</strong> procesador Intel P<strong>en</strong>tium M a 2.26 GHz y 1 GB <strong>de</strong> Ram.<br />

Las tablas 2 a 5 muestran los resultados para cada <strong>con</strong>junto <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes.<br />

La primera columna <strong>de</strong> cada tabla indica el número <strong>de</strong> nuevos c<strong>en</strong>tros <strong>de</strong> la firma <strong>en</strong>trante, la proporción<br />

se refleja <strong>en</strong> la segunda columna, el número <strong>de</strong> candidatos vi<strong>en</strong>e dado <strong>en</strong> la tercera columna, las<br />

localizaciones óptimas se repres<strong>en</strong>tan <strong>en</strong> la cuarta columna, don<strong>de</strong> un número <strong>en</strong>tre 1 y 45 correspon<strong>de</strong> a<br />

nodo población, <strong>en</strong>tre 46 y 84 correspon<strong>de</strong> a un nodo cruce y un número mayor a un punto isodistante.<br />

En la quinta columna se muestra la cuota <strong>de</strong> mercada capturada, los porc<strong>en</strong>tajes <strong>de</strong> cuota capturados total<br />

y parcialm<strong>en</strong>te están expresados <strong>en</strong> la sexta y séptima columnas, respectivam<strong>en</strong>te. Finalm<strong>en</strong>te, el tiempo<br />

<strong>en</strong> segundos para g<strong>en</strong>erar la lista completa <strong>de</strong> candidatos y <strong>en</strong> obt<strong>en</strong>er las soluciones mediante el<br />

optimizador Xpress-MP, se muestran <strong>en</strong> las columnas ocho y nueve respectivam<strong>en</strong>te.<br />

46<br />

4


Figura 9.<br />

Tabla 1: Nodos Poblacionales<br />

1<br />

A modo <strong>de</strong> ilustración, la Figura 10 muestra las localizaciones óptimas <strong>de</strong> 4 nuevos c<strong>en</strong>tros para q i =<br />

4<br />

, si<strong>en</strong>do Murcia, Cartag<strong>en</strong>a y Lorca, los c<strong>en</strong>tros preexist<strong>en</strong>tes. Los nuevos c<strong>en</strong>tros están localizados <strong>en</strong> el<br />

nodo población Murcia y los puntos isodistantes 194, 240 y 247. La cuota <strong>de</strong> mercado capturada es <strong>de</strong><br />

47


663,924.75. Esta cantidad correspon<strong>de</strong> a un 39.91 % <strong>de</strong> la <strong>de</strong>manda total, la cual es totalm<strong>en</strong>te capturada<br />

<strong>en</strong> los correspondi<strong>en</strong>tes nodos por los c<strong>en</strong>tros localizados <strong>en</strong> los puntos isodistantes 194, 240 y 247, y a<br />

un 7.77% <strong>de</strong> la <strong>de</strong>manda total, la cual es parcialm<strong>en</strong>te capturada <strong>en</strong> los correspondi<strong>en</strong>tes nodos por el<br />

c<strong>en</strong>tro localizado <strong>en</strong> Murcia. El tiempo <strong>de</strong> cómputo es <strong>de</strong> 8.266 segundos para g<strong>en</strong>erar la lista completa<br />

<strong>de</strong> candidatos y el optimizador Xpress-MP utilizó 20.124 segundos para <strong>en</strong><strong>con</strong>trar las soluciones óptimas<br />

<strong>de</strong>l mo<strong>de</strong>lo.<br />

1<br />

Tabla 2: F = { Murcia}<br />

2<br />

Tabla 3: F =<br />

{ Murcia, Cartag<strong>en</strong>a}<br />

48


3<br />

Tabla 4: F = { Murcia, Catag<strong>en</strong>a, Lorca}<br />

4<br />

Tabla 5: F = { Murcia, Cartag<strong>en</strong>a, Lorca, Molina <strong>de</strong> Segura}<br />

8. ANÁLISIS DE SENSIBILIDAD<br />

Respecto al número <strong>de</strong> candidatos <strong>en</strong> los difer<strong>en</strong>tes esc<strong>en</strong>arios, se observa lo sigui<strong>en</strong>te: El número <strong>de</strong><br />

candidatos a <strong>localización</strong> óptima <strong>de</strong>crece al aum<strong>en</strong>tar el número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes para cada valor<br />

<strong>de</strong> , Para <br />

el número <strong>de</strong> candidatos es mayor que para 0 <br />

, , pero esa difer<strong>en</strong>cia <strong>en</strong><br />

<br />

número va <strong>de</strong>creci<strong>en</strong>do porc<strong>en</strong>tualm<strong>en</strong>te al aum<strong>en</strong>tar el número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes. Así, <strong>en</strong> la Tabla<br />

767 - 446<br />

548 - 349<br />

2 la difer<strong>en</strong>cia porc<strong>en</strong>tual es ´ 100 = 71.97% y <strong>en</strong> la Tabla 5 es ´ 100 =<br />

446<br />

349<br />

66.56% . Para 1, el número <strong>de</strong> candidatos también es m<strong>en</strong>or que para 0 <br />

, , pero ahora la<br />

49


difer<strong>en</strong>cia <strong>en</strong> número va creci<strong>en</strong>do porc<strong>en</strong>tualm<strong>en</strong>te al aum<strong>en</strong>tar el número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes. En<br />

446 - 378<br />

329 - 274<br />

la Tabla 2 la difer<strong>en</strong>cia es ´ 100 = 15.24 % y <strong>en</strong> la Tabla 5 es ´ 100 = 16.71<br />

446<br />

329<br />

%.<br />

Figura 10.<br />

En cuanto a las localizaciones óptimas, resulta que son bastante estables <strong>en</strong> relación a la variación <strong>de</strong> .<br />

Los cambios <strong>de</strong> <strong>localización</strong> se produc<strong>en</strong> al aum<strong>en</strong>tar significativam<strong>en</strong>te, y alg<strong>una</strong>s localizaciones<br />

permanec<strong>en</strong> óptimas para un gran rango <strong>de</strong> valores. Para valores <strong>de</strong> elevados las ciuda<strong>de</strong>s más<br />

pobladas suel<strong>en</strong> ser las localizaciones óptimas y esté hecho se <strong>con</strong>stata más cuanto mayor es el número<br />

<strong>de</strong> competidores.<br />

Como era <strong>de</strong> esperar, la cuota <strong>de</strong> mercado disminuye <strong>con</strong>forme al número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes es<br />

mayor. Sin embargo, se observa que esta disminución es mayor para valores <strong>de</strong> pequeños, sobre todo<br />

si el número <strong>de</strong> c<strong>en</strong>tros localizados es bajo. Si los valores <strong>de</strong> son altos y se localizan muchos c<strong>en</strong>tros,<br />

se observa que la cuota <strong>de</strong> mercado <strong>de</strong>crece muy l<strong>en</strong>tam<strong>en</strong>te <strong>con</strong> el número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes. En<br />

particular, para 1 y cuatro nuevos c<strong>en</strong>tros la cuota <strong>de</strong> mercado obt<strong>en</strong>ida por la firma <strong>en</strong>trante no se<br />

ve afectada al aum<strong>en</strong>tar el número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes.<br />

A<strong>de</strong>más, <strong>en</strong> cada uno <strong>de</strong> los 16 esc<strong>en</strong>arios, se observa que para valores pequeños <strong>de</strong> sólo <strong>en</strong> muy<br />

pocos nodos la <strong>de</strong>manda es compartida por la firma <strong>en</strong>trante y sus competidores. Para valores gran<strong>de</strong>s <strong>de</strong><br />

ocurre lo <strong>con</strong>trario.<br />

El tiempo máximo <strong>de</strong> cómputo para g<strong>en</strong>erar la lista <strong>de</strong> candidatos y <strong>en</strong><strong>con</strong>trar las localizaciones óptimas<br />

se alcanza para <br />

. El tiempo para g<strong>en</strong>erar la lista <strong>de</strong> candidatos se increm<strong>en</strong>ta al aum<strong>en</strong>tar el número<br />

<br />

<strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes para cada valor pero ocurre lo <strong>con</strong>trario para <strong>en</strong><strong>con</strong>trar las localizaciones<br />

óptimas. Los tiempos que se requier<strong>en</strong> para <strong>en</strong><strong>con</strong>trar las localizaciones óptimas disminuy<strong>en</strong><br />

significativam<strong>en</strong>te, mi<strong>en</strong>tras que los tiempos que se requier<strong>en</strong> para g<strong>en</strong>erar la lista <strong>de</strong> candidatos se<br />

increm<strong>en</strong>tan <strong>en</strong> m<strong>en</strong>or proporción <strong>con</strong>forme aum<strong>en</strong>ta el número <strong>de</strong> c<strong>en</strong>tros preexist<strong>en</strong>tes.<br />

50


9. CONCLUSIONES Y FUTURAS INVESTIGACIONES<br />

Se ha visto que bajo ciertas <strong>con</strong>diciones los mo<strong>de</strong>los <strong>de</strong> precios <strong>en</strong> orig<strong>en</strong> y tipo Huff, el número <strong>de</strong><br />

candidatos a <strong>localización</strong> para <strong>una</strong> nueva firma que <strong>de</strong>sea <strong>en</strong>trar <strong>en</strong> compet<strong>en</strong>cia <strong>con</strong> otras ya exist<strong>en</strong>tes<br />

<strong>con</strong> el objetivo <strong>de</strong> maximizar su cuota <strong>de</strong> mercado, se <strong>red</strong>uce a un <strong>con</strong>junto finito <strong>de</strong> puntos, y mediante<br />

su correspondi<strong>en</strong>te formulación como un problema <strong>de</strong> programación lineal <strong>en</strong>tera mixta se pued<strong>en</strong><br />

<strong>en</strong><strong>con</strong>trar sus soluciones óptimas. A<strong>de</strong>más, cuando ocurre que la <strong>de</strong>manda <strong>de</strong> los <strong>con</strong>sumidores se ti<strong>en</strong>e<br />

que repartir <strong>en</strong>tre las firmas que compit<strong>en</strong> por el mercado, un análisis <strong>de</strong>l parámetro <strong>de</strong> la proporción <strong>de</strong><br />

la <strong>de</strong>manda nos permite <strong>de</strong>terminar quiénes serán los candidatos a <strong>localización</strong>.<br />

Como trabajos futuros <strong>de</strong> investigación, será <strong>de</strong> interés analizar bajo qué <strong>con</strong>diciones otros tipo <strong>de</strong><br />

mo<strong>de</strong>los <strong>de</strong> <strong>localización</strong> competitiva, por ejemplo mo<strong>de</strong>los <strong>con</strong> funciones <strong>de</strong> <strong>elección</strong> no <strong>de</strong>terministas, o<br />

mo<strong>de</strong>los don<strong>de</strong> la <strong>de</strong>manda sea variable, la búsqueda <strong>de</strong> los candidatos a <strong>localización</strong> se <strong>red</strong>uzca a un<br />

<strong>con</strong>junto finito <strong>de</strong> puntos <strong>de</strong> la <strong>red</strong>, los cuales se <strong>en</strong>cu<strong>en</strong>tr<strong>en</strong> d<strong>en</strong>tro <strong>de</strong> <strong>una</strong> distancia umbral <strong>de</strong> los<br />

<strong>con</strong>sumidores.<br />

Agra<strong>de</strong>cimi<strong>en</strong>tos: Este trabajo ha sido realizado como parte <strong>de</strong>l proyecto ECO-2008-00667/ECON<br />

financiado por el Ministerio <strong>de</strong> Ci<strong>en</strong>cia e Innovación <strong>de</strong> España y <strong>de</strong>l Fondo Europeo para el Desarrollo<br />

Regional (FEDER).<br />

RECEIVED JANUARY 2010<br />

REVISED JUNE 2010<br />

REFERENCIAS<br />

[1] DREZNER T. and EISELT H.A. (2002): Consumers in competitive location mo<strong>de</strong>ls. In: Facility<br />

Location: Application and Theory (Drezner Z. and Hamacher H. Eds.), 151-178. Springer Verlag,<br />

Berlin.<br />

[2] EISELT H.A. and LAPORTE G. (1998): Demand allocation functions. Location Sci<strong>en</strong>ce, 6, 175-<br />

187.<br />

[3] FERNÁNDEZ P., PELEGRÍN B., GARCÍA M.D. and PEETERS P. (2007): A discrete location-price<br />

problem un<strong>de</strong>r the assumption of discriminatory pricing: Formulations and parametric analysis.<br />

European Journal of Operations Research, 179, 1050-1062.<br />

[4] FRANCIS R.L., LOWE T.J. and TAMIR A. (2002): Demand point aggregation for location mo<strong>de</strong>ls.<br />

In: Facility Location: Application and Theory (Drezner Z. and Hamacher H. Eds.), 207-232. Springer<br />

Verlag, Berlin.<br />

[5] GABSZEWICZ J.J. and THISSE J.F. (1992): Location. In: Handbook of Game Theory with<br />

E<strong>con</strong>omic Applications (Aumann R. and Hart S. Eds.), 281-304. Elsevier Sci<strong>en</strong>ce Publishers,<br />

Amsterdam.<br />

[6] GARCIA M.D. FERNÁNDEZ P. and PELEGRÍN B. (2004): On price competitive in location-price<br />

mo<strong>de</strong>ls with spatially separated markets. TOP, 12, 351-374.<br />

[7] HAKIMI L. (1983): On locating new facilities in a competitive <strong>en</strong>vironm<strong>en</strong>t. European Journal of<br />

Operations Research, 12, 29-35.<br />

[8] NAKANISHI M. and COOPER L.G. (1974): Parameter estimate for multiplicative interactive choice<br />

mo<strong>de</strong>l: Least square approach. Journal of Marketing Research, 11, 303-311.<br />

[9] PELEGRÍN B., FERNÁNDEZ P., SUÁREZ R. and GARCÍA M.D. (2006): Single facility location<br />

on a network un<strong>de</strong>r mill and <strong>de</strong>live<strong>red</strong> pricing. IMA Journal of Managem<strong>en</strong>t Mathematics, 17, 373-<br />

385.<br />

51


[10] PLASTRIA F. (2001): Static competitive facility location: An overview of optimization approaches.<br />

European Journal of Operational Research, 129, 461-470.<br />

[11] SERRA D. and REVELLE C. (1995): Competitive location in discrete space. In: Facility Location:<br />

A Survey of Applications and Methods (Drezner Z. Ed.), 367-386. Springer Verlag, Berlin.<br />

[12] SUÁREZ-VEGA R., SANTOS PEÑATE D.R. and DORTA-GONZÁLEZ P. (1994): Discretization<br />

and Resolution of the ( r | X )-Medianoid Problem Involving Quality Criteria. TOP, 12, 111-133.<br />

p<br />

[13] Xpress-MP. (2007). Dash Optimization.<br />

52

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!