Introducción a Excel y Resolución de Problemas
Excel básico
Pedro Corcuera
Dpto. Matemática Aplicada y
Ciencias de la Computación
Universidad de Cantabria
corcuerp@[Link]
Índice
Excel 2
Objetivos
Excel 3
Técnicas generales de
resolución de problemas
Excel 4
Técnicas generales de
resolución de problemas
• El análisis de ingeniería es un proceso sistemático
para analizar y comprender los problemas que se
encuentran en los diversos campos de la ingeniería.
• Para llevar a cabo este proceso de manera
satisfactoria se debe estar familiarizado con técnicas
generales de resolución de problemas.
Excel 5
Pasos para la resolución de problemas1
Excel 7
Fundamentos de Ingeniería aplicables
Excel 10
Procedimientos de solución matemática
Excel 11
Rol de las hojas de cálculo
Excel 13
Manejo de Excel
Excel 14
Iniciar Excel
Excel 17
Ventana Excel
Excel 18
Componentes de la ventana1
• Barra de título
Excel 19
Componentes de la ventana2
• Cinta de Opciones
Excel 20
Componentes de la ventana3
Excel 21
Componentes de la ventana4
Excel 22
Fundamentos de Excel1
Excel 23
Fundamentos de Excel2
Excel 24
Fundamentos de Excel3
Excel 25
Movimientos por la hoja1
Excel 26
Movimientos por la hoja2
Excel 27
Movimiento en el libro
Excel 28
Introducción de datos
Excel 29
Datos numéricos
Excel 30
Asignación de nombres
Excel 31
Eliminar o borrar contenido de celdas
Excel 32
Inserción de figuras, texto, imágenes y
ecuaciones1
• Para introducir figuras, esquemas, texto artístico e imágenes
se selecciona:
– Formas: Insertar → Ilustraciones → Formas
– Esquemas: Insertar → Ilustraciones → SmartArt
– Imágenes prediseñadas: Insertar → Ilustraciones → Imágenes en
línea
– Imágenes desde archivo: Insertar → Ilustraciones → Imágenes
– Texto artístico: Insertar → Texto → WordArt
– Cuadros de Texto: Insertar → Texto → Cuadro de texto
– Ecuaciones: Insertar → Símbolos → Ecuación
Excel 33
Inserción de figuras, texto, imágenes y
ecuaciones2
Excel 34
Formato F1C1 para celdas
Excel 35
Fórmulas
Excel 37
Operadores1
Excel 38
Operadores2
Excel 39
Operadores3
Excel 40
Operadores4
Excel 41
Precedencia de operadores
Excel 42
Fórmulas con datos en más de una hoja
Excel 43
Fórmulas con datos en más de una hoja
Excel 44
Fórmulas con datos en más de un libro
Excel 45
Fórmulas con datos en más de un libro
Excel 46
Referencias circulares
Excel 47
Formato de celdas1
Excel 48
Formato de celdas2
Excel 49
Formato de celdas3
• Formato de números.
– General: El contenido se presenta como se ha introducido.
– Número: Adecuado para representar números. Se especifica el
número de decimales, separador de miles y números negativos
– Moneda: Se usa para cantidades monetarias. Se especifica el
número de decimales, la moneda y formato de negativos.
– Contabilidad: Igual que el formato moneda, la diferencia es que
alinea los números por la coma decimal y el símbolo de moneda.
Excel 50
Formato de celdas4
Excel 51
Formato de celdas5
• Fecha-Hora
– Fecha: número (parte decimal cero) que indica los días
transcurridos desde el 1/01/1900 hasta la fecha indicada.
– Hora: fracción decimal (parte entera cero) que tiene como unidad el
día (1 equivale a 24 horas).
Excel 52
Formato de celdas6
Excel 53
Formato de celdas7
• Otros Formatos:
– Porcentaje: Multiplica el valor de la celda por 100 y añade el
símbolo porcentual (%).
– Fracción: Muestra los números en forma de fracción.
– Científica: Parte entera y decimal seguido de la letra E y de un
entero que indica el exponente de 10.
– Texto: Se presenta tal como se introduce el texto.
– Especial: Se usa para números que representan determinados
datos (código postal y teléfono).
– Personalizada: Se escribe el formato que se ajusta a nuestra
necesidades adaptando los códigos predefinidos. Códigos #, 0, ?
Excel 54
Formato de celdas8
Excel 55
Formato de celdas9
• Otras opciones:
– Alineación: permite modificar y establecer la Alineación del texto,
Orientación, Control del texto y Dirección del texto.
– Fuente: permite modificar la Fuente (tipo de letra), Estilo, Tamaño,
Subrayado, Color y Efectos.
– Bordes: permite aplicar distintos tipos de bordes a una celda.
– Relleno: permite dar a las celdas distintos tipos de sombreado
(color del fondo) y de trama.
– Proteger: permite bloquear y ocultar celdas. Para que que este tipo
de formato tenga efecto es necesario activar la opción Revisar →
Proteger hoja (asignar contraseña).
Excel 56
Otras opciones de Formato10
Excel 57
Operaciones en rango de celdas
Excel 58
Operaciones en rango de celdas
Excel 59
Nombres de celdas1
Excel 60
Nombres de celdas2
Excel 61
Nombres de constantes o fórmulas
Excel 62
Copiar y pegar celdas
• Celdas
– Seleccionar la celda o bloque de celdas
donde se desea insertar.
– Ejecutar Inicio → Celdas → Insertar.
Se abre la ventana Insertar celdas.
– Seleccionar la opción que interesa. Pulsar
Aceptar.
• Filas (columnas)
– Seleccionar la fila(s) (columna(s)) donde se desea
insertar.
– Ejecutar Insertar → Filas (Columnas).
Excel 65
Eliminar - Deshacer
Excel 67
Serie de datos o fechas1
Excel 68
Serie de datos o fechas2
Excel 69
Serie de datos o fechas3
Excel 70
Listas personalizadas1
Excel 71
Listas personalizadas2
Excel 72
Referencias de celda1
Excel 73
Referencias relativas2
Excel 74
Referencias absolutas3
Excel 75
Referencias mixtas4
Excel 76
Validación de datos1
Excel 77
Validación de datos2
• Ejemplo: Operadores_formatos_graficos_etc.xlsx
Excel 78
Otras opciones del menú Archivo
Excel 79
Otras opciones
Excel 81
Funciones1
Excel 82
Funciones2
Excel 83
Funciones2
Excel 84
Funciones3
Excel 85
Categorías de Funciones
• Funciones matemáticas y trigonométricas
• Funciones lógicas
• Funciones estadísticas
• Funciones financieras
• Funciones de búsqueda y referencia
• Funciones de información
• Funciones de texto
• Funciones de ingeniería
• Funciones de fecha y hora
• Funciones de base de datos
• Funciones de compatibilidad
• Funciones de cubo
• Funciones web
• Funciones definidas por el usuario instaladas con complementos
Excel 86
Funciones Lógicas
Excel 87
Funciones Fecha y Hora
Excel 88
Funciones Búsqueda y Referencia
Excel 89
Funciones Financieras
Excel 90
Funciones Matemáticas y trigonométricas
Excel 92
Funciones Información y Texto
Excel 93
Funciones de Ingeniería
Excel 94
Funciones Base de datos
Excel 95
Gráficos
Excel 96
Gráficos de datos1
Excel 97
Gráficos de datos2
Excel 98
Elementos de los gráficos
1. El área del gráfico.
2. El área de trazado del gráfico.
3. Los puntos de datos de la serie de
datos que se trazan en el gráfico.
4. Los ejes horizontal (categorías) y
vertical (valores) en los que se
trazan los datos del gráfico.
5. La leyenda del gráfico.
6. Un título de eje y de gráfico que
puede agregar al gráfico.
7. Una etiqueta de datos que puede
usar para identificar los detalles de
un punto de datos de una serie de
datos.
Excel 99
Hojas de gráfico y Gráfico Incrustado
Excel 100
Tipos de gráficos1
• Tipos estándar:
– Columna y Barra, adecuados para comparar categorías.
– Línea, apropiado para mostrar la tendencia de una serie de valores
medidos a intervalos regulares de tiempo.
– Circular, usados para representar las distintas partes que
componen un total.
• Anillos, equivalente al gráfico circular, pero adaptado para representar varias
series de datos.
– Área, iguales a los de líneas, pero rellenan los espacios
comprendidos entre las líneas que representan los valores.
– XY (dispersión), adecuado para representar pares de valores.
• Burbujas, similar al de dispersión pero con un valor adicional para tamaño
del marcador.
Excel 101
Tipos de gráficos2
• Tipos estándar:
– Cotizaciones, gráficos específicos para representra cotizaciones
de valores bursátiles.
– Superficie, crea superficies 3D o curvas de nivel en superficies.
– Radial, radial con marcadores en cada valor de datos.
– Cuadro combinado, permite combiner dos tipos diferentes de
gráficos en uno solo.
• Para cada tipo estándar existen subtipos o variants del
mismo.
Excel 102
Ejemplos de gráficos de datos
Excel 103
Ejemplos de gráficos de datos
Excel 104
Ejemplos de gráficos de datos
Excel 105
Gráficos – Ejes múltiples
Excel 107
Gráficos – Ejes múltiples
Excel 108
Gráficos – tipo radial
Excel 110
Gráficos – Superficies 3D
Excel 111
Gráficos – Superficies 3D
Excel 112
Gráficos – Contorno
Excel 113
Gráficos – Contorno
Excel 114
Gráficos – Combinar tipos
Excel 115
Gráficos – Combinar tipos
Excel 116
Gráficos – Anotaciones
Excel 117
Gráficos – Anotaciones
Excel 118
Minigráficos
Excel 119
Funciones matemáticas
Excel 120
Funciones de suma
Excel 121
Otras divisiones y multiplicaciones
Excel 122
Otras divisiones y multiplicaciones
Excel 123
Otras divisiones y multiplicaciones
Excel 124
Funciones exponenciales y logaritmicas
Excel 126
Funciones de redondeo y truncamiento
Excel 127
Funciones de conversión de sistemas
numéricos
• Problema: Se requiere convertir un número de una base
a otra.
• Ejemplo: Funciones_matematicas.xlsx
Excel 128
Funciones de números complejos
COMPLEJO Convierte coeficientes reales e imaginarios en un número complejo.
• Problema: Se requiere realizar [Link] Devuelve el valor absoluto (módulo) de un número complejo.
IMAGINARIO Devuelve el coeficiente imaginario de un número complejo.
cálculos con números [Link] Devuelve el argumento theta, un ángulo expresado en radianes.
complejos.
[Link] Devuelve la conjugada compleja de un número complejo.
[Link] Devuelve el coseno de un número complejo.
[Link] Devuelve el coseno hiperbólico de un número complejo.
• Ejemplo: IMCOT Devuelve la cotangente de un número complejo.
Funciones_matematicas.xlsx [Link]
[Link]
Devuelve la cosecante de un número complejo.
Devuelve la cosecante hiperbólica de un número complejo.
[Link] Devuelve el cociente de dos números complejos.
[Link] Devuelve el valor exponencial de un número complejo.
[Link] Devuelve el logaritmo natural (neperiano) de un número complejo.
IM.LOG10 Devuelve el logaritmo en base 10 de un número complejo.
IM.LOG2 Devuelve el logaritmo en base 2 de un número complejo.
[Link] Devuelve un número complejo elevado a una potencia entera.
[Link] Devuelve el producto de 2 a 255 números complejos.
[Link] Devuelve el coeficiente real de un número complejo.
[Link] Devuelve la secante de un número complejo.
[Link] Devuelve la secante hiperbólica de un número complejo.
[Link] Devuelve el seno de un número complejo.
[Link] Devuelve el seno hiperbólico de un número complejo.
IM.RAIZ2 Devuelve la raíz cuadrada de un número complejo.
[Link] Devuelve la diferencia entre dos números complejos.
[Link] Devuelve la suma de números complejos.
[Link] Devuelve la tangente de un número complejo.
Excel 129
Funciones para cálculos con matrices
Excel 131
Análisis estadístico de datos
Excel 134
Estadística descriptiva
Excel 135
Estadística descriptiva
Excel 136
Resumen de las funciones estadísticas de
Estadística descriptiva
Estadístico Función Excel
Media =PROMEDIO(Datos)
Error típico =Desviación estándar/RAIZ(Cuenta)
Mediana =MEDIANA(Datos)
Moda =MODA(Datos)
Desviación estándar =DESVEST(Data)
Varianza de la muestra =VAR(Datos)
Curtosis =CURTOSIS(Datos)
Coeficiente de asimetría =[Link](Datos)
Rango =Máximo - Mínimo
Mínimo =MIN(Datos)
Máximo =MAX(Datos)
Suma =SUMA(Datos)
Cuenta =CONTAR(Datos)
Mayor (1) =[Link](Datos,1)
Menor(1) =[Link](Datos,1)
Nivel de confianza(95.0%) =[Link](0.05, Desv est.,100)
Excel 137
Distribuciones de frecuencia - Histograma
Excel 138
Histograma – Frecuencia
Excel 140
Histograma – Frecuencia
Excel 141
Análisis de datos – Histograma
Excel 142
Análisis de datos – Histograma
Excel 143
Intervalos de confianza
Excel 144
Intervalos de confianza
Excel 148
Pruebas estadísticas
Excel 149
Analysis de Variance - ANOVA
Excel 152
ANOVA de un factor
>
< α = 0.05
Excel 153
Generación de números aleatorios
Excel 155
Serie de números aleatorios
Excel 156
Datos de muestra
Excel 157
Distribuciones de probabilidad
Excel 158
Aproximación
Excel 159
Ajuste de ecuaciones a datos
• Es muy común en ingeniería intentar hallar la
ecuación (curva) que mejor aproxime un conjunto de
datos.
• Datos: pares de puntos P1=(x1,y1)… Pn=(xn,yn) o
tuplas de variables independientes y dependientes.
• Se trata de pasar una curva a través del conjunto de
datos. Cuando los resultados se usan para hacer
nuevas predicciones de variables dependientes, se
conoce como regresión.
• Se usa el método de mínimos cuadrados: se
fundamenta en la minimización del error ei = yi – f(xi)
obtenido para cada punto.
• Ejemplos: Ajuste_curvas.xlsx
Excel 160
Ajuste lineal por MMC a datos1
SST
n n
SSE = ∑ [ yi − f ( xi )]2 SST = ∑ [ yi − y ]2
i =1 i =1
Excel 161
Ajuste lineal por MMC a datos1
2.0
1.5
1.0
0.5
Y
0.0
estimado
-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0
-0.5
-1.0
-1.5
-2.0
Excel 162
Ajuste lineal por MMC a datos2
• Otro método rápido de obtener un ajuste lineal (y de otro tipo)
a un conjunto tabulado en columnas de datos x (variable
independiente) e y (variable dependiente) es:
• Graficar los datos como tipo de gráfico X-Y (dispersión) como
puntos.
• Pulsar en uno de los puntos dato para seleccionar como
objeto activo el conjunto de datos y pulsar el botón derecho
del ratón para obtener el menú Gráfico.
• Seleccionar Añadir Línea de Tendencia en el menú Gráfico.
Especificar el tipo de curva (Lineal) y llenar las opciones
correspondientes. Conviene seleccionar en Opciones
Presentar ecuación en el gráfico y el valor R (coeficiente de
correlación). Es posible realizar extrapolación.
Excel 163
Ajuste lineal por MMC a datos2
Excel 164
Ajuste lineal por MMC a datos2
10
8
Fuerza, N
0
0 2 4 6 8 10 12 14 16 18
Desplazamiento desde la posición de equilibrio, cm
Excel 165
Ajuste lineal por MMC a datos2
1.5
1.0
0.5
Desviación (grados)
0.0
Lineal (Desviación
-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 (grados))
-0.5
-1.0
-1.5
-2.0
Excel 166
Ajuste lineal de datos3
• Otro método rápido de obtener un ajuste lineal es utilizar la
función [Link] a un conjunto tabulado en
columnas de datos x (variable independiente) e y (variables
dependientes).
• La sintaxis de la función corresponde a una fórmula matriz
(hay que pulsar Ctrl-Mayus-Entrar) cuando se introduce la
fórmula. Ejemplo:
{=[Link](C5:C13, A5:A13, VERDADERO, VERDADERO)}
Y X Para que calcule la
intersección y las
estadísticas ampliadas
Excel 167
Ajuste lineal de datos3
Excel 168
Ajuste multilineal de datos4
• Para hacer un ajuste lineal múltiple del tipo
y = m1x1 + m2x2 + m3x3 + … mnxn + b
también se puede utilizar la función [Link] a
un conjunto tabulado en columnas de datos x (variables
independientes) e y (variables dependientes).
• Es necesario seleccionar una cuadrícula de celdas de
tamaño n+1 columnas, donde n es el número de variables
independientes (x) y 5 filas.
• La sintaxis de la función es de tipo matriz (hay que pulsar
Ctrl-Mayus-Entrar) cuando se introduce la fórmula. Ejemplo:
{=[Link](B12:B27,C12:H27,VERDADERO,VERDADERO)}
Excel 169
Ajuste multilineal de datos4
72000
70000
68000
66000
64000 y
62000 y-est
60000
58000
56000
54000
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Excel 170
Regresión
Excel 171
Regresión
Excel 172
Otros tipos de ajuste
• Exponencial
• Potencial
• Polinómico: es necesario dar el orden del polinomio
– Ejemplos
Excel 173
Selección de la mejor curva de
ajuste a un conjunto de datos
• Método de prueba y error. Primero se grafican los datos como una
línea recta.
• Si no se obtiene un buen ajuste, intentar diferentes tipos de
curvas, usando evaluación visual ayudado por los resultados de la
suma de cuadrados de los errores y el coeficiente de correlación
(r2).
• Si no se obtienen resultados satisfactorios, intentar graficar los
datos de otra manera (y – 1/x, 1/y-x, etc.)
• En algunos casos se consiguen mejores ajustes escalando los
datos (datos de x e y del mismo orden de magnitud).
• Cambio de escala (se obtiene una recta) para el paso 2:
– Exponencial y = a ebx log y vs. x (semi-log)
– Logarítmico y = a ln x + b y vs. log x (semi-log)
– Potencial y = a xb log y vs. log x (log – log)
Excel 174
Ajuste exponencial de datos
10
8
1
Voltios
0 2 4 6 8 10 12
Voltios
6
4
0
2
0
0 2 4 6 8 10 12 0
Tiempo, segundos Tiempo, segundos
Excel 175
Ajuste logarítmico de datos
55 55
50 50
Temperatura ºC
Temperatura ºC
45 45
40 40
35 35
30 30
25 25
20 20
0 500 1000 1500 2000 2500 0.1 1 10 100 1000 10000
Distancia cm. Distancia cm.
Excel 176
Ajuste potencial de datos
2.5 1.0
TR, moles/sec
TR, moles/sec
1 10 100
2.0
1.5 0.1
1.0
0.0
0.5
0.0
0 20 40 60 80 100 0.0
C, moles/cu ft C, moles/cu ft
Excel 177
Ajuste polinomial de datos
Tiempo de aceleración vs [Link].
y = 1E-08x - 5E-06x4 + 0.0008x3 - 0.0525x2 + 1.777x - 20.958
5
2
R = 0.9998
50
45
40
35
Tiempo, seg
30
25
20
15
10
5
0
0 20 40 60 80 100 120 140 160
Velocidad máx. Km/h
Excel 178
Análisis de series de tiempo
Excel 179
Análisis de series de tiempo
Excel 180
Análisis de series de tiempo -
Visualización
• Problema: Graficar un grupo de datos de series de
tiempo para análisis posteriores.
• Se utiliza el asistente para gráficos.
Excel 181
Análisis de series de tiempo -
Visualización
Temperatura anual (Los Angeles)
69.0
68.5
68.0
67.5
Temperatura ºF
67.0
66.5
66.0
65.5
65.0
64.5
64.0
1920 1930 1940 1950 1960 1970 1980 1990 2000 2010
Año
Excel 182
Análisis de series de tiempo - tendencia
Excel 183
Análisis de series de tiempo - tendencia
67.0
66.5
66.0
65.5
65.0
64.5
64.0
1920 1930 1940 1950 1960 1970 1980 1990 2000 2010
Año
Excel 184
Análisis de series de tiempo - tendencia
67.0
66.5
66.0
65.5
65.0
64.5
64.0
1920 1930 1940 1950 1960 1970 1980 1990 2000 2010
Año
Excel 185
Análisis de series de tiempo –
medias móviles
• Problema: Suavizar una serie de tiempo mediante
medias móviles.
• Se puede calcular medias móviles de varias formas:
– usando la función de gráficos Media móvil de la
línea de tendencia
– usando la función Media móvil de las
Datos→Análisis de datos.
Excel 186
Análisis de series de tiempo –
medias móviles línea de tendencia
• Para la Línea de tendencia, se selecciona la serie de
datos y haciendo clic con el botón derecho del ratón
se selecciona Agregar línea de tendencia.
• Se selecciona Media móvil y el período deseado (3).
En opciones se puede escribir un nombre para esta
nueva línea de tendencia.
• Es sencillo pero no genera datos numéricos.
Excel 187
Análisis de series de tiempo –
medias móviles
Temperatura anual (Los Angeles)
69.0
68.5
68.0
67.5
Temperatura ºF
67.0
66.5
66.0
65.5
65.0
64.5
64.0
1920 1930 1940 1950 1960 1970 1980 1990 2000 2010
Año
Excel 188
Análisis de series de tiempo –
medias móviles Análisis de datos
• En Herramientas→Análisis de datos hay una opción
que es Media móvil.
• La ventaja es que genera datos numéricos para la
serie de media móvil.
• En la ventana de diálogo se selecciona el rango de
las celdas que contiene la serie de datos.
• Se introduce el intervalo (3) sobre el que se desea
calcular las medias.
• Introducir la celda donde se desea colocar los
resultados.
Excel 189
Análisis de series de tiempo –
medias móviles Análisis de datos
Media móvil
70.0
69.0
68.0
67.0
Valor
Real
66.0
65.0 Pronóstico
64.0
63.0
62.0
1 6 11 16 21 26 31 36 41 46 51 56 61 66
Punto de datos
Excel 190
Análisis de series de tiempo –
índices estacionales
• Problema: calcular los índices estacionales de una
serie de tiempo que muestra variaciones estacionales.
• Hay varios métodos para calcular el índice estacional
de una serie. Aquí se muestra el método promedio-
porcentaje.
• Primero se calcula el promedio de la variable cada año.
• Despúes se calcula el porcentaje de cada mes
respecto al promedio anual.
• Finalmente se calcula el promedio de los porcentajes
de cada mes para todos los años. Para comprobar el
promedio de los índices debe ser 1.
Excel 191
Análisis de series de tiempo –
índices estacionales
1996 1997 1998 1999
Ene 49.1 49.2 53.5 54.0
Feb 53.0 53.9 53.6 57.9
Mar 54.7 63.8 58.1 58.6
Abr 64.5 62.2 64.7 70.9
May 76.7 72.4 77.0 73.7
Jun 79.1 78.7 83.5 79.9
Jul 82.2 83.1 85.5 82.2
Ago 80.2 81.3 83.8 85.0
Sep 75.8 78.4 80.5 75.8
Oct 66.9 67.1 70.0 66.7 Indices Estacionales
Nov 58.9 55.1 61.5 59.7
Dic 54.2 49.2 53.5 51.0
Average: 66.3 66.2 68.8 68.0
Indice Estacional
1996 1997 1998 1999 Indice
Ene 0.74 0.74 0.78 0.79 0.76
Feb 0.80 0.81 0.78 0.85 0.81
Mar 0.83 0.96 0.84 0.86 0.87
Abr 0.97 0.94 0.94 1.04 0.97
May 1.16 1.09 1.12 1.08 1.11
Jun 1.19 1.19 1.21 1.18 1.19
Jul 1.24 1.26 1.24 1.21 1.24
Ago 1.21 1.23 1.22 1.25 1.23
Sep 1.14 1.18 1.17 1.12 1.15
Oct 1.01 1.01 1.02 0.98 1.01
Nov 0.89 0.83 0.89 0.88 0.87
Dic 0.82 0.74 0.78 0.75 0.77
Sum: 12.00
Average: 1.00
Excel 192
Análisis de series de tiempo –
Transformada discreta de Fourier
• Problema: usar la transformada discreta de Fourier para
analizar un conjunto de datos.
• En Herramientas→Análisis de datos hay una opción que
es Análisis de Fourier que permite realizar
transformaciones discretas de Fourier (DFT) y
transformaciones inversas.
• El tamaño de la serie debe ser potencia de 2 con un
tamaño máximo de 2^12 = 4096.
• Cálculo de la frecuencia f i = i (ns ) en Hz donde i es el
número de la muestra, n el número de muestras y s el
intervalo de muestra.
Excel 193
Análisis de series de tiempo –
Transformada discreta de Fourier
• Cálculo de la frecuencia f i = i n en ciclos por muestra
donde i es el número de la muestra, n el número de
muestras y s el intervalo de muestra.
• La DFT se obtiene seleccionando el rango de las celdas
que contienen la serie.
• Los resultados de la DFT son números complejos que se
escriben como texto. Para manipularlos Excel dispone de
funciones para ellos.
• Cálculo de la potencia en cada banda de frecuencia
hasta la frecuencia de Nyquist (0.5 ciclos/muestra).
[Link](DFT)^2/n^2
Excel 194
Análisis de series de tiempo –
Transformada discreta de Fourier
• Se puede filtrar los datos en el campo de la frecuencia
para aislar un determinado componente.
• Se construye un filtro adecuado.
EXP(-((ABS(frec cs) – fo)/sig)^2))
• Se aplica el filtro a la DFT multiplicándolo (usar funciones
de números complejos).
• Se calcula DFT inversa y se obtiene la serie numérica
utilizando funciones de números complejos.
[Link](InversaDFT)
Excel 195
Análisis de series de tiempo –
Transformada discreta de Fourier
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
-0.5
-1
-1.5
-2
Y(t) Filtrada
Excel 196
Resolviendo ecuaciones
Excel 197
Resolviendo ecuaciones
Excel 198
Resolviendo ecuaciones
Excel 200
Resolviendo ecuaciones –
método gráfico
x f(x)
0 -5.00 Raíz real de un polinomio
0.1 -5.03
0.2 -5.12
10.00
0.3 -5.27
0.4 -5.46 8.00
0.5 -5.69
0.6 -5.92 6.00
0.7 -6.13 4.00
0.8 -6.26
0.9 -6.25 2.00
f(x)
1 -6.00
1.1 -5.41 0.00
1.2 -4.34 1.00 1.10 1.20 1.30 1.40 1.50 1.60
-2.00
1.3 -2.64
1.4 -0.12 -4.00
1.5 3.44
-6.00
1.6 8.29
1.7 14.73 -8.00
1.8 23.07
X
1.9 33.69
2 47.00
Excel 201
Resolviendo ecuaciones –
método gráfico
Raíces de un polinomio cúbico
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
-2
-4
Y
-6
-8
-10
-12
X
Excel 202
Resolviendo ecuaciones –
usando Buscar objetivo
• Se puede obtener una solución rápida de ecuaciones
algebraicas simples usando la opción Buscar Objetivo en
el menú Herramientas.
• Para ello se sigue:
– Escribir un valor inicial de x en una celda.
– Escribir la fórmula de la ecuación en la forma f(x)=0 en otra
celda. Escribir la variable x como referencia a la celda que
contiene el valor inicial.
– Seleccionar Buscar Objetivo en el menú Herramientas.
– En el diálogo escribir la dirección de la celda que contiene la
fórmula, el valor 0 en Valor y la dirección de la celda que
contiene el valor inicial. Pulsar Aceptar.
Excel 203
Resolviendo ecuaciones –
usando Buscar objetivo
• Ejemplo: f(x) = 2*x5 – 3* x2 – 5 = 0
x= 1.40408295
Excel 204
Resolviendo ecuaciones –
usando Solver
• Solver se usa para resolver problemas de más complejidad y
se puede configurar el método y visualización de solución.
• Instalar Solver desde Complementos en el menú
Herramientas.
• Para ello se sigue:
– Escribir un valor inicial de la variable x en una celda
– Escribir la fórmula de la ecuación en la forma f(x)=0 en otra celda.
Escribir la variable x como referencia a la celda que contiene el valor
inicial.
– Seleccionar Solver en el menú Herramientas.
– En el diálogo escribir la dirección de la celda que contiene la fórmula,
el valor 0 en Valor y la dirección de la celda que contiene el valor
inicial.
– Si se desea restringir el rango de x, pulsar en Agregar.
– Pulsar Resolver. Se puede configurar con Opciones.
Excel 205
Resolviendo ecuaciones –
usando Buscar objetivo
• Ejemplo: f(x) = 2*x5 – 3* x2 – 5 = 0
x= 1.40408619
Excel 206
Resolviendo sistemas de ecuaciones
Excel 207
Operaciones matriciales en Excel
[M] = 1 4 3
5 2 6
2 3 -4
[A] = 3 -1 -2 7 12
4 -7 -6 [N] = 11 8
9 10
Det[A] = 82
[M][N] = 78 74
-0.09756098 0.56097561 -0.12195122
111 136
Inverse [A] = 0.12195122 0.04878049 -0.09756098
-0.20731707 0.31707317 -0.13414634
1.90243902 3.56097561 -4.12195122
A + Trans[A] = 3.12195122 -0.95121951 -2.09756098
2 3 4
3.79268293 -6.68292683 -6.13414634
Transpose [A] = 3 -1 -7
-4 -2 -6 0 0 -8
A - Trans[A] = 0 0 5
8 -5 0
Excel 210
Operaciones matriciales en Excel
Excel 211
Operaciones matriciales en Excel
Excel 212
Solución de sistemas de ecuaciones
lineales mediante matrices
• Un método para resolver un sistema de ecuaciones
lineales simultáneas [A][x] = [b] es mediante métodos
matriciales [x] = [A]-1 [ b]
• Para resolver en Excel:
– Escribir los elementos de la matriz A.
– Escribir los elementos del vector b.
– Seleccionar las celdas para la inversa A-1. Escribir la fórmula
en la celda superior izquierda =MINVERSA() y seleccionar
las celdas de A. Pulsar Ctrl-Mayus-Enter simultáneamente.
– Seleccionar las celdas donde se desea aparezca el vector x.
Escribir la fórmula en la celda superior =MMULT() y
seleccionar las celdas de A-1 y B . Pulsar Ctrl-Mayus-Enter
simultáneamente.
Excel 213
Solución de sistemas de ecuaciones
lineales mediante matrices
Solución de Ecuaciones Simultáneas mediante inversión matricial
[A][x] = [b]
9.233
[b] = 8.205
3.934
0.896
[x] = Inv[A] [b] = 0.765
0.614
Excel 214
Solución de sistemas de ecuaciones
usando Solver1
• Solver ofrece un enfoque diferente para resolver
sistemas ecuaciones simultáneas lineales o no lineales.
• Suponiendo que se tiene un sistema de n ecuaciones y n
incógnitas representados mediante las ecuaciones:
f1 ( x1 , x2 , , xn ) = 0
f 2 ( x1 , x2 , , xn ) = 0
f n ( x1 , x2 , , xn ) = 0
• Se desea hallar los valores de x1, x2 …, xn que produce
que cada ecuación sea cero. Una forma para hacer esto
es forzar a que la función (varianza residual):
y =f12 + f22 + …+ fn2 sea cero.
Excel 215
Solución de sistemas de ecuaciones
usando Solver2
• Procedimiento:
– Escribir un valor inicial para cada variable independiente x1, x2 …, xn
en celdas diferentes
– Escribir las ecuaciones f1 , f2 , …, fn e y en celdas diferentes
expresadas como fórmulas dependientes de las celdas donde están
las variables x1, x2 …, xn
– Seleccionar Solver de la Barra de herramientas. Dar la dirección de
la celda que contiene la fórmula de y para Celda Objetivo.
Seleccionar Valores de 0. En Cambiando las celdas dar el rango de
las celdas que contienen los valores iniciales de las variables x1, x2
…, xn.
– Se puede restringir opcionalmente el rango de los valores de las
variables independientes pulsando Agregar.
– Se puede seleccionar la opción de generar Resultados en otra hoja.
Excel 216
Solución de sistemas de ecuaciones
usando Solver3
• Ejemplo 1 (lineal): x1 = 1,428567883
3 x1 + 2 x2 – 2 x3 = 4
x2 = 4,142852593
x3 = 4,28570703
2 x1 – x2 + x3 = 3
x1 + x2 – 2 x3 = – 3 f(x1,x2,x3) = -5,2244E-06
g(x1,x2,x3) = -9,797E-06
h(x1,x2,x3) = 6,41664E-06
f = 3 x1 + 2 x2 – 2 x3 – 4 = 0 y = 1,64449E-10
g = 2 x1 – x2 + x3 – 3 = 0
h = x1 + x2 – 2 x3 + 3 = 0
y = f2 + g2 + h2
Excel 217
Solución de sistemas de ecuaciones
usando Solver4
• Ejemplo 2 (no lineal):
x12 + 2 x22 – 5 x1 + 7 x2 = 40
3 x12 – x22 + 4 x1 + 2 x2 = 28
y= 4,66E-07
Excel 218
Solución de sistemas de ecuaciones
usando Solver4
• Ejemplo 3 (no lineal):
0.5
y = 1 − e− x 9y^2+4x^2=1
0.4
y=1-e^(-x)
9 y + 4x = 1
2 2 0.3
0.2
0.1
0
-0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5
-0.1
7.5386E-07 1.9462E-07
-0.5
Excel 219
Evaluación de Integrales
Excel 220
Evaluación de Integrales
Excel 222
Evaluación de Integrales –
método trapezoidal
• Ejemplo:
– La presión media de una gas cuando la temperatura del gas varía
en el tiempo se calcula por:
100
Excel 223
Evaluación de Integrales –
método trapezoidal
Integración numérica usando la regla trapezoidal con datos espaciados uniformemente
tiempo Temperatura
0 300
5 360 360
10 420 420
15 480 480
20 540 540
25 600 600
30 660 660
35 720 720
40 780 780
45 840 840
50 900 900
55 960 960
60 1020 1020
65 1080 1080
70 1140 1140
75 1200 1200
80 1260 1260
85 1320 1320
90 1380 1380
95 1440 1440
100 1500
Suma = 17100 Integral = 90000 Presión = 73.547
Excel 224
Evaluación de Integrales –
método trapezoidal
• Regla trapezoidal (datos espaciados no uniformes):
– Datos: n pares de puntos (x1,y1), (x2,y2),… (xn,yn) donde x1 = a y x2
=b
– Estos puntos definen n-1 intervalos rectangulares con un ancho
para el i-ésimo intervalo ∆xi = xi+1 – xi
y +y
– La altura de cada intervalo se expresa como: yi = i i +1
2
– Por tanto la integral se aproxima como:
1 n −1
b
I = ∫ y dx = ∑ ( yi +1 + yi )( xi +1 − xi )
a
2 i =1
Excel 225
Evaluación de Integrales –
método trapezoidal
• Ejemplo:
– La corriente por una inductancia se puede obtener con la fórmula
t
1
i = ∫ vdt
L0
donde: i = corriente (amperios), L = inductancia (henrios),
v = voltaje (voltios) y t=tiempo (seg).
Se induce una corriente de 2.15 amperios por un periodo de 500 milisegundos. La
variación del voltaje con el tiempo en este periodo se muestra en la siguiente
tabla:
t (msec) v (volts) t v t v t v
0 0 40 45 90 45 180 27
5 12 50 49 100 42 230 21
10 19 60 50 120 36 280 16
20 30 70 49 140 33 380 9
30 38 80 47 160 30 500 4
Excel 226
Evaluación de Integrales –
método trapezoidal
Integración numérica usando la regla trapezoidal con datos espaciados no uniformes
Excel 227
Evaluación de Integrales –
método Simpson
• Regla de Simpson (número de datos impar – número de
subintervalos par):
– En lugar de considerar rectángulos entre los puntos, se pasa un
polinomio de segundo orden (parábola) a través de tres puntos
adyacentes igualmente espaciados.
– Por tanto la integral se aproxima como:
b
1
I = ∫ y dx = ( y1 + 4 y2 + 2 y3 + 4 y4 + 2 y5 + + 2 yn − 2 + 4 yn −1 + yn ) ∆x
a
3 1
• Ejemplo: Evaluar la integral I = e dx
∫
− x 2
0
en el rango de 0 a 1 con un ∆=0.1 entre puntos.
Excel 228
Evaluación de Integrales –
método Simpson
Integración numérica usando la regla de Simpson 1
x y
I = ∫e − x2
dx
0.0 1.0000 1.0000 0
0.1 0.9900 3.9602
0.2 0.9608 1.9216 Funcion Exp(-x^2)
Y
0.3 0.9139 3.6557
1.2000
0.4 0.8521 1.7043
0.5 0.7788 3.1152 1.0000
0.6 0.6977 1.3954
0.8000
0.7 0.6126 2.4505
0.8 0.5273 1.0546 0.6000
0.9 0.4449 1.7794
1.0 0.3679 0.3679 0.4000
0.2000
SUMA = 22.4047
0.0000 X
INTEGRAL = 0.7468 0.0 0.2 0.4 0.6 0.8 1.0 1.2
Excel 229
Cálculo de la superficie y centroide
mediante Integrales
• Problema: calcular la superficie y centroide de una función
dada como tabla.
• Se aplica una técnica de integración numérica para el cálculo
de la superficie y el primer momento para el centroide.
• El centroide se calcula tomando el primer momento del área
calculado con:
M x = ∫ ydA = ∫ yxdy M y = ∫ xdA = ∫ xydx
My
• De donde se obtiene el centroide: xc = yc =
Mx
A A
Excel 230
Cálculo del segundo momento
de una superficie
• Problema: calcular el segundo momento de una superficie
(momento de inercia).
• Se usa la misma técnica que la sección anterior para el
primer momento, pero usando x2 e y en lugar de x e y. Es
decir para el eje y: I y = x 2 dA = x 2 ydx
∫ ∫
• El momento de inercia para un eje que pasa por el centro de
la superficie se calcula aplicando el teorema del eje paralelo,
I na = I y − Ad 2
donde Ina es el momento de inercia del área sobre el eje
paralelo a y que pasa por el centroide, A es el área y d es la
distancia al eje y.
Excel 231
Cálculo del centroide y segundo
momento de una superficie
Cálculo del Centro y Momento de Inercia de un Area usando integración numérica
x y [Link]
0.000 0.100 1
0.200 0.300 4 1.200
0.400 0.600 2
0.600 0.900 4 1.000
0.800 1.050 2
1.000 1.000 4
1.200 0.700 2 0.800
1.400 0.400 4
1.600 0.200 2 0.600
1.800 0.100 4
2.000 0.050 1
0.400
s= 0.2
0.200
Area = 1.070
Centro del Area
xc = 0.869
yc = 0.382 0.000
0.000 0.500 1.000 1.500 2.000
Iy = 0.97013333
Iyc = 0.16297404
Excel 232
Cálculo de Integrales dobles
Excel 233
Cálculo del volumen bajo una superficie
Cálculo de Volumen mediante Integrales dobles
sx = 0.1
sy = 0.1
y
Coeficientes x 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0.5 0 1.000 0.990 0.961 0.914 0.852 0.779 0.698 0.613 0.527 0.445 0.368
1 0.1 0.990 0.980 0.951 0.905 0.844 0.771 0.691 0.607 0.522 0.440 0.364
1 0.2 0.961 0.951 0.923 0.878 0.819 0.748 0.670 0.589 0.507 0.427 0.353
1 0.3 0.914 0.905 0.878 0.835 0.779 0.712 0.638 0.560 0.482 0.407 0.336
1 0.4 0.852 0.844 0.819 0.779 0.726 0.664 0.595 0.522 0.449 0.379 0.313
1 0.5 0.779 0.771 0.748 0.712 0.664 0.607 0.543 0.477 0.411 0.346 0.287
1 0.6 0.698 0.691 0.670 0.638 0.595 0.543 0.487 0.427 0.368 0.310 0.257
1 0.7 0.613 0.607 0.589 0.560 0.522 0.477 0.427 0.375 0.323 0.273 0.225
1 0.8 0.527 0.522 0.507 0.482 0.449 0.411 0.368 0.323 0.278 0.235 0.194
1 0.9 0.445 0.440 0.427 0.407 0.379 0.346 0.310 0.273 0.235 0.198 0.164
0.5 1 0.368 0.364 0.353 0.336 0.313 0.287 0.257 0.225 0.194 0.164 0.135
Areas: 0.746 0.739 0.717 0.682 0.636 0.581 0.521 0.457 0.393 0.332 0.275
Coeficientes: 0.5 1 1 1 1 1 1 1 1 1 0.5
Volumen = 0.55683
Curva de Areas
1.000
0.900 0.800
0.800 0.900-1.000 0.700
0.700 0.800-0.900
0.600
0.700-0.800
0.600
0.600-0.700 0.500
z 0.500 Area
0.500-0.600 0.400
0.400 0.400-0.500
0.300 0.300-0.400
0.300
0.200 0.200-0.300 0.200
0 0.100-0.200
0.100 0.100
0.000 0.4 0.000-0.100
x 0.000
0
0.1
0.2
0.8
0.3
y
0.9
Y
1
Excel 234
Optimización
Excel 235
Optimización
• Los problemas de ingeniería se modelan mediante
sistemas de ecuaciones para su análisis.
• Muchos problemas de ingeniería requieren la
optimización de un criterio, como el costo, ganancia,
peso, etc, al que se llama función objetivo.
• Adicionalmente hay una serie de condiciones, tales
como leyes de conservación, restricciones de
capacidad u otra restricción técnica, que deben ser
satisfechas. Estas condiciones se llaman
restricciones.
Excel 236
Optimización
• El objetivo de una solución óptima es determinar una
solución que produce que la función objetivo sea
maximizada o minimizada cumpliendo todas las
restricciones.
• Los problemas de este tipo se conocen como
problemas de optimización.
Excel 237
Optimización - Formalización
Excel 238
Optimización usando Solver
• Procedimiento:
– Escribir un valor inicial para cada variable independiente
x1, x2 …, xn en celdas diferentes.
– Escribir la función objetivo como fórmula Excel en una celda.
– Escribir las ecuaciones de cada restricción como fórmulas Excel.
– Seleccionar Herramientas→Solver. Dar la dirección de la celda que
contiene la función objetivo para Celda objetivo. Seleccionar Máximo o
Mínimo en Valor de la celda objetivo. Dar el rango de las celdas que
contienen los valores iniciales de las variables
x1, x2 …, xn en Cambiando las celdas.
– Escribir las celdas que contienen cada restricción, el tipo de restricción
y el valor del lado derecho usando Agregar.
– Si la función objetivo y las restricciones son lineales, pulsar en el botón
Opciones y seleccionar Adoptar Modelo Lineal.
– Pulsar Aceptar y después Resolver. Se puede seleccionar la opción de
generar Resultados en otra hoja.
Excel 239
Optimización - Programación lineal clásica
Excel 240
Optimización - Programación lineal clásica
14
x2
12
R2
Región de soluciones
10
factibles
8
R1
Recta de la Función
6
Objetivo para 455.833
4
0
0 2 4 6 8 10 12 14 16 18 20
x1
Excel 241
Optimización - Programación lineal clásica
Solución (Solver) de problema de programación lineal clá
x1 = 6.6667
x2 = 5.8333
f(x1, x2) = 455.8333 ? Función Objectivo
Valor Límite
Restricción 1 60 60
Restricción 2 50 50
Excel 242
Optimización - Programación lineal clásica
Excel 243
Optimización - Programación lineal clásica
Excel 244
Optimización - no lineal
Excel 246
Optimización - no lineal
Excel 247
Optimización - no lineal
1 1.5
0.5 1
0 0.5
-0.5
0
3
-1 -3 -2 -1 0 1 2 3
1
-0.5
-1.5 -1
-3
-2.25
-1.5
-1
-0.75
0
-3
0.75
1.5
2.25
-1.5
Excel 248
Evaluación económica
Excel 249
Evaluación económica de alternativas
• Una parte importante en la evaluación de proyectos es
la evaluación económica.
• Se basa en el valor del dinero en el tiempo. La
terminología empleada es el principal para indicar la
cantidad prestada y el interés que es el pago adicional
por el uso del dinero.
• Los cálculos de interés se basan en la tasa de interés i.
• Los cálculos económicos se basan en el uso del interés
compuesto. Así para n períodos de interés, la cantidad
total de dinero acumulado al final del último período de
interés es: F = Fn = P(1 + i)n
• Ejemplos: Comparacion_Economica1.xls
Excel 250
Cálculos financieros básicos
• Problema: Calcular el capital acumulado para un depósito a
un interés y período dado.
Interés compuesto Cantidad acumulada
Final de año F Interés del
Acumulación horario = compuesto
interés 5436.55
0 2000.00
P= 2000
1 2100.00
2 2205.00 6000.00
i (anual) = 0.05 3 2315.25
4 2431.01 5000.00
Total acumulado
5 2552.56
n= 20 6 2680.19 4000.00
7 2814.20
3000.00
8 2954.91
9 3102.66
2000.00
10 3257.79
11 3420.68
1000.00
12 3591.71
13 3771.30
0.00
14 3959.86
0 2 4 6 8 10 12 14 16 18 20
15 4157.86
16 4365.75 Final de Año
17 4584.04
18 4813.24
19 5053.90
20 5306.60
Excel 251
Cálculos financieros básicos
Excel 252
Valor presente de un flujo de caja
0 1 2 3 n-1 n
A= -140,000.00 €
i= 0.08
n= 12
P= 1055050.92
Excel 254
Valor presente
A= 140000
i= 0.08
n= 12
P= -1,055,050.92 €
Excel 255
Valor futuro
A= 140000
i= 0.08
n= 12
F= 2,656,797.70 €
Excel 256
Flujos de caja no uniformes
i= 0.08
6 4000000 0
7 5000000 -2000000 0 1 2 3 4 5 6 7 8 9 10 11 12 13
8 6000000 -4000000
9 5000000
-6000000
10 4000000
-8000000
11 3000000
-10000000
12 2000000
13 1000000 -12000000
Final Año
VPN = 2,380,570.73
Excel 257
Comparación de Alternativas
Flujos de caja no uniformes
• Problema: Comparar varias alternativas de flujos de
caja. Se selecciona la de mayor Valor Presente Neto.
Comparación de dos oportunidades de inversión
i= 0.1
Excel 258
Comparación de Alternativas
Tasa interna de retorno (TIR)
• El método de la Tasa Interna de Retorno (TIR) es otro criterio
muy usado para comparar varias alternativas de inversión. A
diferencia del método del Valor Presente no hay necesidad
de especificar una tasa de interés.
• Si dibujamos el valor presente de un flujo de caja en función
de la tasa de interés, la TIR es el punto de cruce, es decir, el
valor de la tasa de interés al cual el valor presente neto se
hace cero.
• Durante la comparación de alternativas mediante la TIR se
escogerá aquella alternativa que tenga la mayor tasa interna
de retorno.
• Excel tiene la función TIR que calcula la tasa interna de
retorno directamente.
Excel 259
VPN - TIR
30,000
i VPN 20,000
0 65,000
10,000
0.03 46,639
0.06 31,057 0
0.09 17,751 -10,000 0 0.05 0.1 0.15 0.2 0.25
0.12 6,322 -20,000
0.15 -3,549
0.18 -12,119 Tasa de Interés
0.21 -19,597
Excel 260
Comparación de Alternativas
Tasa interna de retorno (TIR)
Comparación de dos oportunidades de inversión
i= 0.1
Excel 262
Conversión de Unidades
Excel 263
Conversión de Unidades
Excel 264
Conversiones simples
Excel 265
Conversiones simples
Excel 266
Conversiones simples en Excel
1 (lb f / in 2 )
2
6.3 lb f 4.44822 N 39.37 in
P= × × = 43437 N / m 2
in 2 1lb f 1m
Excel 269
Ejemplo de
Conversiones complejas
• =CONVERTIR(6.3,”lbf”,”N”)*CONVERTIR(1, “m",“in")^2
Convierte 6.3 libras a newtons y se multiplica por la
conversión de metros a pulgadas.
Excel 270
Transferencia de datos
Excel 271
Transferencia de datos - Lectura
• Algunas aplicaciones requieren que sean leídos o importados
ficheros diferentes de Excel.
• Para leer ficheros tipo texto se siguen los siguientes pasos:
– Asegurarse que el fichero es un fichero texto (extensión
típica .txt, .csv, o .prn).
– En Excel seleccionar Archivo→Abrir. Cuando aparece la
ventana de diálogo seleccionar Archivos de texto.
Seleccionar el archivo.
– Aparece el Asistente. Es necesario seleccionar si el
fichero tiene delimitadores entre campos o si son de
ancho fijo.
– Si hay delimitadores, seleccionar el tipo de separador.
– Finalmente se selecciona el formato.
Excel 272
Transferencia de datos - Lectura
Excel 273
Importación de datos desde
páginas Web
• Es posible importar datos desde una página
Web.
– La forma más fácil es utilizar Datos > Obtener
datos externos > Desde Web. Aparece un
navegador donde se puede colocar la URL
deseada. Ejemplo:
[Link]
Excel 274
Transferencia de datos - Lectura
Excel 275
Transferencia de datos - Escritura
Excel 276
Transferencia de datos - Escritura
Excel 277
Organización y Análisis de
datos – Tablas dinámicas
Excel 278
Organización de datos - Listas
Excel 280
Organización de datos - Ordenación
Excel 281
Organización de datos - Ordenación
Excel 282
Organización de datos - Filtrado
Excel 283
Organización de datos - Filtrado
• Ejercicios (Provincias_Españ[Link]):
– Las 10 provincias que tienen mayor densidad de
población.
– Qué provincias tienen superficies que exceden 15,000
km2.
– Qué provincias tienen poblaciones entre 500000 y 1 millón
de habitantes.
Excel 284
Organización de datos - Filtrado
Excel 285
Organización de datos - Filtrado
Excel 286
Organización de datos - Filtrado
Excel 287
Organización de datos - Filtrado
Excel 288
Búsqueda en tablas
Excel 289
Búsqueda en tablas - BUSCAR
Longitud dada = 20
Unidades = ft
Material = Acero
Factor de diseño = 15
Excel 291
Búsqueda en tablas - BUSCARV
Excel 292
Búsqueda en tablas - BUSCARV
1 2 3 4 5 6 7 8 9
Temperatura (C) Densidad (kg/m3)Energía Interna (kJ/kg) Entalpía (kJ/kg) Entropía (J/g*K) Cv (J/g*K) Cp (J/g*K) [Link] (m/s) Viscosidad (Pa*s)
10 999.7 42.018 42.119 0.15108 4.1906 4.1952 1447.3 0.0013059
20 998.21 83.906 84.007 0.29646 4.1567 4.1841 1482.3 0.0010016
30 995.65 125.72 125.82 0.43673 4.1172 4.1798 1509.2 0.00079735
40 992.22 167.51 167.62 0.57237 4.0734 4.1794 1528.9 0.00065298
50 988.04 209.32 209.42 0.70377 4.0262 4.1813 1542.6 0.00054685
60 983.2 251.15 251.25 0.83125 3.9765 4.185 1551 0.0004664
70 977.76 293.02 293.12 0.95509 3.9251 4.1901 1554.7 0.00040389
80 971.79 334.95 335.06 1.0755 3.8728 4.1968 1554.4 0.00035435
90 965.31 376.96 377.06 1.1928 3.8203 4.2052 1550.5 0.00031441
Excel 293
Búsqueda en tablas - BUSCARH
Excel 294
Búsqueda en tablas - BUSCARH
1 Temperatura (C) 10 20 30 40 50 60 70 80 90
2 Densidad (kg/m3) 999.7 998.21 995.65 992.22 988.04 983.2 977.76 971.79 965.31
3 Energía Interna (kJ/kg) 42.018 83.906 125.72 167.51 209.32 251.15 293.02 334.95 376.96
4 Entalpía (kJ/kg) 42.119 84.007 125.82 167.62 209.42 251.25 293.12 335.06 377.06
5 Entropía (J/g*K) 0.15108 0.29646 0.43673 0.57237 0.70377 0.83125 0.95509 1.0755 1.1928
6 Cv (J/g*K) 4.1906 4.1567 4.1172 4.0734 4.0262 3.9765 3.9251 3.8728 3.8203
7 Cp (J/g*K) 4.1952 4.1841 4.1798 4.1794 4.1813 4.185 4.1901 4.1968 4.2052
8 [Link] (m/s) 1447.3 1482.3 1509.2 1528.9 1542.6 1551 1554.7 1554.4 1550.5
9 Viscosidad (Pa*s) 0.001306 0.001002 0.000797 0.000653 0.000547 0.000466 0.000404 0.000354 0.000314
Excel 295
Búsqueda en tablas - COINCIDIR
Excel 296
Búsqueda en tablas - INDICE
Columna 2 = 0.9333076
0.0390457
0.4968178
0.6043112
Excel 297
Tablas pivot o dinámicas
Excel 298
Tablas pivot o dinámicas
• Ejemplo: (datos_pobl_usa.xls)
– Primero asegurarse que la lista está formada por bloques de
celdas contiguas con un encabezado en cada columna.
– Se selecciona cualquier celda dentro de la lista y seleccionar
en el menú Insertar – Tablas la opción Tabla dinámica.
– Aparece el Asistente para crear tabla dinámica con el rango de
datos incluyendo los encabezados. Seleccionar Nueva hoja de
cálculo. Pulsar el botón Aceptar.
– Aparece una hoja de trabajo que incluye una ventana de
Campos de tabla dinámica.
– Seleccionar los campos según las filas y columnas que se
requieran. También los campos que se quiere mostrar en
valores y los filtros.
Excel 299
Tablas pivot o dinámicas
Excel 300
Tablas pivot o dinámicas
Excel 301
Tablas pivot o dinámicas
– Herramientas – Gráfico dinámico
Excel 302
Series
Excel 303
Series de números
sn = sn −1 + xn
• Evaluación de funciones y expansión de Taylor:
– Del Cálculo se sabe que cualquier función que tiene n+1
derivadas en un punto a tiene una expansion polinómica
nth de Taylor Polynomial centrada en a y un error.
Excel 305
Series de números – Constantes Array
• Ejemplo: [Link]
• Se puede usar en Excel constantes Array para crear
fórmulas de series.
– Una constante array es un array de valores separados por
comas y encerrados entre llaves, usado como argumento
de una función. Ejemplo de array literal: {40,21,300,10}
– Se puede usar una constante array para hacer la
evaluación de una fórmula de serie más compacta y
precisa. Por ejemplo para evaluar:
∞
1
e = 1+ ∑ = 1 +SUMA( 1 /FACT({1,2,3,4,5,6,7,8,9,10}))
k =1 k!
Excel 306
Series de números – Función FILA
• Ejemplo: [Link]
• Se puede usar la función Excel FILA para generar
series de números.
– Si se introduce en una celda =FILA(1:100), se selecciona
y se pulsa la tecla F9 se obtiene:
={1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21;22;23;24;25;26;27;28;29;30;31;
32;33;34;35;36;37;38;39;40;41;42;43;44;45;46;47;48;49;50;51;52;53;54;55;56;57;58;59;
60;61;62;63;64;65;66;67;68;69;70;71;72;73;74;75;76;77;78;79;80;81;82;83;84;85;86;87;
88;89;90;91;92;93;94;95;96;97;98;99;100}
• Ejemplo: [Link]
• Se puede usar la función Excel INDIRECTO para crar
una referencia especificada por una cadena de texto.
– Si se introduce en una celda =INDIRECTO(“A1) crea una
referencia a la celda A1 y devuelve el valor contenida en
esa celda.
– Se puede usar este método junto con FILA para evaluar
fórmulas de series. Por ejemplo para calcular e:
{=1+SUMA(1/FACT(FILA(INDIRECTO("1:20"))))} o
{=1+SUMA(1/FACT(FILA(INDIRECTO("1:"&A1))))} donde el valor
en A1 especifica el número de términos a evaluar.
Excel 308
Series de Taylor
1! 2! n!
f (n +1) (c )
Rn ( x ) = (x − a )n+1
(n + 1)!
– El valor f (k)(a) es la kth derivada evaluada en a. La función
Rn(x) representa el error donde c es un valor entre x y a.
Excel 309
Series de Taylor
• Ejemplo: [Link]
Excel 310
Series de Taylor - Ejemplos
1 − x k =0 (1 − c )n+2
∞ k 2 3 n c
x x x x x e
ex = ∑ = 1+ + + + + + x n +1
k = 0 k! 1! 2! 3! n! (n + 1)!
sin x = ∑
n
(− 1) x 2 k +1
k
x3 x5
= x − + ± ±
x 2 n +1
+
± sin c 2 n + 2
x
k = 0 (2 k + 1)! 3! 5! (2n + 1)! (2n + 2)!
cos x = ∑
n
( − 1) x 2 k
k
x2 x4
= 1− + ± ±
x 2n
+
± cos c 2 n + 2
x
k =0 (2k )! 2! 4! (2n )! (2n + 1)!
Excel 311
Series de Taylor - Ejemplos
– Serie de∞ Taylor para la función f(x) = arctan(x).
1
= ∑ x k = 1 + x + x 2 + x3 +
1 − x k =0
∞
= ∑ (− x ) = 1 − x + x 2 − x 3 ±
1 k
1 + x k =0
∞
1
= ∑ (− 1) x = 1 − x2 + x4 − x6 ±
k 2k
1+ x 2
k =0
arctan x = ∑
∞
(− 1)k x 2 k +1 = x − x 3 + x 5 − x 7 ±
2k + 1 3 5 7
( )
k =0
– Serie de Taylor para la función g ( x ) = x 3 cosh x
( )
∞ ∞
x 2k x2 x4 x6 xk x x 2 x3
cosh x = ∑ = 1+ + + + cosh x =∑ = 1+ + + +
k =0 (2 k )! 2! 4! 6! k =0 (2 k )! 2! 4! 6!
( ) x k +3
∞
x 4 x5 x6
3
x cosh x =∑ = x + + + +
3
k =0 (2 k )! 2! 4! 6!
Excel 312
Series de Taylor - Observaciones
Excel 314
Interpolación
Excel 317
Interpolación lineal mediante VBA
Excel 321
Interpolación numérica
Excel 322
VBA (Visual Basic for
Applications) en Excel
Programación en Excel
Excel 323
Introducción a VBA
Excel 324
El editor de Visual Basic
Excel 325
Ventanas del editor de Visual Basic
Barra de Menús
Barra de Herramientas
Explorador de
Proyectos: diagrama
de árbol que contiene
cada hoja de trabajo.
Para abrir Control+R
Ventana código. Cada elemento
de un proyecto tiene asociada
una ventana de código.
Ventana de Propiedades
Ventana inmediato. Es útil para ejecutar
instrucciones de VBA directamente. Para
abrirla Control+G.
Excel 326
Gestión de módulos en VBA
• En la ventana del Explorador de proyectos se gestionan los
módulos.
• Los módulos pueden ser de cuatro tipos:
• Procedimientos Sub. Conjuntos de instrucciones que ejecutan
alguna acción.
• Procedimientos Function. Es un conjunto de instrucciones que
devuelven un solo valor.
• Procedimientos Property. Son procedimientos especiales que se
usan en módulos de clase.
• Declaraciones. Es información acerca de una variable que se le
proporciona a VBA.
• Un solo módulo de VBA puede guardar cualquier cantidad de
procedimientos Sub, procedimientos Function y
declaraciones.
Excel 327
Objetos
• Excel incluye cerca de 200 objetos, que representan rangos
de celdas, gráficos, hojas de cálculo, libros y la propia
aplicación de Excel.
• Cada objeto tiene propiedades (que permiten acceder y
controlar sus atributos) y métodos (funcionalidades).
• El examinador de objetos es una herramienta que permite
navegar por los objetos para explorar sus propiedades y
métodos.
• Para abrir el examinador de objetos en VBA pulsar F2 o
seleccionar: Ver → Examinador de Objetos
Excel 328
Objetos
Excel 329
Examinador de Objetos
Excel 330
Aplicación Excel
• Excel es una aplicación con un modelo de tres niveles:
• El primer nivel es el de servicios de cliente, que es la interfaz que
permite a los usuarios manejar la aplicación.
• El segundo nivel es el modelo de objetos de Excel, que es el que
se utiliza para realizar las operaciones en el libro de cálculo
(Workbook) o en las hojas de cálculo (Worksheets). Cada comando
de Excel se puede manejar mediante el modelo de objetos.
• El tercer nivel es el de servicios de datos, que es el que mantiene
los datos en las hojas de cálculo que son modificados por los
comandos del modelo de objetos de Excel.
Excel 331
Modelos de Objetos
• El modelo de objetos de Excel contiene una gran cantidad de
elementos ordenados en forma jerárquica. Algunos son:
• Application: Es el objeto que se encuentra en la base de la
jerarquía del modelo de objetos de Excel y representa a la
aplicación en sí.
• Workbooks: Objetos que representan los libros de cálculo o
archivos de Excel. Se encuentra debajo del objeto application
en la jerarquía.
• Worksheets: Objetos que representan las hojas de cálculo de
Excel. Este objeto pertenece al objeto workbook.
• Ranges: Objeto que representa un rango de celdas. Este
objeto pertenece al objeto worksheet.
• Charts: Objetos que representan gráficos.
• Pivot Tables: Objetos que representan tablas dinámicas.
Excel 332
Objeto Application
• El objeto Application representa el programa Excel. Entrega
acceso a las opciones y otras funcionalidades de Excel.
• La propiedad ActiveSheet se refiere a la hoja de cálculo
activa. Ejemplo:
[Link](1, 2) = time
• Le dice a Excel que coloque el valor de time en la celda que
está en la fila 1 y columna 2.
• La propiedad ScreenUpdating le indica a Excel si debe
refrescar la pantalla cuando se ejecuta código.
[Link] = False
Excel 333
Objeto Workbook
• El objeto Workbook representa un archivo Excel.
• El objeto ActiveWorkbook pertenece al objeto Application, y
entrega el objeto Workbook activo. Ejemplo:
[Link]
• El objeto ActiveSheet pertenece al objeto Workbook y se
refiere a la hoja de cálculo activa.
[Link]
• La propiedad Names entrega la lista de nombres que se han
definido en ese Workbook.
• La propiedad Path se refiere al directorio donde se encuentra
el Workbook. Ejemplo:
directorio = [Link]
Excel 334
Colección Workbook
• La colección Workbooks agrupa a todos los archivos de
Excel que se encuentran abiertos.
• El método Open, Save y SaveAs le indican a Excel si debe
abrir, guardar o guardar como el workbook correspondiente.
Ejemplos:
[Link](“ClaseIndustrial”).Save
[Link](“C:\[Link]”)
Workbooks(“Libro1”).SaveAs(“[Link]”,,”clavesecreta”)
• Se pueden entregar los parámetros por nombre a los
métodos. Ejemplos:
[Link] FileName :=“C:\[Link]”, _
ReadOnly:=True, Password:=“clavesecreta”
[Link](“ClaseIndustrial”).Save
Excel 335
Objeto Worksheet
• El objeto Worksheet representa una hoja de cálculo Excel. El
objeto ActiveSheet es un subobjeto del objeto Workbook que
entrega el Worksheet activo.
• Se puede copiar, pegar, imprimir, guardar, activar y borrar la
hoja de cálculo. Ejemplo:
With [Link](“ClaseIndustrial”)
[Link]
[Link]
[Link]
[Link]
[Link]
[Link]
End With
Excel 336
Colección Worksheet
• La colección Worksheets contiene a todas las hojas de cálculo que
pertenecen a algún workbook.
• Se le puede dar un nombre a un worksheet en particular para referirse a él.
Ejemplo:
Dim w As Workbook, s As Worksheet
Set w = Workbooks(“Libro1”)
Set s = [Link](“Hoja1”)
MsgBox [Link](“a1”).Value
• Se pueden nombrar todas las hojas de un archivo usando el comando For
Each … Next Loop.
Sub MuestraNombres()
Dim w As Worksheet
For Each w In Worksheets
MsgBox [Link]
Next
End Sub
Excel 337
Objeto WorksheetFunction
Excel 339
Objeto Range
• El objeto Range representa rangos de celdas. También es posible
acceder a las celdas usando la propiedad Cells de ActiveSheet.
• Ejemplos:
Set notas = Worksheets(“Funciones”).Range(“F2:F13”)
prom = [Link](notas)
Worksheets(“Funciones”).Range(“F14”).Value = prom
Worksheets(“Funciones”).Range(“F15”).Formula = “=average(F2:F13)”
Worksheets(“Funciones”).Cells(2, 1).Select
Workbooks(“Libro1”).Worksheets(“Hoja1”).Range(“A1).Value = 10
Workbooks(“Libro1”).Worksheets(“Hoja1”).Range(“A2.A10”).Value = 5
Workbooks(“Libro1”).Worksheets(“Hoja1”).Range(“A2:A10”).Value = 5
Workbooks(“Libro1”).Worksheets(“Hoja1”).Range(“A2”, ”A10”).Value = 5
Excel 340
Objeto Range
Excel 341
Módulos VBA
• Un módulo VBA se compone de procedimientos que son
códigos de ordenador que realizan alguna acción sobre los
objetos o con ellos.
Sub Prueba()
Sum= 1+1
MSGBox “La respuesta es” & Sum
End Sub
Excel 342
Introducir código VBA
Sub Hola()
Msg = “Su nombre es “ & [Link] & “?”
Ans = MsgBox(Msg, vbYesNo)
If Ans = VbNo Then
MsgBox “No se preocupe”
Else
MsgBox “Debo ser adivino!”
End If
End Sub
Excel 343
Ejecutar código VBA
Excel 344
Subrutinas
Excel 346
Subrutinas
• Las subrutinas se pueden llamar desde otras partes del
código usando su nombre y agregando los parámetros que
necesita.
• Para llamar a una subrutina llamada MiSub se puede usar:
MiSub 4, 2.87
Call MiSub(4, 2.87)
• También se puede agregar el nombre de la subrutina a
botones u otros controles de VBA.
Excel 347
Funciones
• Las funciones son similares a las subrutinas con la diferencia
que se usa Function en vez de Sub y que retornan un valor
después de ejecutarse.
Public Function Calc_q(y1 As Double, y3 As Double) As Double
Calc_q = 1 / ((Abs(y3 ‐ y1)) ^ 0.74)
End Function
• Se pueden usar como cualquier función de Excel.
Public Function MiFactorial(N As Integer) As Integer
‘Funcion que calcula el factorial de un numero N
MiFactorial = 1
For i% = 1 To N
MiFactorial = i * MiFactorial
Next
End Function
Excel 348
Conceptos Básicos del Lenguaje
• Para comentar el código se usa ‘ o Rem
‘Declaración de variables
Dim y As Double
Dim x As Double
Rem Declaración de Matrices
Dim M(1 To 8, 1 To 8) As Double
Dim N(8, 8) As Double
• Para separar múltiples líneas se usa un guión bajo (_):
K2(1) = dt * dy1dt(y(1) + k1(1) / 2#, y(2) + _
k1(2) / 2#, y(3) + k1(3) / 2#, y(4) + _
k1(4) /2#)
Tiene que haber un espacio antes del underscore.
Excel 349
Variables y Tipos de Datos
• Los datos manipulados en VBA residen en objetos (p.e. rangos de
hojas de cálculo) o en variables que se crean.
• Una variable es una localización de almacenamiento con nombre,
dentro de la memoria del ordenador. Debe tener asociado un tipo de
dato.
• Las reglas para nombrar las variables son:
• Se pueden usar caracteres alfabéticos, números y algún carácter
de puntuación, pero el primero de los caracteres debe ser
alfabético
• VBA no distingue entre mayúsculas y minúsculas
• No se pueden usar espacios ni puntos
• No se pueden incrustar en el nombre de una variable los
siguientes símbolos: #, $, %, !
• La longitud del nombre puede tener hasta 254 caracteres
Excel 350
Tipos de Datos en VBA
Tipo de dato Bytes Rango de valores
Byte 1 0 a 255
Boolean 2 True o False
Integer 2 -32768 a 32767
Long 4 - 2147483648 y 2147483647
Currency 8 -922337203685477.5808 a 922337203685477.5807
Single 4 -3.402823E38 a 3.402823E38
Double 8 -1.79769313486231E308 a 1.79769313486232E308
Date 8 1-1-100 al 31-12-9999 y horarios de 0:00:00 a 23:59:59
String longitud variable (2^31 caracteres). longitud fija (2^16)
Object 4
Variant cualquier clase de datos excepto cadena de longitud fija
Excel 351
Definición de Variables
• Con Dim o Public se declaran las variables:
Dim b As Double, a As Double
Dim n, m As Integer
Dim InerestRate As Single
Dim TodaysDate As Date
Dim UserName As String * 20
Dim x As Integer, y As Integer, z As Integer
Si una variable no se declara se asume de tipo
Variant (tipo genérico).
• En general debe ser:
Dim NombreVariable As DataType
Excel 352
Ámbito de las variables
• El ámbito de una variable determina el módulo y el procedimiento
en el que se puede usar una variable.
Cómo se declara una variable en este
Ámbito
ámbito
Excel 354
Arrays
• Por defecto los subíndices de los arrays de VBA empiezan en 0. Si
deseamos que comience en 1 en vez de en 0, incluiremos antes del
primer array y antes del primer procedimiento la expresión:
Option Base 1 o explícitamente el rango de elementos
• Para acceder a los elementos del array:
y(3) = 2.983
M(1, 2) = 4.321
MiArray(1) = 20
MiMatriz(1,2) = 20
• Si no se sabe el tamaño, se puede usar ReDim:
Dim Matriz() As Double
ReDim Matriz(10)
ReDim Preserve Matriz(12) ‘Mantiene lo que estaba
Excel 355
Definición de Constantes
• Con Const se declaran las constantes:
Const MiConstante As Integer = 14
Const MiConstante2 As Double = 1.025
Const NumTrim As Integer = 4
Const Interés = 0.05, Periodo = 12
Const Nombre Mod as String = “Macros Presupuestos”
Public Const NombreApp As String = “Aplicación Presupuestos”
• Las constantes también poseen un ámbito:
– Si se declaran después de Sub o Function es local.
– Si se declara al inicio de un módulo está disponible para
todo el módulo.
– Si se declara con Public al inicio de un módulo está
disponible para todos los módulos de una hoja de trabajo.
Excel 356
Constantes y Cadenas
• Constantes predeterminadas, que se pueden usar sin
necesidad de declararlas.
Sub CalcManual()
[Link] = xlManual
End Sub
• Cadenas, hay dos tipos de cadenas en VBA:
• De longitud fija, que se declaran con un número específico
de caracteres. La máxima longitud es de 65.536
caracteres.
• De longitud variable, que teóricamente pueden tener hasta
2.000 millones de caracteres.
Dim MiCadena As String * 50
Dim SuCadena As String
Excel 357
Fechas y Expresiones
• Trabajar con Fechas
Dim Hoy As Date
Dim HoraInicio As Date
Const PrimerDía As Date = #1/1/2001#
Const MedioDía As date = #12:00:00#
• Expresiones de asignación, expresión que realiza evaluaciones
matemáticas y asigna el resultado a una variable o a un objeto.
Se usa el signo igual “=“ como operador de asignación.
x=1
x=x+1
x = (y * 2) / (z * 2)
FileOpen = true
Range(“Año”). Value = 1995
Excel 358
Operadores
• OPERADORES ARITMÉTICOS
+ Suma, - Resta, * Multiplicación, / División, \ División entera,
Mod Resto, ^ exponencial, & Concatenación
• OPERADORES COMPARATIVOS
= Igual, < Menor, > Mayor, <= Menor o igual, >= Mayor o igual,
<> Distinto
• OPERADORES LÓGICOS
Not (negación lógica, And (conjunción lógica), Or (disyunción
lógica), XoR (exclusión lógica), Eqv (equivalencia en dos
expresiones), Imp (implicación lógica)
Excel 359
Estructuras WITH...END WITH
• VBA ofrece dos estructuras que simplifican el trabajo con objetos y
colecciones.
• Con WITH...END WITH se permite realizar múltiples operaciones en un
solo objeto.
Sub CambiarFuente()
With [Link]
.Name = “Times New Roman”
.FontStyle = “Bold Italic”
.Size = 12
.Underline = xlSingle
.ColorIndex = 5
End With
End Sub
Excel 360
Estructuras FOR EACH...NEXT
• Para una colección no es necesario saber la cantidad de elementos
que existen en ella para usar la estructura For Each...Next.
Sub ContarHojas()
‘Muestra el nombres de las hojas del libro de trabajo activo
Dim Item As Worksheet
For Each Item In [Link]
MsgBox [Link]
Next Item
End Sub
Sub VentanasAbiertas()
‘Cuenta el número de ventanas abiertas
Suma = 0
For Each Item In Windows
Suma = Suma + 1
Next Item
MsgBox “Total de ventanas abiertas”, & Suma
End Sub
Excel 361
Condicionales
• Los tests lógicos en VBA tienen la siguiente sintaxis:
If (time = 32000) Then
MsgBox “time vale 32000”
End If
If (MiCondicion = True) Then
MsgBox “Mi Condición es Verdad”
Else
MsgBox “Mi Condición No es Verdad”
End If
If (contador < 10) Then
MsgBox “El Contador es menor a 10”
ElseIf (contador < 20) Then
MsgBox “El Contador es mayor que 10 y menor que 20”
ElseIf (contador < 30) Then
MsgBox “El Contador es mayor que 20 y menor que 30”
End If
Excel 362
Estructuras Select Case
• La estructura Select Case es útil para elegir entre tres o más
opciones
Sub Positivos_Negativos_Cero()
a = InputBox("Ingrese un número")
Select Case a
Case Is > 0
Msg = "Número Positivo"
Case Is < 0
Msg = "Número negativo"
Case Else
Msg = "Cero"
End Select
MsgBox Msg
End Sub
Excel 363
Bucles For…Next
• Esta sentencia de iteración se ejecuta un número determinado de
veces. Su sintaxis es:
For contador = empezar To finalizar [Step valorincremento]
[Instrucciones]
[Exit For]
[instrucciones]
Next [contador]
______________________________________________________________________________
Sub SumaNúmeros
Sum = 0
For Count = 0 To 10
Sum = Sum + Count
Next Count
MsgBox Sum
End Sub
Excel 364
Bucles For…Next
For i = 1 To n
‘Código
Next i
_________________________________
For i = 1 To n Step 2
‘Código
Next i
_________________________________
For i = 1 To n
‘Código
If tiempo >10 Then
Exit For
End If
‘Más Código
Next i
Excel 365
Bucles Do…While, Do…Until
• El bucle se ejecuta hasta que la condición llegue a ser verdadera. Do
Until tiene la sintaxis.
Do Until [condicion]
[instrucciones]
[Exit Do]
[instrucciones]
Loop
_______________________________________________________________________________________________________
Sub DoUntilDemo()
Do
[Link] = 0
[Link](1, 0).Select
Loop Until Not IsEmpty(ActiveCell)
End Sub
Excel 366
Bucles Do
Do While (tiempo < 10)
‘Código
Loop
______________________________________________________________________________________
Do
‘Código
Loop While (tiempo < 10)
_______________________________________________________________________
Do
‘Más Código
Loop Until (tiempo > 10)
Excel 367
Funciones para cálculos con vectores
x = [Link](1).Value
y = [Link](2).Value
z = [Link](3).Value
v_Mag = Sqr(x ^ 2 + y ^ 2 + z ^ 2)
End Function
Excel 369
Funciones para cálculos con vectores
Excel 371
Código para la función v_CrossProduct
vx = [Link](1).Value
vy = [Link](2).Value
vz = [Link](3).Value
Excel 372
Depuración
Excel 373
Depuración
Excel 374
Formularios
Excel 376
Formularios
Excel 378
Formularios MODIFICAR!!!
Excel 380
Controles - Diseño
Excel 381
Tipos de Controles
Excel 382
Macros
Excel 384
Macros - Diseño
Excel 385
Escribir Macros
• El Editor de Visual Basic es una herramienta para
escribir y modificar código escrito en VBA
• Para abrir el Editor de Visual Basic: En el menú
Herramientas → Macro → Editor de Visual Basic o
Alt+F11.
• Las macros se almacenan en módulos de un libro de
trabajo.
• Los módulos se agregan en el Editor de Visual Basic
seleccionando Módulo en el menú Insertar del editor.
• Debe aparecer una ventana de módulo vacía dentro
de la ventana principal del Editor de Visual Basic.
Excel 386
Macros - Editor VB
Excel 387
Asignar nombre a la Macro
Excel 388
Asignar código a la Macro
• Si se desea mostrar un mensaje simple escribir MsgBox “Mi
primera macro”.
• MsgBox es la palabra que VBA utiliza para los cuadros de
mensaje.
• Si se ejecuta la macro, Excel mostraría un mensaje con el texto
Mi primera macro y un botón Aceptar para cerrar el mensaje.
Excel 389
Macros de Bucle
• Las macros de bucle funcionan recorriendo los datos de
celdas para realizar acciones automáticamente de
manera repetida.
• Hay varias instrucciones que permiten crear este tipo de
macros:
– For Each…Next
– For ... Next
– For ... Next Loop With Step
– Do While ... Loop
– Do Until ... Loop
– Do ... Loop While
– Do ... Loop Until
Excel 390
Macro de Bucle For Each…Next
Excel 391
Propiedad Cells y Range
Excel 392
Ejemplos de Macros
Excel 393
Ejemplos de Macros
Libro Libro
Pelicula Pelicula
Revista Revista
Lee Libro Lee Libro
Ver pelicula Ver pelicula
Vino Vino
Texto Texto
Libro texto Libro texto
Excel 394
Ejemplos de Macros
Excel 395
Evaluación de derivadas
Excel 396
Diferenciación de funciones continuas
Diferenciación 397
Derivadas a partir de datos
x x+Δx
Graphical Representation of forward difference approximation of first derivative.
Diferenciación 399
Aproximación por diferencia en atraso
lim f ( x + Δx ) − f ( x )
• Sabemos que f ′( x ) =
Δx → 0 Δx
f (x + ∆x ) − f ( x )
• Para un finite ' Δx' , f ′( x ) ≈
∆x
Diferenciación 400
Aproximación por diferencia en atraso
where
x
x-Δx x
Diferenciación 401
Obtención de la adad a partir de las series
de Taylor
• Taylor’s theorem says that if you know the value of a
function f at a point xi and all its derivatives at that
point, provided the derivatives are continuous
between xi and xi +1 , then
f ′′( xi )
f ( xi +1 ) = f ( xi ) + f ′( xi )( xi +1 − xi ) + (xi +1 − xi )2 +
2!
Substituting for convenience Δx = xi +1 − xi
f ′′( xi )
f (xi +1 ) = f (xi ) + f ′(xi )Δx + (Δx )2 +
2!
f ( xi +1 ) − f ( xi ) f ′′(xi )
f ( xi ) =
′ − (∆x ) +
∆x 2!
f ( xi +1 ) − f ( xi )
f ′( xi ) = + O(∆x )
∆x
Diferenciación 402
Obtención de la adad a partir de las series
de Taylor
• The O(∆x ) term shows that the error in the approxima-
tion is of the order of ∆x . It is easy to derive from
Taylor series the formula for backward divided
difference approximation of the first derivative.
• As shown above, both forward and backward divided
difference approximation of the first derivative are
accurate on the order of O(∆x ) .
• Can we get better approximations? Yes, another
method is called the Central difference approxima-
tion of the first derivative.
Diferenciación 403
Obtención de la adc a partir de las series
de Taylor
• From Taylor series
f ′′( xi ) f ′′′(xi )
f ( xi +1 ) = f ( xi ) + f ′( xi )Δx + (Δx ) +
2
(Δx )3 + (1)
2! 3!
f ′′( xi ) ′′′( )
f ( xi −1 ) = f ( xi ) − f ′( xi )Δx + (Δx )2 − f xi (Δx )3 + (2)
2! 3!
f (xi +1 ) − f (xi −1 )
f ′( xi ) = + O(∆x )
2
2∆x
Diferenciación 404
Obtención de la adc a partir de las series
de Taylor
• Hence showing that we have obtained a more
accurate formula as the error is of the order of O(∆x )2
f(x)
x
x-Δx x x+Δx
Diferenciación 405
Fórmula de 5 puntos
• Ejemplo: [Link]
Diferenciación 406
Fórmulas de diferencia finita en adelanto
Diferenciación 407
Fórmulas de diferencia finita en atraso
Diferenciación 408
Fórmulas de diferencia finita centrada
Diferenciación 409
Ecuaciones diferenciales
ordinarias
Excel 410
Ecuaciones diferenciales de primer orden
Excel 411
Ecuaciones diferenciales de primer orden
y valor inicial
• Problema: se requiere hallar la solución de una
ecuación diferencial de primer orden de la forma:
𝑑𝑑𝑑𝑑
= 𝑓𝑓 𝑥𝑥, 𝑦𝑦
𝑑𝑑𝑑𝑑 y = e^x - x - 1
𝑦𝑦(0) = 𝑎𝑎 0.8
0.7
𝑑𝑑𝑑𝑑
• Ejemplo: = 𝑥𝑥 + 𝑦𝑦
0.6
𝑑𝑑𝑑𝑑 0.5
0.4
𝑦𝑦 0 = 0 0.3
0.2
Solución: 𝑦𝑦 = 𝑒𝑒 𝑥𝑥 − 𝑥𝑥 − 1 0.1
0
0 0.2 0.4 0.6 0.8 1
Excel 413
Método de Euler
Excel 414
Método de Euler
Excel 415
Método de Euler
• Código:
Public Sub DoEuler1stOrder() For i = 1 To n
Dim yn As Double yn1 = yn + (xn + yn) * dx
Dim yn1 As Double xn = xn + dx
Dim xn As Double yn = yn1
Dim dx As Double [Link](i + 1, 1) = xn
Dim n As Integer [Link](i + 1, 2) = yn
yn = 0 Next i
xn = 0 End Sub
dx = 0.001
n = 1000
Excel 419
Método Runge-Kutta aplicado a
problemas de valor inicial de 2do orden
• Considerar la siguiente ecuación y condiciones
iniciales: 2
𝑑𝑑 𝑠𝑠 𝑑𝑑𝑠𝑠
𝑚𝑚 2 + 𝐶𝐶𝑑𝑑 = 𝑇𝑇
𝑑𝑑𝑡𝑡 𝑑𝑑𝑡𝑡
𝑠𝑠 0 = 0
𝑑𝑑𝑠𝑠
0 =0
𝑑𝑑𝑡𝑡
• Físicamente representa la ecuación del movimiento
de un objeto sujeto a un empuje T. m es la masa, Cd
un factor de rozamiento y s la posición del objeto.
Excel 420
Método Runge-Kutta aplicado a
problemas de valor inicial de 2do orden
• Para resolver la ecuación de movimiento se reescribe
para obtener dos ecuaciones de primer orden:
𝑑𝑑𝑑𝑑
si hacemos: 𝑣𝑣 =
𝑑𝑑𝑑𝑑
𝑑𝑑𝑣𝑣
𝑚𝑚 = 𝑇𝑇 − 𝐶𝐶𝑑𝑑 𝑣𝑣
𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑
= 𝑣𝑣
𝑑𝑑𝑑𝑑
𝑠𝑠𝑡𝑡=0 = 0
𝑣𝑣𝑡𝑡=0 = 0
• Se obtiene dos ecuaciones de primer orden
acopladas, a las que se aplican técnicas numéricas.
Excel 421
Método Runge-Kutta aplicado a
problemas de valor inicial de 2do orden
• El método de Runge Kutta se basa en tomar más
términos de la serie de Taylor de la función, que se
traduce en expandir más series de Taylor para
estimar las derivadas de mayor orden.
• El enfoque RK reduce el error de truncamiento del
orden de (dt)5 en oposición a (dt)2 del método de
Euler, con lo que se puede aumentar el paso
manteniendo la precisión.
• El compromiso es que hay que realizer más cálculos
en cada paso.
Excel 422
Método Runge-Kutta aplicado a
problemas de valor inicial de 2do orden
• Las ecuaciones generales de Runge Kutta para la
integración son:
𝑘𝑘1 = 𝑦𝑦 ′ (𝑥𝑥, 𝑦𝑦)(∆𝑥𝑥)
′
∆𝑥𝑥 𝑘𝑘1
𝑘𝑘2 = 𝑦𝑦 (𝑥𝑥 + , 𝑦𝑦 + )(∆𝑥𝑥)
2 2
′
∆𝑥𝑥 𝑘𝑘2
𝑘𝑘3 = 𝑦𝑦 (𝑥𝑥 + , 𝑦𝑦 + )(∆𝑥𝑥)
2 2
𝑘𝑘4 = 𝑦𝑦 ′ (𝑥𝑥 + ∆𝑥𝑥, 𝑦𝑦 + 𝑘𝑘3 )(∆𝑥𝑥)
(𝑘𝑘1 +2𝑘𝑘2 +2𝑘𝑘3 +𝑘𝑘4 )
𝑦𝑦 𝑥𝑥 + ∆𝑥𝑥 = 𝑦𝑦 𝑥𝑥 +
6
donde: 𝑦𝑦 ′ representa 𝑑𝑑𝑑𝑑⁄𝑑𝑑𝑑𝑑
Excel 423
Método Runge-Kutta aplicado a
problemas de valor inicial de 2do orden
For i = 1 To n ' Start iterations
Public Sub DoRK2ndOrder() F = (t - (Cd * Vn)) ' Compute k1
Dim t, Cd, M, dt As Double ' Thrust, Drag coefficient, Mass A =F/M
Dim dt, F, A As Double ' Time step size, Force, Acceleration k1 = dt * A
Dim Vn As Double ' Velocity at time t F = (t - (Cd * (Vn + k1 / 2))) ' Compute k2
Dim Vn1 As Double ' Velocity at time t + dt A =F/M
Dim Sn As Double ' Displacement at time t k2 = dt * A
Dim Sn1 As Double ' Displacement at time t + dt F = (t - (Cd * (Vn + k2 / 2))) ' Compute k3
Dim time As Double ' Total time A=F/M
Dim k1, k2, k3, k4 As Double ' RK k1, RK k2, RK k3, RK k4 k3 = dt * A
Dim n As Integer ' Counter controlling total number of time steps F = (t - (Cd * (Vn + k3))) ' Compute k4
Dim C As Integer ' Counter controlling output of results to spreadsheet A=F/M
Dim k As Integer ' Counter controlling output row k4 = dt * A
Dim r As Integer ' Number of output rows Vn1 = Vn + (k1 + 2 * k2 + 2 * k3 + k4) / 6 ' Compute velocity at t + dt
With ActiveSheet ' Extract given data from the active spreadsheet: Sn1 = Sn + Vn1 * dt ' Compute displacement at t + dt using Euler
dt = .Range("dt") time = time + dt ' Update variables
t = .Range("T") Vn = Vn1
M = .Range("M") Sn = Sn1
Cd = .Range("Cd") If C >= n / r Then ' Output results to the active spreadsheet
n = .Range("n") [Link](k + 1, 1) = time
r = .Range("r_") [Link](k + 1, 2) = Sn
End With [Link](k + 1, 3) = Vn
k=1 ' Initialize variables k=k+1
time = 0 C=0
C=n/r Else
Vn = 0 C=C+1
Sn = 0 End If
Next i
End Sub
Excel 424
Ecuaciones diferenciales con condiciones
de contorno o de frontera
• Hay problemas que se modelizan mediante una
ecuación diferencial de segundo orden con
condiciones en sus dos extremos [a, b], que se
denomina como ecuación diferencial ordinaria con
valores en la frontera o contorno. Se formula como:
𝑦𝑦 ′′ = 𝑓𝑓 𝑡𝑡, 𝑦𝑦, 𝑦𝑦 ′ , 𝑦𝑦 𝑎𝑎 = 𝛼𝛼, 𝑦𝑦 𝑏𝑏 = 𝛽𝛽
Excel 428
Sistemas de ecuaciones
lineales
[A]{x} = {b}
[A] = [x1 x3 ]
−1
x2
Sistemas de ec. lineales 434
Matriz inversa y sistemas estímulo -
respuesta
• Recall that LU factorization can be used to efficiently
evaluate a system for multiple right-hand-side vectors
- thus, it is ideal for evaluating the multiple unit vectors
needed to compute the inverse.
• Many systems can be modeled as a linear
combination of equations, and thus written as a matrix
equation:
[Interactions]{response} = {stimuli}
• The system response can thus be found using the
matrix inverse.
Sistemas de ec. lineales 435
Normas vectoriales y matriciales
= ∑ x i
p
X
i=1
p
i =1
p = ∞ : maximum − magnitude X ∞
= max xi
1≤i ≤ n
i=1 j=1
n
row - sum norm A ∞ = max ∑ aij
1≤i≤n
j=1
A 2 = (µ max )
1/2
spectral norm (2 norm)
a11
b2 − a21 x1j − a23 x3j−1
x =
j
2
a22
b3 − a31 x1j − a32 x2j
x =
j
3
a33
Sistemas de ec. lineales 440
Iteración de Jacobi
Excel 445
Ecuaciones en derivadas parciales
∂ 2u ∂ 2u ∂ 2u
A 2 +B +C 2 + D = 0
∂x ∂x∂y ∂y
• where A, B, and C are functions of x and y , and
.D is a function of
∂u ∂u
x, y, u and , .
∂x ∂y
• can be:
Elliptic if B2 – 4AC < 0
Parabolic if B2 – 4AC = 0
Hyperbolic if B2 – 4AC > 0
Ec. derivadas parciales 447
Ejemplos de EDPs de 2 orden
• Elliptic A = 1, B = 0, C = 1
∂ 2T ∂ 2T
+ 2 =0 Laplace equation
∂x 2
∂y
• Parabolic A = k , B = 0, C = 0
∂T ∂ T 2
=k 2 Heat equation
∂t ∂x
• Hyperbolic
1
A = 1, B = 0, C = −
c2
∂2 y 1 ∂2 y
= 2 2 Wave equation
∂x 2
c ∂t
Tt
W Tl Tr
x
Tb
L
x (i, j − 1)
(0,0)
Tb (m,0)
∂ 2T T ( x + ∆x, y ) − 2T ( x, y ) + T ( x − ∆x, y )
( x, y ) ≅
∂x 2
(∆x )2
∂ 2T T ( x, y + ∆y ) − 2T ( x, y ) + T ( x, y − ∆y )
( x , y ) ≅
∂y 2 (∆y )2
Ec. derivadas parciales 450
Discretizando la PDE elíptica
y
Tt
(0, n)
∆x (i, j + 1)
∆y ∆x
∆y Tr
Tl (i, j ) (i − 1, j ) (i, j ) (i + 1, j )
x (i, j − 1)
(0,0)
Tb (m,0)
∂ 2T T ( x, y + ∆y ) − 2T ( x, y ) + T ( x, y − ∆y ) ∂ 2T Ti , j +1 − 2Ti , j + Ti , j −1
( x , y ) ≅ ≅
∂y 2 (∆y )2 ∂y 2 i, j (∆y )2
Ec. derivadas parciales 451
Discretizando la PDE elíptica
∂ 2T ∂ 2T
+ =0
∂x 2
∂y 2
• if, ∆x = ∆y
• the Laplace equation can be rewritten as
Ti +1, j + Ti −1, j + Ti , j +1 + Ti , j −1 − 4Ti , j = 0 (Eq. 1)
• there are several numerical methods that can be used to solve the
problem:
Direct Method
Gauss-Seidel Method
Lieberman Method
Ec. derivadas parciales 452
Ejemplo 1: Método directo
• Consider a plate 2.4 m × 3.0 m that is subjected to the boundary
conditions shown below. Find the temperature at the interior nodes
using a square grid with a length of 0.6 m by using the direct method.
y
300 °C
75 °C
W 3.0 m 100 °C
x
50 °C
2.4 m
L
Ec. derivadas parciales 453
Ejemplo 1: Método directo
• We discretize the plate by taking, ∆x = ∆y = 0.6m
L W
m= = 4 n= =5
∆x ∆y
T0, 4
T1, 4 T2, 4 T3, 4 T4, 4
T0,3
T1,3 T2,3 T3,3 T4,3
T0, 2
T1, 2 T2, 2 T3, 2 T4, 2
T1,1 T1,2 T1,3 T1,4 T2,1 T2,2 T2,3 T2,4 T3,1 T3,2 T3,3 T3,4 RHE
-4 1 0 0 1 0 0 0 0 0 0 0 -125 T1,1 74.8719
1 -4 1 0 0 1 0 0 0 0 0 0 -75 T1,2 95.8959 300.0 300.0 300.0
0 1 -4 1 0 0 1 0 0 0 0 0 -75 T1,3 127.8036
75.0 196.9 206.3 185.0 100.0
0 0 1 -4 1 0 0 1 0 0 0 0 -375 T1,4 196.9288
1 0 0 0 -4 1 0 0 1 0 0 0 -50 T2,1 78.5917 75.0 127.8 143.4 133.5 100.0
0 1 0 0 1 -4 1 0 0 1 0 0 0 T2,2 105.9082 75.0 95.9 105.9 105.8 100.0
0 0 1 0 0 1 -4 1 0 0 1 0 0 T2,3 143.3896
75.0 74.9 78.6 83.6 100.0
0 0 0 1 0 0 1 -4 0 0 0 1 -300 T2,4 206.3200
0 0 0 0 1 0 0 0 -4 1 0 0 -150 T3,1 83.5868
50.0 50.0 50.0
0 0 0 0 0 1 0 0 1 -4 1 0 -100 T3,2 105.7554
0 0 0 0 0 0 1 0 0 1 -4 1 -100 T3,3 133.5267
0 0 0 0 0 0 0 1 0 0 1 -4 -400 T3,4 184.9617
75 °C
W 3.0 m 100 °C
x
50 °C
2.4 m
L
460
Ejemplo 2: Método Gauss-Seidel
• Discretizing the plate by taking, ∆x = ∆y = 0.6m
L W
m= = 4 n= =5
∆x ∆y
a i, j
Ti ,present
j
i=1, j=1 T1,1 = 42.9688 º C ε a 1,1 = 27.27% i=2, j=3 T2,3 = 56.4881 º C ε a 2,3 = 83.58%
i=1, j=2 T1, 2 = 38.7596 º C ε a 1, 2 = 31.49% i=2, j=4 T2, 4 = 156.150 º C ε a 2, 4 = 34.46%
i=1, j=3 T1,3 = 55.7862 º C ε a 1,3 = 54.49% i=3, j=1 T3,1 = 56.3477 º C ε a 3,1 = 24.44%
i=1, j=4 T1, 4 = 133.283 º C ε a 1, 4 = 24.90% i=3, j=2 T3, 2 = 56.0425 º C ε a 3, 2 = 31.70%
i=2, j=1 T2,1 = 36.8164 º C ε a 2,1 = 44.83% i=3, j=3 T3,3 = 86.8394 º C ε a 3,3 = 57.44%
i=2, j=2 T2, 2 = 30.8594 º C ε a = 62.03% i=3, j=4 T3, 4 = 160.747 º C ε a 3, 4 = 16.12%
2, 2
Residuals-squared
8
=(-4*D5+D4+D6+C5+E5)^2
7 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
6 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
5 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
4 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
3 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
2 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
1 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0
x\y 0 1 2 3 4 5 6 7 8
75 °C
3.0 m Insulated
x
50 °C
2.4 m
x
50 °C
2.4 m
• We can then substitute this into the original equation gives us,
∂T
2Tm −1, j + 2(∆x) + Tm , j −1 + Tm , j +1 − 4Tm , j = 0
∂x m, j
x ∆x ∆x
i −1 i i +1
Schematic diagram showing interior nodes
Ec. derivadas parciales 476
Solución EDP Parabólica:
Método explícito
• If we define ∆x =
L
we can then write the finite central divided
n
difference approximation of the left hand side at a general interior node
( i ) as ∂ T ≅ Ti +1 − 2Ti + Ti −1 where ( j ) is the node number along
2 j j j
∂x 2 (∆x )2
the time.
i, j
i −1 i i +1
Schematic diagram showing interior nodes
Ec. derivadas parciales 477
Solución EDP Parabólica:
Método explícito
• Substituting these approximations into the governing equation yields
Ti +j1 − 2Ti j + Ti −j1 Ti j +1 − Ti j
α =
(∆x )2 ∆t
• Solving for the temp at the time node j +1 gives
∆t
Ti j +1
= Ti + α
j
(∆x) 2
T j
i +1 −(2Ti
j
+ T j
i −1 )
• choosing, λ = α ∆t
2
(∆x)
• we can write the equation as, Ti = Ti + λ Ti +1 − 2Ti + Ti −1
j +1 j j j j
( )
• we can be solved explicitly: for each internal location node of the rod for
time node j + 1 in terms of the temperature at time node j. If we know
the temperature at node j = 0 , and the boundary temperatures, we can
find the temperature at the next time step. We continue the process until
we reach the time at which we are interested in finding the temperature.
Ec. derivadas parciales 478
Ejemplo 1 EDP Parabólica:
Método explícito
• Consider a steel rod that is subjected to a temperature of 100°C on the
left end and 25°C on the right end. If the rod is of length 0.05m ,use the
explicit method to find the temperature distribution in the rod from t = 0
and t = 9 seconds. Use ∆x = 0.01m , ∆t = 3s .
W kg
• Given: k = 54 , ρ = 7800 3 , C = 490
J
m−K m kg − K
• The initial temperature of the rod is 20°C.
i=0 1 2 3 4 5
T =100 ° C T = 25 °C
0.01m
• All internal nodes are at 20°C for t = 0 sec: Ti 0 = 20°C , for all i = 1,2,3,4
T00 = 100°C
T10 = 20°C We can now calculate the temperature at each node
T20 = 20°C
explicitly using the equation formulated earlier,
( )
Interior nodes
T30 = 20°C
Ti j +1 = Ti j + λ Ti +j1 − 2Ti j + Ti −j1
T40 = 20°C
T50 = 25°C
• Rearranging yields
− λTi −j1+1 + (1 + 2λ )Ti j +1 − λTi +j1+1 = Ti j
given that
∆t
λ =α
(∆x )2
• The rearranged equation can be written for every node during each
time step. These equations can then be solved as a simultaneous
system of linear equations to find the nodal temperatures at a particular
time.
Ec. derivadas parciales 483
Ejemplo 2 EDP Parabólica:
Método implícito
• Consider a steel rod that is subjected to a temperature of 100°C on the
left end and 25°C on the right end. If the rod is of length 0.05m ,use the
implicit method to find the temperature distribution in the rod from t = 0
and t = 9 seconds. Use ∆x = 0.01m , ∆t = 3s .
W kg
• Given: k = 54 , ρ = 7800 3 , C = 490
J
m−K m kg − K
• The initial temperature of the rod is 20°C.
i=0 1 2 3 4 5
T =100 ° C T = 25 °C
0.01m
• All internal nodes are at 20°C for t = 0 sec: Ti 0 = 20°C , for all i = 1,2,3,4
T00 = 100°C
T10 = 20°C
We can now form the system of equations for the
T20 = 20°C
first time step by writing the approximated heat
T30 = 20°C
Interior nodes conduction equation for each node
T40 = 20°C − λTi −j1+1 + (1 + 2λ )Ti j +1 − λTi +j1+1 = Ti j
T50 = 25°C
1
0 0 − 0 . 4239 1 . 8478 T4 30.598
• The above coefficient matrix is tri-diagonal, so special algorithms
([Link]’ algorithm) can be used to solve. The solution is given by
T01 100
1
T11 39.451 T1 39.451
1
T2 = 24.792 T21 24.792
T31 21.438 1 =
1 T3 21.438
T4 21.477 T 1 21.477
41
T5 25
Ec. derivadas parciales 486
Ejemplo 2 EDP Parabólica:
Método implícito
• Nodal temperatures when: t = 3 sec t = 6 sec t = 9 sec
T01 100 T02 100 T03 100
1 2 3
T
39 . 451 T
51. 326
1 1 T
1
59 . 043
T2 24.792
1
T2 30.669
2
T23 36.292
1 = 2 = 3 =
3
T 21 . 438 3
T 23 . 876 T3 26.809
T 21.477
1 T 22.836
2 T 3 24.243
41 42 43
5
T 25 T
5 25
5
T 25
∂x 2
O(∆x) 2, while our approximation of ∂T was of O(∆t ) accuracy.
∂t
• One can achieve similar orders of accuracy by approximating the
second derivative, on the left hand side of the heat equation, at the
midpoint of the time step. Doing so yields
∂ 2T α Ti +j1 − 2Ti j + Ti −j1 Ti +j1+1 − 2Ti j +1 + Ti −j1+1
≈ +
∂x 2 i, j
2 (∆x )2 (∆x )2
• The first derivative, on the right hand side of the heat equation, is
approximated using the forward divided difference method at time level
j +1, j +1
∂T j
Ti − Ti
≈
∂t i, j ∆t
• giving
− λTi −j1+1 + 2(1 + λ )Ti j +1 − λTi +j1+1 = λTi −j1 + 2(1 − λ )Ti j + λTi +j1
• where
∆t
λ =α
(∆x )2
• Having rewritten the equation in this form allows us to discretize the
physical problem. We then solve a system of simultaneous linear
equations to find the temperature at every node at any point in time.
i=0 1 2 3 4 5
T =100 ° C T = 25 °C
0.01m
• All internal nodes are at 20°C for t = 0 sec: Ti 0 = 20°C , for all i = 1,2,3,4
T00 = 100°C
T10 = 20°C
We can now form the system of equations for the
T20 = 20°C
first time step by writing the approximated heat
Interior nodes conduction equation for each node
T30 = 20°C
T40 = 20°C − λTi −j1+1 + 2(1 + λ )Ti j +1 − λTi +j1+1 = λTi −j1 + 2(1 − λ )Ti j + λTi +j1
T50 = 25°C
1
0 0 − 0. 4239 2. 8478 4
T 52 . 718
• The above coefficient matrix is tri-diagonal, so special algorithms
([Link]’ algorithm) can be used to solve. The solution is given by
T01 100
1
T11 44.372 T1 44.372
1
T
= 23 .746
T21 23.746
1 =
2
T31 20.797
1 3
T 20 . 797
T4 21.607 T 21.607
1
41
T5 25
1 = 2 = 3 =
3
T 20. 797 3
T 23 .174 3
T 26 . 562
T 21.607
1 T 22.730
2 T 24.042
3
41 42 43
5
T 25
T
5 25 T5 25









