0% encontró este documento útil (0 votos)
125 vistas17 páginas

Práctica MATLAB: Sumas de Riemann y Trapezoidal

Este documento presenta una práctica de laboratorio sobre métodos numéricos para aproximar el cálculo de integrales. Se explican conceptos básicos de integración y se proponen ejercicios para implementar la suma de Riemann y la suma trapezoidal usando Matlab, comparando los resultados con el valor exacto.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
125 vistas17 páginas

Práctica MATLAB: Sumas de Riemann y Trapezoidal

Este documento presenta una práctica de laboratorio sobre métodos numéricos para aproximar el cálculo de integrales. Se explican conceptos básicos de integración y se proponen ejercicios para implementar la suma de Riemann y la suma trapezoidal usando Matlab, comparando los resultados con el valor exacto.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

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

También podría gustarte