Trabajo Practico Fourier

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 12

Matemática C – 1er semestre 2022 - Trabajo Práctico

Alumno=Francisco Lerra

Dni=44339490

Numero de alumno =74163/5

SERIES DE FOURIER

1_Jean-Baptiste-Joseph Fourier

(Auxerre, Francia, 1768 - París, 1830) Ingeniero y matemático francés. Era hijo de un sastre, y
fue educado por los benedictinos. Los puestos en el cuerpo científico del ejército estaban
reservados para familias de estatus reconocido, así que aceptó una cátedra militar de
matemáticas. Las principales aportaciones de Fourier a la Física y las Matemáticas fueron el
desarrollo de la teoría del calor, ecuaciones diferenciales, ecuaciones algebraicas, series
trigonométricas, estadística matemática y teoría de las probabilidades; siendo su obra
monumental la Théorie analytique de la chaleur; en la que desarrolla ecuaciones para explicar
la conducción térmica en metales.

Fuentes:

https://www.biografiasyvidas.com/

https://principia.io/

2_Una función es periódica cuando su valor se repite cada vez que la variable independiente
recorre un cierto intervalo. El valor de este intervalo se llama periodo. Se llama período un
lapso de tiempo determinado en el cual se desenvuelve una acción, un fenómeno o una
sucesión de eventos concatenados.

f(x)=Sen(x)
3_Una función es par si f(-x) = f(x) para toda x en el dominio de f. La gráfica de una función par
es simétrica con respecto al eje y.

F(x)=cos(x)

4_Una función es impar si f(-x) = -f(x) para toda x en el dominio de f. La gráfica de una función
impar es simétrica con respecto al origen.

F(x)=sen(x)
5_

a) El producto de dos funciones pares es una función par.


Sea F y G funciones pares
F(x)=F(-x) G(x)=G(-x)
H(x)=F(x).G(x)
H(-x)=F(-x).G(-x)
H(x)=H(-x)

Ejemplo =
f(x)=x2 g(x)= cos(x)

f(x) X g(x) =x2.Cos(x) =h(x)


b) El producto de dos funciones impares es una función par.
Sea F y G funciones impares

F(-x)=-F(x) G(-x)=-G(X)

H(-x)=F(-x).G(-x)

H(x)=-F(x).(-G(x))

H(x)=H(-x) (función par)

Ejemplo =

f(x)=x3 g(x)=sen(x)

f(x) X g(x)=x3.sen(x)=h(x)

c) El producto de una función par y una función impar es una función impar.

Sea F (x) una función par y G (x) una función impar

F(x)=F(-x) G(-x)=-G(X)

H(x)=F(x).(-G(x))

H(-x)=F(-x).G(-x)

H(-x)=-H(x)(función impar
f(x)=x2 g(x)=sen(x)

f(x) X g(x)=x2.sen(x)=h(x)

d) Sea 𝑓 integrable en el intervalo cerrado [–L, L].


𝐿
Si 𝑓 es una función impar, entonces:∫−𝐿 𝑓(𝑥)𝑑𝑥 = 0
Supongamos que f(x) es un impar, se cumple que f(-x)=-f(x)
𝐿 0 𝐿
∫ 𝑓(𝑥) 𝑑𝑥 = ∫ 𝑓(𝑥) 𝑑𝑥 + ∫ 𝑓(𝑥) 𝑑𝑥
−𝐿 −𝐿 0

Con -x=u , -dx=du, planteamos lo siguiente


0 0 𝐿 𝐿 𝐿
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓(𝑥)𝑑𝑥 + ∫ 𝑓 (𝑥)𝑑𝑥 = − ∫ 𝑓 (𝑢)𝑑𝑢 + ∫ 𝑓 (𝑥)𝑑𝑥 = 0
−𝐿 −𝑙 0 0 0

Ejemplo
π
Sea 𝑓(x)=sen(x) ∫−π 𝑠𝑒𝑛(𝑥) 𝑑𝑥 = 0

e) Sea 𝑓 integrable en el intervalo cerrado [–L, L].


𝐿 𝐿
Si 𝑓 es una función par, entonces: ∫−𝐿 𝑓 (𝑥)𝑑𝑥 = 2 ∫0 𝑓(𝑥)𝑑𝑥

𝐿 0 𝐿
∫ 𝑓 (𝑥)𝑑𝑥 = ∫ 𝑓(𝑥)𝑑𝑥 + ∫ 𝑓 (𝑥)𝑑𝑥
−𝐿 −𝐿 0

Con x=-u dx=-du


0 0 𝐿 𝐿
∫ 𝑓 (𝑥)𝑑𝑥 = − ∫ 𝑓 (−𝑢)𝑑𝑢 = ∫ 𝑓 (−𝑢)𝑑𝑢 = ∫ 𝑓 (𝑢)𝑑𝑢
−𝐿 𝐿 0 0

Replazando
𝐿 0 𝐿 𝐿 𝐿 𝐿
∫ 𝑓 (𝑥)𝑑𝑥 = ∫ 𝑓 (𝑥)𝑑𝑥 + ∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓(𝑢)𝑑𝑢 + ∫ 𝑓(𝑥)𝑑𝑥 = 2 ∫ 𝑓 (𝑥)𝑑𝑥
−𝐿 −𝐿 0 0 0 0
Ejemplo

π/2 𝜋/2 𝜋/2


Sea 𝑓(x)=cos(x) ∫−π/2 𝑓 (𝑥)𝑑𝑥 = ∫−𝜋/2 cos(x) dx = 2 ∫0 cos(𝑥)𝑑𝑥 = 2

6_

La función no tiene paridad ya que g(-x)≠g(x) y g(-x)≠-g(x)


Serie de Fourier:
∞ ∞
𝑎0 𝑛𝜋𝑥 𝑛𝜋𝑥
𝑓 (𝑥) = + ∑ 𝑎𝑛 cos ( ) + ∑ 𝑏𝑛 𝑠𝑒𝑛( )
2 𝐿 𝐿
𝑛=1 𝑛=1

1 𝐿
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥
𝐿 −𝐿

1 𝐿 𝑛𝜋𝑥
𝑏𝑛 = ∫ 𝑓 (𝑥)𝑠𝑒𝑛( )𝑑𝑥
𝐿 −𝐿 𝐿

1 𝐿 𝑛𝜋𝑥
𝑎𝑛 = ∫ 𝑓 (𝑥) cos ( ) 𝑑𝑥
𝐿 −𝐿 𝐿

F(x)=1-4x en [-1,1]

1
𝑎0 = 1 ∫ 1 − 4𝑥 𝑑𝑥 = 2
−1

1
𝑛𝜋𝑥
𝑎𝑛 = 1 ∫ 1 − 4𝑥 cos ( ) 𝑑𝑥 = 0
−1 1
1
𝑛𝜋𝑥 8(−1)𝑛
𝑏𝑛 = 1 ∫ 1 − 4𝑥 𝑠𝑒𝑛 ( ) 𝑑𝑥 =
−1 1 𝜋𝑛
8(−1)𝑛
𝑓 (𝑥 ) = 1 + ∑ ∞
𝑛=1 𝜋𝑛
𝑠𝑒𝑛(𝑛𝜋𝑥)}

8
𝑆1 (𝑥) = 1 − 𝑠𝑒𝑛(𝜋𝑥)
𝜋
8 4
𝑆2 (𝑥) = 1 − 𝑠𝑒𝑛(𝜋𝑥) + 𝑠𝑒𝑛(2𝜋𝑥)
𝜋 𝜋

8 4 8
𝑆3 (𝑥) = 1 − 𝑠𝑒𝑛(𝜋𝑥) + 𝑠𝑒𝑛(2𝜋𝑥 ) − 𝑠𝑒𝑛(3𝜋𝑥)
𝜋 𝜋 3𝜋

A la hora de resolver las integrales use Symbolab


Primer suma

Segunda suma
Tercer suma
Derivas parciales

1_ La ecuación de difusión es
∂U(x, t) ∂2 U(x, t)
( )−α( )=0
∂t ∂x 2

Esta ecuación, de segundo orden en x y primer orden en t, permite determinar, por


ejemplo, la temperatura U(x, t) de una barra conductora de calor de longitud L en función
de la posición x y el tiempo t, si se conocen la distribución inicial de temperatura U(x, 0) y
las temperaturas en los bordes U(0, t) y U(L, t) (se asume que la barra esta térmicamente
aislada salvo en los extremos). En este caso α es una constante positiva que depende de la
conductividad térmica, el calor especifico y la densidad del material.

2𝑡
𝑈(𝑥, 𝑡) = 𝐴𝑠𝑒𝑛(𝜆𝑥 + 𝜑)𝑒 −𝛼𝜆
∂U(x, t) ∂2 U(x, t)
( ) = α( )
∂t ∂x 2

2 2
∂(𝐴𝑠𝑒𝑛(𝜆𝑥 + 𝜑)𝑒 −𝛼𝜆 𝑡 ) ∂2 (𝐴𝑠𝑒𝑛(𝜆𝑥 + 𝜑)𝑒 −𝛼𝜆 𝑡 )
( ) = α( )
∂t ∂x 2
2 𝛼𝑡 2 𝛼𝑡
−𝐴𝜆2 𝑒 −𝜆 𝛼𝑠𝑒𝑛(𝜆𝑥 + 𝜑)= −𝐴𝜆2 𝑒 −𝜆 𝛼𝑠𝑒𝑛(𝜆𝑥 + 𝜑)

Se puede ver que para cualquier valor de A y 𝜑 contantes,


2
𝑈(𝑥, 𝑡) = 𝐴𝑠𝑒𝑛(𝜆𝑥 + 𝜑)𝑒 −𝛼𝜆 𝑡
Es solución de la ecuación de difusión.

2_
El método de separación de variables consiste en proponer; 𝑈(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
X(x) depende solamente de x para 0 < x < 1
T(t) depende solamente de t para t > 0
∂U(x,t) ∂2U(x,t)
( )=( ) 𝛼=1
∂t ∂x2
𝑑𝑈(𝑥,𝑡)
𝑑𝑡
=X(x)T'(t)

𝑑2 𝑈(𝑥,𝑡)
𝑑2𝑥
= 𝑋′′(𝑥)T(t)

X(x)T'(t)= X''(x)T(t)
T′(t) X′′(x)
=
T(t) X(x)

Como el primer cociente depende solo de x y el segundo solo de t, siendo x y t variables


Independientes, la única posibilidad es que ambos cocientes sean constantes:

T′(t) X′′(x)
= = −𝑘 2
T(t) X(x)

La constante k2 se denomina constante de separación. Se obtienen as ́ı dos ecuaciones

Diferenciales ordinarias, una de segundo orden para X y una de primer orden para T:
X''(x)+k2X(x)=0
T'(t)+k2T(t)=0
Si k=0
X(x)=A+Bx , T(t)=C
Solución estacionaria = X(x)T(t)=A+Bx (independiente de t )

Si k ≠0
X(x)=Acos(kx)+Bsen(kx)
2𝑡
T(t)=C𝑒 −𝑘
2𝑡
Soluciones no estacionarias = X(t)T(t)=(Acos(kx)+Bsen(kx))𝑒 −𝑘
Las condiciones de contorno implican ahora ,
U(0,t)=X(0)T(t)=0 , U(L,t)=X(L)T(t)=0

∀ t > 0. Por lo tanto, si T(t) ≠ 0, debe cumplirse X(0) = 0, X(L) = 0.


2
Para k =0, esto implica A = B = 0,ya que , X(x)=Acos(kx)+Bsen(kx) ,T(t)=C𝑒 −𝑘 𝑡 . El valor
k = 0 conduce en este caso solo a la solución trivial. En cambio, para k ≠ 0,implica A = 0,
2
X(t)T(t)=(Acos(kx)+Bsen(kx))𝑒 −𝑘 𝑡 , y entonces
B sen(kL) = 0
Con B≠0 esto conduce a Kl=nπ, con n entero

K=nπ/l , con n =1,2,3,……

Por la propiedad de superposición proponemos la solución general ,


𝑛𝜋 2
−( ) 𝑡 𝑛𝜋𝑥
∑∞
𝑛=1 𝐵𝑛 𝑒 𝐿 𝑠𝑒𝑛( )
𝐿

Con
1 𝐿 𝑛𝜋𝑥
𝐵𝑛 = 𝐿 ∫−𝐿 𝑓 (𝑥)𝑠𝑒𝑛( 𝐿
) 𝑑𝑥 n=1,2,3,……

1
Bn=1 ∫−1((2𝑠𝑒𝑛(2𝜋𝑥) − 3𝑠𝑒𝑛(4𝜋𝑥)) 𝑠𝑒𝑛(𝑛𝜋𝑥))𝑑𝑥

Llegamos a que Bn es igual a la resta de dos integrales


1 1
Bn=2 ∫−1 𝑠𝑒𝑛(2𝜋𝑥 )𝑠𝑒𝑛(𝑛𝜋𝑥)𝑑𝑥 -3 ∫−1 𝑠𝑒𝑛(4𝜋𝑥 )𝑠𝑒𝑛( 𝑛𝜋𝑥)𝑑𝑥

Como tenemos un conjunto de funciones ortogonales usaremos la siguiente propiedad ,

Para la primer integral tenemos que para n=2 la integral vale 1 y para n≠2 vale 0
Para la segunda integral tenemos que para n=4 vale 1y para n≠4 vale 0
𝑛𝜋 2
−( ) 𝑡 𝑛𝜋𝑥
La solucion de la ecuación quedaría U(x,t)=∑∞
𝑛=1 𝐵𝑛 𝑒 𝐿 𝑠𝑒𝑛( 𝑙
)

2 2
U(x,t)=2𝑒 −4𝜋 𝑡 𝑠𝑒𝑛(2𝜋𝑥) − 3𝑒 −16𝜋 𝑡 𝑠𝑒𝑛(4𝜋𝑥)

2
El termino dominante de la suma es el primero, 2𝑒 −4𝜋 𝑡 𝑠𝑒𝑛(2𝜋𝑥 )

También podría gustarte