viernes, 20 de mayo de 2011

FUNCIONES CONVEXAS Y CÓNCAVAS.

El concepto de convexidad es fundamental en el análisis y resolución de los problemas de optimización.




Sea S Rn, un conjunto convexo y no vacío, y sea f: S R f es una función convexa

en S, si y solo si:

f [ λ x1 + (1 – λ) x2 ] ≤ λ f (x1) + (1 – λ) f (x2)


 
Sea S Rn, un conjunto convexo y no vacío, y sea f: S R, f es una función



estrictamente convexa en S, si:



f[λx1 + (1 – λ) x2 ] < λ f(x1) + (1 – λ) f(x2)



λ ] 0,1 [     y              ∀ x1, x2 ∈ S. con x1 ≠ x2



Sea S ⊆ Rn, un conjunto convexo y no vacío, y sea f: S → R f es una función cóncava en S, si y solo si:


f [ λ x1 + (1-λ) x2 ] ≥ λ f(x1) + (1-λ) f(x2)

∀ λ ∈[0,1] y ∀ x1,x2 ∈ S.


 
Una función es estrictamente cóncava si la desigualdad se verifica en sentido estricto, es decir:
f [ λ x1 + (1-λ) x2 ] > λ f(x1) + (1-λ) f(x2)
λ (0,1) y x1,x2 S. con x1 ≠ x2
Es importante hacer notar que las definiciones que hemos dado con anterioridad no exigen ni la continuidad ni la diferenciabilidad de la función.
Si f(x) es una función convexa en S (convexo y no vacío), entonces la función [-f(x)] es una función cóncava en S.

Prueba:

Si f(x) es una función convexa verifica que:

λ [0,1] y x1,x2 S.
f [ λ x1 + (1-λ) x2 ] ≤ λ f(x1) + (1-λ) f(x2)

si multiplicamos esta expresión por (-1), tenemos

- f [ λ x1 + (1-λ) x2 ] ≥ - λ f(x1) - (1-λ) f(x2)

o lo que es lo mismo:

(-f) [ λ x1 + (1-λ) x2 ] ≥ λ (-f)(x1) + (1-λ) (-f)(x2)

con lo cual (-f) es una función cóncava.

Así por ejemplo, las funciones lineales son cóncavas y convexas a la vez, dado que cumplen la definición de función cóncava y convexa como una igualdad entre los dos miembros de la definición, pero precisamente por este motivo no pueden ser ni estrictamente cóncavas ni convexas.

Por el contrario, la función coseno ( cos(x) ) no es cóncava ni convexa sobre todo su dominio ( R ), pero sin embargo, sobre ciertos subdominios si tiene algunas de estas propiedades. Así, en el dominio [ π/2, 3π/2 ] es un función convexa, mientras que en el dominio [ 3π/2, 5π/2 ] se trata de una función cóncava, y además lo es estrictamente en ambos casos.



 
PROPIEDADES DE LAS FUNCIONES CONVEXAS.


Toda combinación lineal con coeficientes positivos de funciones convexas es una función convexa.


Sea S ⊆ Rn un conjunto convexo y no vacío, y sea f: S → R una función convexa. Entonces el conjunto de nivel inferior Sα = { x ∈ S / f(x) ≤ α }, es un conjunto convexo.


Prueba:


Sean x1 , x2 ∈ Sα , lo que significa que:

f(x1) ≤ α

f(x2) ≤ α


lo que tenemos que probar es que: ∀ λ ∈ [0,1] se verifica que: λ x1 + (1-λ) x2 ∈ Sα , o lo que es lo mismo que : f[λ x1 + (1-λ) x2] ≤ α.


Como ayuda para hacer más compresible esta prueba, definimos:


x° = λ x1 + (1-λ) x2


por lo que


f( x° ) = f[ λ x1 + (1-λ) x2]


entonces:


f[ λ x1 + (1-λ) x2] ≤λ f(x1) + (1-λ) f(x2) ≤ λ α + (1-λ) α = α


lo que significa que


λ x1 + (1-λ) x2 ∈ Sα


que es lo que queríamos probar, que el conjunto Sα es un conjunto convexo.



De igual manera tenemos la siguiente propiedad:


Si f es un función cóncava el conjunto de nivel superior Sα ={x ∈ S/f(x)≥ α}, es un conjunto convexo.


El reciproco de estas dos propiedades no es cierto, es decir, que el conjunto de nivel sea sea convexo, no implica que la función sea convexa (cóncava), aunque esta propiedad se cumple para las funciones cuasiconcavas y cuasiconvexas.




 
CARACTERIZACIONES DE LAS FUNCIONES CONVEXAS.


La aplicación de la definición de convexidad o concavidad a una función puede, en muchas ocasiones, resultar complicado, por lo que se recurre a las caracterizaciones, es decir, a ciertas condiciones que pueden verificar las funciones y que nos permiten clasificar a las funciones en convexas o cóncavas.


Caracterización de funciones de clase C2. Para las funciones que admiten derivadas continuas hasta el segundo orden, podemos utilizar una caracterización basada en el hessiano de la función. En este caso podemos acudir a la siguiente proposición.

Dada una función f: S ⊂ Rn → R, donde S es un conjunto convexo y no vacío, y f’ ∈ C2(S)-función con segunda derivada continua en S-, entonces se cumple que: 


a) f es convexa en S sii se cumple Hf(x) es semidefinida positiva en S.

b) f es cóncava en S sii se cumple que Hf(x) es semidefinida negativa en S.

c) f es estrictamente convexa solamente si Hf(x) es definida positiva en S.

d) f es estrictamente cóncava solamente si Hf(x) es definida negativa en S.



Ejemplo:



F(x) = x2



En primer lugar, y por representación gráfica de la función (una parábola con centro el origen) podremos aplicarle la definición de función convexa:



f [ λ x1 + (1-λ) x2 ] ≤ λ f(x1) + (1-λ) f(x2)



∀ λ ∈[0,1] y ∀ x1,x2 ∈ S



Para probar desigualdades recurrimos a probar su diferencia respecto de cero, es decir, poner todos los componentes en un único miembro de la desigualdad, es decir:


0≤ λ f(x1) + (1-λ) f(x2) -f [ λ x1 + (1-λ) x2 ]



Sustituyendo se tiene:



λ x21 + (1-λ) x22 - (λx1 + (1-λ) x2)2



operando



λ x21 + (1-λ) x22 - 2 λ(1- λ) x1 x2 =



 x21[λ – λ2] + x22 [(1 – λ)2] – 2 λ (1 – λ) x1 x2 =



λ (1 – λ) x12 + (1 – λ) λ x22 - 2 λ (1 – λ) x1 x2 =



λ (1 – λ) [x12 + x22 – 2 x1x2] = λ(1- λ) (x1 – x2)2 ≥ 0



También recurriendo a la segunda derivada:



F'(x) = 2 x

F''(x) = 2 > 0 Definida positiva, por tanto F(x) es convexa.

Incluso se podría decir que es estrictamente convexa.




FUENTE: Universidad de Valencia. http://www.uv.es/~sala/clase03.pdf







jueves, 7 de abril de 2011

TALLER 1. - PROGRAMACIÓN LINEAL

TALLER 1.

2. Se cuentan con 210.000 euros para invertir en valores. Nos recomiendan dos tipos de acciones. Las del tipo 1, que dejan utilidad del 10 y las del tipo 2, que dejan utilidad del 8%. El monto máximo a invertir es de 130.000 euros en las del tipo 1 y como mínimo 60.000 en las del tipo 2. Se espera que la inversión en las del tipo 1 sea menor que el doble de la inversión en 2. ¿Cómo se debe distribuir la inversión para obtener el máximo de utilidad?


Utilidad
Monto a Invertir
Acción tipo 1
10%
≤ 130.000
Acción tipo 2
8%
≥ 60.000

Restricciones

X1 < 2X2
X1 ≤ 130.000
X2 ≥ 60.000
X1 + X2 ≤ 210.000
X1 ; X2 ≥ 0

Función Objetivo

            Z = 0.1X1 + 0.08X2


3. En Postres y Ponques se hacen dos tipos de tortas: Genovesa y Tropical. Cada torta genovesa necesita un cuarto de relleno por cada Kg de bizcocho y produce una utilidad de 250 pesos, mientras que una torta tropical necesita medio kg de relleno por cada kg de bizcocho y produce 400 pesos de utilidad. En Postres y Ponqués se pueden hacer diariamente hasta 150 kg de bizcocho y 50 kg de relleno, aunque por problemas de maquinaria no pueden hacer más de 125 tortas de cada tipo. ¿Cuántas tortas Genovesa y cuantas de Tropical deben vender al día para que sea máxima la utilidad?



Relleno
Bizcocho
Utilidad
Máximo a producir
Genovesa. X1
¼ kg
1
$ 250
125 Und
Torta Tropical. X2
½ kg
1
$ 400
125 Und
Disponibilidad
50 kg
150kg




Restricciones

X1 ≤ 125
X2 ≤ 125
X1(1/4) + X2 (1/2) ≤ 150
X1 + X2  ≤ 50
X1 ; X2 ≥ 0

Función Objetivo
Z = X1(250) + X2(400)

4. Se planea una excursión estudiantil para 400 alumnos. La empresa de transporte tiene 8 buses de 40 puestos y 10 buses de 50 puestos, pero solo dispone de 9 conductores. El alquiler de un bus grande cuesta 80 dólares y el de uno pequeño, 60 dólares. ¿Cuantos buses de cada tipo hay que utilizar para que la excursión resulte lo más económica posible para la escuela?


Disponibilidad
Costo
Bus, 40 puestos X1
40 und
$60/und
Bus, 50 puestos. X2
10 und
$80/und

Restricciones
X1 ≤ 40
X2 ≤ 10
X1 + X2 ≤ 9
X1(40) + X2(50) ≥ 400
X1 ; X2 ≥ 0

Función Objetivo
Z = X1(60) + X2(8)

5. En un taller de vehículos van a trabajar latoneros y mecánicos. Por necesidades de mercado, es necesario que haya mayor o igual número de latoneros que de mecánicos y que  el número de latoneros no supere al doble que el de mecánicos. En total hay disponibles 30 latoneros y 20 mecánicos. La utilidad de la empresa por jornada es de 250 dólares por latonero y 200 dólares por mecánico. ¿Cuántos trabajadores de cada clase deben elegirse para obtener el máximo beneficio y cual es este?


Disponibilidad
Utilidad/jornada
Latoneros
30
250
Mecánicos
20
200

Restricciones
            X1 ≥ X2
            2X1 X2
X1 ; X2 ≥ 0

Función Objetivo
            Z = X1(250) + X2(200)

6. Avianca desea ofertar a lo sumo, 5.000 puestos de dos tipos: T (clase turista) y P (primera clase). La ganancia correspondiente a cada plaza de  tipo T es de 30.000 pesos, mientras que la ganancia del tupo P es de 40.000 pesos. El número de puestos tipo T no puede exceder de 4.500 y el del tipo P, debe ser, como máximo, la tercera parte de los del tipo T que se oferten. Calcular cuántos tienen que ofertarse de cada clase para que las ganancias sean máximas.


Ganancia
Disponibilidad
X1. C.Tipo T
30.000
≤4.500
X2. C.Tipo P
40.000
≤ 1/3 X1

Restricciones
X1≤ 4.500
X21/3 X1
X1 + X2 ≤ 5.000
X1 ; X2 ≥ 0

Función Objetivo
            Z = X1(30.000) + X2(40.000)