Universidad Laica Eloy Alfaro de
Manabí
Carrera de Análisis Numérico
Ingeniería Eléctrica
Practica #1
Paralelo A
Nombre del profesor Prof. MSc. Marcos Ponce Jara
Nombre del alumno Arroyo Cerezeda Dilia Daniela
Semestre 4
Fecha 20/09/2023
Fecha presentación 10/10/2023
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
Universidad Laica Eloy Alfaro de
Manabí
PRACTICAS DE LABORATORIO
Facultad de ingeniería Análisis Numérico - MATLAB
Ingeniería eléctrica
Practica Nº1
Suma de Riemann y Suma trapezoidal
Prof. MSc. Marcos Ponce Jara
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
ÍNDICE
1.1. Comandos básicos Matlab para integración .........................................................................4
2. Introducción a la Practica .............................................................................................................5
3. Objetivos de la practica.................................................................................................................5
4. Programas utilizados para las practicas ........................................................................................5
5. Ejercicios de evaluación practica ..................................................................................................6
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
1.1. Comandos básicos Matlab para integración
El comando “int” calcula de manera simbólica la integral de la función f.
Ejemplo: syms x
int(x^2/(x^6 - 8))
Para aproximar la integral de una función f mediante sumas de Riemann y la respectiva
representación de rectángulos se utiliza el comando “rsum”
Ejemplo: syms x
rsums (exp(-2*x))
Se utiliza el comando “trapz(x,y)” para calcular la integral de una función aproximando el
valor de la integral de Y con respecto a X
Ejemplo: X=0:0.001:1
Y = (exp(-2*x))
trapz(X, Y)
Para calcular las primitivas de una función se debe declara las variables como simbólicas y
posteriormente realizar los cálculos pertinentes:
Ejemplo: syms x a b
f = sin(a*x)+cos (b*x);
integral=diff (f, x);
pretty(integral)
Si queremos evaluar una integral entre un intervalo determinado, una de las soluciones es
utilizar la siguiente función “integral (función, intervalo1, intervalor2)”
Ejemplo: fun = @(x) exp(-x^2)
(aquí introducimos el intervalo de integración) q
= integral(fun,0,Inf)
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
2. Introducción a la Practica
Esta segunda práctica nos introduce a los métodos numéricos como forma de calcular de
forma aproximada el área bajo una curva utilizando métodos iterativos de cálculo. En este
caso se desea programar un script en Matlab para graficar y calcular áreas mediante el
método de Suman de Riemann y el método de Suma Trapezoidal. Se comparará cada método
respecto el resultado exacto, obteniendo de esta manera los errores cometidos en cada caso.
3. Objetivos de la practica
• Entender con profundidad el concepto de integración
• Utilización de métodos iterativos para el cálculo de áreas
• Familiarizarse con los distintos comandos para calcular integrales definidas e
indefinidas de forma exacta
• Graficar el método de la Suma de Riemann y Suma Trapezoidal
• Cálculo de errores
• Utilización del bucle FOR como método iterativo para el cálculo de áreas
4. Programas utilizados para las practicas
Se utilizará el software Matlab rb2013 o en su defecto la versión más reciente instalada
en los laboratorios de la facultad de ingeniería. Adicionalmente se utilizará el programa
Wolfram CDF (versión de prueba) para comprobar que nuestros resultados son correctos.
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
5. Ejercicios de evaluación practica
I. Suma de Riemann (sumatoria puntos externos derechos):
Sea A el área de la región que está bajo la gráfica de f(x) = 𝒆−𝟐𝒙, entre x = 0 y x = 2.
Se pide:
a) Explique el concepto de integral con sus propias palabras
La integral para calcular el área bajo la curva de una función o la acumulación de una
cantidad a lo largo de una variable, nos permite encontrar la suma infinitesimal.
b) Encuentre una expresión para A como un límite (no evalué ese límite)
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
c) Estime el área tomando los puntos muestra 2, 4, 8 y 10 (evalúelos manualmente).
d) Estime el área creando un programa para calcular áreas mediante el método
numérico de Riemann. Compruebe los resultados calculados manualmente y calcule
el área tomando 100, 500 subintervalos. El programa tiene que generar una gráfica
como la que se muestra a continuación:
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
e) Dentro del programa, calcule la integral exacta de f(x) en el intervalo [0, 2] y muestre
el error cometido entre la aproximación y el cálculo exacto.
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
Nota: Para comprobar que están realizando el programa de forma correcta comprobar
los resultados con los mostrados en la siguiente aplicación con la opción “right sum”:
http://demonstrations.wolfram.com/ComparingBasicNumericalIntegrationMethods/
II. Suma de Riemann (sumatoria puntos externos izquierdos):
f) Realice exactamente lo mismo que el ejercicio
anterior pero tomando en cuenta la sumatoria con los
puntos externos de la izquierda de los rectángulos
(apartados c, d y e)
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
Nota: comprobar los resultados con la opción “right sum”:
http://demonstrations.wolfram.com/ComparingBasicNumericalIntegrationMethods/
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
III. Suma de Riemann (Sumatorios con punto medio)
g) Realice exactamente lo mismo que el ejercicio
anterior pero tomando en cuenta el punto medio de los
rectángulos (apartados c, d y e)
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
Nota: comprobar los resultados con la opción “midpoint”:
http://demonstrations.wolfram.com/ComparingBasicNumericalIntegrationMethods/
IV. Compare los resultados y exponga cual es la mejor aproximación y explique porque.
Observando las graficas y los resultados de los puntos derecho, izquierdo y medio se puede
observar que la mejor aproximación es la del punto medio ya que se acerca a la área exacta
y también tiene menos error
V. Suma de trapecios Tomando la misma función de los apartados anteriores se pide
a) Encuentre una expresión para A como un límite (no evalué ese límite)
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
b) Estime el área tomando los puntos muestra 2, 4, 8 y 10 (evalúelos manualmente).
c) Estime el área creando un programa para calcular áreas mediante el método numérico
de Riemann. Compruebe los resultados calculados manualmente y calcule el área tomando
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
100, 500 subintervalos. El programa tiene que generar una gráfica como la que se muestra a
continuación:
Prof. MSc. Ing. Marcos Ponce Jara
Universidad Laica
Análisis Numérico Ingeniería Electrica
Eloy Alfaro de Manabí
c) ¿Introduce alguna mejora esta nueva aproximación? ¿por qué?
Si, ya que con el método el trapecio se puede observar que el valor de error es mínimo
y también se acerca más al área exacta
Nota: comprobar los resultados con la opción “trapez”:
http://demonstrations.wolfram.com/ComparingBasicNumericalIntegrationMethods/
Para los dos ejercicios se tiene que presentar e código que se creó, explicando paso a paso los
comandos utilizados en Matlab.
Prof. MSc. Ing. Marcos Ponce Jara