| Autor |
Mensaje |
Basterman
Nivel 9
Edad: 34
Registrado: 28 Nov 2008
Mensajes: 2329
Carrera: Mecánica
|
|
Tengo la siguiente formula
![[tex]\int_{-1}^{1}f(x)*(1+x^{2})dx=af(-1)+bf(0)+cf(1)[/tex]](images/latex/2d1b7027d4420457a80355da24d9c7f61dd83f68_0.png "[tex]\int_{-1}^{1}f(x)*(1+x^{2})dx=af(-1)+bf(0)+cf(1)[/tex]")
Me piden encontrar a b c tal que la solucion sea exacta para polinomios mayor grado posible, en este caso seria 5.
Plantee el sistema de ecuaciones y suponiendo que conozco x1 x2 y x3 (-1 0 y 1) despeje a b y c. Aca tengo mi primer problema, ya que tenia entendido que las raices de este metodo y las constantes ya estan normalizadas, pero partiendo desde las x's que uso, nada me da lo mismo.
Asi que hice de cuenta que no sabia ninguna raiz ni nada y llegue a lo que tenia que llegar, hasta ahi todo piola.
El punto (b) se trata de aproximar ![[tex]\int_{0}^{1}x^{4}[/tex]](images/latex/d408be3940412c19948928b36b0a8aa78047041b_0.png "[tex]\int_{0}^{1}x^{4}[/tex]") usando la formula que halle en la primera parte. Usando mis datos, y despues de hacer el cambio de variables, no me da ni cerca de la solucion exacta, y usando las raices y coeficientes normalizados de este metodo, tampoco me da.
Cualquier ayuda es bienvenida, rindo hoy.
|
|
|
|
|
|
   |
    |
|
Basterman
Nivel 9
Edad: 34
Registrado: 28 Nov 2008
Mensajes: 2329
Carrera: Mecánica
|
|
Que foro choto, nadie sabe un carajo(?)
Bue, contestenla para Febrero al menos, no para mi porque hoy apruebo de alguna manera, pero a alguien mas le tiene que servir.
|
|
|
|
|
|
| |
|
|
violethill
Nivel 5
Edad: 35
Registrado: 27 Ago 2009
Mensajes: 152
Carrera: Informática
|
|
| no me acuerdo ni en pedo, baster... =(
|
|
|
|
_________________
And as i turned to you, you smiled at me, how could we say no?
|
|
   |
 |
|
Lautaz
Nivel 8
Registrado: 05 Sep 2008
Mensajes: 550
Carrera: Informática y Sistemas
|
|
Suena a cuadratura de Gauss pero la verdad ni idea jeje
¿En qué curso estás?
|
|
|
|
_________________
61.7
Death ... By exile
|
|
| |
|
|
Spike Spiegel
Nivel 9
Edad: 36
Registrado: 10 Ago 2007
Mensajes: 1507
Carrera: Informática
|
|
Ese problema es una cuadratura de Gauss PERO NO de Legendre, ahí está el punto. La idea es resolver con la definición pura de cuadratura de Gauss, donde ese (1+x^2) es la función de ponderación (busquen en Wikipedia).
El problema por el cual el resultado de aproximar x^4 no es exacto a pesar de tener grado menor a 5 lo voy a aclarar más abajo.
Vos ya tenés los puntos en los cuales evaluar para que sea una cuadratura de Gauss, así que solo te queda plantear como estimo que habrás hecho para sacar a, b y c.
Lo que hay que notar acá, es que esos valores en los que evaluás f(x) son las raíces de x(1-x^2), el cual es un polinomio ortogonal a 1+x^2 para el intervalo [-1,1] (la integral del producto de los polinomios da 0, remember producto interno en espacios de polinomios en Álgebra II). Es un detalle, pero vale la pena saberlo al que le interese entender bien el tema.
Ahora vamos al punto (b):
Vos tenés la integral entre 0 y 1 de x^4, que es lo mismo que 1/2 integral
entre -1 y 1 por ser función impar.
Lo que hay que hacer es llevar eso a la forma de nuestra cuadratura. Dividís y multiplicás por la función de peso 1+x^2, por lo que ahora la f(x) sobre la que vas a evaluar es x^4 / 1 +x^2.
Con eso ya tenés tu aproximación (te falta dividir por 2, obviamente).
Ahora, la pregunta que se hacen todos: ¿por qué no da exacto, si hasta grado 5 se banca?
Porque lo que tiene que ser un polinomio es esa f(x)con la que aproximo
En el caso de la cuadratura Gauss-Legendre, como la función de ponderación es w(x) = 1, todo mi integrando es igual a la f(x) con la que aproximo. Sin embargo, esta que tenemos acá es otra cuadratura inventada y nuestra f(x) usada para aproximar x^4 no es un polinomio sino una homográfica.
No sé si se entendió del todo, cualquier cosa preguntá de nuevo (vos o cualquiera que esté leyendo esto).
|
|
|
|
_________________
![[tex]\mbox{Detrás de todo 'tengo hambre' hay un gran 'comete esta'}[/tex]](images/latex/8abc5691357a123d82fd07cae688fd7767bf6633_0.png "[tex]\mbox{Detrás de todo 'tengo hambre' hay un gran 'comete esta'}[/tex]")
|
|
  |
    |
|
Basterman
Nivel 9
Edad: 34
Registrado: 28 Nov 2008
Mensajes: 2329
Carrera: Mecánica
|
|
Lo lei 2 meses despues casi, se entendio barbaro, dejo una pregunta nomas para cerrar este tema:
| Cita:
|
|
El problema por el cual el resultado de aproximar x^4 no es exacto a pesar de tener grado menor a 5 lo voy a aclarar más abajo.
|
O sea que esta bien que no me haya dado exacto? Ponele que me dio 0.18 y el valor exacto es 0.2.
Gracias.
|
|
|
|
|
|
| |
|
|
Spike Spiegel
Nivel 9
Edad: 36
Registrado: 10 Ago 2007
Mensajes: 1507
Carrera: Informática
|
|
| Claro, es totalmente lógico.
|
|
|
|
_________________
|
|
| |
|
|
AlanB
Nivel 9
Edad: 33
Registrado: 08 Mar 2010
Mensajes: 977
Ubicación: Quilmes
Carrera: Mecánica
|
|
| Espero que mañana no caigan con algo de esto que casi ni lo vimos este cuatrimestre. Como sacaron las constantes? Perdon pero no la cazo y el cuaderno no ayuda
|
|
|
|
_________________
|
|
 |
|
|
AlanB
Nivel 9
Edad: 33
Registrado: 08 Mar 2010
Mensajes: 977
Ubicación: Quilmes
Carrera: Mecánica
|
|
Calculé las constantes con la fórmula de la integral de la productoria pero no sé si no hice cualquiera. Alguien que me diga como lo sacan y cuánto les da? 
|
|
|
|
_________________
|
|
| |
|
|
Franzl
Nivel 7
Edad: 33
Registrado: 23 Ago 2011
Mensajes: 384
Carrera: Mecánica
|
|
Cuadratura de Gauß:
se parte de
![[tex]\int_{-1}^{1} f(x)w(x)dx[/tex]](images/latex/219581eb6e4cc4d231826dee606ab3ae37f417ad_0.png "[tex]\int_{-1}^{1} f(x)w(x)dx[/tex]") , y ![[tex]w(x) \geq 0[/tex]](images/latex/0f027575c8fcf518350a63a43ed4ffe1e39c3c90_0.png "[tex]w(x) \geq 0[/tex]")
se aproxima a la ![[tex]f(x)[/tex]](images/latex/07e15811bb43b993b6af429e34e77ed7c4129232_0.png "[tex]f(x)[/tex]") por un polinomio que será la sumatoria de la función evaluada en los ![[tex]x_i[/tex]](images/latex/1d6b6e3424f1379720112f7094b0d4dc7ab71b5e_0.png "[tex]x_i[/tex]") por la función loca (que me parece que es de Lagrange) ![[tex]\delta_i (x)[/tex]](images/latex/088c3a28d37b0141c8b9a414fdc155f897d261f9_0.png "[tex]\delta_i (x)[/tex]")
![[tex]\int_{-1}^{1} f(x)w(x)dx= \int_{-1}^{1} p(x)w(x)dx= \int_{-1}^{1} \sum f(x_i) \delta_i (x) dx = \sum f(x_i) \int_{-1}^{1} \delta_i (x) dx[/tex]](images/latex/58ac7fa11d3c42bebb8b01e8f59d221ed1f52b44_0.png "[tex]\int_{-1}^{1} f(x)w(x)dx= \int_{-1}^{1} p(x)w(x)dx= \int_{-1}^{1} \sum f(x_i) \delta_i (x) dx = \sum f(x_i) \int_{-1}^{1} \delta_i (x) dx[/tex]")
![[tex]\sum f(x_i) \int_{-1}^{1} \delta_i (x) dx = \sum f(x_i)A_i[/tex]](images/latex/83dbc5e18288b9a62277ddb32d63b82023f7f62c_0.png "[tex]\sum f(x_i) \int_{-1}^{1} \delta_i (x) dx = \sum f(x_i)A_i[/tex]")
con
![[tex]A_i= \int_{-1}^{1} w(x) \delta_i (x) dx[/tex]](images/latex/b4687d17dfaf455cefed2583a5e03508fbd977cb_0.png "[tex]A_i= \int_{-1}^{1} w(x) \delta_i (x) dx[/tex]")
![[tex] \delta_i (x)= \frac{ \prod_{j \neq i} (x-x_j)} { \prod_{j \neq i} (x_i - x_j)} [/tex]](images/latex/7e97fac72af8342277e45e4f90bcaa60887c7517_0.png "[tex] \delta_i (x)= \frac{ \prod_{j \neq i} (x-x_j)} { \prod_{j \neq i} (x_i - x_j)} [/tex]")
Fijate que el caso particular ![[tex] w(x)=1 [/tex]](images/latex/7307ac33bf78b16cfb6386fa651b7a44117c3765_0.png "[tex] w(x)=1 [/tex]") es el caso de Gauß-Legendre
Se entendió? Tengo el ejercicio hecho, si querés te lo escaneo.
|
|
|
|
|
|
 |
|
|
AlanB
Nivel 9
Edad: 33
Registrado: 08 Mar 2010
Mensajes: 977
Ubicación: Quilmes
Carrera: Mecánica
|
|
Ahi esta la integral+productoria que quise decir para los coeficientes. Dale copate y pasamelo por fb. Gracias querido.
Este cuatrimestre, por lo que tengo, Griggio lo tiro en los ultimos 10 minutos del cuatrimestre integrar por estos metodos 
|
|
|
|
_________________
|
|
| |
|
|
Franzl
Nivel 7
Edad: 33
Registrado: 23 Ago 2011
Mensajes: 384
Carrera: Mecánica
|
|
Mirá la verificación de abajo jaajajajaj
|
|
|
|
|
|
| |
|
|
Spike Spiegel
Nivel 9
Edad: 36
Registrado: 10 Ago 2007
Mensajes: 1507
Carrera: Informática
|
|
Alan, la cuadratura tiene 3 puntos de evaluación, por lo que debería integrar con exactitud polinomios de hasta grado 5 (acompañados por la función de peso).
Hacés 3 integrales (una por cada punto de evaluación) de polinomios boludos de distintos grados (1, x, x^2, por ej) (no te olvides de la función de peso) y te armás un sistema de 3x3 boludísimo de resolver. Listo, ya encontraste a, b y c.
|
|
|
|
_________________
|
|
| |
|
|
AlanB
Nivel 9
Edad: 33
Registrado: 08 Mar 2010
Mensajes: 977
Ubicación: Quilmes
Carrera: Mecánica
|
|
Ah cierto, en el Burden vi algo así. Gracias Spike maestro, ahí lo hice de esa forma y me dio igual que usando:
![[tex]A_i= \int_{-1}^{1} w(x) \frac{ \prod_{j \neq i} (x-x_j)} { \prod_{j \neq i} (x_i - x_j)} dx[/tex]](images/latex/3e934ced26b1f60fb6307773bbcae6c50672abc2_0.png "[tex]A_i= \int_{-1}^{1} w(x) \frac{ \prod_{j \neq i} (x-x_j)} { \prod_{j \neq i} (x_i - x_j)} dx[/tex]")
Ya que estamos, por si a alguien le interesa verificar o lo que sea, me dio:
a = 8/15
b = 8/5
c = 8/15
|
|
|
|
_________________
|
|
| |
|
|
Basterman
Nivel 9
Edad: 34
Registrado: 28 Nov 2008
Mensajes: 2329
Carrera: Mecánica
|
|
| Eso que dice Spike hice yo, y me dio lo mismo en solo 3 renglones, una fiesta.
|
|
|
|
|
|
| |
|
|
|
|
