9.3: Serie de Fourier
( \newcommand{\kernel}{\mathrm{null}\,}\)
Nuestra solución de las ecuaciones de difusión y onda requerirá el uso de una serie de Fourier. Una función periódicaf(x) con punto2L, se puede representar como una serie de Fourier en la forma f(x)=a02+∞∑n=1(ancosnπxL+bnsinnπxL).
Determinación de los coeficientesa0,a1,a2,… yb1,b2,b3,… hace uso de relaciones de ortogonalidad para seno y coseno. Primero definimos el delta de Kronecker ampliamente utilizadoδnm comoδnm={1if n=m;0otherwise.
Las relaciones de ortogonalidad paran y los enterosm positivos se dan entonces con notación compacta como las fórmulas de integración ∫L−Lcos(mπxL)cos(nπxL)dx=Lδnm,∫L−Lsin(mπxL)sin(nπxL)dx=Lδnm,∫L−Lcos(mπxL)sin(nπxL)dx=0.
Para ilustrar la técnica de integración utilizada para obtener estos resultados, derivamos(???) asumiendo quen ym son enteros positivos conn≠m. Cambiando variables aξ=πx/L, obtenemos∫L−Lsin(mπxL)sin(nπxL)dx=Lπ∫π−πsin(mξ)sin(nξ)dξ=L2π∫π−π[cos((m−n)ξ)−cos((m+n)ξ)]dξ=L2π[1m−nsin((m−n)ξ)−1m+nsin((m+n)ξ)]π−π=0.
Puesm=n, sin embargo,∫L−Lsin2(nπxL)dx=Lπ∫π−πsin2(nξ)dξ=L2π∫π−π(1−cos(2nξ))dξ=L2π[ξ−12nsin2nξ]π−π=L.
Las fórmulas de integración dadas por(???) y(???) pueden derivarse de manera similar.
Para determinar el coeficientean, multiplicamos ambos lados de(???) porcos(nπx/L) conn un entero no negativo, y cambiamos la variable de suma ficticio den am. Integrando sobrex de−L aL y asumiendo que la integración se puede hacer término por término en la suma infinita, obtenemos∫L−Lf(x)cosnπxLdx=a02∫L−LcosnπxLdx+∞∑m=1{am∫L−LcosnπxLcosmπxLdx+bm∫L−LcosnπxLsinmπxLdx}.
Sin=0, entonces la segunda y tercera integrales en el lado derecho son cero y la primera integral es2L para que el lado derecho se conviertaLa0. Sin es un entero positivo, entonces la primera y tercera integrales en el lado derecho son cero, y la segunda integral esLδnm. Para este caso, tenemos∫L−Lf(x)cosnπxLdx=∞∑m=1Lamδnm=Lan, donde todos los términos en la suma exceptom=n son cero en virtud del delta del Kronecker. Por lo tanto, obtenemos paran=0,1,2,… an=1L∫L−Lf(x)cosnπxLdx.
Para determinar los coeficientesb1,b2,b3,…, multiplicamos ambos lados de(???) porsin(nπx/L), conn un entero positivo, y nuevamente cambiamos la variable de suma ficticio den am. Integrando, obtenemos∫L−Lf(x)sinnπxLdx=a02∫L−LsinnπxLdx+∞∑m=1{am∫L−LsinnπxLcosmπxLdx+bm∫L−LsinnπxLsinmπxLdx}.
Aquí, la primera y segunda integrales en el lado derecho son cero, y la tercera integral esLδnm para que∫L−Lf(x)sinnπxLdx=∞∑m=1Lbmδnm=Lbn.
Por lo tanto, paran=1,2,3,…, bn=1L∫L−Lf(x)sinnπxLdx.
Nuestros resultados para la serie de Fourier de una funciónf(x) con periodo2L son dados por(???),(???) y(???).