48 votos

Un integral de Fresnel integrales 0((2 S(x)1)2+(2 C(x)1)2)2dxdex0((2 S(x)1)2+(2 C(x)1)2)2dxdex

Necesito calcular la siguiente integral: 0((2 S(x)1)2+(2 C(x)1)2)2x dx,0((2 S(x)1)2+(2 C(x)1)2)2x dx, donde S(x)=x0sinπz22dz,S(x)=x0sinπz22dz, C(x)=x0cosπz22dzC(x)=x0cosπz22dz son las integrales de Fresnel.

Integración numérica da un resultado aproximado de 0.31311841522422385...0.31311841522422385... que es cercano a los 16log28π216log28π2, por lo que podría ser la respuesta.

60voto

psychotik Puntos 171

Paso 1. La reducción de la integral

Deje que II denotar la integral en cuestión. Con el cambio de variable v=πx22v=πx22, tenemos

I=1π0{(12C(x))2+(12S(x))2}2dvI=1π0{(12C(x))2+(12S(x))2}2dv

donde x=2v/πx=2v/π es entendida como una función de la vv. Señalando que

12S(x)=2πv\pecadouuduy12C(x)=2πvcosuudu12S(x)=2πv\pecadouuduy12C(x)=2πvcosuudu

podemos escribir I $$ como

I=4π30|A(v)|4 dv\etiqueta1I=4π30|A(v)|4 dv\etiqueta1

donde A(v)A(v) denota la función definida por

A(v)=veuiudu.A(v)=veuiudu.


Paso 2. Simplificación de |A(v)|2|A(v)|2.

Ahora queremos simplificar |A(v)|2|A(v)|2. Para este fin, se nota que por u>0,

1u=1Γ(12)Γ(12)u1/2=1π0euxxdx=2π0eux2dx\etiqueta2

El uso de esta identidad,

Una(v)=2πveui0eux2dxdu=2π0vex2i)ududx=2eivπ0evx20ex2i)ududx=2eivπ0evx2x2idx.

Así, por medio de las coordenadas polares de cambio (x,y)(r,θ) seguido por las sustituciones r2=s y tanθ=t, obtenemos

|A(v)|2Y=A(v)¯A(v)=4π00ev(x2+y2)(x2i)(y2+i)dxdy=4π0π20revr2(r2cos2θi)(r2sin2θ+i)dθdr=2π0π20evs(scos2θi)(ssin2θ+i)dθds=2π0evssπ20(1scos2θi+1ssin2θ+i)dθds=2π0evss0(1si(t2+1)+1st2+i(t2+1))dtds.

Evaluación de la interna de la integral es fácil, y obtenemos

|A(v)|2Y=20evss(i1+)ds.

La aplicación de (2) de nuevo, nos encontramos con que

|A(v)|2Y=20evss(iπ0e(1+)uudu)ds=2π0evss0eusin(su)ududs=2π0euu0sin(su)sevsdsdu=2π0euuarctan(uv)du=4vπ0evx2arctan(x2)dx(u=vx2)\etiqueta3

Aquí, hemos explotado la identidad

0sinxxesxdx=arctan(1s),

lo que puede ser demostrado por la diferenciación de ambos lados con respecto a los s.


Paso 3. Evaluación de la I.

El taponamiento de (3) (1) y la aplicación de la coordenada polar cambio, I se reduce a

Yo=64π4000vev(x2+y2)arctan(x2)arctan(x2)dxdydv=64π400arctan(x2)arctan(x2)(x2+y2)2dxdy=64π4π200arctan(r2cos2θ)arctan(r2sin2θ)r3drdθ=32π4π200arctan(scos2θ)arctan(ssin2θ)s2dsdθ.(s=r2)\etiqueta4

Ahora vamos a denotar

J(u,v)=0arctan(nosotros)arctan(vs)s2ds.

Un simple cálculo muestra que

2JuvJ(u,v)=0ds(u2s2+1)(v2s2+1)=π2(u+v).

De hecho, tanto el contorno método de integración global o parcial de la fracción de descomposición método de trabajo aquí. La integración, hemos

J(u,v)=π2{(u+v)log(u+v)uloguvlogv}.

Conectando de a (4), se sigue que

Yoy=64π3π20sin2θlogsinθdθ=16π3βz(32,12)

donde β(z,w) es la función beta, la satisfacción de las siguientes beta de la función identidad

β(z,w)=20sin2z1θcos2w1θdθ=Γ(z)Γ(w)Γ(z+w).

Por lo tanto, tenemos

\begin{align*} Yo &= \frac{16}{\pi^{3}} \frac{\Gamma\left(\frac{3}{2}\right)\Gamma\left(\frac{1}{2}\right)}{\Gamma(2)} \left\{ \psi_{0} (2) - \psi_{0} \left(\tfrac{3}{2} \right) \right\} = \frac{8}{\pi^2} \int_{0}^{1} \frac{\sqrt{x} de {- x}{1 - x} \, dx = \frac{8 (2 \log 2 - 1)}{\pi^2}, \end{align*}

donde ψ0(z)=Γ(z)Γ(z) es la función digamma, que satisface la siguiente identidad

ψ0(z+1)=γ+101xz1xdx.

25voto

Graham Hesketh Puntos 4746

A raíz de Ron...

Hacer las sustituciones: x=y, p=1/2+1/21+4r, para obtener:

64π20x(1\pecado(2πx2p(p1))pdp)2dx,

=32π20(02sin(2πy r)1+4r(1+1+4r)dr)2dy, consulte el Apéndice,

=16π2(02sin(2πy\r)1+4r(1+1+4r)dr)2dy,

=16π2(H(r)sin(2πy\r)1+4|r|(1+1+4|r|)dr)2dy : H \left( r \right) =\casos{1&$0\leq x$\cr -1&$x<0$\cr}.

Nota entonces que:

(H(r)sin(2πy\r)1+4|r|(1+1+4|r|)dr)2=|i2H(r)ei2πyr1+4|r|(1+1+4|r|)dr+i2H(r)ei2πyr1+4|r|(1+1+4|r|)dr|2

=|(H(r)H(r))ei2πyr21+4|r|(1+1+4|r|)dr|2=|ei2πyr1+4|r|(1+1+4|r|)dr|2.

Ahora podemos usar el teorema de Plancherel que establece:

\displaystyle \int _{-\infty }^{\infty }\! \left| f \left y \right) 
 \right| ^{2}{dy}=\int _{-\infty }^{\infty }\! 
 \left| F \left( r \right) \right| ^{2}{dr}
: \displaystyle F \left( r \right) =\int _{-\infty }^{\infty }\!f \left y \right) {
e^{-i2\pi r}}{dy}
,

y por lo tanto:

16π2|ei2πyr1+4|r|(1+1+4|r|)dr|2dy
=16π2|11+4|r|(1+1+4|r|)|2dr,

=16π202(1+4\r)(1+1+4r)2dr.

Ahora si queremos deshacer la sustitución realizada por dejar r=p2p:

16π202(1+4\r)(1+1+4r)2dr=16π2112p2(2p1)dp,

y, a continuación, hacer un último sustitución p=1q+2, se obtiene:

16π2112p2(2p1)dp=16π21212+1qdq=16ln(2)8π2.

Apéndice:

Tenga en cuenta que la integral: 02sin(2πy\r)1+4r(1+1+4r)dr

converge por la Chartier-Dirichlet de la prueba debido a que: f(r)=21+4r(1+1+4r)

es monótona y continua en R+ y f(r)0 como r, y porque: |b0sin(2πy\r)dr| es delimitada como b.

15voto

Ron Gordon Puntos 96158

Esta no es una solución, pero creo que es un paso hacia adelante.

Definir

f(x)=(12S(x))cos(π2x2)(12C(x))sin(π2x2)

g(x)=(12C(x))cos(π2x2)+(12S(x))sin(π2x2)

A continuación, se trata de un sencillo ejercicio para demostrar que

g(x)+if(x)=eiπx2/2xdteiπt2/2

y que

(2S(x)1)2+(2C(x)1)2=4[g(x)2+f(x)2]=4x21du1dveiπx2(u2v2)/2

Puedo convertir la integral doble para una sola integral por cambio de coordenadas a p=u+v, q=uv, p[2,), p[(p2),p2]. El Jacobiano es de 1/2, y tenemos

(2S(x)1)2+(2C(x)1)2=2x22dpp2(p2)dqeiπx2pq/2

que después de la evaluación de la interna integral y algunos reescalado, obtenemos

(2S(x)1)2+(2C(x)1)2=8π1dpsin[2πx2p(p1)]p

No estoy seguro de cómo evaluar esta integral analíticamente, ni estoy seguro de que esa sería la mejor mover aquí. El deseado integral es, por tanto,

0dxx[(2S(x)1)2+(2C(x)1)2]2=64π20dxx1dpsin[2πx2p(p1)]p1dqsin[2πx2q(q1)]q

A partir de aquí, no estoy muy seguro de qué hacer. Naturalmente, me gustaría invertir el orden de integración de manera que la integral de más de x es interior, pero no estoy seguro de cómo justificar que, dado que las integrales no son absolutamente convergentes. Además, incluso si puedo justificar que, creo que el resultado de la integral es algunas funciones delta de que no se ven prometedores. Sigo a ver esto, pero me imagino que tal vez alguien más puede contribuir también desde aquí, o sólo dime si estoy fuera de la marca.

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