Apunts d’Anàlisi Real
Anà lisi real - Departament de Matemà tica Aplicada IV - UPC
Anà lisi real - Departament de Matemà tica Aplicada IV - UPC
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
4.5 El fenomen de Gibbs 95<br />
S’observa la presència de pics a prop de les discontinuïtats de la funció, la magnitud dels<br />
quals no disminieix amb N, encara que s’acosten a les discontinuïtats al crèixer N. Per a cada<br />
punt hi ha convergència puntual, però per a una N donada sempre hi ha punts a una distància<br />
fixada del valor cap al qual tendeixen eventualment. Això és el fenomen de Gibbs.<br />
Anem a calcular el valor del màxim (x ∗ , S N (x ∗ )) més proper a zero, i veurem que no decreix<br />
quan aumenta N. Tenim<br />
S ′ N(x) = 4 π<br />
N∑<br />
cos(2k + 1)x = 2 π<br />
k=0<br />
sin(2N + 2)x<br />
.<br />
sinx<br />
Cal observar que S<br />
N ′ (x) > 0 en un entorn prou petit (que disminueix amb N) de zero. Imposant<br />
S<br />
N ′ (x) = 0 tindrem (2N + 2)x = kπ, k ∈ ¤,<br />
de manera que el primer màxim relatiu a la dreta del zero (és màxim per la positivitat de la<br />
derivada en un entorn de zero) és<br />
x ∗ π<br />
=<br />
2N + 2 .<br />
Llavors<br />
S N (x ∗ ) = 4 π<br />
N∑<br />
k=0<br />
1 (2k + 1)π<br />
sin<br />
2k + 1 2N + 2 .<br />
Això es pot interpretar com una suma de Riemann corresponent a la integral de la funció<br />
sinc(x) = sin x<br />
x<br />
. En efecte, dividint [0, π] en N + 1 intervals<br />
[n π π<br />
, (n + 1)<br />
N + 1 N + 1 ],<br />
n = 0, . . .,N,<br />
i escollint el punt central de cada interval per fer la suma. Efectivament<br />
on<br />
Per tant<br />
S N (x ∗ ) =<br />
4 π<br />
= 2 π<br />
x k =<br />
N∑<br />
k=0<br />
N∑<br />
k=0<br />
sin (2k+1)π<br />
2N+2<br />
(2k+1)π<br />
2N+2<br />
sin (2k+1)π<br />
2N+2<br />
(2k+1)π<br />
2N+2<br />
} {{ }<br />
sin x k<br />
x k<br />
(2k + 1)π<br />
2N + 2 .<br />
∫ π<br />
π<br />
2N + 2 =<br />
π<br />
N + 1<br />
} {{ }<br />
∆x k<br />
lim S N(x ∗ ) =<br />
2 siny<br />
dy =<br />
N→∞ π 0 y<br />
( ∫ 1 π<br />
siny<br />
= 1 + 2 dy − 1 )<br />
π 0 y 2<br />
∼ 1 + 2 × 0.0895,<br />
on hem utilitzat<br />
∫ π<br />
0<br />
siny<br />
y<br />
dy ∼ 1.852.<br />
Carles Batlle i Enric Fossas 2002