Simetrías de Ondas Periódicas Las simetrías permiten calcular más fácilmente los coeficientes de Fourier de una señal periódica x(t). Tipos de Simetrías: Combinaciones de estas simetrías: a) Par + Media Onda = Cuarto de Onda Par b) Impar + Media Onda = Cuarto de Onda Impar
1) Par 2) Impar 3) Media Onda
4) Escondida : Mediante desplazamientos en el eje de las ordenadas y/o de las absisas permiten transformar la señal periódica dada en otra serie que tenga una de las simetrías anteriores, de este modo calcular más fácilmente los coeficientes de Fourier para luego volver a la señal original y así obtener su serie. Antes de desarrollar las propiedades de las simetrías veamos algunas de las propiedades de las funciones pares e impares. Regla de los signos: fpar = (+) ; fimpar = (-)
fpar fpar fimpar fimpar
. fpar . fimpar . fpar . fimpar
= = = =
(+ ) . (+ ) = (+ ) ……. fpar (+ ) . ( - ) = ( - ) ……. fimpar ( - ) . (+ ) = ( - ) ……. fimpar ( - ) . ( - ) = (+ ) ……. fpar
Área de funciones pares e impares en intervalos simétricos X(t) par T /2
∫
-T/2
t
T/2 X(t) impar
−T / 2
T /2
x ( t )dt = 2. ∫ x (t )dt 0
T /2
∫ x(t )dt = 0
-T/2 T/2
−T / 2
t
Descomposición de una función f(t) en una función par más una función impar Toda función f(t) se puede descomponer en la suma de una función par más una función impar. f(t) = fp(t) + fi(t)
(1)
fpar: fp(t) = fp(-t)
fimpar: fi(t) = - fi(-t)
Si evaluamos la expresión anterior en –t, obtenemos: f(-t) = fp(-t) + fi(-t)
(2)
Pero como fp(-t) = fp(t) y fi(-t) = - fi(t) , de (1) y (2) sumando y restando las ecuaciones miembro a miembro obtenemos: Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 1 de 9
fp(t) = [ f(t) + f(-t) ] / 2
(3)
fi(t) = [ f(t) - f(-t) ] / 2
(4)
Ej: Obtener las funciones par e impar de la sig. función f(t) f(-t) 3 2
f(t) 3 2
Espejando la función obtenemos f(-t)
-1
-2
1
-1
1
2
-2 -1
2
-1
-2
-2
Ahora podemos usar las expresiones (3) y (4) anteriores y sumando y restando gráficamente obtener fi(t) y fp(t) . f(t)
f(t)
3 2
3 2
-2 -1
-1
-2
1
-1
2
1
-1
2
-2
-2
- f(-t)
f(-t) 3 2
2 1 1
2
-2 -1 1 2
-2 -1 -1
-2
-2 -3 4
2
2fi(t) = f(t) - f(-t)
2fp(t) = f(t) + f(-t)
-4 fi(t) fp(t) 1
2
-2
Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 2 de 9
Desarrollo de las propiedades de cada simetría: 1) Simetría Par: x(t) es una función par y el coseno también, por lo tanto el an queda:
an =
2 T
T /2
∫ x(t ). cos(nωot )dt
=
−T / 2
4 T
T /2
∫ x(t ). cos(nω t )dt o
0
En el caso de los bn tenemos la multiplicación de una función par ,x(t) por una función impar integrada en intervalo simétrico queda:
bn =
2 T
T /2
∫ x(t ).sen(nω t )dt o
= 0
−T / 2
De este modo el coeficiente de la serie compleja de Fourier es real puro
an − jbn 2
cn =
=
an 2
Y la serie de Fourier toma la sig. forma:
x(t ) =
ao ∞ + ∑ an cos( nω0t ) 2 n =1
+∞
o su forma exponencial
x(t ) =
an n = −∞ 2
∑
e jnωot
Notar que: 1) Una función par x(t) se descompone en una suma de funciones par (los cosenos y el ao/2) 2) En la forma exponencial por cada par de exponenciales se genera un coseno, por ejemplo se puede tomar n=3 y n=-3 y formar el coseno correspondiente. 3) an = a-n ya que el coseno es par respecto de “n”. (ver definición de an ) 2) Simetría Impar:
x(t) es una función impar y el seno también, por lo tanto el bn queda:
bn =
2 T
T /2
∫ x(t ).sen(nω t )dt o
−T / 2
=
4 T
T /2
∫ x(t ).sen(nω t )dt o
0
En el caso de los an tenemos la multiplicación de una función impar ,x(t), por una función par integrada en intervalo simétrico queda:
an =
2 T
T /2
∫ x(t ). cos(nω t )dt o
= 0
−T / 2
Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 3 de 9
De este modo el coeficiente de la serie compleja de Fourier es imaginario puro
an − jbn 2
cn =
= −j
bn 2
Y la serie de Fourier toma la sig. forma: ∞
x( t ) = ∑ bn sen( nω0t )
+∞
x(t ) =
o su forma exponencial
∑
j
n = −∞
n =1
bn 2
e jnωot
Notar que: 1) Una función impar x(t) se descompone en una suma de funciones impares (los senos) 2) En la forma exponencial por cada par de exponenciales se genera un seno. 3) bn = -b-n ya que el seno es impar respecto de “n”. (ver definición de bn ) 3) Simetría de Media Onda: En este caso x(t) cumple con lo siguiente: x(t) = - x(t + T/2) Para darnos cuenta si una función periódica tiene simetría de media onda hacemos lo sig.:
-T/2
-T/2
-T/2
T/2
T/2
T/2
3) Desplace el tramo seleccionado medio período a la izquierda o a la derecha. Si coincide con lo que está dibujado entonces tiene simetría de media onda.
2) Multiplique el tramo seleccionado por -1
1) Seleccione medio período de la señal, puede ser en cualquier lugar de x(t), pero tome uno que sea de fácil identificación.
Las señales x(t) que tienen simetría de media onda tienen an y bn distintos de cero para n impares y valen 0 para los n pares. Vamos a demostrar esto para los an ya que para los bn es similar. Partimos de la definición del an:
an =
2 T
T /2
∫ x(t ). cos(nωot )dt
−T / 2
=
2 T
0
∫ x(t ). cos(nωot )dt
+
−T / 2
2 T
T /2
∫ x(t ). cos(nω t )dt o
0
Reemplazando en la primera integral x(t) por - x(t + T/2) obtenemos: Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 4 de 9
an =
T /2
2
∫ x(t ). cos(nω t )dt o
T
=
−
−T / 2
Haciendo la sustitución
2
an = −
2 T
0
∫ x(t + T / 2). cos(nω t )dt
T
T /2
∫ x(t ). cos(nω t )dt o
0
u = t + T/2 tenemos:
T /2
T
+
o
−T / 2
2
∫ x(u ). cos[nωo (u − T / 2)]du
+
0
2 T
T /2
∫ x(t ). cos(nω t )dt o
0
0 Reemplazando el
cos[ n ω o ( u − T
teniendo en cuenta que ωo =
2π T
/ 2 )]
= cos( n ω o u ). cos n π + sen ( n ω o u ). sen ( n π )
, obtenemos la expresión definitiva para el an. 4 T
T /2
2
an = [1 − cos(nπ )] ∫ x( t ). cos(nωot )dt = T
0
4 T
T /2
bn = [1 − cos(nπ )] ∫ x(t ).sen( nωot )dt = T
T /2
∫ x(t ). cos(nω t )dt o
0
si n es impar
0
0
Análogamente: 2
y además
si
n es par
T /2
∫ x(t ).sen(nω t )dt o
si n es impar
0
0
si
n es par
3a) Simetría de cuarto de Onda Par: Se denomina así a una señal x(t) que tenga simetría de media onda y además sea par. Las expresiones para an y bn son:
an =
8 T
T /4
∫ x(t ). cos(nω t )dt o
si
n es impar ; an = 0
si
n es par
;
bn = 0
0
3b) Simetría de cuarto de Onda Impar: Se denomina así a una señal x(t) que tenga simetría de media onda y además sea impar. Las expresiones para an y bn son:
bn =
8 T
T /4
∫ x(t ).sen(nω t )dt o
si n es impar ; bn = 0 si n es par
;
an = 0
0
Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 5 de 9
Ejemplos: x2(t)
x1(t)
3
2 -1 1
t
-2 -1
1
-2
t
2
-3
x1(t) con simetría de cuarto de onda impar T/2 = 1
x2(t) con simetría de cuarto de onda par T/2 = 2
Porqué se llama simetría de cuarto de onda par o impar? Tomemos por ejemplo la señal x1(t) dibujada anteriormente, como es de cuarto de onda impar, significa que tiene simetría de media onda y es impar, por lo tanto sólo tendrá bn distintos de 0 con n impares
bn =
2 T
T /2
∫ x (t ).sen(nω t )dt o
1
,
para la señal x1(t) el período es T=2 ,
ωo = π
−T / 2
Como vemos de la definición del bn este coeficiente es 2 / T multiplicado por una integral, pero esa integral representa un área, dibujemos entonces esa área que es el producto de x1(t) con el sen(nπ t), y para hacer las cosas más fáciles graficaremos el seno para n=1 como si estuviesemos calculando el coef. b1, para los demás coeficientes el análisis es similar. x1(t)
Entonces el coeficiente b1 sería el área graficada en el último gráfico integrada entre -1 y 1. Pero si observamos detenidamente el gráfico vemos que esa área está compuesta de 4 áreas exactamente iguales
2 -1 1
t
-2
x1(t) sen(nπ t) 2
A1
sen(nπ t) 1 -1
-1 1 -1
T/4
x1(t) sen(nπ t) 2
-1
1/2
1
1
t
t
t
b 1=
2 2
T/2
1
∫ x (t ).sen(π t )dt = 1
−1
4. A = 1
8 2
1/ 2
∫ x (t ).sen(π t )dt 1
0
Es decir que toda la información que se necesita para calcular el coeficiente b1 está contenida en sólo ¼ del período. Por eso se llama simetría de cuarto de onda !
Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 6 de 9
4) Simetría escondida: Mediante desplazamientos en el eje de las ordenadas y/o de las absisas se define una nueva señal y(t) a la cuál es más fácil calcularle los coeficientes de Fourier. Una vez hecho esto se vuelve a la señal original x(t). x(t) y(t)
xo
y ( t ) = x ( t + to ) − x o t
to
t
0
Primero hacemos el desplazamiento en el eje de las ordenadas, quedando y(t) = x(t) – xo , luego la desplazamos en el eje de las absisas, obteniendo finalmente y ( t ) = x ( t + to ) − x o
Luego calculamos la serie de y(t), obteniendo:
y (t ) =
a0 ∞ + ∑ an cos( nω0t ) + bn sen( nω0t ) 2 n =1
Por último volvemos a la serie original despejando x(t), para ello, resto to en el eje ‘t’ en toda la ecuación: y (t − to) = x (t + to − to) − xo ,
quedando
x(t ) = y (t − to) + xo
Luego reemplazamos la serie calculada y ya tenemos la serie para x(t) que era lo que buscábamos.
x(t ) =
∞ a0 + xo + ∑ an cos[ nω0 (t − to)] + bn sen[nω0 (t − to)] 2 n =1
Veamos un ejemplo para que esto quede más claro: Obtener la serie de Fourier de la sig. señal x(t) Queremos que este pto, el (1,5) de x(t) valla a parar al (0,0) del gráfico de y(t), de modo de transformar esta señal x(t) en una y(t) con simetría impar.
x(t) 6 5 4 -3
1
5
9
t
De la expresión para y(t) deducida anteriormente tenemos: y (t ) = x (t + to) − xo = x( t + 1) − 5
Es decir a los valores de x(t) en el eje de las ordenadas le restamos 5 para lograr el y=0 y desplazamos en el eje de las absisas a la izquierda en 1. Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 7 de 9
y(t) 1 -1
-4
t
8
4
Luego si miramos y(t) vemos que tiene simetría de media onda y también es impar, por lo tanto tendrá bn con armónicos impares como hemos visto antes, luego: ∞
y (t ) = ∑ bn sen( nω o t )
T = 8 , ω0 =
n =1
bn =
bn =
8 T
8
T /4
∫ x(t ).sen(nω t )dt o
si n es impar ; bn = 0 si n es par
0
8/4
∫ 8 0
bn =
2π π = 8 4
π
− cos(n t ) 2 4 1 . sen (n t )dt = nπ 4 0 4
4 nπ
π
=
4 π [1 − cos(n )] si n es impar nπ 2
n impar
Reemplazando lo hallado en la ecuación de y(t) obtenemos: ∞
4 sen( nω 0t ) , n impar n =1 nπ
y (t ) = ∑
Y usando la ecuación que da x(t) tenemos: ∞
x( t ) = 5 + y (t − 1) = 5 + ∑ n =1
4 π seno [ n (t − 1)] nπ 4
, n impar
Por último podríamos reemplazar
π
π
π
π
π
sen [n (t − 1)] = cos(n ) . sen ( n t ) − sen( n ) . cos ( n t ) 4 4 4 4 4 Simplemente para verificar que la señal original x(t) al no ser ni par ni impar tiene en su serie términos con senos y cosenos ∞
x(t ) = 5 + ∑ n =1
[ −nπ4 sen(n π4 )] . cos (n π4 t ) + [ n4π cos(n π4 )] . sen (n π4 t ) an de x(t)
, n impar
bn de x(t)
ao/2 de x(t) Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 8 de 9
Apéndice: Otra forma de expresar una serie trigonométrica de Fourier Partimos de la expresión ya conocida de la serie de Fourier
x(t ) =
a0 ∞ + ∑ an cos( nω0t ) + bn sen( nω0t ) 2 n =1
an2 + bn2 obteniendo:
Multiplicamos y dividimos dentro de la sumatoria por
x(t ) =
a0 ∞ +∑ 2 n =1
an2 + bn2
[
an a +b 2 n
2 n
bn
cos( nω0t ) +
a + bn2 2 n
sen( nω 0t )
]
Si nos fijamos los cocientes que acompañan al coseno y seno se pueden representar en un triángulo rectángulo.
an
an2 + bn2
φ
θ an
a +b 2 n
bn
2 n
bn a +b 2 n
2 n
bn
= cosθ
a + bn2 2 n
an
= cos φ
a + bn2 2 n
= senθ
= senφ
Reemplazando en x(t) para las dos versiones, en función de θ o de Φ nos queda:
x(t ) =
a0 ∞ +∑ 2 n =1
an2 + bn2
[cosθ . cos(nω t ) + senθ . sen(nω t )] ............(1)
a0 ∞ +∑ 2 n =1
an2 + bn2
[senφ . cos(nω t ) + cos φ . sen(nω t )] ............(2)
0
0
o
x (t ) =
0
0
Y finalmente si reemplazamos cos( nω0 t − θ ) =
cos θ . cos( nω0t ) + senθ . sen( nω0t )
en la ecuación (1)
sen( nω0 t + φ ) =
sen φ . cos( nω0t ) + cos φ . sen(nω0t )
en la ecuación (2)
Obtenemos:
x(t ) =
∞ a0 +∑ 2 n =1
an2 + bn2
cos( nω0t − θ )
donde θ = arctg
bn an
sen( nω0 t + φ )
donde φ = arctg
an bn
o la expresión equivalente:
x (t ) =
∞ a0 +∑ 2 n =1
an2 + bn2
Autor: Ing. Jorge Calcagno - JTP Matemática Avanzada Facultad de Ingeniería – Universidad Nacional de Mar del Plata - Argentina
Página 9 de 9