Series de Fourier
Series de Fourier
Series de Fourier
CALDAS
CÁLCULO INTEGRAL
Hayden Julian Bejarano Gonzalez
Johan Sebastian Vargas Agudelo
Febrero 2019
Demuestre que:
Proof.
∞
4 X sen(2k + 1)x sen3x sen5x
f (x) ≈ = senx + + + ······
n 2k + 1 3 5
k=0
.
Demostración.
1 0 1 π
Z Z
a0 = −1dx + 1dx
π −π π 0
Z 0 Z π
1
= − dx + dx
π −π 0
1
= (−x|π0 + x|π0 ) = 0
π
1
ii.Ahora calculesmos el coeficiente an
Z 0 Z π
1 1
an = −cos(nx)dx + cos(nx)dx
π −π π 0
11 0 1 1 π
=− sen(nx)|−π + sen(nx)|0
πn π n
1 1 1 1 1 1 1
− sen(0) − sen(−πn) + sen(nπ) − − sen(0) = 0
π n n π n n n
∴ an = 0
1 0 1 π
Z Z
bn = − sen(nx)dx + sen(nx)dx
π −π π 0
1 1 0 1 1 π
=− − cos(nx)|−π + − cos(nx)|0
π n π n
1 1 1 1
=− − (cos(0) − cos(−nπ)) + − (cos(nπ) − cos(0))
π n π n
1 1 1 1
=− − (1 − cos(nπ) + − (cos(nπ) − 1)
π n π n
1 1 cos(nπ 1 1 1
=− − + + − cos(nπ) +
π n n π n n
2 − 2cos(nπ)
∴ bn =
nπ
Ahora sabemos que la forma general para la serie de Fourier esta dada por:
∞
a0 X
f (x) ≈ + [an cos(nx) + bn sen(nx)]
2 n=1
∞
X 2 − 2cos(nπ)
f (x) ≈ sen(nx)
n=1
nπ
2
sabemos que cos(nπ) = (−1)n ).
∞
2 X sen(nx)(1 − (−1)n )
f (x) ≈
π n=1 n
∞
4 X sen(2n + 1)
f (x) ≈
π 0 2n + 1
Ahora miraremos mediante aproximaciones que esta serie se acerca mucho a la funcion original.
1. f (x)
2. S1
3
3. S2
4. S4
5. S10
4
6. S15
7. S20
8. S30
5
9. S50
10. S75
11. S100
De esta manera podemos ver como la Serie de Fourier se empieza a aproximar a la función original de una
manera muy exacta.
6
1 1
2x, si − 2 ≤x≤ 2
2.Si f (x) = −2(x − 1), si 12 ≤ x ≤ 1 es una función periódica de periodo 2.
−2(x + 1), si − 1 ≤ x ≤ − 21
Proof. Demostración.
Z − 12 Z 1
2
Z 1
a0 = (−2(x + 1))(cos(0πx))dx + (2x)(cos(0πx))dx + (−2(x − 1))(cos(0πx))dx
−1 − 21 1
2
Z − 12 Z 1
2
Z 1
= −2 (x + 1)(cos(0))dx + 2 x(cos(0))dx − 2 (x − 1)(cos(0))dx
−1 − 12 1
2
Z − 21 Z 1
2
Z 1
= −2 (x + 1)dx + 2 xdx − 2 (x − 1)dx
−1 − 12 1
2
x2
2
2
− 12 x 12 x 1
= −2 + x|−1 + 2 | 1 −2 − x| 1
2 2 −2 2 2
1 1 1 1 1 1 1 1
= −2 − − −1 +2 − −2 −1 − −
8 2 2 8 8 2 8 2
1 1 1 1
= −2 + 2 (0) − 2 − =− +0+ =0
8 8 4 4
Z − 12 Z 1
2
Z 1
an = (−2(x + 1))(cos(nπx))dx + (2x)(cos(nπx))dx + (−2(x − 1))(cos(nπx))dx
−1 − 12 1
2
Z − 21 Z 1
2
Z 1
= −2 (x + 1)(cos(nπx))dx + 2 x(cos(nπx))dx − 2 (x − 1)(cos(nπx))dx
−1 − 21 1
2
7
Calculamos:
Z − 21
(x + 1)(cos(nπx))dx
−1
DERIVADA INTEGRAL
x+1 cos(nπx)
1
1 nπ (sen(nπx))
0 − n21π2 (cos(nπx))
x+1 1 − 12
= (sen(nπx)) + 2 2 (cos(nπx)) |−1
nπ n π
nπ
1 nπ 1 1
= sen − + 2 2 cos − 0 + 2 2 cos(−nπ)
2nπ 2 n π 2 n π
1 nπ 1 nπ 1
=− sen + 2 2 cos − − 2 2 (cos(nπ))
2nπ 2 n π 2 n π
Ahora:
Z 1
2
x(cos(nπx))dx
− 21
DERIVADA INTEGRAL
x cos(nπx)
1
1 nπ (sen(nπx))
0 − n21π2 (cos(nπx))
x 1 1
= (sen(nπx)) + 2 2 (cos(nπx)) |− 2
1
nπ n π 2
nπ
1 nπ 1 nπ 1 nπ 1
= sen + 2 2 cos − − sen − + 2 2 cos −
2nπ 2 n π 2 2nπ 2 n π 2
1 nπ 1 nπ 1 nπ 1 nπ
= sen + 2 2 cos −( sen − 2 2 cos =0
2nπ 2 n π 2 2nπ 2 n π 2
8
Ahora: Z 1
(x − 1)(cos(nπx))dx
1
2
DERIVADA INTEGRAL
x−1 cos(nπx)
1
1 nπ (sen(nπx))
0 − n21π2 (cos(nπx))
x−1 1 1
= (sen(nπx)) + 2 2 (cos(nπx)) | 1
nπ n π 2
nπ
1 1 nπ 1
= 0 + 2 2 (cos(nπ)) − − sen + 2 2 cos
n π 2nπ 2 n π 2
1 1 nπ 1 nπ
= 2 2 (cos(nπ)) + sen − 2 2 cos
n π 2nπ 2 n π 2
Ası́:
1 nπ 1 nπ 1 1 1 nπ
an = −2 − sen + 2 2 cos − 2 2 (cos(nπ)) +2(0)−2 (cos(nπ)) + sen −
2nπ 2 n π 2 n π n2 π 2 2nπ 2 n
1 nπ 2 nπ 2 2 1 nπ 2 nπ
= sen − 2 2 cos + 2 2 (cos(nπ))− 2 2 (cos(nπ))− sen + 2 2 cos =0
nπ 2 n π 2 n π n π nπ 2 n π 2
∴ an = 0
Z − 21 Z 1
2
Z 1
bn = (−2(x + 1))(sen(nπx))dx + (2x)(sen(nπx))dx + (−2(x − 1))(sen(nπx))dx
−1 − 12 1
2
Z − 12 Z 1
2
Z 1
= −2 (x + 1)(sen(nπx))dx + 2 x(sen(nπx))dx − 2 (x − 1)(sen(nπx))dx
−1 − 21 1
2
9
Calculamos:
Z − 21
(x + 1)(sen(nπx))dx
−1
DERIVADA INTEGRAL
x+1 sen(nπx)
1
1 − nπ (cos(nπx))
0 − n21π2 (sen(nπx))
x+1 1 − 12
= − (cos (nπx)) + 2 2 (sen (nπx)) |−1
nπ n π
nπ
1 nπ 1
= − cos − + 2 2 sen − − (0 + 0)
2nπ 2 n π 2
1 nπ 1 nπ
=− cos − + 2 2 sen −
2nπ 2 n π 2
1 nπ 1 nπ
=− cos − 2 2 sen
2nπ 2 n π 2
Ahora: Z 1
2
x(sen(nπx))dx
− 12
DERIVADA INTEGRAL
x sen(nπx)
1
1 − nπ (cos(nπx))
0 − n21π2 (sen(nπx))
x 1 1
= − (cos (nπx)) + 2 2 (sen (nπx)) |− 1 2
nπ n π 2
nπ
1 nπ 1 nπ 1 nπ 1
= − cos + 2 2 sen − cos − + 2 2 sen −
2nπ 2 n π 2 2nπ 2 n π 2
nπ
1 nπ 1 nπ 1 nπ 1
= − cos + 2 2 sen − cos − 2 2 sen
2nπ 2 n π 2 2nπ 2 n π 2
1 nπ 1 nπ 1 nπ 1 nπ
=− cos + 2 2 sen − cos + 2 2 sen
2nπ 2 n π 2 2nπ 2 n π 2
2 nπ 2 nπ 1 nπ 2 nπ
=− cos + 2 2 sen =− cos + 2 2 sen
2nπ 2 n π 2 nπ 2 n π 2
10
Ahora: Z 1
(x − 1)(sen(nπx))dx
1
2
DERIVADA INTEGRAL
x−1 sen(nπx)
1
1 − nπ (cos(nπx))
0 − n21π2 (sen(nπx))
x−1 1 1
= − (cos (nπx)) + 2 2 (sen (nπx)) | 1
nπ n π 2
nπ
1 nπ 1
= (0 + 0) − cos + 2 2 sen
2nπ 2 n π 2
1 nπ 1 nπ
=− cos − 2 2 sen
2nπ 2 n π 2
Ası́:
nπ nπ
1 nπ 1 1 nπ 2
bn = −2 − cos − 2 2 sen +2 − cos + 2 2 sen ···
2nπ 2 n π 2 nπ 2 n π 2
nπ
1 nπ 1
··· − 2 − cos − 2 2 sen
2nπ 2 n π 2
nπ 1
1 nπ 2 nπ 2 nπ 4 nπ
= cos + 2 2 sen + − cos + 2 2 sen + cos ···
nπ 2 n π 2 nπ 2 n π 2 nπ 2
2 nπ
· · · + 2 2 sen
n π 2
1 nπ 2 nπ 2 nπ 4 nπ 1 nπ 2 nπ
= cos + 2 2 sen − cos + 2 2 sen + cos + 2 2 sen
nπ 2 n π 2 nπ 2 n π 2 nπ 2 n π 2
8 nπ 8 nπ
= 0 + 2 2 sen = 2 2 sen
n π 2 n π 2
Ahora sabemos que la forma general para la serie de Fourier esta dada por:
∞
a0 X
f (x) ≈ + [an cos(nx) + bn sen(nx)]
2 n=1
∞
X 8 nπ
f (x) ≈ sen sen(nx)
n=1
n2 π 2 2
11
Ahora miraremos mediante aproximaciones que esta serie se acerca mucho a la funcion original.
1. f (x)
2. S1
3. S3
12
4. S5
5. S10
6. S20
13
7. S50
8. S100
De esta manera podemos ver como la Serie de Fourier se empieza a aproximar a la función original de una
manera muy exacta.
3.La base de un tetraedro es un triángulo equilátero de 3cm de lado. El cuarto vértice está directamente
encima de uno de los vértices de la base, de modo que toda sección recta paralela a la base es un triángulo
equilátero. Si la altura del tetraedro es de 2cm, hallar su volumen.
Proof. Para resolver esta ejercicio se hará primero una generalización en el caso de que los lados de los
triángulo sean de lados de lado l y altura h.
Ahora trabajando con una de las caras del tetraedro y hallando la ecuación de la recta L con puntos A = (0, 0)
y B = (h, 2l ), si 0 ≤ x ≤ h, como se observa en la siguiente figura.
14
Sea la pendiente:
l
2 l lx
m= = ⇒y=
h 2h 2h
l
Proof. Sea 2 la base de medio triangulo, ahora por teorema de pitagoras:
2
l l2
l2 = + h2 ⇒ h2 = l2 −
2 4
√
3l
⇒h=
2
√
3l2
A4 =
4
15
Dado que los lados l del triángulo depende del valor de la función hallada en (1), por lo tanto sabemos
que el volumen de la sección transversal es:
Z h
√ 2 !
3 l
V = x dx
0 4 h
√ h
√
3l2 3l2
Z
1 3h
V = x2 dx = h |
4h2 0 4h2 3 0
3 2
∴V = l h
12
√ √
3 2 3 3
⇒V = (3) (2) =
12 2
4.Calcule mediante el metodo de Simpson o el metodo de los trapecios la siguiente integral definida.
Z 1
2
ex dx
0
16
1. MÉTODO DE LOS TRAPECIOS
El método del trapecio es un método de integración numérica para integrales definidas las cuales no cuentan
con una primitiva, por lo que se hace una aproximación del área bajo la curva mediante la construcción de
trapecios si subdividimos el intervalo [a, b] en subintervalos de la misma longitud h donde:
b−a
h=
n
Z b Z x1 Z x2 Z b
f (x)dx = f (x)dx + f (x)dx + · · · + f (x)dx
a a x1 xn−1
Z b
f (a) + f (x1 ) f (xn−1 ) + f (b)
f (x)dx ≈ (x1 − a) + · · · + (b − xn−1 )
a 2 2
Z b
h
[f (a) + 2f (x1 ) + · · · + 2f (xn ) + f (b))]
f (x)dx ≈
a 2
Z b " Pn−1 #
f (a) + 2 i=1 f (xi ) + f (b)
f (x)dx ≈ (b − a)
a 2n
Ahora aplicándolo a: Z 1
2
ex dx = 1, 46
0
17
Z 1
x2 f (0) + f (1) 1+e
e dx ≈ (1 − 0) = = 1.85914
0 2 2
18
Z 1
x2 f (0) + 2f (0.2) + 2f (0.4) + 2f (0.6) + 2f (0.8) + f (1)
e dx ≈ (1 − 0)
0 2(5)
" 2 2 2 2
#
1 + e(0.2) + e(0.4) + e(0.6) + e(0.8) + e
= = 1, 48065
10
De esta manera podemos ver que escogiendo una particion mas fina se va acercando al valor del área bajo
la curva.
1. MÉTODO DE SIMPSON
Este método es muy parecido al de los trapecios pero se tiene en cuenta el punto medio de cada subintervalo
que se tome y esta definido por la siguiente formula:
" Pn Pn−1 #
Z b
f (a) + 4 i=1 f (xmi ) + 2 i=1 f (xi ) + f (b)
f (x)dx ≈ (b − a)
a 6n
Z 1
x2 f (0) + 4 [f (0.1) + f (0.3) + f (0.5) + f (0.7) + f (0.9)] + 2 [f (0.2) + f (0.4) + f (0.6) + f (0.8)] + f (1)
e dx ≈ (1−0)
0 6(5)
1 h h 2 2 2 2 2
i h 2 2 2 2
i i
= 1 + 4 e(0.1) + e(0.3) + e(0.5) + e(0.7) + e(0.9) + 2 e(0.2) + e(0.4) + e(0.6) + e(0.8) + e = 1, 4626
30
Esta es una aproximación muy cercana al valor real de esta área.Por lo que podemos apreciar que los métodos
antes mencionados son muy eficaces para la resolución de integrales definidas que no tengan primitiva.
19