Autor Tema: Aproximar las integrales

0 Usuarios y 2 Visitantes están viendo este tema.

27 Octubre, 2023, 12:09 pm
Respuesta #10

ancape

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,081
  • País: es
  • Karma: +0/-5
  • Sexo: Masculino
Hola

Creo que el problema está puesto para hacer notar que es muy importante examinar la convergencia de una integral antes de aplicar un método numérico. En el caso a) no tendremos ningún problema pues el valor exacto de la integral es \( pi/4 \) En el caso b), el método del trapecio pone de manifiesto el problema de dividir por cero cuando va a aplicarse, pero en de Newton-Côtes abierto, no considera los extremos del intervalo de integración y da como resultado un número finito incompatible con la divergencia de la integral.

Saludos

29 Octubre, 2023, 04:22 pm
Respuesta #11

esmeraldabrown

  • $$\Large \color{#5e8d56}\pi\,\pi\,\pi$$
  • Mensajes: 130
  • País: ve
  • Karma: +0/-1
Me pueden explicar cómo debo comparar las aproximaciones para saber si son compatibles con la fórmula de error?

29 Octubre, 2023, 08:12 pm
Respuesta #12

esmeraldabrown

  • $$\Large \color{#5e8d56}\pi\,\pi\,\pi$$
  • Mensajes: 130
  • País: ve
  • Karma: +0/-1
Para calcular la integral del
\( \displaystyle\int_{1}^{10}\frac{1}{x} dx \)

Por la regla del trapecio

\( \frac{h}{2}[f(x_0)+f(x_1)] \)

Luego
\( \frac{10-1}{2}[\frac{1}{1}+\frac{1}{10}] \)
\( = \frac{9}{2}[\frac{11}{10}]= 4.95 \)

\( E =|\frac{h^3}{12}f"(\xi)| \)

\( E =|\frac{(10-1)^3}{12} \frac{2}{t^3}|≤0.1215 \)

Es correcto? Cómo hago para saber si las aproximaciones son compatibles con la fórmula de error?

Ahora cómo puedo aplicar el método de Newton cotes?

30 Octubre, 2023, 10:56 am
Respuesta #13

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 58,907
  • País: es
  • Karma: +0/-0
Hola

Para calcular la integral del
\( \displaystyle\int_{1}^{10}\frac{1}{x} dx \)

Por la regla del trapecio

\( \frac{h}{2}[f(x_0)+f(x_1)] \)

Luego
\( \frac{10-1}{2}[\frac{1}{1}+\frac{1}{10}] \)
\( = \frac{9}{2}[\frac{11}{10}]= 4.95 \)

\( E =|\frac{h^3}{12}f"(\xi)| \)

\( E =|\frac{(10-1)^3}{12} \color{red}\frac{2}{t^3}\color{black}|≤0.1215 \)

No está bien. Tienes que acotar \( 2/t^3 \) cuando \( t\in [1,10] \). No se de de donde sale el \( 0.1215 \).

\( E =|\frac{(10-1)^3}{12} \color{red}\frac{2}{t^3}\color{black}|\leq \dfrac{9^3\cdot 2}{12\cdot 1^3} \)

Citar
Es correcto? Cómo hago para saber si las aproximaciones son compatibles con la fórmula de error?

Compara la estimación del error con el error real.

Citar
Ahora cómo puedo aplicar el método de Newton cotes?

¿Cuál es la duda? ¡Aplica la fórmula!.

\( h=(10-1)/5 \)
\( x_0=2.8,\qquad x_1=4.6,\quad x_2=6.4,\quad x_3=8.2 \)


Saludos.

CORREGIDO (no me di cuenta de que eran abiertas)

30 Octubre, 2023, 12:48 pm
Respuesta #14

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 58,907
  • País: es
  • Karma: +0/-0
Hola

Luis desde donde parte la fórmula del error?

\( \displaystyle\int_{1}^{10} \frac{1}{x}dx= \frac{5(3)}{24}[11(1)+0.25+0
.14+11(0.1)]= 7.80625 \)

Corregí los valores; no me di cuenta de que era Newton-Cotes ABIERTA:

¿Cuál es la duda? ¡Aplica la fórmula!.

\( h=(10-1)/5 \)
\( x_0=2.8,\qquad x_1=4.6,\quad x_2=6.4,\quad x_3=8.2 \)


Saludos.

CORREGIDO (no me di cuenta de que eran abiertas)

Citar
Luis desde donde parte la fórmula del error?

¿A qué te refieres con esa pregunta?. No entiendo.

Saludos.

30 Octubre, 2023, 12:53 pm
Respuesta #15

esmeraldabrown

  • $$\Large \color{#5e8d56}\pi\,\pi\,\pi$$
  • Mensajes: 130
  • País: ve
  • Karma: +0/-1
Es que he venido calculando la integral y la fórmula de error aparte, como en el método del trapecio, no sé si me entiendas. En el caso de la primera integral como quedan los valores de x? De dónde sacas los cálculos,? Porque en esa debo aplicar también esta fórmula

30 Octubre, 2023, 06:49 pm
Respuesta #16

esmeraldabrown

  • $$\Large \color{#5e8d56}\pi\,\pi\,\pi$$
  • Mensajes: 130
  • País: ve
  • Karma: +0/-1
Luis no sé que hacer porque he hecho los cálculos y la integral me da 3.45625 y el error 3790,8 definitivamente creo que lo he hecho mal

30 Octubre, 2023, 07:03 pm
Respuesta #17

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 58,907
  • País: es
  • Karma: +0/-0
Hola

 Con Newton Cotes abierto para \( n=3 \) tienes:

\( \displaystyle\int_{a}^{b}f(x)dx\approx \dfrac{5h}{24}(11f(x_0)+f(x_1)+f(x_2)+11f(x_3)) \)

 con error \( \dfrac{95}{144}h^5\cdot f''''(t) \) con \( t\in (a,b) \), donde:

 \( h=\dfrac{b-a}{5},\quad x_i=a+(i+1)\cdot h \)

 Para \( \displaystyle\int_{0}^{\pi/2}sin^2(x)dx \) tienes:

\(  h=\pi/10 \)
\( x_0=\pi/10,\quad x_1=2\pi/10,\quad x_3=3\pi/10,\quad x_4=4\pi/10 \)

 La aproximación queda:

\(  I=0.785398163397448 \)

 Para cota de error \( f''''(x)=-8cos(2x) \) y así \( |f''''(x)|\leq 8 \). Queda

\( error\leq \dfrac{95}{144}h^5\cdot 8=0.0161510389192232 \)

  Para \( \displaystyle\int_{1}^{10}\dfrac{1}{x}dx \) tienes:

\(  h=(10-1)/2=9/5=1.8 \)
\( x_0=2.8,\quad x_1=4.6,\quad x_3=6.4,\quad x_4=8.2 \)

 La aproximación queda:

\(  I=2.11637855533253 \)

 Para cota de error \( f''''(x)=25/x^5 \) y así \( |f''''(x)|\leq 25 \) cuando \( x\in [1,10] \). Queda

\( error\leq \dfrac{95}{144}h^5\cdot 25=311.6475 \)

Saludos.

30 Octubre, 2023, 07:50 pm
Respuesta #18

esmeraldabrown

  • $$\Large \color{#5e8d56}\pi\,\pi\,\pi$$
  • Mensajes: 130
  • País: ve
  • Karma: +0/-1
Luis una última pregunta respecto a este tema, ¿cómo sé si las aproximaciones son compatibles con la fórmula de error? Y respecto a la integral b) si en realidad se debe calcular la integral de \( 0 \) hasta \( 10 \) de \( 1/x \) entonces no tiene solución por ninguno de los dos métodos?

31 Octubre, 2023, 11:05 am
Respuesta #19

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 58,907
  • País: es
  • Karma: +0/-0
Hola

Luis una última pregunta respecto a este tema, ¿cómo sé si las aproximaciones son compatibles con la fórmula de error?

Comprobando si:

\( |\text{valor aproximado}-\text{valor real}|\leq \text{error estimado} \)

Citar
Y respecto a la integral b) si en realidad se debe calcular la integral de \( 0 \) hasta \( 10 \) de \( 1/x \) entonces no tiene solución por ninguno de los dos métodos?

Sería bueno que más allá de las cuentas no pierdas de vista el significado de lo que estás haciendo.

Los métodos que estás estudiando son métodos para aproximar el valor de unas ciertas integrales definidas.

Si consideras la integral:

\( \displaystyle\int_{0}^{10}\dfrac{1}{x}dx \)

puedes comprobar que no converge; es una integral de una función no acotada que no existe (o si prefieres el resultado es infinito), dado que:

\( \displaystyle\lim_{a \to 0^+}{}\displaystyle\int_{a}^{10}\dfrac{1}{x}dx=+\infty \)

Entonces cualquier método que te de un valor concreto aproximado para esa integral estaría fallando, en el sentido de que pretendería aproximar una integral cuyo valor es infinito.

No obstante un podría intentar aplicarlos.

Ahora, el primer método (Trapecio):

\( \displaystyle\int_{a}^{b}f(x)dx=\dfrac{h}{2}(f(b)\color{red}+\color{black}f(a))-\underbrace{\dfrac{\color{blue}h^3\color{black}}{12}f''(t)}_{error} \)

NO puedes usarlo, porque en este caso \( a=0 \) y \( f(a)=f(0)=1/0 \) que NO existe.

Sin embargo Newton Cotes abierto para \( n=3 \) SI puedes aplicarlo:

Con Newton Cotes abierto para \( n=3 \) tienes:

\( \displaystyle\int_{a}^{b}f(x)dx\approx \dfrac{5h}{24}(11f(x_0)+f(x_1)+f(x_2)+11f(x_3)) \)

 con error \( \dfrac{95}{144}h^5\cdot f''''(t) \) con \( t\in (a,b) \), donde:

 \( h=\dfrac{b-a}{5},\quad x_i=a+(i+1)\cdot h \)

 En tu caso tendrías \( h=(10-0)/5=2 \) y:

\( x_0=2,\quad x_1=4,\quad x_2=6,\quad x_3=8 \).

Saludos.