Trabajo Colaborativo Probabilidad

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 18

FACULTAD DE INGENIERÍA Y

CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

PROBABILIDAD
DISTRIBUCIÓN DE PROBABILIDAD DE WEIBULL

JORGE EDUARDO GIRALDO MEDINA COD. 1721020459


SANDRA ANYUL YARURO MARIN COD. 1911022569
LUIS OSWALDO QUINTERO VELASQUEZ COD.1811024524
MANUEL DAVID GUERRERO RAMIREZ COD. 1911022142

POLITECNICO GRANCOLOMBIANO
FACULTAD DE INGENIERIA DE SOFTWARE
2020
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

TABLA DE CONTENIDOS

INTRODUCCIÓN ..................................................................................................................................3
OBJETIVOS ..........................................................................................................................................4
SOLUCIÓN DE LOS EJERCICIOS............................................................................................................5
CONCLUSIONES ................................................................................................................................17
BIBLIOGRAFÍA ...................................................................................................................................18
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

INTRODUCCIÓN

En el desarrollo del presente trabajo se explica paso a paso las soluciones para
cada uno de los ejercicios planteados, para lo cual, a través de la distribución de
probabilidad Weibull se pusieron en práctica los conceptos vistos a lo largo de la
materia, de esta manera se evidencia la aplicación de las respectivas formulas y se
presentan las gráficas que ayudan a entender más fácilmente los datos entregados.
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

OBJETIVOS

• Realizar cada uno de los ejercicios indicando los pasos efectuados en sus respectivas
soluciones.
• Entender claramente los datos proporcionados de cada departamento para realizar una
adecuada solución de los planteamientos
• Identificar las respectivas variables indicadas en cada ecuación con el fin de reemplazar
los valores adecuadamente.
• Generar un debate constructivo con los integrantes del equipo con el fin de plantear
diferentes soluciones y corregir los posibles errores.
• Poner en práctica los conceptos de probabilidad vistos en clase.
.
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

SOLUCIÓN DE LOS EJERCICIOS


PUNTO 1
El grupo debe seleccionar 2 departamentos y realizar un histograma de la variable velocidad
del viento y temperatura. ¿Cuál presenta mayor variabilidad? justificar.

Los departamentos seleccionados son Boyacá y Amazonas, se estiman resultados


favorables en Boyacá por su altitud, y en Amazonas por las lluvias, por lo tanto, para ambos
departamentos se extrajeron los datos del Excel del IDEAM [1].
Se analizaron las medidas de tendencia central, media, moda, mediana, mínimos y
máximos, la varianza, la desviación estándar, los rangos, los intervalos y amplitud de dichos
intervalos para posterior mente agrupar los datos.

Amazonas:
Medidas de tendencia central
Para realizar el cálculo de las medidas de tendencia central se filtró la búsqueda de los
datos en el archivo del IDEAM, luego se copiaron en Excel, y se reescribieron con la opción
Ctrl+B, la mayoría de las variables, tienen una función en Excel, por lo cual fue bastante útil
esta herramienta a continuación unas de ellas.
moda (número1: número2)
MEDIANA (NÚMERO1: NÚMERO2)
DESVEST . M (NÚMERO1: NÚMERO2)
NÚMERO DE INTERVALOS 1+3,322*LOG (NÚMERO DE DATOS)

T ABLA 1

Variable Temperatura °C Velocidad del viento m/s


Media 13,22327408 1,795245969
Moda 10,1 1,9
Mediana 12,5 1,8
Mínimo 2,4 0,1
Máximo 36,8 6,3
Varianza 16,62401695 0,485930681
Desviación estándar 4,077256057 0,697087284
Rango 34,4 6,2
Intervalos (Regla de Sturges) 15 15
Amplitud 2,3 0,42
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

Agrupación de datos para Temperatura Amazonas

T ABLA 2

AMAZONAS Temperatura
Marca de
numero Intervalo Frecuencia Frecuencia relativa
Clase
1 19,5 20,9 20,2 45 3,466872111
2 21 22,4 21,7 306 23,57473035
3 22,5 23,9 23,2 234 18,02773498
4 24 25,4 24,7 179 13,79044684
5 25,5 26,9 26,2 144 11,09399076
6 27 28,4 27,7 144 11,09399076
7 28,5 29,9 29,2 103 7,935285054
8 30 31,4 30,7 65 5,00770416
9 31,5 32,9 32,2 38 2,927580894
10 33 34,4 33,7 27 2,080123267
11 34,5 35,9 35,2 13 1,001540832

Histograma Temperatura Amazonas


25 23 4/7

20 18

15 13 4/5
11 11
10 8
5
5 3 1/2 3
2
1
0
20,2 21,7 23,2 24,7 26,2 27,7 29,2 30,7 32,2 33,7 35,2

FIGURA 1 H ISTOGRAMA TEMPERATURA AMAZONAS


FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

Agrupación de datos para Velocidad del viento en amazonas


T ABLA 3

AMAZONAS Velocidad del viento


Marca de
numero Intervalo Frecuencia Frecuencia relativa
Clase
1 0,2 0,35 0,275 14 1,078582435
2 0,36 0,51 0,435 112 8,628659476
3 0,52 0,67 0,595 125 9,630200308
4 0,68 0,83 0,755 297 22,88135593
5 0,84 0,99 0,915 170 13,09707242
6 1 1,15 1,075 324 24,9614792
7 1,16 1,31 1,235 175 13,48228043
8 1,32 1,47 1,395 28 2,157164869
9 1,48 1,63 1,555 37 2,850539291
10 1,64 1,79 1,715 10 0,770416025
11 1,8 1,95 1,875 6 0,462249615

Histograma Velocidad del viento Amazonas


30
25
25 22 7/8

20

15 13 13 1/2

9 5/8
10 8 5/8

5 2 6/7
2 1/6
1 7/9 1/2
0
0,275 0,435 0,595 0,755 0,915 1,075 1,235 1,395 1,555 1,715 1,875

FIGURA 2 H ISTOGRAMA V ELOCIDAD DEL VIENTO EN AMAZONAS


FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

Boyacá:
Medidas de tendencia central
Para realizar el cálculo de las medidas de tendencia central se filtró la búsqueda de los
datos en el archivo del IDEAM, luego se copiaron en Excel, y se reescribieron con la
opción Ctrl+B, la mayoría de las variables, tienen una función en Excel, por lo cual fue
bastante útil esta herramienta a continuación unas de ellas.
moda (número1: número2)
mediana (número1: número2)
desvest.m (número1: número2)
número de intervalos 1+3,322*LOG (número de datos)
T ABLA 4

Variable Temperatura °C Velocidad del viento m/s


Media 13,22327408 1,795245969
Moda 10,1 1,9
Mediana 12,5 1,8
Mínimo 2,4 0,1
Máximo 36,8 6,3
Varianza 16,62401695 0,485930681
Desviación estándar 4,077256057 0,697087284
Rango 34,4 6,2
Intervalos (Regla de Sturges)
15 15
Amplitud 2,3 0,42
Agrupación de datos para Temperatura en Boyacá
T ABLA 5

BOYACA Temperatura
numero Intervalo Marca de Clase Frecuencia Frecuencia relativa
1 0,1 2,4 1,25 41 0,282485876
2 2,4 4,7 3,55 165 1,136833402
3 4,7 7 5,85 1761 12,13311286
4 7 9,3 8,15 3864 26,62257131
5 9,3 11,6 10,45 3161 21,77897203
6 11,6 13,9 12,75 2560 17,63814248
7 13,9 16,2 15,05 1471 10,13504203
8 16,2 18,5 17,35 750 5,167424556
9 18,5 20,8 19,65 393 2,707730467
10 20,8 23,1 21,95 198 1,364200083
11 23,1 25,4 24,25 88 0,606311148
12 25,4 27,7 26,55 28 0,192917183
13 27,7 30 28,85 19 0,130908089
14 30 32,3 31,15 9 0,062009095
15 32,3 34,6 33,45 6 0,041339396
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

Histograma Temperatura Boyacá


30
26 5/8

25
21 7/9

20 17 2/3

15
12 1/8
10 1/7
10
5 1/6
5 2 5/7
1 1/7 1 3/8 3/5
2/7 1/5 1/8 0 0
0
1,25 3,55 5,85 8,15 10,45 12,75 15,05 17,35 19,65 21,95 24,25 26,55 28,85 31,15 33,45

FIGURA 3 H ISTOGRAMA TEMPERATURA EN B OYACÁ


Agrupación de datos Velocidad del viento para Boyacá
T ABLA 6

BOYACA Velocidad del viento


numero Intervalo Marca de Clase Frecuencia Frecuencia relativa
1 2,4 4,7 3,55 220 1,51577787
2 4,7 7 5,85 957 6,593633733
3 7 9,3 8,15 1952 13,44908364
4 9,3 11,6 10,45 3938 27,13242387
5 11,6 13,9 12,75 4647 32,01736255
6 13,9 16,2 15,05 1771 12,20201185
7 16,2 18,5 17,35 450 3,100454733
8 18,5 20,8 19,65 151 1,040374811
9 20,8 23,1 21,95 115 0,792338843
10 23,1 25,4 24,25 98 0,675210142
11 25,4 27,7 26,55 81 0,558081852
12 27,7 30 28,85 61 0,420283864
13 30 32,3 31,15 44 0,303155574
14 32,3 34,6 33,45 22 0,151577787
15 34,6 36,9 35,75 7 0,048229296
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

Histograma Velocidad del viento Boyacá


35
32

30
27 1/8

25

20

15 13 4/9
12 1/5

10
6 3/5

5 3 1/9
1 1/2 1 4/5 2/3 5/9 3/7 1/3 1/7 0
0
3,55 5,85 8,15 10,45 12,75 15,05 17,35 19,65 21,95 24,25 26,55 28,85 31,15 33,45 35,75

FIGURA 4 H ISTOGRAMA VELOCIDAD DEL VIENTO

OBSERVACIÓN:
Gracias a los histogramas se puede ver una variabilidad en ambos departamentos, si
comparamos la temperatura con el viento en cada uno de los departamentos, quien
presenta mayor variabilidad es Amazonas, si comparamos ambos departamentos en
temperatura son muy parecidos los histogramas, tienen una campana con bastante
distribución, y si comparamos los histogramas de velocidad del viento el que presenta
mayor variabilidad es Boyacá.
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

PUNTO 2
Calcular para cada ciudad los parámetros k y c de la distribución de Weibull, para ello se
usaron las ecuaciones 3 y 4 y se sustituyeron en la función de probabilidad f(v) (ecuación
1).
𝜎 −1.09
𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 3. 𝐾=[ ]
𝑣

Por lo tanto, sigma es la desviación estándar y v es el promedio de la velocidad del viento,


los cuales fueron calculados anteriormente.
𝑣
𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 4. 𝑐= 1
𝑟(1 + )
𝑘
Donde r corresponde a la función Gamma para el valor dentro del paréntesis.
Se calcula el valor de k y c para ambos departamentos:

AMAZONAS:
Valores de velocidad promedio y de sigma

𝜎 = 0,28713216
𝑣 = 0,915177196
Valor de k

0,28713216 −1.09
𝐾=[ ]
0,915177196

𝑘 = 3,537785142

Valor de c

0,915177196
𝑐= 1
𝑟 (1 + )
3,537785142

0,915177196
𝑐=
𝑟(1,282662728)

0,915177196
𝑐=
0,90026152

𝑐 = 1,016568163
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

BOYACÁ:
Valores de Velocidad media y de sigma

𝜎 = 0,697087284
𝑣 = 1,795245969
Valor de k

0,697087284 −1.09
𝐾=[ ]
1,795245969

𝑘 = 2,804220189
Valor de c

1,795245969
𝑐= 1
𝑟 (1 + )
2,804220189

1,795245969
𝑐=
𝑟(1,356605378)

1,795245969
𝑐=
0,890501949

𝑐 = 2,015993307

El valor de GAMMA se halla con Excel teniendo en cuenta el valor de (1+1/K) para cada
departamento con la siguiente fórmula

T ABLA 7

GAMA PARA BOYACÁ 0,890501949


FORMULA GAMMA(1,356605378)
GAMMA PARA AMAZONAS 0,900261517
FORMULA GAMMA(1,282662728)
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

Las funciones quedarían expresadas de la siguiente forma

𝑓 (𝑣 )𝐴𝑀𝐴𝑍𝑂𝑁𝐴𝑆
3,537785142−1 3,537785142
3,537785142 𝑣 [−[
𝑣
] ]
=[ ]∗[ ] 𝑒 1,016558163
1,016558163 1,016558163

𝑓 (𝑣 )𝐵𝑂𝑌𝐴𝐶Á
2,804220191−1 2,804220191
2,804220191 𝑣 [−[
𝑣
] ]
=[ ]∗[ ] 𝑒 2,015993301
2,015993301 2,015993301

PUNTO 3
Graficar la distribución de probabilidad f(v) obtenida en el ítem anterior y
compararla con el histograma obtenido en el punto 1) de la velocidad del viento
para los departamentos o ciudades seleccionados. ¿Qué puede concluir de las
gráficas en relación con el comportamiento de las variables?
Para Graficar la distribución de probabilidad se hace uso de la plataforma
DESMOS que es una herramienta graficadora de libre uso en internet,

AMAZONAS

FIGURA 5 GRAFICA DE LA FUNCIÓN

BOYACÁ

FIGURA 6 GRAFICA DE LA FUNCIÓN


FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

OBSERVACIÓN:
Son similares, muestran una campana donde la variabilidad de los datos de la
velocidad del viento es clara, más sin embargo a grandes rasgos se puede estimar
que AMAZONAS se descartaría para una posible intervención y se realizaría la
instalación de las turbinas eólicas en BOYACÁ.

PUNTO 4
Para cada departamento o ciudad seleccionado, obtener el valor de velocidad del
viento más probable y el valor de la velocidad del viento que entregaría la máxima
energía eólica, use las ecuaciones 5 y 6. Al comparar los valores para los dos
departamentos o ciudades seleccionados, ¿cuál de ellos tiene más probabilidad
de generar mayor energía eólica?

𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 5. 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑣𝑖𝑒𝑛𝑡𝑜 𝑚á𝑠 𝑝𝑟𝑜𝑏𝑎𝑏𝑙𝑒


1
𝑘−1 𝑘
𝑣𝑚𝑝 = 𝑐 [ ]
𝑘

𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 6. 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑣𝑖𝑒𝑛𝑡𝑜 𝑞𝑢𝑒 𝑒𝑛𝑡𝑟𝑒𝑔𝑎 𝑙𝑎 𝑚á𝑧𝑖𝑚𝑎 𝑒𝑛𝑒𝑟𝑔í𝑎 𝑒ó𝑙𝑖𝑐𝑎


1
𝑘+2 𝑘
𝑣𝑚á𝑥𝐸 = 𝑐 [ ]
𝑘

Para el departamento de Amazonas los datos son

𝑐 = 1,016568163
𝑘 = 3,537785142

Velocidad del viento más probable

1
3,537785142 − 1 3,537785142
𝑣𝑚𝑝 = 1,016568163 [ ]
3,537785142
𝑣𝑚𝑝 = 0,9254541047

Velocidad del viento que entrega la máxima energía eólica

1
3,537785142 + 2 3,537785142
𝑣𝑚á𝑥𝐸 = 1,016568163 [ ]
3,537785142
𝑣𝑚á𝑥𝐸 = 1,1538355735
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

Para el departamento de Boyacá los datos son

𝑐 = 2,015993307
𝑘 = 2,804220189

Velocidad del viento más probable

1
2,804220189 − 1 2,804220189
𝑣𝑚𝑝 = 2,015993307 [ ]
2,804220189
𝑣𝑚𝑝 = 1,72262643

Velocidad del viento que entrega la máxima energía eólica

1
2,804220189 + 2 2,804220189
𝑣𝑚á𝑥𝐸 = 2,015993307 [ ]
2,804220189
𝑣𝑚á𝑥𝐸 = 2,442683911

OBSERVACIÓN:
Lo que se observa es que Boyacá tiene gran probabilidad de generar mayor
energía eólica que Amazonas, con una diferencia pequeña pero significante, con
esto se corrobora lo dicho en el punto anterior.

PUNTO 5

¿Cuál es la probabilidad de que la velocidad del viento en uno de los


departamentos o ciudades seleccionadas sobrepase percentil 75? (haga uso de la
función f(v) obtenida en el punto 4 o de F(v) ecuación 2)

𝑘 (𝑛 + 1)
𝑃𝑜𝑠𝑖𝑐𝑖ó𝑛 =
100

75(14514 + 1)
𝑃75 = = 10886,25
100

La probabilidad de que la velocidad del viento supere el percentil 75 es de 8.164,5


FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

PUNTO 6
¿Cuál es la probabilidad que en uno de los departamentos o ciudades
seleccionadas se registre velocidades entre los percentiles 50 a 70? (haga uso de
la función f(v) obtenida en el punto 4)

Se escoge el departamento de Boyacá ya que los gráficos nos demuestran que es


más probable que pueda generar mayor energía eólica que Amazonas, los Percentil
que se toman son k=50 y 70 respectivamente, y n=número total de datos.
Para Boyacá se tienen N=14514 Datos

𝑘(𝑛 + 1)
𝑃𝑜𝑠𝑖𝑐𝑖ó𝑛 =
100
50(14514 + 1)
𝑃50 = = 7257.5
100

70(14514 + 1)
𝑃70 = = 10160.5
100

Teniendo en cuenta estas posiciones, reducimos nuestra tabla de valores para la


velocidad del viento a (7257,10161) y pasamos a un n=2895 de datos.
Boyacá
ṽ = 1,7201
𝜎 = 1,3117

1,3117 −1,09
𝑘=[ ]
1,7201
𝑘 = 1,3436

1,7201
𝑐= 1
𝑟 (1 + 1,3436)
1,7201
𝑐=
0,917771538

𝑐 = 1,8742

Con estos nuevos datos ya se puede calcular la probabilidad

1,3436 1,7201 (1,3436−1) [−(1,7201)]1,3436


𝑓(𝑣) = ( )∗( ) 𝑒 1,8742
1,8742 1,8742

𝑓(𝑣) = 2,1381
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

CONCLUSIONES

• Los modelos de probabilidad y estadística son de gran utilidad al momento de


tomar decisiones en todos los ámbitos de la vida laboral y personal ya que ofrecen
una visión más amplia de la información.
• El uso de la herramienta de Excel permitió un manejo de los datos eficiente, así
mismo como las fórmulas que este programa tiene preestablecidas, además la
visualización de gráficos permitió aclarar dudas y comparar respuestas.
• A través del desarrollo de este taller se afianzaron los conocimientos y conceptos
de probabilidad vistos en clase.
• Se pudo apreciar por medios de los modelos de distribución de datos que es más
fácil entender la información por medio de graficas que los representen
correctamente.
FACULTAD DE INGENIERÍA Y
CIENCIAS BÁSICAS

PROYECTO GRUPAL Institución Universitaria Politécnico Grancolombiano

BIBLIOGRAFÍA

[1]https://docs.google.com/spreadsheets/d/1PndlDQJ5esXXtaxlvJWLWzGPiztOh9NzgCk5hX9rKQk/
edit#gid=50990622

https://es.wikipedia.org/wiki/Distribuci%C3%B3n_de_Weibull

https://www.youtube.com/watch?v=Xsg7nR-GvyE

https://www.youtube.com/watch?v=CuKr7GzohbI&t=443s

https://www.youtube.com/watch?v=Z7Cd4OdFWUg&t=229s

https://www.youtube.com/watch?v=XznHWLqwor8&t=1s

https://www.youtube.com/watch?v=lbYZW92Hh34&t=2s

https://www.youtube.com/watch?v=9piGuCi1R-w

También podría gustarte