Gli integrali di Fresnel ,
S
(
x
)
{\displaystyle S(x)}
e
C
(
x
)
{\displaystyle C(x)}
, sono due funzioni speciali trascendenti introdotte in ottica dall'ingegnere francese Augustin-Jean Fresnel per studiare i fenomeni della diffrazione .
Grafico degli integrali di Fresnel normalizzati:
S
(
x
)
=
∫
0
x
sin
(
π
2
t
2
)
d
t
{\displaystyle S(x)=\int _{0}^{x}\sin \left({\frac {\pi }{2}}t^{2}\right)\,dt}
e
C
(
x
)
=
∫
0
x
cos
(
π
2
t
2
)
d
t
{\displaystyle C(x)=\int _{0}^{x}\cos \left({\frac {\pi }{2}}t^{2}\right)\,dt}
.
Esse sono definite attraverso le seguenti rappresentazioni:
Grafico degli stessi integrali non normalizzati:
S
(
x
)
=
∫
0
x
sin
(
t
2
)
d
t
{\displaystyle S(x)=\int _{0}^{x}\sin \left(t^{2}\right)\,dt}
e
C
(
x
)
=
∫
0
x
cos
(
t
2
)
d
t
{\displaystyle C(x)=\int _{0}^{x}\cos \left(t^{2}\right)\,dt}
.
S
(
x
)
:=
∫
0
x
sin
(
π
2
t
2
)
d
t
{\displaystyle S(x):=\int _{0}^{x}\sin \left({\frac {\pi }{2}}t^{2}\right)\,dt}
C
(
x
)
:=
∫
0
x
cos
(
π
2
t
2
)
d
t
{\displaystyle C(x):=\int _{0}^{x}\cos \left({\frac {\pi }{2}}t^{2}\right)\,dt}
anche se altri autori preferiscono definirle senza il
π
2
{\displaystyle {\frac {\pi }{2}}}
nell'argomento di seno e coseno .
S
(
x
)
{\displaystyle S(x)}
e
C
(
x
)
{\displaystyle C(x)}
sono funzioni dispari .
C
(
i
z
)
=
i
C
(
z
)
{\displaystyle C(iz)=iC(z)}
S
(
i
z
)
=
−
i
S
(
z
)
{\displaystyle S(iz)=-iS(z)}
Gli integrali di Fresnel non possono essere calcolati in forma chiusa in termini di funzioni elementari , salvo casi particolari. Infatti essi convergono all'infinito e si ha:
lim
x
→
+
∞
S
(
x
)
=
lim
x
→
+
∞
C
(
x
)
=
1
2
.
{\displaystyle \lim _{x\to +\infty }S(x)=\lim _{x\to +\infty }C(x)={\frac {1}{2}}\,.}
Poiché gli integrali di Fresnel non possono essere calcolati coi metodi tradizionali, una possibile dimostrazione di
∫
0
+
∞
cos
(
x
2
)
d
x
=
π
8
{\displaystyle \int _{0}^{+\infty }\cos {(x^{2})}\,dx={\sqrt {\frac {\pi }{8}}}}
sfrutta l'analisi complessa e il risultato dell'integrale di Gauss
∫
0
+
∞
e
−
x
2
d
x
=
π
2
{\displaystyle \int _{0}^{+\infty }e^{-x^{2}}dx={\frac {\sqrt {\pi }}{2}}}
. L'integrale di partenza può essere scritto come parte reale di un numero complesso secondo quella che è la forma polare di un numero complesso:
∫
0
+
∞
cos
(
x
2
)
d
x
=
ℜ
(
∫
0
+
∞
e
i
x
2
d
x
)
.
{\displaystyle \int _{0}^{+\infty }\cos {(x^{2})}\,dx=\Re {\left(\int _{0}^{+\infty }e^{ix^{2}}dx\right)}.}
Curva semplice chiusa
γ
{\displaystyle \gamma }
nel piano complesso, suddivisa in
γ
1
{\displaystyle \gamma _{1}}
,
γ
2
{\displaystyle \gamma _{2}}
e
γ
3
{\displaystyle \gamma _{3}}
.
Per calcolare il secondo integrale si sfrutta il teorema di Cauchy-Goursat scegliendo come cammino chiuso di integrazione la curva chiusa
γ
{\displaystyle \gamma }
suddivisibile nei tre tratti
γ
1
{\displaystyle \gamma _{1}}
,
γ
2
{\displaystyle \gamma _{2}}
e
γ
3
{\displaystyle \gamma _{3}}
come in figura:
∮
γ
e
i
z
2
d
z
=
∫
γ
1
e
i
z
2
d
z
+
∫
γ
2
e
i
z
2
d
z
+
∫
γ
3
e
i
z
2
d
z
=
0
.
{\displaystyle \oint _{\gamma }e^{iz^{2}}dz=\int _{\gamma _{1}}e^{iz^{2}}dz+\int _{\gamma _{2}}e^{iz^{2}}dz+\int _{\gamma _{3}}e^{iz^{2}}dz=0\,.}
Questa operazione si può fare perché la funzione
e
i
z
2
{\displaystyle e^{iz^{2}}}
è analitica in
C
{\displaystyle \mathbb {C} }
, che è semplicemente connesso .
Nel piano complesso
γ
3
{\displaystyle \gamma _{3}}
ha equazione
z
=
r
e
i
ϑ
{\displaystyle z=re^{i\vartheta }}
, con
r
{\displaystyle r}
variabile; per ricondursi all'integrale della gaussiana si impone che l'inclinazione di tale retta sia tale che
i
z
2
=
i
(
r
e
i
ϑ
)
2
=
−
r
2
{\displaystyle iz^{2}=i{\left(re^{i\vartheta }\right)}^{2}=-r^{2}}
, ovvero
ϑ
=
π
4
{\displaystyle \vartheta ={\frac {\pi }{4}}}
. Il terzo integrale diventa quindi
∫
γ
3
e
i
z
2
d
z
=
∫
R
0
e
−
r
2
e
i
π
4
d
r
,
{\displaystyle \int _{\gamma _{3}}e^{iz^{2}}dz=\int _{R}^{0}e^{-r^{2}}e^{i{\frac {\pi }{4}}}\,dr\,,}
che per
x
{\displaystyle x}
, ovvero
R
→
+
∞
{\displaystyle R\rightarrow +\infty }
, vale
lim
R
→
+
∞
e
i
π
4
∫
R
0
e
−
r
2
d
r
=
−
e
i
π
4
π
2
.
{\displaystyle \lim _{R\rightarrow +\infty }e^{i{\frac {\pi }{4}}}\int _{R}^{0}e^{-r^{2}}dr=-e^{i{\frac {\pi }{4}}}{\frac {\sqrt {\pi }}{2}}\,.}
La curva
γ
2
{\displaystyle \gamma _{2}}
può essere parametrizzata come
z
=
R
e
i
ϑ
{\displaystyle z=Re^{i\vartheta }}
, questa volta con
ϑ
{\displaystyle \vartheta }
variabile. Il secondo integrale diventa
∫
γ
2
e
i
z
2
d
z
=
∫
0
π
4
e
i
R
2
e
2
i
ϑ
(
i
R
)
e
i
ϑ
d
ϑ
=
i
R
∫
0
π
4
e
i
R
2
cos
(
2
ϑ
)
+
i
ϑ
e
−
R
2
sin
(
2
ϑ
)
d
ϑ
.
{\displaystyle \int _{\gamma _{2}}e^{iz^{2}}dz=\int _{0}^{\frac {\pi }{4}}e^{iR^{2}e^{2i\vartheta }}(iR)\,e^{i\vartheta }\,d\vartheta =iR\int _{0}^{\frac {\pi }{4}}e^{iR^{2}\cos {(2\vartheta )}+i\vartheta }e^{-R^{2}\sin {(2\vartheta )}}\,d\vartheta \,.}
Per
0
≤
ϑ
≤
π
4
{\displaystyle 0\leq \vartheta \leq {\frac {\pi }{4}}}
,
sin
(
2
ϑ
)
≥
0
{\displaystyle \sin {(2\vartheta )}\geq 0}
e
cos
(
2
ϑ
)
≥
0
{\displaystyle \cos {(2\vartheta )}\geq 0}
, e vale la disuguaglianza
sin
ϑ
>
2
π
ϑ
{\displaystyle \sin {\vartheta }>{\frac {2}{\pi }}\vartheta }
. Ponendo
2
ϑ
=
φ
{\displaystyle 2\vartheta =\varphi }
, è possibile fare la seguente maggiorazione:
0
≤
|
R
2
∫
0
π
2
e
i
R
2
cos
φ
+
i
φ
2
e
−
R
2
sin
φ
d
φ
|
≤
R
2
∫
0
π
2
|
e
i
R
2
cos
φ
+
i
φ
2
|
|
e
−
R
2
sin
φ
|
d
φ
=
R
2
∫
0
π
2
e
−
R
2
sin
φ
d
φ
≤
R
2
∫
0
π
2
e
−
R
2
2
π
φ
d
φ
=
−
π
4
R
(
e
−
R
2
−
1
)
,
{\displaystyle 0\leq \left|{\frac {R}{2}}\int _{0}^{\frac {\pi }{2}}e^{iR^{2}\cos {\varphi }+i{\frac {\varphi }{2}}}e^{-R^{2}\sin {\varphi }}\,d\varphi \right|\leq {\frac {R}{2}}\int _{0}^{\frac {\pi }{2}}\left|e^{iR^{2}\cos {\varphi }+i{\frac {\varphi }{2}}}\right|\,\left|e^{-R^{2}\sin {\varphi }}\right|d\varphi ={\frac {R}{2}}\int _{0}^{\frac {\pi }{2}}e^{-R^{2}\sin {\varphi }}\,d\varphi \leq {\frac {R}{2}}\int _{0}^{\frac {\pi }{2}}e^{-R^{2}{\frac {2}{\pi }}\varphi }\,d\varphi =-{\frac {\pi }{4R}}\left(e^{-R^{2}}-1\right),}
e dal teorema del confronto , segue che per
R
→
+
∞
{\displaystyle R\rightarrow +\infty }
il secondo integrale vale
0
{\displaystyle 0}
.
La curva
γ
1
{\displaystyle \gamma _{1}}
, infine, può essere parametrizzata come
z
=
x
=
R
{\displaystyle z=x=R}
. Dal teorema di Cauchy-Goursat
∫
0
+
∞
e
i
x
2
d
x
=
∫
γ
1
e
i
z
2
d
z
=
−
∫
γ
3
e
i
z
2
d
z
=
e
i
π
4
π
2
.
{\displaystyle \int _{0}^{+\infty }e^{ix^{2}}dx=\int _{\gamma _{1}}e^{iz^{2}}dz=-\int _{\gamma _{3}}e^{iz^{2}}dz=e^{i{\frac {\pi }{4}}}{\frac {\sqrt {\pi }}{2}}\,.}
L'integrale di Fresnel cercato diventa perciò
∫
0
+
∞
cos
(
x
2
)
d
x
=
ℜ
(
∫
0
+
∞
e
i
x
2
d
x
)
=
ℜ
(
e
i
π
4
π
2
)
=
π
8
,
{\displaystyle \int _{0}^{+\infty }\cos {(x^{2})}\,dx=\Re {\left(\int _{0}^{+\infty }e^{ix^{2}}\,dx\right)}=\Re {\left(e^{i{\frac {\pi }{4}}}{\frac {\sqrt {\pi }}{2}}\right)}={\sqrt {\frac {\pi }{8}}}\,,}
come volevasi dimostrare .
C
(
z
)
+
i
S
(
z
)
=
z
M
(
1
2
,
3
2
,
i
π
2
z
2
)
,
{\displaystyle C(z)+iS(z)=zM\left({\frac {1}{2}},{\frac {3}{2}},i{\frac {\pi }{2}}z^{2}\right),}
dove
M
{\displaystyle M}
denota una funzione ipergeometrica confluente .
La relazione con la funzione degli errori è:
C
(
z
)
+
i
S
(
z
)
=
1
+
i
2
e
r
f
[
π
2
(
1
−
i
)
z
]
.
{\displaystyle C(z)+iS(z)={\frac {1+i}{2}}\mathrm {erf} \left[{\frac {\sqrt {\pi }}{2}}(1-i)z\right].}