Articulo
Articulo
Resumen Abstract
En este documento se encuentran detalladas las Redactar aquí el resumen en inglés con las mismas
instrucciones para que los autores preparen un especificaciones del formato descrito en español.
artículo en el formato requerido de la Revista ……………..
INGENIUS y además puede emplearse como plantilla.
Se presenta el formato de publicación, tamaños y
tipos de fuente, contiene además las normas para
presentar ecuaciones, figuras, tablas y referencias.
Los autores deben seguir las instrucciones para
mantener el estándar de publicación.
Esta primera sección es para generar un resumen del
contenido del artículo dando una clara indicación de
la justificación del tema, objetivos, metodología, y
principales resultados para que los lectores puedan
determinar si el texto completo será de su particular
interés. Debe contener un máximo de 230 palabras,
no debe incluir ecuaciones o referencias y debe estar
escrito de manera impersonal. Use la fuente Times
New Roman en tamaño 12. El contenido del
resumen debe estar completamente justificado.
1
En esta sección escriba la adscripción institucional actual en el siguiente orden: Dependencia a la que pertenece en su institución,
Institución a la que pertenece, país, ORCID; como se presenta en el ejemplo siguiente.
2
Laboratorio de Investigación en Sistemas Informáticos e Inteligencia Artificial , Universidad Politécnica de Valencia – España, ORCID:
0000-0002-8235-4886
3
……….
Autor para correspondencia: correoelectrónico@autorparacorrespondencia
1
Artículo Científico / Scientific Paper
Para la simulación del método de Monte Carlo se otras funciones de densidad de probabilidad son:
necesita aplicar una distribución de probabilidad, en este permitir estimaciones satisfactorias de la asimetría de la
caso utilizamos la probabilidad Normal paralelamente la distribución de densidad de probabilidad [14]. El
distribución de Weibull. nombre del método se lo lleva por su creador Weibull.
Para ajustar la distribución de Weibull se emplea el
método de máxima verosimilitud. Cabe recalcar que
2.2.1. Distribución de probabilidad Normal para encontrar el valor de k se debe utilizar un algoritmo
La distribución de probabilidad Normal es la más de iteración numérica tal como el método de Newton-
importante en toda la probabilidad y estadística, Raphson para obtener una solución aproximada. Una
conocida como gaussiana o acampanada es simétrica, la vez hallado el parámetro de forma (k) se puede hallar el
media, moda y mediana coinciden, y es descrita parámetro de escala (c). La fórmula que se utiliza para
completamente por sus dos parámetros mu ( μ) y sigma ( la distribución de Weibull es la siguiente: [15].
( ( ))
k
x
σ ). El primero en detallar dicho paradigma fue Abraham
()
k−1 −
k x c
(2)
de Moivre, no obstante, fue Gauss quien colocó al f ( x )= e
c c
alcance de las personas [12]. La fórmula general de la
distribución de probabilidad normal es la siguiente:
10-4
2
6
− ( X −μ )
1 2σ
2
(1)
f ( X ; μ , σ )= ∙e
√2 π σ 5
4
10-4
6
3
5
2
4
1
3
0
500 1000 1500 2000 2500 3000 3500 4000 4500
2
Figura 2. Gráfica de la densidad de una Distribución de Weibull. [4]
Power [kW]
muestran a continuación:
2500
2000
La Demanda Eléctrica se relaciona con la cantidad de
energía consumida cada hora durante el mes de marzo,
como resultado una matriz de datos que contiene 24 1500
registros por cada uno de los 31 días del mes, los datos
1
2
3
4
5
6
7
8
9
10
12
13
14
15
16
17
18
19
20
21
22
23
24
11
Time [H]
utilizados hacen énfasis a un promedio (media) mensual
de cada hora respecto al mes de marzo como se Figura 4. Promedio mensual de la Demanda Eléctrica.
representa en la Figura 4.
La Radiación Solar y la Velocidad del Viento
reflejan mediciones tomadas cada 5 minutos durante el
mes de agosto, dándonos una matriz de 288 datos por 30
días, al igual que la Demanda Eléctrica se usó los datos
promedio de cada 5 minutos como se representa en la
Figura 5 y la Figura 6.
Al analizar la Demanda Eléctrica, se observó que
los datos históricos presentan una variabilidad moderada
con respecto al consumo eléctrico. Esto sugiere que a lo
largo del mes no se registraron cambios significativos en 600 Average solar radiation
el consumo eléctrico. Por otro lado, en el caso de la
Radiación Solar, se identificó una menor variabilidad en 500
la intensidad de la radiación, ya que tiende a mantener
Solar Radiation [W/m 2]
100
0
10 0
11 0
12 0
12 0
13 0
13 0
14 0
40
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10 0
60
2
8
4
0
6
2
8
0
12
18
24
30
36
42
48
54
60
66
72
78
84
90
96
Time [min]
4
Apellido Autor et al / Titulo del Articulo
4
Distribución de Probabilidad Normal, que calcula tanto
3.5 el valor promedio (media) como la desviación estándar
3 y otro basado en la Distribución de Weibull, que calcula
los parámetros de forma y escala. Estos valores son
2.5 esenciales para comprender el comportamiento de los
2 datos y permiten generar valores aleatorios.
Posteriormente, se aplica una restricción de error, que
1.5 compara los valores generados con los valores reales
0
60
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10 0
10 0
11 0
12 0
12 0
13 0
13 0
14 0
40
para llegar al resultado deseado, como se detalla más
12
18
24
30
36
42
48
54
60
66
72
78
84
90
96
2
8
4
0
6
2
8
Time [min]
adelante.
Figura 6. Promedio mensual de la Velocidad del Viento. En la Figura 9, Figura 10 y Figura 11 se muestra la
generación de series sintéticas utilizando la Distribución
de Probabilidad Normal. En la Figura 12, Figura 13 y
B. Simulación del Método de Monte Carlo Figura 14 se presenta la generación de series sintéticas
El método de Montecarlo puede usarse para abordar utilizando la Distribución de Weibull. Cada imagen
problemas complejos a través de un proceso iterativo. muestra una representación visual de los datos, que
Este proceso utiliza datos históricos para generar un incluye el vector original y 10 series sintéticas
conjunto de valores aleatorios que se combinan para generadas con el método mencionado anteriormente.
obtener una estimación del resultado. Se puede
visualizar en el diagrama de flujo de la Figura 7 y
Figura 8.
5
Figura 8. Diagrama de flujo del método de Montecarlo aplicada la
INICIO
distribución de Weibull.
4500
Ingresar datos:
Typical daily demand
Vector datos, n Synthetic Series
4000
Calcular:
m = length (Vector datos) 3500
mu = mean (Vector datos)
Power [kW]
sigma = std (Vector datos)
3000
ǡ ൌ
ൌ ͳ
2500
i, k n, m 2000
1500
Calcular:
1
2
3
4
5
6
7
8
9
10
12
13
14
15
16
17
18
19
20
21
22
23
24
11
Si Time [H]
sinteticas (k, i) = normrnd (mu, sigma, 1, 1)
No No
Figura 9. Series sintéticas generadas de Demanda Eléctrica con el
método Montecarlo (Distribución de probabilidad normal).
Restricción:
ER ͲെͳͲͲΨ
FIN
Figura 7. Diagrama de flujo del método de Montecarlo aplicada la
Distribución de probabilidad normal.
400
Calcular:
m = length (Vector datos)
mu = mean (Vector datos) 300
weibull = wblfit (Vector datos)
200
i, k == 1
100
i, k n, m
0
10 0
11 0
12 0
12 0
13 0
13 0
14 0
40
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10 0
60
2
8
4
0
6
2
8
Calcular:
0
12
18
24
30
36
42
48
54
60
66
72
78
84
90
96
Si
sinteticas (k, i) = wblrnd (weibull (k, 1), weibull (k, 2), 1, 1) Time [min]
No No
FIN
6
Apellido Autor et al / Titulo del Articulo
4 500
3
300
2.5
200
2
100
1.5
0
60
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10 0
10 0
11 0
12 0
12 0
13 0
13 0
14 0
40
12
18
24
30
36
42
48
54
60
66
72
78
84
90
96
2
8
4
0
6
2
8
0
Time [min]
10 0
11 0
12 0
12 0
13 0
13 0
14 0
40
Figura 11. Series sintéticas generadas de Velocidad del Viento con el
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10 0
60
2
8
4
0
6
2
8
0
12
18
24
30
36
42
48
54
60
66
72
78
84
90
96
método Montecarlo (Distribución de probabilidad normal). Time [min]
5.5
Typical daily wind speed
5 Synthetic Series
4.5
Wind Speed [m/s]
3.5
2.5
4500 1.5
Typical daily demand
0
60
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10 0
10 0
11 0
12 0
12 0
13 0
13 0
14 0
40
12
18
24
30
36
42
48
54
60
66
72
78
84
90
96
2
8
4
0
6
2
8
Synthetic Series
Time [min]
4000
3500 Figura 14. Series sintéticas generadas de Velocidad del Viento con el
método Montecarlo (Distribución de Weibull).
Potencia [kW]
3000
C. Validación
2500 El análisis de error significativo, requiere que la medida
de error propuesta se relacione con los datos reales
2000 obtenidos en los experimentos que se van a comparar,
sosteniendo dos requisitos necesarios para una buena
1500 medida de error: (a) una medida de error es una
cuantificación de las diferencias entre el observable
10
12
13
14
15
16
17
18
19
20
21
22
23
24
11
1
2
3
4
5
6
7
8
9
Time [H]
objetivo y el observable aproximado que se está
7
midiendo; en particular, debería indicar correctamente En donde:
los casos en los que el objetivo y los observables
V E=valor exacto
aproximados concuerdan y en los que no; (b) el error
puede estimarse a partir de los datos obtenidos en el V A =valor aproximado
experimento en cuestión y una medición de referencia
ideal del objetivo observable [23]. E A =error absoluto
Es fundamental tener en cuenta la importancia del
análisis de error para evaluar el proceso mencionado
anteriormente. Este análisis nos permite controlar que 2.2. Error relativo
los datos generados a partir de una distribución se El error relativo es una medida de la precisión relativa
ajusten a los datos reales, lo que implica una menor de una medición o cálculo en relación con el valor real.
dispersión en las gráficas sintéticas. A continuación, los Se calcula como el cociente del error absoluto y el valor
tipos de errores mencionados se aplicaron real (o una estimación de la magnitud real). El error
posteriormente tanto a los datos reales como a cada serie relativo se expresa típicamente como un porcentaje o
sintética creada. una fracción [26]. Fórmula del error relativo:
Para garantizar la generación de números aleatorios
controlados se estableció errores relativos de forma EA
Erelativo = ( 4)
independiente para cada contexto en función de las VE
características específicas de cada conjunto de datos
ingresados. Cada error relativo se determinó de manera En donde:
diferente de acuerdo al comportamiento de los datos, V E=valor exacto
con el objetivo de lograr un margen de Error Cuadrático
Medio Normalizado (ECMN) como se observa en el E A =error absoluto
diagrama de flujo de la Figura 15 que fuera menor al
5%. Esto asegura que el (ECMN) se mantenga dentro de 2.3. Error cuadrático medio (ECM)
un rango óptimo y demuestre la efectividad del método
El error cuadrático medio (ECM) es una medida
de Montecarlo como una estrategia confiable para
comúnmente utilizada para evaluar la precisión de un
generar series sintéticas precisas. En este proceso, se
modelo de predicción o estimación en comparación con
evaluó el error relativo entre el valor generado y el valor
los valores reales. Representa el promedio de los errores
real. Si el error relativo es menor que la restricción
al cuadrado entre las predicciones del modelo y los
establecida, se guarda el valor generado. En caso
valores observados [11]. La fórmula del error cuadrático
contrario, se procede a generar un nuevo valor hasta que
medio es:
se cumpla la condición preestablecida.
n
1
MSE=ECM = ∑ ( Y i−Y i ) (5)
2
8
Apellido Autor et al / Titulo del Articulo
young and elderly men” Gait Pos vol. 20, pp. 291- Tesis:
298. 2004 [11] J. Basantes, F. Torres, “Desarrollo de un
[3] J. Riess, J. J. Abbas, “Adaptive control of cyclic Sistema de Control para un Brazo Robótico
movements as muscles fatigue using functional mediante Adquisición y Procesamiento de
neuromuscular stimulation”. IEEE Trans. Neural Imágenes” Proyecto de titulación, Escuela
Syst. Rehabil. Eng vol. 9, pp.326–330, 2001. Politécnica Nacional, Quito, Ecuador. Feb. 2009.
[12] S.M. Newman, “Active damping control of a
Libros: flexible space structure using piezoelectric sensors
[4] L. Stein, “Random patterns,” in Computers and and actuators” Master Thesis, U.S. Naval
You, J. S. Brake, Ed. New York: Wiley, 2004, pp. Postgraduate School, 1992.
55-70. [13] L.M. Moreno, "Computación paralela y
[5] A. Barrientos, L. Peñin, et. Al, “Coordenadas” en entornos heterogéneos," Tesis doctoral, Dep.
Fundamentos de Robótica, 2da ed., vol. 2, Ed. Estadística, Investigación Operativa y
McGraw-Hill, España, 2007, pp. 217–29. Computación, Universidad de La Laguna, La
[6] G. O. Young, “Synthetic structure of industrial Laguna, 2005.
plastics” in Plastics, 2nd ed., vol. 3, J. Peters, Ed.
New York: McGraw-Hill, 1964, pp. 15-64. Estándares:
[14] IEEE Guide for Application of Power
Reportes Técnicos: Apparatus Bushings, IEEE Standard C57.19.100-
[7] E. E. Reber, R. L. Mitchell, y C. J. Carter, 1995, Aug. 1995.
"Oxygen absorption in the Earth's atmosphere,"
Aerospace Corp., Los Angeles, CA, Tech. Rep. Patentes:
TR-0200 (4230-46)-3, Nov. 1968. [15] G. Brandli, M. Dick, "Alternating current fed
[8] M. A. Brusberg and E. N. Clark, “Installation, power supply," U.S. Patente 4 084 217, Nov. 4,
operation, and data evaluation of an oblique- 1978.
incidence ionosphere sounder system,” in “Radio [16] J. P. Wilkinson, “Nonlinear resonant circuit
Propagation Characteristics of the Washington- devices,” U.S. Patent 3 624 125, July 16, 1990.
Honolulu Path,” Stanford Res. Inst., Stanford, CA,
Contract NOBSR-87615, Final Rep., Feb. 1995, Manuales:
vol. 1 [17] Motorola Semiconductor Data Manual,
Motorola Semiconductor Products Inc., Phoenix,
Artículos presentados en conferencias (No AZ, 1989.
publicados):
[9] Vázquez, Rolando, Presentación curso “Realidad Recursos de Internet:
Virtual”. National Instruments. Colombia, 2009. [18] J. Lau, “Directrices internacionales para la
alfabetización informativa” [online]. México:
Artículos de Memorias de Conferencias Universidad Veracruzana, 2004 Disponible
(Publicados): en:[Link]
[10] L. I. Ruiz, A. García, J. García, G. Taboada. [Link]
“Criterios para la optimización de sistemas [19] E. H. Miller, “A note on reflector arrays”
eléctricos en refinerías de la industria petrolera: [online] IEEE Trans. Antennas Propag., vol. 53,
influencia y análisis en el equipo eléctrico,” IEEE pp. 475, 2005. Disponible en:
CONCAPAN XXVIII, Guatemala 2008. [Link]
arnumber=1549967&tag=1
11