Ejercicio 1: Filtro Analógico

Cada estudiante de manera individual debe resolver el siguiente ejercicio. Un filtro

analógico se encuentra representado por la siguiente función de transferencia:

60∗a
2
H ( s )= 2
s + a∗103.7 s+ 40000

La constante “a” y la constante “b”, serán asignadas en el foro por su tutor al comienzo de
la actividad.

a=12

Determine:

a) Los valores de la frecuencia natural del sistema w 0 y el coeficiente de

amortiguamiento δ .

La función de transferencia de un filtro pasa-bajo de segundo orden es:

k s2
H (s)=
w
s 2 + 0 + w02
Q

Haciendo uso de la constante, se tiene que la función de transferencia es:

360
H (s)= 2
s +1244.4 s +40000

La frecuencia natural del sistema w 0 se encuentra de la siguiente manera:

w 02=40000
w 0=√ 40000
w 0=200 rad / s

Para encontrar el coeficiente de amortiguamiento δ , se debe hallar primero el

factor de calidad dado que:

1
2δ=
Q
 w0
=1244.4
Q
w0
Q=
1244.4

200
Q=
1244.4

Q=0.160

Al conocer el factor de calidad ya se puede hallar el coeficiente de

amortiguamiento δ .

1
2δ=
Q

1
δ=
2Q

1
δ=
2 ( 0.160 )

δ =3.125

b) La ganancia máxima del filtro en dB.

Dada la función de transferencia

360
H (s)= 2
s +1244.4 s +40000

Se debe hacer el uso del teorema del valor final, para esto se debe evaluar H ( s)
en cero, al inicio del filtro antes de hacer cualquier atenuación:

360
H ( 0 )=G=
40000

Ahora, se debe encontrar la ganancia en dB, haciendo uso de:

GdB=20∗log 10 G

360
G dB=20∗log 10 ( 40000 )
 G dB=−40.915 dB

c) La frecuencia en la cual se presenta la ganancia del sistema.

Para esto se debe encontrar la frecuencia de corte f 0.

Se entiende que:

w 0=200 rad / s

Además, se tiene la relación que w 0=2 π f 0

w0
f 0=
2π

200
f 0=
2π

f 0=31.83 Hz

d) El ancho de banda.

Como se tiene un filtro pasa-bajo de segundo orden se tiene que la banda de paso
va desde la frecuencia cero hasta la frecuencia de corte. Entonces el ancho de
anda es:

Bw=f 0
Bw=31.83 Hz
 Al observar la figura anterior es evidente que inicia la gráfica en G dB=−40.915 dB
e) La respuesta al impulso del filtro en el dominio del tiempo h ( t ).

Dada la función de transferencia

360
H (s)= 2
s +1244.4 s +40000

Se debe hallar h ( t ). Para esto se debe hacer uso de las propiedades de la

Transformada de Laplace y su Inversa. Por otra parte:
 Y ( s)
H (s)=
X (s )

Luego, Y ( s )=H ( s )∗X ( s )

Donde X ( s )=L { δ ( t ) } =1

Entonces reemplazando en la ecuación se tiene lo siguiente:

360
Y ( s )= 2
∗1
s +1244.4 s +40000

Ahora, se debe aplicar a ambos lados de la igualdad la Transformada de Laplace.

360
L−1 { Y ( s ) }=L−1
{ 2
s +1244.4 s +40000 }
Luego

360
y ( t ) =h ( t ) =L−1
{ 2
s +1244.4 s+ 40000 }
Para resolver la parte derecha de la expresión anterior, encontrar las raíces del
polinomio del denominador con el fin de aplicar fracciones parciales.

s2 +1244.4 s+ 40000=0

Entonces las soluciones de la ecuación cuadrática están dados por:

2
−1244.4 ± √ ( 1244.4 ) −4 ( 1 )∗40000
s=
2

s=−655.22± 1833.57 i

Para este caso las raíces son complejas conjugadas, entonces:

A A¿
Y ( s )= +
s−(−655.22+1833.57 i ) s− (−655.22−1833.57 i )

A A¿
Y ( s )= +
s +655.22−1833.57 i s+655.22+1833.57 i
 Ahora se debe encontrar las constantes A y A¿ haciendo uso de:
k i=( s− pi ) H ( s ) ⋮ s− p i

Luego;

360
A=( s +655.22−1833.57 i ) ∨S=(−655.22+1833.57 i)
( s +655.22−1833.57 i ) ( s+ 655.22+1833.57 i )

360
A=
(−655.22+1833.57 i ) +655.22+1833.57 i

360
A=
3667.14 i

360 i
A=
3667.14 i i

−360 i
A=
3667.14

A=−0.0981i

Como A es el conjugado de A , se tiene que A =−0.0981 i. Ahora reemplazando

¿ ¿

los valores encontrados se cumple:

−0.0981 i −0.0981 i
y ( t ) =h ( t ) =L−1 { s+ 655.22−1833.57 i s +655.22+ 1833.57i }
+

Usando las propiedades de linealidad de la Inversa de Laplace

−0.0981i +−0.0981 i
h ( t )=L−1 {s +655.22−1833.57 }i L { s+ 655.22+1833.57
−1
i}

Aplicando la información en la tabla de la Inversa de Laplace se tiene:

h ( t )=2 e−655.22t [ 0 cos ( 1833.57 t )+ 0.0981sin ( 1833.57 t ) ] u ( t )

h ( t )=2 e−655.22t [ 0.0981 sin (1833.57 t ) ] u ( t )

h ( t )=0.1962 e−655.22t sin ( 1833.57 t ) u ( t )

f) La respuesta en estado estable para una entrada x 1 ( t )=4 cos ( 0.02 t+ a ). (a está en
grados y se debe pasar a radianes).
De esta manera x 1 ( t )=4 cos ( 0.02 t+12 )

Para calcular la respuesta de estado permanente de un sistema LTI a una entrada

senoidal de la forma x ( t )= A cos ( ω0 t+θ ).

Y la función de transferencia

360
H (s)= 2
s +1244.4 s +40000

Donde s= j ώ

Como primera medida se debe pasar el desfase a radianes haciendo uso de una
regla de tres simple directa. Para este caso θ=12.

Por lo tanto, el desfase dado en radianes es:

π
12 ≈ 0.2094 rad
180

Además:

ώ=0.02 rad /s

De esta manera la señal de entrada es de la forma x 1 ( t )=4 cos ( 0.02 t+ 0.2094 )

Ahora haciendo uso de las propiedades de Fourier Y 1 ( ώ )=H ( ώ ) X 1 ( ώ )

Donde:

360
H ( ώ)= 2
( j ώ ) + 1244.4 ( j ώ ) + 40000

ώ=0.02 rad /s
j 2 =−1

360
H (0.02)= 2
(−0.02 ) +1244.4 ( j 0.02 ) +40000

360
H (0.02)= 2
(−0.02 ) +24.88 j+40000

Se asume que 40000 ≈ 40000−0.022.

Luego
 360
H (0.02)=
24.88 j+40000

Ahora se multiplica por el conjugado a la expresión anterior.

360
∗40000−24.88 j
24.88 j+ 40000
H ( 0.02 )=
40000−24.88 j

14400000−2488.8
H ( 0.02 )=
1.6 ×109

1.44 ×107 2488.8

H ( 0.02 )= −
1.6 ×109 1.6 ×109

H ( 0.02 )=9 ×10−3−1.55 ×10−6 j

La expresión anterior tiene una forma compleja. Por lo tanto es posible encontrar
un módulo así:

2 2
√
H ( 0.02 )= ( 9 ×10−3 ) + ( 1.55× 10−6 )

H ( 0.02 )=0.009

Además,

−1.55× 10−6
¿ H ( 0.02 )=tan−1
( 9 ×10−3 )
¿ H ( 0.02 )=−1.722 ×10−4

De esta manera

K=0.009
∅=−1.722×10−4

Ahora, se reemplaza H ( 0.02 ) en Y 1 ( ώ )=H ( ώ ) X 1 ( ώ ), se tiene lo siguiente:

Y 1 ( ώ )=( 9 ×10−3 −1.55× 10−6 j ) X 1 ( ώ )

Como el objetivo es encontrar y 1 (t ), se debe hacer uso de la Transformada

Inversa de Fourier
 F
−1
{Y 1 ( ώ ) }=F−1 {( 9 ×10−3−1.55 ×10−6 j ) X 1 ( ώ ) }
Usando las propiedades de la linealidad

F−1 {Y 1 ( ώ ) }=( 9× 10−3−1.55 × 10−6 j ) F−1 { X 1 ( ώ ) }

y 1 (t )=( 9 ×10−3 −1.55× 10−6 j ) x 1 ( t )

Donde x 1 ( t )=4 cos ( 0.02 t+ 0.2094 ) de esta manera la solución es:

y 1 (t )=( 9 ×10−3 −1.55× 10−6 j ) 4 cos ( 0.02 t+0.2094 )

Por otra parte, usando la magnitud y el ángulo

K=0.009
∅=−1.722×10−4

y 1 (t )=0.009 ( 4 ) cos ( 0.02t +0.2094+ ( −1.722× 10−4 ) )

y 1 (t )=0.036 cos ( 0.02t +0.2092 )

g) La respuesta en estado estable para una entrada x 2 ( t )=4 cos ( 300 t+ a ). (a está en
grados y se debe pasar a radianes).

De esta manera x 2 ( t )=4 cos ( 300 t+12 )

Por lo tanto, el desfase dado en radianes es:

π
12 ≈ 0.2094 rad
180

De esta manera la señal de entrada es de la forma x 2 ( t )=4 cos ( 300 t+ 0.2094 )

Ahora haciendo uso de las propiedades de Fourier Y 2 ( ώ )=H ( ώ ) X 2 ( ώ )

Donde:

360
H ( ώ)= 2
( j ώ ) + 1244.4 ( j ώ ) + 40000

ώ=3 00 rad /s
j 2 =−1
 360
H (3 00)= 2
(−3 00 ) +1244.4 ( 3 00 j ) +40000

360
H (3 00)=
323320 j

Ahora se multiplica por el conjugado a la expresión anterior.

360
∗−323320 j
323320 j
H (3 00)=
−323320 j

116395200
H (3 00)=
( 323320 )2

H ( 300 )=1.11 ×10−3 j

La expresión anterior tiene una forma compleja. Por lo tanto es posible encontrar
un módulo así:

2
√
H ( 300 )= ( 1.11 ×10−3 )

H ( 300 )=1.11 ×10−3 j

Además,

π
¿ H ( 300 )= rad
2

De esta manera

K=1.11× 10−3 j

π
∅= rad
2

Ahora, se reemplaza H ( 300 ) en Y 2 ( ώ )=H ( ώ ) X 2 ( ώ ), se tiene lo siguiente:

Y 2 ( ώ )= ( 1.11 ×10−3 j ) X 2 ( ώ )

Como el objetivo es encontrar y 2 ( t ) , se debe hacer uso de la Transformada

Inversa de Fourier
 F
−1
{Y 2 ( ώ ) }=F−1 {( 1.11 ×10−3 j ) X 2 ( ώ) }
Usando las propiedades de la linealidad

F−1 {Y 2 ( ώ ) }=( 1.11 × 10−3 j ) F−1 { X 2 ( ώ ) }

y 1 (t )=( 1.11× 10−3 j ) x 2 ( t )

Donde x 2 ( t )=4 cos ( 300 t+ 0.2094 ) de esta manera la solución es:

y 2 ( t )= ( 1.11 ×10−3 j ) 4 cos ( 300 t+0.2094 )

y 2 ( t )=0.0044 cos ( 300 t+ 0.2094 )

Por otra parte, usando la magnitud y el ángulo

K=1.11× 10−3 j

π
∅= rad
2

π
(
y 2 ( t )= ( 1.11 ×10−3 j ) 4 cos 300t +0.2094+
2 )
π
(
y 2 ( t )=0.0044 cos 300t +0.2094+
2 )
h) La respuesta en estado estable para una entrada x 3 ( t )=4 cos ( 9 00 t +a ) . (a está en
grados y se debe pasar a radianes).

De esta manera x 3 ( t )=4 cos ( 9 00 t +12 )

Por lo tanto, el desfase dado en radianes es:

π
12 ≈ 0.2094 rad
180

De esta manera la señal de entrada es de la forma x 3 ( t )=4 cos ( 9 00 t +0.2094 )

Ahora haciendo uso de las propiedades de Fourier Y 3 ( ώ ) =H ( ώ ) X 3 ( ώ )

Donde:

360
H ( ώ)= 2
( j ώ ) + 1244.4 ( j ώ ) + 40000

ώ=9 00 rad / s
j 2 =−1

360
H (9 00)= 2
(−9 00 ) +1244.4 ( 900 j )+ 40000
360
H (9 00)=
−850000+1119960 j

Ahora se multiplica por el conjugado a la expresión anterior.

360
∗−850000+1119960 j
−850000+1119960 j
H (9 00)=
−850000+1119960 j

−306000000−403185600
H (9 00)=
9 .520 ×1011

3 . 06 ×10 8 4 . 03 ×108
H (9 00)= −
9.520 ×10 11 9.520× 1011

H ( 9 00 )=−3.21 ×10−4 −4 . 23 ×10−4 j

La expresión anterior tiene una forma compleja. Por lo tanto es posible encontrar
un módulo así:

2 2
√
H (9 00)= ( 3.21 ×10−4 ) + ( 4.23× 10−4 )

H (9 00)=5 . 3 1× 10−4

Además,

−4.23× 10−4
¿ H (9 00)=tan −1 ( −3.21 ×10−4 )
¿ H (9 00)=0.9216

De esta manera
 K=5.31×10−4
∅=0.9216

Ahora, se reemplaza H (9 00) en Y 3 ( ώ ) =H ( ώ ) X 3 ( ώ ), se tiene lo siguiente:

Y 3 ( ώ ) =(−3.21× 10−4 −4.23 × 10−4 j ) X 3 ( ώ )

Como el objetivo es encontrar y 3 ( t ) , se debe hacer uso de la Transformada

Inversa de Fourier

F−1 {Y 3 ( ώ ) }=F−1 {(−3.21 ×10−4−4.23 ×10−4 j ) X 3 ( ώ ) }

Usando las propiedades de la linealidad

F−1 {Y 3 ( ώ ) }=( −3.21× 10−4 −4.23 ×10−4 j ) F−1 { X 3 ( ώ ) }

y 3 ( t ) =(−3.21× 10−4 −4.23 × 10−4 j ) x 3 ( t )

Donde x 3 ( t )=4 cos ( 900 t +0.2094 ) de esta manera la solución es:

y 3 ( t ) =(−3.21× 10−4 −4.23 × 10−4 j ) 4 cos ( 900 t+ 0.2094 )

Por otra parte, usando la magnitud y el ángulo

K=5.31×10−4
∅=0.9216

y 1 (t )=5.31 ×10−4 ( 4 ) cos ( 900 t+0.2094 +0.9216 )

y 1 (t )=0.002 cos ( 900 t +0.2094+ 0.9216 )

i) Utilice software para verificar sus resultados (se debe realizar el diagrama de
bode).