EJEMPLO PROGRAMACION NO LINEAL SIN RESTRICIONES

CON VARIAS VARIABLE METODO ASCENSO ACELERADO

1. Buscar el punto Óptimo de

𝑴𝒊𝒏 𝒁 = (𝒙𝟏 − √𝟓 )𝟐 + ( 𝒙𝟐 − 𝝅 )𝟐 + 𝟏𝟎 𝝐 = 𝟎. 𝟓 𝑿𝟎 = [ 𝟗. 𝟐𝟓𝟎 ; 𝟕. 𝟓𝟒𝟔 ]

La función original está en 𝑴𝒊𝒏𝒊𝒎𝒊𝒛𝒂𝒄𝒊𝒐𝒏 el Método aplicado es Ascenso Acelerado, por lo tanto la función debe pasars
𝑴𝒂𝒙𝒊𝒎𝒊𝒛𝒂𝒄𝒊𝒐𝒏 multiplicando por -1

𝑀𝑎𝑥 𝑍 = 2
−1( (𝑥1 − √5 ) + ( 𝑥2 − 𝜋 )2 + 10 )

𝟐
𝑴𝒂𝒙 𝒁 = − (𝒙𝟏 − √𝟓 ) − ( 𝒙𝟐 − 𝝅 )𝟐 − 𝟏𝟎

Primer Paso:

Calculamos el punto 𝑿𝟎 en la función original 𝑴𝒂𝒙

𝑋0 = [9.250 ; 7.546] 𝑋𝑇 9.250|

0 =|
7.546

𝒙𝟏 𝒙𝟐

𝒁 = − (𝒙 − √ 𝟐 ( 𝒙 − 𝝅 )𝟐 − 𝟏𝟎
𝟏 𝟓) − 𝟐

𝑍𝑋0 2
= − (9.250 − √5 ) ( 7.546 − 𝜋 )2 − 10
−

𝒁𝑿𝟎 = −𝟕𝟖. 𝟓𝟗𝟒

 Segundo Paso:
Hallamos el Gradiente 𝛁𝒁𝒙 de la función original

𝜕𝑍𝜕𝑍
∇𝑍𝑥 = [𝜕𝑥 ; 𝜕𝑥 ]
1 2

𝜕𝑍 𝜕𝑍
= −2( 𝑥1 − √5 ) = −2( 𝑥2 − 𝜋 )
𝜕𝑥1 𝜕𝑥2

∇𝑍𝑥 = [ −2( 𝑥1 − √5 ) ; −2( 𝑥2 − 𝜋 ) ] −2( 𝑥1 − √5 )

∇𝑍𝑥𝑇 = | |
−2( 𝑥2 − 𝜋

Hallamos el Gradiente 𝛁𝒁𝒙 de la función en el punto 𝑿𝟎

−2( 9.250 − √5 )
∇𝑍𝑇 = | | ∇𝑍𝑇𝑋0 = |−14.027|
𝑋0 −8.808
−2( 7.546 − 𝜋

Tercer paso:

Primer Iteración:

Utilizamos la siguiente fórmula para las iteraciones

1𝑋𝑇 = [ 𝑋𝑇0 + 𝜆∗ ∇0𝑋

𝑍𝑇0 ]
 9.250 −14.027 9.250 −14.027 𝜆0∗
𝑋1𝑇 = [ | | + 𝜆∗ 0 | |] 𝑋1𝑇 = [ |] 7.546
|+ | |
7.546 −8.808 −8.808 𝜆∗ 0

9.250 − 14.027 𝜆0∗ 𝑋1 = [ 𝟗. 𝟐𝟓𝟎 − 𝟏𝟒. 𝟎𝟐𝟕 𝝀∗ ; 𝟕. 𝟓𝟒𝟔 − 𝟖. 𝟖𝟎𝟖 𝝀∗ ]

𝑋1𝑇 = [ ] 𝟎 𝟎
7.546 − 8.808 𝜆0∗

𝒙𝟏 𝒙𝟐

Remplazamos los valores de 𝑿𝟏 en la función original 𝒁

𝟐
𝑴𝒂𝒙 𝒁 = − (𝒙𝟏 − √𝟓 ) − ( 𝒙𝟐 − 𝝅 )𝟐 − 𝟏𝟎

𝒁𝑿 𝟐
𝟏
= − (9.250 − 14.027 𝜆∗ − 𝟓) ( 7.546 − 8.808 𝜆∗ − 𝝅 )𝟐 − 𝟏𝟎
√ 0 − 0

𝒁𝑿 = − (7.013 − 14.027 𝜆∗ )𝟐 − ( 4.404 − 8.808 𝜆∗ )𝟐 − 10

𝟏 0 0

𝒁𝑿 = −[(7.013)2 − 2 ∗ 7.013 ∗ 14.027𝜆∗ + (14.027 𝜆∗ )2]

𝟏 0 0
− [(4.404)2 − 2 ∗ 4.404 ∗ 8.808𝜆∗ + (8.808𝜆∗ )2] − 10
0 0

𝒁𝑿 = −[49.182 − 196.742𝜆∗ + 196.756𝜆∗ 2] − [19.395 − 77.580𝜆∗ + 77.580𝜆∗ 2] − 10

𝟏 0 0 0 0

𝒁𝑿 = −49.182 + 196.742𝜆 − 196.756𝜆

∗ ∗ 2
− 19.395 + 77.580𝜆 − 77.580𝜆 ∗ ∗ 2
− 10
𝟏 0 0 0 0

𝒁𝑿 = −78.577 + 274.322𝜆∗ − 274.336𝜆∗ 2

𝟏 0 0

Derivamos la función 𝒁𝑿con

𝟏
respecto a 𝝀∗𝟎

𝜕 𝑍𝑋1 274.322 − 548.672𝜆∗ = 0

= 274.322 − 548.672𝜆 ∗
0
𝜕𝜆∗0 0
 274.322 ∗
548.672 = 𝜆0 𝜆∗
0 = 0.499

Calculamos 𝑿𝑻𝟏remplazando 𝝀∗𝟎

∗
1 𝑋𝑇 = [ 9.250
7.546−−14.027 𝜆 0]
8.808 𝜆∗0

𝑇 9.250 − 14.027 ∗ 0.499

1𝑋 =[
7.546 − 8.808 ∗ 0.499 ]

𝑇 2.250 𝑋1 = [ 2.250 ;3.150 ]

1𝑋 =[ ]
3.150

𝒙𝟏 𝒙𝟐

Calculamos la función 𝒁 con el punto𝑿𝟏

𝒁 = − (𝒙 − √ 𝟐 ( 𝒙 − 𝝅 )𝟐 − 𝟏𝟎
𝟏 𝟓) − 𝟐

𝑍𝑋1 2
= − (2.250 − √5 ) ( 3.150 − 𝜋 )2 − 10
−

𝒁𝑿𝟏 = −𝟏𝟎

𝑿𝟏 se le resta al hallado en 𝒁𝑿𝟎 si el resultado es menor que la tolerancia de error 𝝐 = 𝟎. 𝟓el ejercicio termina y el resulta
 𝒁𝑿𝟏 = −𝟏𝟎 𝒁𝑿𝟎 = −𝟕𝟖. 𝟓𝟗𝟒

𝒁𝑿𝟏 − 𝒁𝑿𝟎 = −𝟏𝟎 − (−𝟕𝟖. 𝟓𝟗𝟒) = 𝟔𝟖. 𝟓𝟗𝟒

Resultado mayor a la tolerancia continuamos……..

Segunda Iteración:

Hallamos el Gradiente 𝛁𝒁𝒙 de la función en el punto 𝑿𝟏

−2( 𝑥1 − √5 ) 2.250
∇𝑍𝑥𝑇 = | −2( 𝑥 − 𝜋 | 𝑋1𝑇 = [ 3.150 ]
2

−2( 2.250 − √5 ) −0.0278

𝑋1𝑇 = |
∇𝑍 ∇𝑍𝑋𝑇1 = | −0.0168|
−2( 3.150 − 𝜋 |

Utilizamos la siguiente fórmula para las iteraciones

2𝑋𝑇 = [ 𝑋𝑇1 + 𝜆∗ ∇1𝑋

𝑍𝑇1 ]

2.250 −0.0278 2.250 −0.0278 𝜆1∗

𝑇
𝑋2 = [ | +| 𝜆∗ 1 | |] 𝑋2𝑇 = [ | |+ | |
3.150 −0.0168 ] 3.150 −0.0168 𝜆∗ 1
 ∗ 𝑋2 = [ 𝟐. 𝟐𝟓𝟎 − 𝟎. 𝟎𝟐𝟕𝟖 𝝀∗ ; 𝟑. 𝟏𝟓𝟎 − 𝟎. 𝟎𝟏𝟔𝟖 𝝀∗ ]
1] 𝟏 𝟏
2 𝑋𝑇 = [ 2.250
3.150 −− 0.0278
0.0168 𝜆𝜆∗ 1

𝒙𝟏 𝒙𝟐

Remplazamos los valores de 𝑿𝟐 en la función original 𝒁

𝟐
𝑴𝒂𝒙 𝒁 = − (𝒙𝟏 − √𝟓 ) − ( 𝒙𝟐 − 𝝅 )𝟐 − 𝟏𝟎

𝒁𝑿 𝟐
𝟐
= − (2.250 − 0.0278 𝜆∗ − 𝟓) ( 3.150 − 0.0168 𝜆∗ − 𝝅 )𝟐 − 𝟏𝟎
√ 1 − 1

𝒁𝑿 = − (0.0139 − 0.0278 𝜆∗ )𝟐 − ( 8.407X10−3 − 0.0168 𝜆∗ )𝟐 − 10

𝟐 1 1

𝒁𝑿 = −[(0.0139) − 2 ∗ 0.0139 ∗ 0.0278𝜆 + (0.0278 𝜆 )2]

2 ∗ ∗

𝟐 1 1
− [(8.407X10−3)2 − 2 ∗ 8.407X10−3 ∗ 0.0168 𝜆∗ + (0.0168 𝜆∗ )2] − 10
1 1

2
𝒁𝑿 = − [1.932X10−4 − 7.728𝑋10−4𝜆1∗ + 7.728X10−4𝜆1∗ ]
𝟐

− [7.067X10−5
− 2.824X10−4 𝜆∗ + 2.822X10−4 𝜆∗ ] − 10
1 21

2
𝒁 = −1.932X10−4 + 7.728𝑋10−4𝜆∗ − 7.728X10−4𝜆∗ − 7.067X10−5 + 2.824X10−4𝜆∗
𝑿𝟐 1 1 1
2
− 2.822X10 𝜆 1 − 10
−4 ∗

2
𝒁 = −10 + 1.055𝑋10−3𝜆∗ − 1.055X10−3𝜆∗
𝑿𝟐 1 1

Derivamos la función 𝒁𝑿con

𝟐
respecto a 𝝀∗𝟏

𝜕𝑍𝑋2 1.055𝑋10−3 − 2.11𝑋10−3𝜆∗ = 0

= 1.055𝑋10−3 − 2.11𝑋10−3𝜆∗1 1
𝜕𝜆∗1
 1.055𝑋10−3 𝜆∗
1 = 0.5
2.11𝑋10−3 = 𝜆1
∗

Calculamos 𝑿𝑻𝟐remplazando 𝝀∗𝟏

∗
1]
2 𝑋𝑇 = [ 2.250
3.150 −− 0.0278
0.0168 𝜆𝜆∗1

𝑇 2.250 − 0.0278 ∗ 0.5

2𝑋 =[
3.150 − 0.0168 ∗ 0.5 ]

𝑇 2.236 𝑋2 = [ 2.236 ;3.141 ]

2𝑋 =[ ]
3.141

𝒙𝟏 𝒙𝟐

Calculamos la función 𝒁 con el punto𝑿𝟐

𝒁 = − (𝒙 − √ 𝟐 ( 𝒙 − 𝝅 )𝟐 − 𝟏𝟎
𝟏 𝟓) − 𝟐

𝑍𝑋2 2
= − (2.236 − √5 ) ( 3.141 − 𝜋 )2 − 10
−

𝒁𝑿𝟐 = −𝟏𝟎

𝑿𝟐 se le resta al hallado en 𝒁𝑿𝟏 si el resultado es menor que la tolerancia de error 𝝐 = 𝟎. 𝟓el ejercicio termina y el resulta
 𝒁𝑿𝟐 = −𝟏𝟎 𝒁𝑿𝟏 = −𝟏𝟎

𝒁𝑿𝟐 − 𝒁𝑿𝟏 = −𝟏𝟎 − (−𝟏𝟎) = 𝟎

Resultado:

𝑿𝟐 = [ 𝟐. 𝟐𝟑𝟔 ;𝟑. 𝟏𝟒𝟏 ]

𝒁𝑿𝟐 = −𝟏𝟎