Loading [MathJax]/jax/element/mml/optable/MathOperators.js

4 votos

Integral de exponencial con x(1x) término

¿La siguiente integral tiene una solución de forma cerrada?

y0exp(x(1x))dx

¿O debo conformarme con una aproximación?


Editar: La forma actual de la integral tiene una α : y0exp(αx(1x))dx Necesito que la solución (o solución aproximada) funcione con grandes α (y si es posible para pequeños α también


2ª edición: ¿Existe una expresión analítica para la integral definida que contiene el α ? 10exp(αx(1x))dx


3ª Edición: En realidad es α no α

7voto

JohnDoe Puntos 16

Possible direction

x(1x)=14(x12)2=1214(x12)2. haciendo así un cambio de variable cos(t)=2(x12),sin(t)dt=2dx. podemos reescribir la integral como 120πeα2sin(t)sin(t)dt=12π0eα2sin(t)sin(t)dt

por lo que la integral original se ve como 12π0eα2sin(s)sin(s)ds=12π[L1(α2)+I1(α2)]1.48983 update: Parece haber una discrepancia entre esta "respuesta" y el resultado numérico obtenido en los comentarios anteriores. Si pudieran por favor abstenerse de votar hasta que las diferencias sean explicadas, eso sería muy grande. ¡Salud!

update 2 La última ecuación ha sido modificada de mi respuesta anterior debido a mi tonto error de convertir los límites de la integración! De todos modos, usando la matemática (disculpas) la integral final es de la forma función Struve modificada L1(x) . La respuesta es aproximadamente 1.48983 para α=1 y limita de 0 a 1, lo que corresponde al comentario anterior.

update 3 Respuesta modificada como se menciona en los comentarios a la pregunta original. Ahora trata de encontrar una manera de representar las funciones especiales con argumentos negativos.

5voto

Claude Leibovici Puntos 54392

Para la gama 0y1 puedes tener una estimación muy precisa de la integral expandiendo primero la integral como una serie de Taylor construida en x=0 . Esto da ex(1x)=1+x+x2x3/2311x22411x5/23059x372013x7/2630+1513x440320311x9/222680+14761x5362880031417x11/21247400594659x64790016001877123x13/2972972008409491x787178291200+O(x15/2) Entonces, para la integral, tienes y0ex(1x)dx=y+2y3/23+y242y5/21511y37211y7/210559y4288013y9/22835+1513y5201600311y11/2124740+14761y62177280031417y13/28108100594659y733530112001877123y15/27297290008409491y8697426329600+O(y17/2) lo cual es bastante exacto.

Añadido más tarde después de la solicitud de la OP

Si consideramos el caso de y0eax(1x)dx asumiendo que a es pequeño, un procedimiento idéntico primero lleva a eax(1x)=1+ax+a2x2+16a(a23)x3/2+124a2(a212)x2+1120a(a430a215)x5/2+1720a4(a260)x3+a(a6105a4+315a2315)x7/25040+a4(a4168a2+1680)x440320+a(a8252a6+5670a4+3780a214175)x9/2362880+a6(a4360a2+15120)x53628800+a(a10495a8+34650a6103950a4+155925a21091475)x11/239916800+a6(a6660a4+71280a2665280)x6479001600+a(a12858a10+135135a82702700a62027025a4+12162150a2127702575)x13/26227020800+a8(a61092a4+240240a28648640)x787178291200+O(x15/2) y, luego, para la integral y0eax(1x)dx=y+23ay3/2+a2y24+115a(a23)y5/2+172a2(a212)y3+1420a(a430a215)y7/2+a4(a260)y42880+a(a6105a4+315a2315)y9/222680+a4(a4168a2+1680)y5201600+a(a8252a6+5670a4+3780a214175)y11/21995840+a6(a4360a2+15120)y621772800+a(a10495a8+34650a6103950a4+155925a21091475)y13/2259459200+a6(a6660a4+71280a2665280)y73353011200+a(a12858a10+135135a82702700a62027025a4+12162150a2127702575)y15/246702656000+a8(a61092a4+240240a28648640)y8697426329600+O(y17/2)

1voto

stealth_angoid Puntos 429

Siento no poder escribir un buen contenido matemático en este sitio todavía, me acabo de inscribir.

Creo que es posible encontrar una solución analítica a esto. Primero la función : x>x(1x) es bijectivo cuando x está entre 0 y 12 . Así que podemos asumir que x verifica eso, o separamos la integral en dos, una parte de 0 a 12 y el otro de 12 sabiendo que el método para calcularlo es el mismo.

Primero cambiamos las variables con u2=x(1x) que da x=(114u22) (dos soluciones pero sólo una de ellas es más pequeña que 12 ))

reemplazarlo en el integral y utilizar una integración por partes en 4u(14u2) que debería aparecer después de la correcta reescritura de la integral tras el cambio de variable.

Entonces, para simplificar la notación, use otra: v=2u ( podríamos hacerlo antes de la integración por partes)

Si estoy en lo cierto, deberías obtener una cierta función de y (la variable principal) que está hecha de raíces y exponenciales por lo que se conoce) y (1v2)ev2 delimitado por los valores adecuados (aquí 0 y 12y(1y) ). Luego usas un cambio de pecado de la variable : v=sin(t) y te da :

el mismo término función de y, y cos(t)2esin(t)2 A partir de aquí, se utiliza la integración de dos partes para clasificar la cosa y debería darte la solución analítica, por fea que sea :) . Puede que me haya equivocado en algún momento, así que por favor no dude en decirme si resulta ser así. Gracias a los que leen esto hasta el final :)

1voto

fcop Puntos 2891

y0eαx(1x) dx

=y0eα(x2x) dx

=y0eα(x2x+1414) dx

=y0eα14(x12)2 dx

=y1212eα14x2 dx

=cos1(2y1)πeα14(cosx2)2 d(cosx2)

=12πcos1(2y1)eαsinx2sinx dx

=πcos1(2y1)n=0α2nsin2n+1x22n+1(2n)!dxπcos1(2y1)n=0α2n+1sin2n+2x4n+1(2n+1)!dx

Para n es cualquier número entero no negativo,

sin2n+2x dx=(2n+2)!x4n+1((n+1)!)2nk=0(2n+2)!(k!)2sin2k+1xcosx4nk+1((n+1)!)2(2k+1)!+C

Este resultado puede lograrse mediante la integración sucesiva por partes.

sin2n+1x dx

=sin2nx d(cosx)

=(1cos2x)n d(cosx)

=nk=0Cnk(1)kcos2kx d(cosx)

=nk=0(1)k+1n!cos2k+1xk!(nk)!(2k+1)+C

= \left [ \sum\limits_ {n=0}^ \infty\sum\limits_ {k=0}^n \dfrac {(-1)^{k+1}n! \alpha ^{2n} \cos ^{2k+1}x}{2^{2n+1}(2n)!k!(n-k)!(2k+1)} \right ]_{ \cos ^{-1}(2y-1)}^ \pi - \left [ \sum\limits_ {n=0}^ \infty\dfrac { \alpha ^{2n+1}x}{2^{4n+3}n!(n+1)!}- \sum\limits_ {n=0}^ \infty\sum\limits_ {k=0}^n \dfrac {(k!)^2 \alpha ^{2n+1} \sin ^{2k+1}x \cos x}{2^{4n-2k+3}n!(n+1)!(2k+1)!} \right ]_{ \cos ^{-1}(2y-1)}^ \pi

= \sum\limits_ {n=0}^ \infty\sum\limits_ {k=0}^n \dfrac {(-1)^kn! \alpha ^{2n}((2y-1)^{2k+1}+1)}{2^{2n+1}(2n)!k!(n-k)!(2k+1)}+ \sum\limits_ {n=0}^ \infty\dfrac { \alpha ^{2n+1}( \cos ^{-1}(2y-1)- \pi )}{2^{4n+3}n!(n+1)!}+ \sum\limits_ {n=0}^ \infty\sum\limits_ {k=0}^n \dfrac {(k!)^2 \alpha ^{2n+1}(2y-1)(1-(2y-1)^2)^{k+ \frac {1}{2}}}{2^{4n-2k+3}n!(n+1)!(2k+1)!}

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X