23.11.2013 Views

respuesta en frecuencia de sistemas lineales, invariantes en el ...

respuesta en frecuencia de sistemas lineales, invariantes en el ...

respuesta en frecuencia de sistemas lineales, invariantes en el ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Universidad Nacional <strong>de</strong> San Juan<br />

Facultad <strong>de</strong> Ing<strong>en</strong>iería<br />

Departam<strong>en</strong>to <strong>de</strong> Electrónica<br />

y Automática<br />

“RESPUESTA EN FRECUENCIA DE SISTEMAS LINEALES,<br />

INVARIANTES EN EL TIEMPO.”<br />

Cátedra: Control I.<br />

Carreras: Ing<strong>en</strong>iería Electrónica y Bioing<strong>en</strong>iería.<br />

Autores:<br />

Ing. Mario Alberto Perez.<br />

Ing. Analía Perez Hidalgo.<br />

Dra. Bioing. Elisa Perez Ber<strong>en</strong>guer.


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

RESPUESTA FRECUENCIAL<br />

1- Introducción.<br />

Ya se ha investigado la <strong>respuesta</strong> <strong>de</strong> compon<strong>en</strong>tes y <strong>sistemas</strong> a varios tipos <strong>de</strong><br />

<strong>en</strong>tradas <strong>en</strong> <strong>el</strong> dominio temporal. Se vio que la función <strong>respuesta</strong> c (t)<br />

conti<strong>en</strong>e<br />

dos términos, un término transitorio (la solución complem<strong>en</strong>taria) y un<br />

término <strong>de</strong> estado estacionario o constante (la solución particular), obt<strong>en</strong>idos<br />

ambos por la solución <strong>de</strong> la ecuación <strong>de</strong>l sistema, cuando es aplicada una<br />

excitación <strong>en</strong> la <strong>en</strong>trada.<br />

El pres<strong>en</strong>te capítulo se <strong>de</strong>dicará al estudio <strong>de</strong> la <strong>respuesta</strong> <strong>en</strong> estado estable<br />

<strong>de</strong> compon<strong>en</strong>tes y <strong>sistemas</strong> cuando sean excitados por una señal s<strong>en</strong>oidal <strong>de</strong><br />

amplitud fija pero con una frecu<strong>en</strong>cia que varía <strong>en</strong> un cierto rango.<br />

Este concepto se ilustra <strong>en</strong> la figura 1, <strong>en</strong> la cual un sistema lineal es excitado<br />

por una señal s<strong>en</strong> t<br />

b ( ω ) s<strong>en</strong> ωt<br />

+ φ(<br />

ω)<br />

.<br />

a ω , la <strong>respuesta</strong> es [ ]<br />

Figura 1.<br />

La forma <strong>de</strong> las ondas <strong>de</strong> <strong>en</strong>trada y salida se ilustran <strong>en</strong> la figura 2.<br />

Figura 2.<br />

Este resultado obt<strong>en</strong>ido concuerda totalm<strong>en</strong>te con lo ya visto como solución<br />

particular <strong>de</strong> un sistema cuando era excitado por una señal armónica <strong>de</strong> la<br />

forma r( t)<br />

= H s<strong>en</strong>ωt<br />

, don<strong>de</strong> H es la amplitud constante y ω la frecu<strong>en</strong>cia<br />

angular <strong>de</strong> <strong>en</strong>trada, estando <strong>el</strong> sistema <strong>en</strong> estado estacionario.<br />

Para <strong>el</strong> caso <strong>de</strong> un sistema <strong>de</strong> primer ord<strong>en</strong> la solución particular era:<br />

- 1 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

K H<br />

c ( t)<br />

= s<strong>en</strong> t<br />

φ = arctg ω<br />

1<br />

2 2<br />

( 1+<br />

ω T )<br />

( ω − φ)<br />

Como se ve, tanto la amplitud <strong>de</strong> la <strong>respuesta</strong> c<br />

1(<br />

t ) como la fase son ambas<br />

funciones <strong>de</strong> la frecu<strong>en</strong>cia ω <strong>de</strong> la <strong>en</strong>trada, lo que corrobora lo dicho<br />

anteriorm<strong>en</strong>te.<br />

El sigui<strong>en</strong>te gráfico ilustra lo que se acaba <strong>de</strong> m<strong>en</strong>cionar:<br />

- 2 -<br />

Figura 3.<br />

Es común <strong>en</strong> <strong>el</strong> análisis frecu<strong>en</strong>cial que <strong>el</strong> interés se c<strong>en</strong>tre <strong>en</strong> <strong>el</strong> estudio <strong>de</strong> las<br />

sigui<strong>en</strong>tes r<strong>el</strong>aciones:<br />

a) La r<strong>el</strong>ación <strong>de</strong> amplitud b a<br />

la <strong>de</strong>signa como M (ω)<br />

.<br />

, que se la d<strong>en</strong>omina r<strong>el</strong>ación <strong>de</strong> magnitud y se<br />

b) El ángulo <strong>de</strong> fase φ (ω)<br />

. Un ángulo <strong>de</strong> fase negativo recibe <strong>el</strong> nombre <strong>de</strong><br />

retardo <strong>de</strong> fase, y un ángulo <strong>de</strong> fase positiva es d<strong>en</strong>ominado a<strong>de</strong>lanto <strong>de</strong> fase.<br />

Se tratará ahora la <strong>de</strong>terminación <strong>de</strong> información sobre la <strong>respuesta</strong> a la frecu<strong>en</strong>cia,<br />

<strong>de</strong> un modo analítico, aunque tales datos se pued<strong>en</strong> obt<strong>en</strong>er<br />

experim<strong>en</strong>talm<strong>en</strong>te si <strong>el</strong> sistema existe.<br />

Las mediciones <strong>de</strong> <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia <strong>en</strong> g<strong>en</strong>eral son simples y pued<strong>en</strong><br />

ser efectuadas con exactitud usando g<strong>en</strong>eradores <strong>de</strong> señal s<strong>en</strong>oidales<br />

fácilm<strong>en</strong>te obt<strong>en</strong>ibles y equipos <strong>de</strong> medición precisos.<br />

Frecu<strong>en</strong>tem<strong>en</strong>te se pued<strong>en</strong> <strong>de</strong>terminar experim<strong>en</strong>talm<strong>en</strong>te, las funciones<br />

transfer<strong>en</strong>cia <strong>de</strong> compon<strong>en</strong>tes complicados <strong>en</strong> prueba <strong>de</strong> <strong>respuesta</strong> <strong>de</strong><br />

frecu<strong>en</strong>cia. A<strong>de</strong>más, <strong>el</strong> método <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia ti<strong>en</strong>e la v<strong>en</strong>taja <strong>de</strong><br />

que se pue<strong>de</strong> diseñar un sistema <strong>de</strong> manera que los efectos <strong>de</strong>l ruido<br />

in<strong>de</strong>seable sean <strong>de</strong>spreciables, y <strong>de</strong> que ese análisis y diseño pueda<br />

ext<strong>en</strong><strong>de</strong>rse a ciertos <strong>sistemas</strong> <strong>de</strong> control no <strong>lineales</strong>.<br />

Obt<strong>en</strong>er la <strong>respuesta</strong> frecu<strong>en</strong>cial es importante puesto que proporciona medios<br />

conv<strong>en</strong>i<strong>en</strong>tes para obt<strong>en</strong>er la <strong>respuesta</strong> <strong>en</strong> <strong>el</strong> estado estable pera cualquier<br />

sistema lineal sujeto a una señal s<strong>en</strong>oidal. También está r<strong>el</strong>acionado<br />

íntimam<strong>en</strong>te con <strong>el</strong> método <strong>de</strong> análisis frecu<strong>en</strong>cial que se verá más a<strong>de</strong>lante.<br />

El procedimi<strong>en</strong>to simple para obt<strong>en</strong>er la <strong>respuesta</strong> a la frecu<strong>en</strong>cia, cont<strong>en</strong>ido<br />

<strong>en</strong> cuatro pasos, es <strong>el</strong> sigui<strong>en</strong>te:<br />

1 - Se obti<strong>en</strong>e la función <strong>de</strong> transfer<strong>en</strong>cia para <strong>el</strong> compon<strong>en</strong>te o <strong>sistemas</strong> a<br />

analizar. Es <strong>de</strong>cir:


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

C(<br />

s)<br />

F ( s)<br />

=<br />

R(<br />

s)<br />

Don<strong>de</strong> C (s)<br />

es la transformada <strong>de</strong> la salida y R (s)<br />

la transformada <strong>de</strong> la<br />

<strong>en</strong>trada, y don<strong>de</strong> se han <strong>de</strong>spreciado todas las condiciones iniciales porque se<br />

vio no afectada la <strong>respuesta</strong> <strong>en</strong> estado estable.<br />

2 - En la función <strong>de</strong> transfer<strong>en</strong>cia se sustituye s por j ω . La justificación <strong>de</strong><br />

esta sustitución se hará más a<strong>de</strong>lante.<br />

3 - Para varios valores <strong>de</strong> la frecu<strong>en</strong>cia ω , se <strong>de</strong>termina la r<strong>el</strong>ación <strong>de</strong><br />

magnitud M (ω)<br />

y <strong>el</strong> ángulo <strong>de</strong> fase φ (ω)<br />

.<br />

4 - Se grafican los resultados <strong>de</strong> 3 <strong>en</strong> coord<strong>en</strong>adas polares o rectangulares.<br />

Estas gráficas no solam<strong>en</strong>te san medios conv<strong>en</strong>i<strong>en</strong>tes para pres<strong>en</strong>tar los datos<br />

<strong>de</strong> <strong>respuesta</strong> a la frecu<strong>en</strong>cia, sino que también son la base para los métodos<br />

<strong>de</strong> análisis y diseño que se verán <strong>en</strong> capítulos posteriores.<br />

Aunque la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> un sistema <strong>de</strong> control <strong>de</strong> una imag<strong>en</strong><br />

cualitativa <strong>de</strong> la <strong>respuesta</strong> transitoria, la corr<strong>el</strong>ación <strong>en</strong>tre frecu<strong>en</strong>cia y<br />

<strong>respuesta</strong>s transitorias, es indirecta, excepto <strong>en</strong> <strong>el</strong> caso <strong>de</strong> <strong>sistemas</strong> <strong>de</strong><br />

segundo ord<strong>en</strong>. Al proyectar un sistema <strong>de</strong> lazo cerrado, se pue<strong>de</strong> ajustar la<br />

característica <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia, usando diversos criterios <strong>de</strong> diseño<br />

para obt<strong>en</strong>er características <strong>de</strong> <strong>respuesta</strong> transitoria aceptables.<br />

Una vez <strong>en</strong>t<strong>en</strong>dida la corr<strong>el</strong>ación indirecta <strong>en</strong>tre diversas mediciones <strong>de</strong> la<br />

<strong>respuesta</strong> transitoria y la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia, pue<strong>de</strong> utilizarse<br />

v<strong>en</strong>tajosam<strong>en</strong>te <strong>el</strong> método <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia.<br />

El diseño <strong>de</strong> un sistema <strong>de</strong> control basado <strong>en</strong> este procedimi<strong>en</strong>to, se funda <strong>en</strong><br />

la interpretación <strong>de</strong> las características <strong>de</strong> <strong>respuesta</strong> a la frecu<strong>en</strong>cia. Este<br />

análisis <strong>de</strong> un sistema <strong>de</strong> control, indica gráficam<strong>en</strong>te qué modificaciones hay<br />

que hacer <strong>en</strong> la función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> lazo abierto para obt<strong>en</strong>er las<br />

características <strong>de</strong>seadas <strong>de</strong> <strong>respuesta</strong> transitoria.<br />

Lógicam<strong>en</strong>te se podría hacer la pregunta <strong>de</strong> por qué es tan importante <strong>el</strong><br />

análisis <strong>de</strong> la <strong>respuesta</strong> <strong>de</strong> <strong>sistemas</strong> a una señal s<strong>en</strong>oidal, cuando <strong>en</strong> la<br />

práctica, los <strong>sistemas</strong> <strong>de</strong> control raram<strong>en</strong>te están expuestos a señales<br />

armónicas. La <strong>respuesta</strong> es que la información obt<strong>en</strong>ida por <strong>el</strong> análisis<br />

s<strong>en</strong>oidal pue<strong>de</strong> usarse para establecer la naturaleza <strong>de</strong> la <strong>respuesta</strong> a una<br />

gran variedad <strong>de</strong> señales. A<strong>de</strong>más, <strong>el</strong> análisis es conv<strong>en</strong>i<strong>en</strong>te para manejarlo<br />

matemática y experim<strong>en</strong>talm<strong>en</strong>te.<br />

- 3 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

2- Justificación <strong>de</strong> la Sustitución <strong>de</strong> s por jω.<br />

Ya se ha visto que para realizar <strong>el</strong> análisis frecu<strong>en</strong>cial es necesario sustituir s<br />

por j ω <strong>en</strong> la función <strong>de</strong> transfer<strong>en</strong>cia. Esta sustitución será ahora justificada.<br />

El procedimi<strong>en</strong>to consistirá <strong>en</strong> trabajar con una función <strong>de</strong> transfer<strong>en</strong>cia<br />

g<strong>en</strong>eral, obt<strong>en</strong>i<strong>en</strong>do primero la <strong>respuesta</strong> <strong>en</strong> <strong>el</strong> estado estable a una señal<br />

s<strong>en</strong>oidal usando la transformada <strong>de</strong> Laplace, y <strong>de</strong>spués haci<strong>en</strong>do la sustitución<br />

<strong>de</strong> s por j ω .<br />

Si ambas soluciones resultan idénticas la sustitución es válida. Se supone una<br />

función <strong>de</strong> transfer<strong>en</strong>cia g<strong>en</strong>eral con un numerador N (s)<br />

y un d<strong>en</strong>ominador<br />

D (s) .<br />

C(<br />

s)<br />

N(<br />

s)<br />

F ( s)<br />

= =<br />

(1)<br />

R(<br />

s)<br />

D(<br />

s)<br />

La señal <strong>de</strong> <strong>en</strong>trada será una s<strong>en</strong>oi<strong>de</strong> con amplitud igual a uno ( s<strong>en</strong>ω t ). La<br />

transformada <strong>de</strong> R (s)<br />

es:<br />

s<br />

2<br />

ω<br />

2<br />

+ ω<br />

; por lo tanto:<br />

N(<br />

s)<br />

ω N(<br />

s)<br />

ω N(<br />

s)<br />

C( s)<br />

= R(<br />

s)<br />

=<br />

=<br />

(2)<br />

2 2<br />

D(<br />

s)<br />

s + ω D(<br />

s)<br />

( s + jω)( s − jω) D(<br />

s)<br />

Hay que recordar que todas las condiciones iniciales pued<strong>en</strong> <strong>de</strong>spreciarse porque<br />

éstas no afectan la <strong>respuesta</strong> <strong>en</strong> <strong>el</strong> estado estable.<br />

Desarrollando <strong>en</strong> fracciones parciales:<br />

C(<br />

s)<br />

=<br />

C<br />

+<br />

C<br />

+<br />

C<br />

1 2<br />

3<br />

2<br />

.......<br />

( s + jω) ( s + jω) ( s + r ) ( s + r ) ( s + r )<br />

1<br />

+<br />

C<br />

2<br />

+<br />

+<br />

C<br />

n<br />

n−2<br />

(3)<br />

s + , …., son factores <strong>de</strong> D (s)<br />

. Para un sistema estable los<br />

C3<br />

transitorios <strong>de</strong>saparec<strong>en</strong> (o sea las transformadas inversas <strong>de</strong> … se<br />

( s + r1)<br />

anulan conforme t → ∞ ) y la <strong>respuesta</strong> <strong>en</strong> <strong>el</strong> estado estable es:<br />

Don<strong>de</strong> ( s + r 1<br />

), ( r 2<br />

)<br />

C<br />

− jω<br />

t<br />

SS<br />

( t)<br />

= C1e<br />

+<br />

C<br />

2<br />

e<br />

jω<br />

t<br />

Las constantes C<br />

1<br />

y C<br />

2<br />

se <strong>de</strong>terminan por <strong>el</strong> método:<br />

(4)<br />

C<br />

n<br />

⎡ N(<br />

s)<br />

= lim<br />

S→<br />

r<br />

⎢<br />

n<br />

⎣ D(<br />

s)<br />

⎤<br />

( s + r ) ⎥⎦<br />

n<br />

o sea para este caso:<br />

C<br />

1<br />

=<br />

N(<br />

s)<br />

D(<br />

s)<br />

( s + r )<br />

1<br />

C<br />

N(<br />

s)<br />

=<br />

ω<br />

ω<br />

s + jω<br />

N(<br />

− jω)<br />

=<br />

( )( ) ( ) s + j s − jω<br />

D(<br />

− jω)<br />

− jω<br />

1<br />

D(<br />

s)<br />

2<br />

S =−r<br />

S =− jω<br />

1<br />

∴<br />

ω<br />

- 4 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

1 N(<br />

− jω)<br />

∴ C1<br />

= −<br />

(5)<br />

2 j D(<br />

− jω)<br />

N(<br />

s)<br />

ω<br />

C2 =<br />

s + jω<br />

=<br />

D(<br />

s)<br />

N(<br />

jω)<br />

ω<br />

( )( ) ( ) s + jω<br />

s − jω<br />

D(<br />

jω)<br />

jω<br />

2<br />

S = + jω<br />

1 N(<br />

jω)<br />

∴ C<br />

1<br />

=<br />

(6)<br />

2 j D(<br />

jω)<br />

Como alternativa, la fracción<br />

N(<br />

jω)<br />

D(<br />

jω)<br />

pue<strong>de</strong> expresarse <strong>de</strong> la sigui<strong>en</strong>te forma:<br />

N(<br />

jω)<br />

= A(<br />

ω)<br />

+ jB(<br />

ω)<br />

(7)<br />

D(<br />

jω)<br />

N(<br />

− jω)<br />

= A(<br />

ω)<br />

−<br />

D(<br />

− jω)<br />

jB(<br />

ω)<br />

(8)<br />

Don<strong>de</strong> A (ω)<br />

y B (ω)<br />

son números reales y son funciones <strong>de</strong> la frecu<strong>en</strong>cia. Por lo<br />

tanto:<br />

1<br />

1<br />

C1 = − [ A(<br />

ω)<br />

− jB(<br />

ω)<br />

] y C<br />

2<br />

= [ A(<br />

ω)<br />

+ jB(<br />

ω)<br />

]<br />

2 j<br />

2 j<br />

Sustituy<strong>en</strong>do estos valores <strong>en</strong> la ecuación (4):<br />

C<br />

SS<br />

1<br />

− jω<br />

t 1<br />

+ jω<br />

t<br />

( t)<br />

= − [ A(<br />

ω)<br />

− jB(<br />

ω)<br />

] e + [ A(<br />

ω)<br />

+ jB(<br />

ω)<br />

] e<br />

(9)<br />

2 j<br />

2 j<br />

C<br />

A(<br />

ω)<br />

2 j<br />

B(<br />

ω)<br />

2 j<br />

A(<br />

ω)<br />

2 j<br />

− jω<br />

t<br />

− jω<br />

t<br />

jω<br />

t<br />

SS<br />

( t)<br />

= − e + j e + e +<br />

B(<br />

ω)<br />

j e<br />

2 j<br />

jω<br />

t<br />

(10)<br />

C<br />

SS<br />

( t)<br />

=<br />

⎡e<br />

A(<br />

ω ) ⎢<br />

⎣<br />

jω<br />

t<br />

− e<br />

2 j<br />

− jω<br />

t<br />

jω<br />

t<br />

⎤ ⎡ e + e<br />

⎥ + B(<br />

ω)<br />

⎢<br />

⎦ ⎣ 2<br />

− jω<br />

t<br />

⎤<br />

⎥<br />

⎦<br />

∴ ( t)<br />

= A(<br />

ω)<br />

s<strong>en</strong>ω<br />

t + B(<br />

ω)<br />

cosω<br />

t<br />

(11)<br />

C SS<br />

La <strong>respuesta</strong> <strong>en</strong> estado estable se ve que está compuesta por una s<strong>en</strong>oi<strong>de</strong><br />

más una cos<strong>en</strong>oi<strong>de</strong>. La suma <strong>de</strong> estas dos ondas se lleva a cabo por adición<br />

vectorial como se muestra <strong>en</strong> la figura 4.<br />

- 5 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Figura 4.<br />

Aquí un vector o fasor <strong>de</strong> magnitud A (ω)<br />

se consi<strong>de</strong>ra que esta girando <strong>en</strong><br />

dirección contraria a las manecillas <strong>de</strong>l r<strong>el</strong>oj con una frecu<strong>en</strong>cia ω . La<br />

proyección <strong>de</strong> este vector sobre <strong>el</strong> eje imaginario produce la s<strong>en</strong>oi<strong>de</strong><br />

requerida. Un vector <strong>de</strong> magnitud B (ω)<br />

se muestra a<strong>de</strong>lantado 90° con<br />

respecto al primer vector, su proyección sobré <strong>el</strong> eje imaginario produce la<br />

cos<strong>en</strong>oi<strong>de</strong> requerida, Los dos vectores pued<strong>en</strong> reemplazarse par un solo<br />

vector <strong>de</strong> magnitud:<br />

2<br />

A ( ω)<br />

+ B<br />

2<br />

( ω)<br />

Y ángulo <strong>de</strong> fase:<br />

B(<br />

ω)<br />

φ(<br />

ω)<br />

= arc.<br />

tg<br />

A(<br />

ω)<br />

Por lo tanta la (10) pue<strong>de</strong> escribirse como:<br />

C SS<br />

2<br />

2<br />

( t)<br />

= A ( ω ) + B ( ω)<br />

s<strong>en</strong> t +<br />

( ω φ(<br />

ω)<br />

)<br />

Don<strong>de</strong>:<br />

B(<br />

ω)<br />

φ(<br />

ω)<br />

= arc.<br />

tg<br />

(12)<br />

A(<br />

ω)<br />

Hay que recordar que la señal <strong>de</strong> <strong>en</strong>trada era una s<strong>en</strong>oi<strong>de</strong> can amplitud<br />

unitaria Por tanto, la r<strong>el</strong>ación <strong>de</strong> magnitud, es <strong>de</strong>cir <strong>el</strong> coci<strong>en</strong>te <strong>en</strong>tre las<br />

amplitu<strong>de</strong>s <strong>de</strong> salida y <strong>en</strong>trada es:<br />

2<br />

2<br />

M ( ω)<br />

= A ( ω)<br />

+ B ( ω)<br />

(13)<br />

Y <strong>el</strong> ángulo <strong>de</strong> fase:<br />

B(<br />

ω)<br />

φ(<br />

ω)<br />

= arc.<br />

tg<br />

(14)<br />

A(<br />

ω)<br />

Debe reconocerse que <strong>el</strong> ángulo <strong>de</strong> fase también pue<strong>de</strong> ser negativo.<br />

- 6 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Habi<strong>en</strong>do completado la solución usando transformada <strong>de</strong> Laplace, falta<br />

<strong>de</strong>mostrar que <strong>el</strong> mismo resultado pue<strong>de</strong> también pue<strong>de</strong> obt<strong>en</strong>erse<br />

directam<strong>en</strong>te haci<strong>en</strong>do la sustitución <strong>de</strong> s por jω <strong>en</strong>:<br />

Cs () Ns () C( j ) N( j )<br />

A( ) jB( )<br />

Rs () = Ds () → ω ω<br />

ω ω<br />

R( jω) = D( jω)<br />

= + (15)<br />

C( jω)<br />

2 2<br />

= ⎡A ( ω) + B ( ω) ⎤ s<strong>en</strong> ωt+<br />

φ( ω)<br />

R( jω)<br />

⎣ ⎦<br />

( )<br />

2 2<br />

∴ M( ω) = ⎡<br />

⎣A ( ω) + B ( ω)<br />

⎤<br />

⎦ (16)<br />

⎡B( ω)<br />

⎤<br />

φω ( ) = arctag ⎢ A ( ω)<br />

⎥<br />

⎣ ⎦<br />

(17)<br />

Que es <strong>el</strong> mismo resultado obt<strong>en</strong>ido anteriorm<strong>en</strong>te, ahora con m<strong>en</strong>os<br />

esfuerzo.<br />

- 7 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

3- Respuesta <strong>de</strong> Frecu<strong>en</strong>cia a partir <strong>de</strong> los Diagramas <strong>de</strong> Polos y<br />

Ceros.<br />

Se pue<strong>de</strong> <strong>de</strong>terminar gráficam<strong>en</strong>te la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia a partir <strong>de</strong> los<br />

diagramas <strong>de</strong> polos y ceros <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia. Sea la sigui<strong>en</strong>te<br />

función <strong>de</strong> transfer<strong>en</strong>cia:<br />

Gs () =<br />

Ks ( + z)<br />

ss ( + p)<br />

(18)<br />

Don<strong>de</strong> p y z son reales. Se pue<strong>de</strong> obt<strong>en</strong>er la <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia <strong>de</strong> esta<br />

función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> la r<strong>el</strong>ación:<br />

G( jω)<br />

=<br />

K( jω<br />

+ z)<br />

s( jω<br />

+ p)<br />

(19)<br />

Los factores jω + z , jω + p son magnitu<strong>de</strong>s complejas como pue<strong>de</strong> verse <strong>en</strong> la<br />

figura 5.<br />

Figura 5: Determinación <strong>de</strong> la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>en</strong> <strong>el</strong> plano complejo.<br />

La amplitud <strong>de</strong> G( jω ) es:<br />

G( jω)<br />

=<br />

K jω<br />

+ z<br />

jω<br />

jω+<br />

p<br />

(20)<br />

KAP<br />

G( jω ) = (21)<br />

OP BP<br />

- 8 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Y <strong>el</strong> ángulo <strong>de</strong> fase <strong>de</strong> G( jω ) es:<br />

G( jω)<br />

= jω+ z− jω− jω+ p<br />

(22)<br />

⎡ω<br />

⎤<br />

⎡ω<br />

⎤<br />

G( jω) = arctag<br />

⎢<br />

−90º<br />

−arctag<br />

z ⎥<br />

⎢<br />

p<br />

⎥<br />

⎣ ⎦ ⎣ ⎦<br />

(23)<br />

G( jω)<br />

= φ−φ − φ<br />

(24)<br />

1 2<br />

Don<strong>de</strong> los ángulos φ , φ 1<br />

y φ 2<br />

están <strong>de</strong>finidos <strong>en</strong> la figura 5. Se hace notar que<br />

se <strong>de</strong>fine como s<strong>en</strong>tido positivo para la medición <strong>de</strong> ángulo a la rotación anti<br />

horaria.<br />

D<strong>el</strong> análisis <strong>de</strong> la <strong>respuesta</strong> transitoria <strong>de</strong> los <strong>sistemas</strong> <strong>de</strong> lazo cerrado, se<br />

sabe que un par <strong>de</strong> polos complejos conjugados cercanos al eje jω produce<br />

un modo altam<strong>en</strong>te oscilatorio <strong>de</strong> la <strong>respuesta</strong> transitoria. En <strong>el</strong> caso <strong>de</strong> la<br />

<strong>respuesta</strong> frecu<strong>en</strong>cial, un par <strong>de</strong> polos así ubicados han <strong>de</strong> producir una<br />

<strong>respuesta</strong> <strong>de</strong> pico <strong>el</strong>evado.<br />

Sea, por ejemplo, la sigui<strong>en</strong>te función <strong>de</strong> transfer<strong>en</strong>cia:<br />

K<br />

Gs () =<br />

( s+ p ) + ( s+<br />

p )<br />

1 2<br />

(25)<br />

Don<strong>de</strong> p<br />

1<br />

y p<br />

2<br />

son complejos conjugados como se ve <strong>en</strong> la figura 6. Se pue<strong>de</strong><br />

hallar la <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia <strong>de</strong> esta función <strong>de</strong> transfer<strong>en</strong>cia.<br />

Figura 6: Determinación <strong>de</strong> la <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia <strong>en</strong> <strong>el</strong> plano complejo.<br />

G( jω)<br />

=<br />

K<br />

jω+ p jω+<br />

p<br />

1 2<br />

(26)<br />

- 9 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

K<br />

G( jω ) = (27)<br />

AP BP<br />

G( jω)<br />

=−φ − φ<br />

(28)<br />

1 2<br />

Don<strong>de</strong> los ángulos φ 1<br />

y φ 2<br />

están <strong>de</strong>finidos <strong>en</strong> la figura 6. Como AP BP es muy<br />

pequeño <strong>en</strong> la cercanía <strong>de</strong> ω = ω1<br />

, G( jω<br />

1)<br />

es muy gran<strong>de</strong>.<br />

De manera que un par <strong>de</strong> polos complejos conjugados cerca <strong>de</strong>l eje jω ha <strong>de</strong><br />

producir una <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> pico <strong>el</strong>evado. Inversam<strong>en</strong>te, si la<br />

<strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia no es <strong>de</strong> picos <strong>el</strong>evados, la función <strong>de</strong> transfer<strong>en</strong>cia<br />

no ha <strong>de</strong> t<strong>en</strong>er polos complejos conjugados cerca <strong>de</strong>l eje jω . Una función <strong>de</strong><br />

transfer<strong>en</strong>cia así no ha <strong>de</strong> pres<strong>en</strong>tar <strong>respuesta</strong> transitoria altam<strong>en</strong>te<br />

oscilatoria.<br />

Como la <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia indirectam<strong>en</strong>te <strong>de</strong>scribe la posición <strong>de</strong> los<br />

polos y ceros <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia, se pued<strong>en</strong> estimar las<br />

características <strong>de</strong> <strong>respuesta</strong> transitoria <strong>de</strong>l conocimi<strong>en</strong>to <strong>de</strong> sus características<br />

<strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia (esto se <strong>en</strong>t<strong>en</strong><strong>de</strong>rá más claram<strong>en</strong>te cuando se<br />

estudie <strong>el</strong> concepto <strong>de</strong> estabilidad r<strong>el</strong>ativa).<br />

- 10 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

4- Especificaciones <strong>en</strong> <strong>el</strong> Dominio Frecu<strong>en</strong>cial.<br />

Se ha puesto <strong>de</strong> manifiesto anteriorm<strong>en</strong>te que la información que se busca <strong>en</strong><br />

<strong>el</strong> análisis <strong>de</strong> un sistema <strong>de</strong> control es normalm<strong>en</strong>te la <strong>respuesta</strong> temporal.<br />

Sin embargo, <strong>en</strong> g<strong>en</strong>eral, la <strong>respuesta</strong> temporal es difícil <strong>de</strong> obt<strong>en</strong>er<br />

analíticam<strong>en</strong>te a causa <strong>de</strong> la cantidad <strong>de</strong> cálculos que implica. Por<br />

consigui<strong>en</strong>te, la <strong>respuesta</strong> frecu<strong>en</strong>cial <strong>de</strong> un sistema <strong>de</strong> control se obti<strong>en</strong>e a<br />

m<strong>en</strong>udo por medio <strong>de</strong> métodos gráficos (repres<strong>en</strong>taciones polares y<br />

rectangulares tales como los diagramas <strong>de</strong> Nyquist y Bo<strong>de</strong> que se estudian<br />

más a<strong>de</strong>lante) y luego, la interpretación <strong>de</strong>l comportami<strong>en</strong>to <strong>de</strong>l sistema <strong>en</strong> <strong>el</strong><br />

dominio temporal se basa <strong>en</strong> las r<strong>el</strong>aciones <strong>en</strong>tre ambos dominios, temporal y<br />

frecu<strong>en</strong>cial.<br />

El punto <strong>de</strong> partida para <strong>el</strong> análisis <strong>en</strong> <strong>el</strong> dominio frecu<strong>en</strong>cial es la función <strong>de</strong><br />

transfer<strong>en</strong>cia. Para un sistema <strong>de</strong> control con realim<strong>en</strong>tación unitaria, la<br />

función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> lazo cerrado es:<br />

C(<br />

s)<br />

R(<br />

s)<br />

G(<br />

s)<br />

= (29)<br />

1+<br />

G(<br />

s)<br />

En condiciones <strong>de</strong> régim<strong>en</strong> sinusoidal,<br />

s =<br />

jω<br />

, la ecuación 29 se convierte <strong>en</strong>:<br />

C(<br />

jω)<br />

G(<br />

jω)<br />

= M ( jω)<br />

=<br />

(30)<br />

R(<br />

jω)<br />

1+<br />

G(<br />

jω)<br />

Cuando M ( jω)<br />

se escribe <strong>en</strong> forma <strong>de</strong> amplitud y fase, se ti<strong>en</strong>e:<br />

M( jω) = M( ω) φ ( ω)<br />

(31)<br />

Don<strong>de</strong>:<br />

m<br />

G(<br />

jω)<br />

M ( ω)<br />

= (32)<br />

1+<br />

G(<br />

jω)<br />

Y<br />

φm ( ω) = G( jω) − 1 + G( jω)<br />

(33)<br />

El significado <strong>de</strong> M ( jω)<br />

<strong>en</strong> un sistema <strong>de</strong> control es similar a la ganancia o<br />

amplificación <strong>de</strong> un amplificador <strong>el</strong>ectrónico. En un amplificador <strong>de</strong><br />

audiofrecu<strong>en</strong>cia, por ejemplo, <strong>el</strong> criterio i<strong>de</strong>al <strong>de</strong> proyecto es que <strong>el</strong><br />

amplificador t<strong>en</strong>ga una curva <strong>de</strong> ganancia plana para las audiofrecu<strong>en</strong>cias. En<br />

los <strong>sistemas</strong> <strong>de</strong> control, sin embargo. La situación i<strong>de</strong>al <strong>en</strong> algunas ocasiones<br />

es que la salida siga a la <strong>en</strong>trada <strong>en</strong> todo instante o, simplem<strong>en</strong>te, que <strong>el</strong><br />

módulo <strong>de</strong> M ( jω)<br />

sea igual a la unidad para todas las frecu<strong>en</strong>cias. Pero <strong>en</strong> la<br />

expresión <strong>de</strong> la ecuación 32 se ve que M ( jω)<br />

solo pe<strong>de</strong> ser igual a la unidad<br />

cuando G ( jω)<br />

es infinito o, <strong>en</strong> otras palabras, la ganancia <strong>de</strong>l sistema <strong>de</strong>be ser<br />

infinita para todas las frecu<strong>en</strong>cias. Esto es imposible <strong>de</strong> conseguir <strong>en</strong> la<br />

práctica y a<strong>de</strong>más no es conv<strong>en</strong>i<strong>en</strong>te puesto que la mayoría <strong>de</strong> los <strong>sistemas</strong><br />

<strong>de</strong> control resultan inestables para valores <strong>el</strong>evados <strong>de</strong> ganancia. A<strong>de</strong>más,<br />

- 11 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

todos los <strong>sistemas</strong> <strong>de</strong> control están sujetos a ruidos, es <strong>de</strong>cir, a<strong>de</strong>más <strong>de</strong><br />

respon<strong>de</strong>r a la señal <strong>de</strong> <strong>en</strong>trada los <strong>sistemas</strong> <strong>de</strong>b<strong>en</strong> ser capaces <strong>de</strong> rechazar y<br />

suprimir los ruidos y las señales involuntarias. Esto significa que, <strong>en</strong> g<strong>en</strong>eral,<br />

la <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia <strong>de</strong> un sistema <strong>de</strong> control <strong>de</strong>be <strong>de</strong> t<strong>en</strong>er un punto <strong>de</strong><br />

corte característico y, algunas veces una banda pasante o no pasante<br />

característica.<br />

Figura 7: Comparación <strong>de</strong> las características <strong>de</strong> magnitud y fase <strong>de</strong> un filtro pasa bajo i<strong>de</strong>al<br />

(a) y <strong>de</strong> un sistema <strong>de</strong> control típico (b)<br />

La característica <strong>de</strong> fase <strong>de</strong> la <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia ti<strong>en</strong>e también su<br />

importancia. La situación i<strong>de</strong>al es que la fase sea una función lineal <strong>de</strong> la<br />

frecu<strong>en</strong>cia d<strong>en</strong>tro <strong>de</strong> la banda <strong>de</strong> la señal <strong>de</strong> <strong>en</strong>trada. La figura 7 repres<strong>en</strong>ta<br />

las características <strong>de</strong> ganancia y <strong>de</strong> fase <strong>de</strong> un filtro pasa bajo i<strong>de</strong>al, imposible<br />

<strong>de</strong> realizar físicam<strong>en</strong>te. Las características típicas <strong>de</strong> amplitud y fase <strong>de</strong> un<br />

sistema <strong>de</strong> control están dibujadas <strong>en</strong> la figura 7b. Se ve que la ganancia<br />

disminuye al crecer la frecu<strong>en</strong>cia. Ello es <strong>de</strong>bido a los efectos <strong>de</strong> las inercias e<br />

inductancias <strong>de</strong> los <strong>sistemas</strong> físicos, por lo que toda <strong>respuesta</strong> cesa cuando la<br />

frecu<strong>en</strong>cia ti<strong>en</strong><strong>de</strong> a infinito.<br />

Las especificaciones más usadas <strong>de</strong> las características <strong>de</strong> un sistema <strong>de</strong><br />

control <strong>en</strong> <strong>el</strong> dominio frecu<strong>en</strong>cial son las sigui<strong>en</strong>tes;<br />

1- Ancho <strong>de</strong> banda. El ancho <strong>de</strong> banda, A. B., se <strong>de</strong>fine como la frecu<strong>en</strong>cia a<br />

la cual <strong>el</strong> módulo <strong>de</strong> ( ) M jω vale <strong>el</strong> 70.7 % <strong>de</strong>l niv<strong>el</strong> a frecu<strong>en</strong>cia cero o 3 dB<br />

por <strong>de</strong>bajo <strong>de</strong>l niv<strong>el</strong> <strong>de</strong> frecu<strong>en</strong>cia nula (figura 8). En g<strong>en</strong>eral, <strong>el</strong> ancho <strong>de</strong><br />

banda indica les características <strong>de</strong> filtraje <strong>de</strong> ruido <strong>de</strong>l sistema. El ancho <strong>de</strong><br />

banda da, también, una mecida <strong>de</strong> las propieda<strong>de</strong>s <strong>de</strong> la <strong>respuesta</strong> transitoria.<br />

Un ancho <strong>de</strong> banda largo indica, normalm<strong>en</strong>te, que las señales <strong>de</strong> alta<br />

frecu<strong>en</strong>cia pasarán a la salida. Es <strong>de</strong>cir, la <strong>respuesta</strong> transitoria <strong>de</strong>be <strong>de</strong> t<strong>en</strong>er<br />

un tiempo <strong>de</strong> subida rápido acompañado <strong>de</strong> un amplio rebase.<br />

- 12 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Por <strong>el</strong> contrario, si <strong>el</strong> ancho <strong>de</strong> banda es chico, solo pasarán las señales <strong>de</strong><br />

baja frecu<strong>en</strong>cia; por consigui<strong>en</strong>te, la <strong>respuesta</strong> temporal será l<strong>en</strong>ta.<br />

2- Factor <strong>de</strong> resonancia ( M<br />

r<br />

). Se <strong>de</strong>fine como <strong>el</strong> valor máximo <strong>de</strong> M ( ω ),<br />

que proporciona a<strong>de</strong>más indicación <strong>de</strong> la estabilidad r<strong>el</strong>ativa <strong>de</strong>l sistema. Si se<br />

recuerda la figura 10 y como se verá más a<strong>de</strong>lante, cuando se estudió <strong>el</strong><br />

diagrama <strong>de</strong> Bo<strong>de</strong> <strong>de</strong> un sistema <strong>de</strong> segundo ord<strong>en</strong> se obti<strong>en</strong><strong>en</strong> distintas curvas<br />

<strong>de</strong> <strong>respuesta</strong> para distintos valores <strong>de</strong> δ .Es evid<strong>en</strong>te que a valores<br />

<strong>el</strong>evados <strong>de</strong> M correspond<strong>en</strong> amplios rebases <strong>de</strong> la <strong>respuesta</strong> temporal.<br />

r<br />

Cuando se proyecta se admite, g<strong>en</strong>eralm<strong>en</strong>te, que <strong>el</strong> valor oprimo <strong>de</strong><br />

estar compr<strong>en</strong>dido <strong>en</strong>tre 1.1 y 1.5.<br />

M<br />

r<br />

<strong>de</strong>be<br />

Figura 8: Curva <strong>de</strong> amplitud <strong>de</strong> un sistema.<br />

3 — Frecu<strong>en</strong>cia <strong>de</strong> resonancia ( ω r<br />

). Es la frecu<strong>en</strong>cia para la que se produce<br />

<strong>el</strong> factor <strong>de</strong> resonancia.<br />

4 - Razón <strong>de</strong> corte. A m<strong>en</strong>udo para frecu<strong>en</strong>cias <strong>el</strong>evadas es importante la<br />

razón <strong>de</strong> corte <strong>de</strong> la curva <strong>de</strong> <strong>respuesta</strong>, <strong>en</strong> frecu<strong>en</strong>cia, puesto que indica la<br />

capacidad <strong>de</strong>l sistema para distinguir la señal <strong>de</strong> ruido. Sin embargo, las<br />

características <strong>de</strong> corte agudo s<strong>en</strong> acompañadas normalm<strong>en</strong>te <strong>de</strong> valores<br />

<strong>el</strong>evados <strong>de</strong> M , lo que significa que <strong>el</strong> sistema es poco estable.<br />

r<br />

5 - Marg<strong>en</strong> <strong>de</strong> amplitud y Marg<strong>en</strong> <strong>de</strong> fase. Estas dos cantida<strong>de</strong>s que son<br />

una medida <strong>de</strong> la estabilidad r<strong>el</strong>ativa <strong>de</strong> un sistema. Serán <strong>de</strong>finidos más<br />

a<strong>de</strong>lante cuando se estudi<strong>en</strong> los diagramas <strong>de</strong> Bo<strong>de</strong> y Nyquist.<br />

- 13 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

4.1- Especificaciones En El Dominio Frecu<strong>en</strong>cial Para Un Sistema De<br />

Segundo Ord<strong>en</strong><br />

En un sistema <strong>de</strong> segundo ord<strong>en</strong>, <strong>el</strong> factor <strong>de</strong> resonancia M<br />

r<br />

y la frecu<strong>en</strong>cia<br />

<strong>de</strong> resonancia ω <strong>de</strong>p<strong>en</strong><strong>de</strong>r, únicam<strong>en</strong>te ri<strong>el</strong> coefici<strong>en</strong>te <strong>de</strong> amortiguami<strong>en</strong>to δ<br />

r<br />

y <strong>de</strong> la frecu<strong>en</strong>cia natural sin amortiguami<strong>en</strong>to ω<br />

n<br />

<strong>de</strong>l sistema. Si se consi<strong>de</strong>ra<br />

la función <strong>de</strong> lazo cerrado <strong>de</strong> segundo ord<strong>en</strong>:<br />

2<br />

ω<br />

n<br />

M ( jω)<br />

= (34)<br />

2<br />

( jω)<br />

+ 2δω<br />

( jω)<br />

+ ω<br />

n<br />

2<br />

n<br />

1<br />

M ( jω)<br />

=<br />

(35)<br />

⎛ ω ⎞ ω<br />

1 + 2 j<br />

⎜ δ −<br />

ω<br />

⎟<br />

⎝ n ⎠ ω<br />

n<br />

Al utilizar la frecu<strong>en</strong>cia reducida<br />

expresión:<br />

ω<br />

u = , <strong>el</strong> módulo <strong>de</strong> M ( jω ) toma la<br />

ω<br />

n<br />

1<br />

M ( u)<br />

= (36a)<br />

2 2<br />

( 1−<br />

u ) + ( 2δu<br />

) 2<br />

Y la fase <strong>de</strong> M ( jω ):<br />

−1<br />

2δu<br />

M ( jω ): φ<br />

m(<br />

u )<br />

= −tg<br />

(36b)<br />

2<br />

1−<br />

u<br />

La frecu<strong>en</strong>cia <strong>de</strong> resonancia se <strong>de</strong>termina <strong>de</strong>rivando M ( u ) con respecto a u e<br />

igualando a cero, es <strong>de</strong>cir:<br />

∂M<br />

u)<br />

1<br />

= −<br />

∂u<br />

2<br />

2<br />

−<br />

2<br />

2<br />

[( ) ( ) ] 3 2 3<br />

1−<br />

u + 2δu<br />

( 4u<br />

− 4u<br />

+ 8uδ<br />

) = 0<br />

( 2<br />

(37)<br />

De don<strong>de</strong>:<br />

3<br />

2<br />

4u − 4u<br />

+ 8uδ<br />

= 0<br />

(38)<br />

Por consigui<strong>en</strong>te:<br />

u = 0<br />

(39)<br />

Y<br />

u =<br />

u r<br />

=<br />

ω<br />

ω<br />

2 r<br />

1 − 2δ =<br />

(40)<br />

n<br />

- 14 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

La solución dada <strong>en</strong> la ecuación 39 indica meram<strong>en</strong>te que la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la<br />

curva M ( ω ) vale cero para ω = 0 0; no es un máximo. De la ecuación 40 se<br />

obti<strong>en</strong>e la frecu<strong>en</strong>cia <strong>de</strong> resonancia que vale:<br />

ω<br />

2<br />

r<br />

= ω n<br />

1−<br />

2δ<br />

(41)<br />

2<br />

Evid<strong>en</strong>tem<strong>en</strong>te, la ecuación 41 es válida solam<strong>en</strong>te para 1 ≥ 2δ<br />

o δ ≤ 0. 707 ,<br />

puesto que <strong>de</strong> otra manera ω r<br />

sería imaginario. Esto significa simplem<strong>en</strong>te<br />

que para todos los valores <strong>de</strong> δ > 0. 707 no hay resonancia (o M<br />

r<br />

= 1) <strong>en</strong> la<br />

curva M (ω)<br />

<strong>en</strong> función <strong>de</strong> ω . La curva M (ω)<br />

es inferior a uno para todos los<br />

valores <strong>de</strong> ω > 0 si <strong>el</strong> coefici<strong>en</strong>te <strong>de</strong> amortiguami<strong>en</strong>to es mayor que 0.707<br />

Sustituy<strong>en</strong>do la ecuación 40 <strong>en</strong> la 36a y simplificando se obti<strong>en</strong>e:<br />

1<br />

M<br />

r<br />

=<br />

(42)<br />

2<br />

2δ 1− δ<br />

Es importante observar que M<br />

r<br />

es función <strong>de</strong> δ solam<strong>en</strong>te, mi<strong>en</strong>tras que ω r<br />

es función <strong>de</strong> ω y δ . Las figuras 9 y 10 repres<strong>en</strong>tan respectivam<strong>en</strong>te las<br />

curvas <strong>de</strong> M<br />

r<br />

y <strong>de</strong> u r<br />

<strong>en</strong> función <strong>de</strong> δ .<br />

Figura 9:<br />

M<br />

r<br />

<strong>en</strong> función <strong>de</strong> δ para un sistema <strong>de</strong> segundo ord<strong>en</strong>.<br />

- 15 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Figura 10:<br />

u<br />

r<br />

<strong>en</strong> función <strong>de</strong> δ para un sistema <strong>de</strong> segundo ord<strong>en</strong>.<br />

En resum<strong>en</strong> se observa que a medida que <strong>el</strong> factor <strong>de</strong> amortiguami<strong>en</strong>to δ<br />

ti<strong>en</strong><strong>de</strong> a cero, la frecu<strong>en</strong>cia <strong>de</strong> resonancia ti<strong>en</strong><strong>de</strong> a ω r<br />

. Para 0 < δ < 0. 707 , la<br />

frecu<strong>en</strong>cia <strong>de</strong> resonancia ω es m<strong>en</strong>or que la frecu<strong>en</strong>cia natural con amortiguami<strong>en</strong>to<br />

ω<br />

n<br />

= ω d<br />

r<br />

2<br />

1−<br />

δ , que aparece <strong>en</strong> la <strong>respuesta</strong> transitoria. Para<br />

δ > 0.707 no hay pico resonante, la amplitud <strong>de</strong> G ( jω)<br />

<strong>de</strong>crece<br />

monótonam<strong>en</strong>te cuando la frecu<strong>en</strong>cia ω crece. Esto significa que no hay pico<br />

<strong>en</strong> la curva <strong>de</strong> <strong>respuesta</strong> para δ > 0. 707 (La amplitud es m<strong>en</strong>or que 0 dB para<br />

todos los valores <strong>de</strong> ω > 0 . Debe recordarse que para 0 .707 < δ < 1, la <strong>respuesta</strong><br />

escalón es oscilatoria, pero las oscilaciones son muy bi<strong>en</strong> amortiguadas y<br />

ap<strong>en</strong>as perceptibles.<br />

Cuando δ ti<strong>en</strong><strong>de</strong> a cero M<br />

r<br />

ti<strong>en</strong><strong>de</strong> a infinito. Esto significa que si <strong>el</strong> sistema<br />

no amortiguado es excitado a su frecu<strong>en</strong>cia natural, la magnitud <strong>de</strong> G ( jω)<br />

se<br />

hace infinito.<br />

Se pue<strong>de</strong> obt<strong>en</strong>er <strong>el</strong> ángulo fase <strong>de</strong> G( jω ), a la frecu<strong>en</strong>cia que se produce <strong>el</strong><br />

pico <strong>de</strong> resonancia.<br />

2<br />

−1 1−<br />

2δ<br />

G( jω)<br />

=− tg<br />

(43)<br />

δ<br />

- 16 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

5 - Gráficas Polares<br />

Como se había m<strong>en</strong>cionado anteriorm<strong>en</strong>te, para varios valores <strong>de</strong> ω se obti<strong>en</strong><strong>en</strong><br />

los valores <strong>de</strong> M (ω)<br />

y φ (ω)<br />

los que podían ser graficados <strong>en</strong> coord<strong>en</strong>adas polares<br />

o rectangulares. Se verá a continuación la primera <strong>de</strong> estas formas <strong>de</strong><br />

repres<strong>en</strong>tación frecu<strong>en</strong>cial.<br />

El diagrama polar <strong>de</strong> una función <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidal G ( jω)<br />

es un diagrama<br />

<strong>de</strong> la amplitud <strong>de</strong> G ( jω)<br />

<strong>en</strong> función <strong>de</strong>l ángulo <strong>de</strong> G ( jω)<br />

<strong>en</strong> coord<strong>en</strong>adas polares al<br />

variar ω <strong>de</strong>s<strong>de</strong> cero a infinito. Se hace notar que <strong>en</strong> los diagramas polares, se<br />

mi<strong>de</strong> un ángulo '<strong>de</strong> fase positivo (negativo) <strong>en</strong> s<strong>en</strong>tido antihorario (<strong>en</strong> s<strong>en</strong>tido<br />

horario) <strong>de</strong>s<strong>de</strong> <strong>el</strong> eje positivo real. El diagrama polar frecu<strong>en</strong>tem<strong>en</strong>te recibe <strong>el</strong><br />

nombre <strong>de</strong> diagrama <strong>de</strong> Nyquist. En la figura 11 hay un ejemplo <strong>de</strong> ese diagrama.<br />

Cada punto <strong>de</strong>l diagrama polar <strong>de</strong> G ( jω)<br />

repres<strong>en</strong>ta <strong>el</strong> punto terminal <strong>de</strong> un<br />

vector pera un valor <strong>de</strong>terminado <strong>de</strong> ω . En <strong>el</strong> diagrama polar es importante<br />

mostrar la graduación <strong>de</strong> frecu<strong>en</strong>cia sobre <strong>el</strong> diagrama. Las proyecciones <strong>de</strong> G ( jω)<br />

<strong>en</strong> los ejes real e imaginario son sus compon<strong>en</strong>tes real e imaginario. Tanto la<br />

amplitud G ( jω)<br />

como <strong>el</strong> ángulo <strong>de</strong> fase G( jω ) <strong>de</strong>b<strong>en</strong> ser calculados directam<strong>en</strong>te<br />

para cada frecu<strong>en</strong>cia ω para pe<strong>de</strong>r construir los diagramas polares. Sin embargo,<br />

como <strong>el</strong> diagrama rectangular logarítmico que se verá posteriorm<strong>en</strong>te es fácil <strong>de</strong><br />

construir, la información necesaria para trazar <strong>el</strong> diagrama polar pue<strong>de</strong> ser<br />

obt<strong>en</strong>ida directam<strong>en</strong>te <strong>de</strong>l diagrama logarítmico si se dibuja previam<strong>en</strong>te aquél y<br />

se convierte <strong>de</strong>cib<strong>el</strong>es <strong>en</strong> una magnitud ordinaria usando la figura 33.<br />

Figura 11.<br />

- 17 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Pera dos <strong>sistemas</strong> conectados <strong>en</strong> cascada la función <strong>de</strong> transfer<strong>en</strong>cia total <strong>de</strong> la<br />

combinación, <strong>en</strong> aus<strong>en</strong>cia <strong>de</strong> efectos <strong>de</strong> carga, es <strong>el</strong> producto <strong>de</strong> las dos<br />

funciones <strong>de</strong> transfer<strong>en</strong>cia individuales. Si se necesita la multiplicación <strong>de</strong> dos<br />

funciones <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidales, esta pue<strong>de</strong> lograrse multiplicando las funciones<br />

<strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidales <strong>en</strong> cada frecu<strong>en</strong>cia realizando la multiplicación<br />

con álgebra compleja.<br />

Es <strong>de</strong>cir, si:<br />

G ( jω)<br />

= G1 ( jω)<br />

G2<br />

( jω)<br />

Entonces:<br />

G( jω) = G( jω) G( jω)<br />

(44)<br />

Don<strong>de</strong>:<br />

G<br />

jω)<br />

= G ( jω)<br />

− G ( j )<br />

(45a)<br />

(<br />

1 2<br />

ω<br />

G( jω) = G ( jω) + G ( jω)<br />

(45b)<br />

1 2<br />

En la figura 12 se muestra <strong>el</strong> producto <strong>de</strong> G jω).<br />

G ( j )<br />

1( 2<br />

ω<br />

Figura 12.<br />

En g<strong>en</strong>eral, si se <strong>de</strong>sea un diagrama polar <strong>de</strong> G<br />

1( jω).<br />

G2<br />

( jω)<br />

es conv<strong>en</strong>i<strong>en</strong>te trazar<br />

un diagrama logarítmico <strong>de</strong> G jω).<br />

G ( j ) y luego convertirla <strong>en</strong> un diagrama polar<br />

1( 2<br />

ω<br />

- 18 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

<strong>en</strong> lugar <strong>de</strong> dibujar los diagramas polares G ( j ) y <strong>de</strong> G ( ) 2<br />

jω multiplicando estos<br />

dos <strong>en</strong> <strong>el</strong> plano complejo.<br />

1<br />

ω<br />

Una v<strong>en</strong>taja <strong>de</strong> usar un diagrama polar es que pres<strong>en</strong>ta las características <strong>de</strong><br />

<strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> un sistema <strong>en</strong> todo <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cias <strong>en</strong> un único<br />

diagrama. Una <strong>de</strong>sv<strong>en</strong>taja es que <strong>el</strong> diagrama no indica claram<strong>en</strong>te las<br />

contribuciones <strong>de</strong> cada uno <strong>de</strong> los factores individuales <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia<br />

<strong>de</strong> lazo abierto.<br />

5.1- Factores Integral Y Derivativo<br />

( ω)<br />

± 1<br />

j .<br />

El diagrama polar <strong>de</strong> G ( jω)<br />

= 1 es <strong>el</strong> eje imaginario negativo, pues:<br />

1 1 1<br />

G( jω) = j 90º<br />

jω =− ω = ω<br />

(46)<br />

El diagrama polar <strong>de</strong><br />

G ( jω)<br />

= jω<br />

es <strong>el</strong> eje positivo imaginario, pues:<br />

G( jω) = jω = ω 90º<br />

(47)<br />

±<br />

5.2- Factores De Primer Ord<strong>en</strong> ( ) 1<br />

1+ jω T<br />

Para la función <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidal:<br />

1 1<br />

1+ jωT<br />

1+<br />

ω T<br />

−1<br />

G( jω)<br />

= = −tg ωT<br />

2 2<br />

(48)<br />

Los valores <strong>de</strong> G ( jω)<br />

<strong>en</strong> ω = 0 y<br />

ω = 1 son respectivam<strong>en</strong>te:<br />

T<br />

G( j 0) = 1 0º<br />

(49)<br />

Y<br />

⎛ 1 ⎞ 1<br />

G⎜<br />

j ⎟ =<br />

⎝ T ⎠ 2<br />

45º<br />

Al t<strong>en</strong><strong>de</strong>r ω a infinito:<br />

(50)<br />

G( j∞ ) = 0 − 90º<br />

(51)<br />

El diagrama polar <strong>de</strong> esta función <strong>de</strong> transfer<strong>en</strong>cia es una semicircunfer<strong>en</strong>cia al<br />

variar la frecu<strong>en</strong>cia <strong>de</strong>s<strong>de</strong> cero a infinito, como se ve <strong>en</strong> la figura 13a. El c<strong>en</strong>tro<br />

está ubicado <strong>en</strong> 0.5 sobre <strong>el</strong> eje real y <strong>el</strong> radio es igual a 0.5.<br />

- 19 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Figura 13. a: Diagrama polar <strong>de</strong><br />

1<br />

(1<br />

+<br />

jωT )<br />

, b: Diagrama <strong>de</strong> G ( jω)<br />

<strong>en</strong> <strong>el</strong> plano x-y.<br />

Para probar que <strong>el</strong> diagrama polar es una semicircunfer<strong>en</strong>cia se <strong>de</strong>fine:<br />

G ( jω ) = X + jY<br />

(52)<br />

Don<strong>de</strong>:<br />

X<br />

1<br />

2<br />

1+<br />

ω T<br />

= 2<br />

=<br />

Parte real <strong>de</strong> G ( jω)<br />

(53)<br />

Y<br />

− ωT<br />

2<br />

1+<br />

ω T<br />

=<br />

2<br />

=<br />

Parte imaginaria <strong>de</strong> G ( jω)<br />

(54)<br />

Entonces se obti<strong>en</strong>e:<br />

⎛<br />

⎜ X<br />

⎝<br />

1 ⎞<br />

− ⎟<br />

2 ⎠<br />

2<br />

+ Y<br />

2<br />

⎛ 1 1−<br />

ω T<br />

=<br />

⎜<br />

2<br />

⎝ 2 1+<br />

ω T<br />

2<br />

2<br />

⎞<br />

⎟<br />

⎠<br />

2<br />

⎛ − ωT<br />

+ ⎜<br />

2<br />

⎝1+<br />

ω T<br />

2<br />

2<br />

⎞<br />

⎟<br />

⎠<br />

⎛ 1 ⎞<br />

= ⎜ ⎟<br />

⎝ 2 ⎠<br />

2<br />

(55)<br />

- 20 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Entonces <strong>en</strong> <strong>el</strong> plano x-y, G ( jω)<br />

es un círculo con c<strong>en</strong>tro <strong>en</strong> X=1/2; Y=0 y<br />

con un radio igual a 1/2 como pue<strong>de</strong> verse <strong>en</strong> la figura 13b. La<br />

semicircunfer<strong>en</strong>cia inferior correspon<strong>de</strong> a 0 ≤ ω ≤ ∞ y la semicircunfer<strong>en</strong>cia<br />

superior correspon<strong>de</strong> a − ∞ ≤ ω ≤ 0 .<br />

El diagrama polar <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia 1 + jωT<br />

es simplem<strong>en</strong>te la<br />

mitad superior <strong>de</strong> la línea recta que pasa por <strong>el</strong> punto (1,0) <strong>en</strong> <strong>el</strong> plano<br />

complejo y paral<strong>el</strong>o al eje imaginario, como pue<strong>de</strong> verse <strong>en</strong> la figura 14. El<br />

Diagrama polar <strong>de</strong> 1 + jωT<br />

ti<strong>en</strong>e un aspecto totalm<strong>en</strong>te difer<strong>en</strong>te al <strong>de</strong><br />

−<br />

( 1+ T ) 1<br />

jω .<br />

Figura 14: Diagrama polar <strong>de</strong><br />

1 + jωT<br />

Si:<br />

2 2 1<br />

G( jω) 1 jωT ω T tg − ωT<br />

= + = 1+ (56)<br />

Para ω = 0 : G( j 0) = 1 0º<br />

(57)<br />

Para<br />

ω = ∞ : G( j∞ ) =∞ 90º<br />

(58)<br />

5.3- Factores Cuadráticos<br />

2<br />

⎡ ⎛ ω ⎞ ⎛ jω<br />

⎞ ⎤<br />

⎢1<br />

+ 2δ<br />

j<br />

⎜<br />

⎟ +<br />

⎜<br />

⎟ ⎥<br />

⎢⎣<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠ ⎥⎦<br />

Las porciones <strong>de</strong> alta y baja frecu<strong>en</strong>cia <strong>de</strong>l diagrama polar <strong>de</strong> la función<br />

<strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidal sigui<strong>en</strong>te:<br />

± 1<br />

- 21 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

1<br />

G( jω) = δ > 0<br />

2<br />

⎛ ω ⎞ ⎛ ω ⎞<br />

1+ 2δ<br />

⎜ j ⎟+<br />

⎜ j ⎟<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠<br />

Están dadas respectivam<strong>en</strong>te por:<br />

(59)<br />

lim G( jω) = 1 0º y lim G( jω) = 0 − 180º<br />

ω→0<br />

ω→∞<br />

El diagrama polar <strong>de</strong> esta función <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidal comi<strong>en</strong>za <strong>en</strong><br />

1 0º y finaliza <strong>en</strong> 0 180º − al aum<strong>en</strong>tar ω <strong>de</strong>s<strong>de</strong> cero a infinito.<br />

Así, la porción alta frecu<strong>en</strong>cia <strong>de</strong> G( jω ) es tang<strong>en</strong>te al eje real negativo.<br />

Figura 15.<br />

Los valores <strong>de</strong> G( jω ) <strong>en</strong> <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cias <strong>de</strong> interés pued<strong>en</strong> ser<br />

catalogados directam<strong>en</strong>te.<br />

En la figura 15, hay ejemplos <strong>de</strong> diagramas polares <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia<br />

recién analizada. La forma exacta <strong>de</strong> un diagrama polar <strong>de</strong>p<strong>en</strong><strong>de</strong> <strong>de</strong>l<br />

valor <strong>de</strong> la r<strong>el</strong>ación <strong>de</strong> amortiguami<strong>en</strong>to δ , pues la forma g<strong>en</strong>eral es la<br />

misma, tanto para <strong>el</strong> caso subamortiguado ( 0< δ < 1) como para <strong>el</strong> caso<br />

sobreamortiguado ( δ > 1).<br />

Para <strong>el</strong> caso subamortiguado <strong>en</strong> ω = ωn<br />

se ti<strong>en</strong>e G( jω)<br />

= 1<br />

j2δ<br />

y <strong>el</strong> ángulo <strong>de</strong><br />

fase es - 90°. Por tanto, se pue<strong>de</strong> ver que la frecu<strong>en</strong>cia a la cual <strong>el</strong> lugar<br />

G( jω ) intersecta al eje imaginario es la frecu<strong>en</strong>cia natural no amortigua da<br />

ω<br />

n<br />

. En <strong>el</strong> diagrama polar, <strong>el</strong> punto <strong>de</strong> frecu<strong>en</strong>cia cuya distancia <strong>de</strong>s<strong>de</strong> <strong>el</strong><br />

- 22 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

orig<strong>en</strong> es máxima, correspon<strong>de</strong> a la frecu<strong>en</strong>cia <strong>de</strong> resonancia ω<br />

r<br />

. Se obti<strong>en</strong>e<br />

<strong>el</strong> valor pico <strong>de</strong> G( jω ) como la r<strong>el</strong>ación <strong>en</strong>tre <strong>el</strong> módulo <strong>de</strong>l vector a la<br />

frecu<strong>en</strong>cia <strong>de</strong> resonancia y <strong>el</strong> módulo <strong>de</strong>l vector <strong>en</strong> ω = 0 . Se indica la<br />

frecu<strong>en</strong>cia <strong>de</strong> resonancia <strong>en</strong> <strong>el</strong> diagrama polar como pue<strong>de</strong> verse <strong>en</strong> la<br />

figura 16.<br />

Figura 16: Diagrama polar que muestra <strong>el</strong> pico <strong>de</strong> resonancia y la frecu<strong>en</strong>cia <strong>de</strong> resonancia<br />

ω<br />

Para <strong>el</strong> caso sobreamortiguado, cuando δ es bastante mayor que la unidad<br />

<strong>el</strong> lugar <strong>de</strong> G( jω ) se aproxima a una semicircunfer<strong>en</strong>cia. Se pue<strong>de</strong> ver esto<br />

<strong>de</strong>l hecho <strong>de</strong> que para un sistema fuertem<strong>en</strong>te amortiguado, las raíces<br />

características son reales y es una mucho más pequeña que la otra. Como<br />

para un valor <strong>de</strong> δ sufici<strong>en</strong>tem<strong>en</strong>te gran<strong>de</strong> <strong>el</strong> efecto <strong>de</strong> la raíz más gran<strong>de</strong><br />

<strong>en</strong> la <strong>respuesta</strong> se hace muy pequeña, <strong>el</strong> sistema se comporta como un<br />

sistema <strong>de</strong> primer ord<strong>en</strong>.<br />

Para la función <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidal:<br />

r<br />

⎛ ω ⎞ ⎛ ω ⎞<br />

G( jω) = 1+ 2δ<br />

⎜ j ⎟+<br />

⎜ j ⎟<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠<br />

2<br />

⎛ ω ⎞ ⎛2δω<br />

⎞<br />

∴ G( jω) = ⎜1− j<br />

2 ⎟+<br />

⎜ ⎟<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠<br />

2<br />

(60a)<br />

(60b)<br />

La porción <strong>de</strong> baja frecu<strong>en</strong>cia <strong>de</strong> la curva es:<br />

lim G( jω) = 1 0º<br />

(61)<br />

ω→0<br />

- 23 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Y la porción <strong>de</strong> alta frecu<strong>en</strong>cia es:<br />

lim G( jω) =∞ 180º<br />

(62)<br />

ω→∞<br />

Como la parte imaginaria <strong>de</strong> G( jω ) es positiva para ω > 0 y es<br />

monótonam<strong>en</strong>te creci<strong>en</strong>te; la parte real <strong>de</strong> G( jω ) es monótonam<strong>en</strong>te<br />

<strong>de</strong>creci<strong>en</strong>te <strong>de</strong>s<strong>de</strong> la unidad, la forma g<strong>en</strong>eral <strong>de</strong>l diagrama polar <strong>de</strong> G( jω )<br />

<strong>de</strong> la ecuación 60a es como pue<strong>de</strong> verse <strong>en</strong> la figura 17.<br />

El ángulo <strong>de</strong> fase está <strong>en</strong>tre 0 o y 180°.<br />

Figura 17: Diagrama polar <strong>de</strong><br />

⎛ ω ⎞ ⎛ ω ⎞<br />

G( jω) = 1+ 2δ<br />

⎜ j ⎟+<br />

⎜ j ⎟<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠<br />

2<br />

Ejemplo 1: Sea la sigui<strong>en</strong>te función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> segundo ord<strong>en</strong>:<br />

1<br />

Gs () =<br />

sTs<br />

( + 1)<br />

Trazar un diagrama polar para esta función <strong>de</strong> transfer<strong>en</strong>cia.<br />

En primer lugar se escribe la función <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidal:<br />

1 T<br />

1<br />

G( jω)<br />

= =− − j<br />

( jω) jωT 1 1 ω T ω 1 ω T<br />

2 2 2 2<br />

( + ) + ( + )<br />

El tramo <strong>de</strong> baja frecu<strong>en</strong>cia <strong>de</strong> la curva es:<br />

lim G( jω) =−T − j∞=∞ − 90º<br />

ω→0<br />

El tramo <strong>de</strong> alta frecu<strong>en</strong>cia <strong>de</strong> la curva es:<br />

lim G( jω) = 0 − j0 = 0 − 180º<br />

ω→∞<br />

- 24 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

En la figura 18 se muestra la forma g<strong>en</strong>eral <strong>de</strong>l diagrama polar G( jω ).El<br />

diagrama <strong>de</strong> G( jω ) es asintótico a la recta vertical que pasa por <strong>el</strong> punto (-<br />

T, 0). Como esta función <strong>de</strong> transfer<strong>en</strong>cia involucra integración ( 1 s<br />

), la<br />

forma <strong>de</strong>l diagrama polar difiere sustancialm<strong>en</strong>te <strong>de</strong> las funciones <strong>de</strong><br />

transfer<strong>en</strong>cia que no ti<strong>en</strong><strong>en</strong> integración.<br />

5.4- Retardo De Transporte.<br />

Figura 18.<br />

Algunos <strong>el</strong>em<strong>en</strong>tos <strong>en</strong> los <strong>sistemas</strong> <strong>de</strong> control se caracterizan por su tiempo<br />

muerto o retraso <strong>de</strong> transporte, durante <strong>el</strong> cual no dan salida a la señal <strong>de</strong><br />

<strong>en</strong>trada que se les haya aplicado; hay un periodo muerto <strong>de</strong> inactividad que<br />

<strong>de</strong>mora la transmisión <strong>de</strong> la señal recibida. La figura 19 muestra <strong>el</strong><br />

diagrama <strong>de</strong> bloques <strong>de</strong>l <strong>el</strong>em<strong>en</strong>to <strong>de</strong>morador.<br />

Figura 19.<br />

La figura 20 indica la salida ct () = rt ( −T). ut ( − T)<br />

que ti<strong>en</strong>e la misma forma<br />

que la <strong>en</strong>trada, pero <strong>de</strong>mora un tiempo T. La <strong>de</strong>mora es una característica<br />

no lineal que, afortunadam<strong>en</strong>te, pue<strong>de</strong> ser repres<strong>en</strong>tada por su<br />

transformada <strong>de</strong> Laplace. Así:<br />

Cs () −st<br />

Gs () = = e<br />

(63)<br />

Rs ()<br />

- 25 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Figura 20.<br />

− j T<br />

∴ G( jω)<br />

= e ω<br />

(64)<br />

Pue<strong>de</strong> ser escrito:<br />

G( jω) = 1 cosωT − js<strong>en</strong>ωT<br />

(65)<br />

Como <strong>el</strong> módulo G( jω ) es simple la unidad y <strong>el</strong> ángulo <strong>de</strong> fase varía<br />

linealm<strong>en</strong>te con ω , <strong>el</strong> diagrama polar <strong>de</strong>l retardo <strong>de</strong> transporte es un círculo<br />

<strong>de</strong> radio unitario, como pue<strong>de</strong> verse <strong>en</strong> la figura 21.<br />

Figura 21: Diagrama polar <strong>de</strong>l retardo <strong>de</strong> transporte.<br />

- 26 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

A frecu<strong>en</strong>cias bajas, <strong>el</strong> retardo <strong>de</strong> transporte<br />

− j T y <strong>el</strong> retardo <strong>de</strong> primer<br />

ord<strong>en</strong> 1(1 + jωT<br />

) se comportan <strong>en</strong> forma similar, lo que pue<strong>de</strong> verse <strong>en</strong> la<br />

figura 22.<br />

e ω<br />

Figura 22: Diagramas polares <strong>de</strong><br />

j T<br />

e ω<br />

− y 1 (1 + jωT<br />

)<br />

j T<br />

e ω<br />

Los diagramas polares <strong>de</strong><br />

− y 1 (1 + jωT<br />

) son tang<strong>en</strong>tes <strong>en</strong>tre sí <strong>en</strong> ω = 0 .<br />

Pue<strong>de</strong> verse esto <strong>de</strong>l hecho que para ω 1 T :<br />

− jωT<br />

1<br />

e = 1− jωT<br />

y 1 jωT<br />

1+<br />

jω<br />

T<br />

= −<br />

Sin embargo, para ω 1 T ; hay una difer<strong>en</strong>cia es<strong>en</strong>cial <strong>en</strong>tre<br />

1(1 + jωT<br />

) como pue<strong>de</strong> verse también <strong>en</strong> la figura 22.<br />

j T<br />

e − ω<br />

Ejemplo 2: Obt<strong>en</strong>er <strong>el</strong> diagrama polar <strong>de</strong> la sigui<strong>en</strong>te función <strong>de</strong><br />

transfer<strong>en</strong>cia<br />

− jωL<br />

e<br />

G( jω)<br />

=<br />

1 + jωT<br />

Se pue<strong>de</strong> escribir:<br />

− jωL⎛ 1 ⎞<br />

G( jω)<br />

= e ⎜ ⎟<br />

⎝1+<br />

jωT<br />

⎠<br />

y<br />

- 27 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

El módulo y <strong>el</strong> ángulo <strong>de</strong> fase son respectivam<strong>en</strong>te:<br />

− jωL<br />

1 1<br />

G( jω) = e . = 1 + jωT<br />

1+<br />

ω T<br />

2 2<br />

Y<br />

1<br />

( ) j ωL<br />

G jω e −<br />

−<br />

= + 11+ jωT = ωL−<br />

tg ωT<br />

Figura 23.<br />

Como <strong>el</strong> módulo disminuye <strong>de</strong>s<strong>de</strong> la unidad <strong>en</strong> forma monótona y <strong>el</strong> ángulo<br />

<strong>de</strong> fase también disminuye monótona e in<strong>de</strong>finidam<strong>en</strong>te, <strong>el</strong> diagrama polar<br />

<strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia dada es una espiral, como pue<strong>de</strong> verse <strong>en</strong> la<br />

figura 23.<br />

5.5- Formas G<strong>en</strong>erales De Los Diagramas Polares<br />

Los diagramas polares <strong>de</strong> una función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> la forma<br />

G( jω)<br />

=<br />

G<br />

K( + jωTa)( + jωTb)<br />

( jω) λ<br />

( + jωT )( + jωT<br />

)<br />

1 1 .....<br />

1 1 ....<br />

1 2<br />

m<br />

m−1<br />

( jω) + b1<br />

( jω)<br />

+<br />

n<br />

n<br />

( jω) + a ( jω) + ...<br />

(66)<br />

0<br />

....<br />

) = b<br />

jω (67)<br />

a<br />

−<br />

(<br />

1<br />

0<br />

1<br />

Don<strong>de</strong> <strong>el</strong> grado <strong>de</strong>l polinomio d<strong>en</strong>ominador es mayor que <strong>el</strong> <strong>de</strong>l<br />

d<strong>en</strong>ominador, t<strong>en</strong>drá las sigui<strong>en</strong>tes formas g<strong>en</strong>erales:<br />

- 28 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

1 - Para λ = 0 o <strong>sistemas</strong> <strong>de</strong> tipo 0: <strong>el</strong> punto <strong>de</strong> iniciación <strong>de</strong>l diagrama<br />

polar (que correspon<strong>de</strong> a ω = 0 ) es finito y está sobre <strong>el</strong> eje positivo real.<br />

La tang<strong>en</strong>te al diagrama polar <strong>en</strong> ω = 0 es perp<strong>en</strong>dicular al eje real. El punto<br />

terminal que correspon<strong>de</strong> a ω =∞, está <strong>en</strong> <strong>el</strong> orig<strong>en</strong> y la curva converge al<br />

orig<strong>en</strong> y es tang<strong>en</strong>te a uno <strong>de</strong> los ejes.<br />

2 - Para λ = 1 o <strong>sistemas</strong> tipo 1: <strong>el</strong> término jω <strong>en</strong> <strong>el</strong> d<strong>en</strong>ominador<br />

contribuye -90° al ángulo <strong>de</strong> fase total <strong>de</strong> G( jω ) para 0 ≤ ω ≤∞. En ω = 0 , <strong>el</strong><br />

módulo <strong>de</strong> G( jω ) es infinito y <strong>el</strong> ángulo <strong>de</strong> fase es -90°. A frecu<strong>en</strong>cias bajas,<br />

<strong>el</strong> diagrama polar es asintótico a una línea paral<strong>el</strong>a al eje imaginario<br />

negativo. En ω =∞, <strong>el</strong> módulo se hace cero y la curva converge al orig<strong>en</strong> y<br />

es tang<strong>en</strong>te a uno <strong>de</strong> los ejes.<br />

3 - Para λ = 2 o <strong>sistemas</strong> tipo 2: <strong>el</strong> término jω <strong>en</strong> <strong>el</strong> d<strong>en</strong>ominador<br />

contribuye -180° al ángulo <strong>de</strong> fase total <strong>de</strong> G( jω ) para 0 ≤ ω ≤∞. En ω = 0 , la<br />

amplitud <strong>de</strong> G( jω ) es infinita y <strong>el</strong> ángulo <strong>de</strong> fase es igual a -180°. En bajas<br />

frecu<strong>en</strong>cias, <strong>el</strong> diagrama polar es asintótico a una línea recta paral<strong>el</strong>a al eje<br />

real negativo. En ω =∞, <strong>el</strong> módulo se hace cero y la curva es tang<strong>en</strong>te a uno<br />

<strong>de</strong> los ejes.<br />

En la figura 24, se v<strong>en</strong> las formas g<strong>en</strong>erales <strong>de</strong> las porciones <strong>de</strong> baja<br />

frecu<strong>en</strong>cia <strong>de</strong> los diagramas polares <strong>de</strong> los <strong>sistemas</strong> tipo 0, tipo 1 y tipo 2.<br />

Se pue<strong>de</strong> ver que si <strong>el</strong> arado <strong>de</strong>l polinomio d<strong>en</strong>ominador <strong>de</strong> G( jω ) es mayor<br />

que <strong>el</strong> <strong>de</strong>l numerador, los lugares <strong>de</strong> G( jω ) converg<strong>en</strong> al orig<strong>en</strong> <strong>en</strong> s<strong>en</strong>tido<br />

horario. En ω =∞, los lugares son tang<strong>en</strong>tes a uno u otro <strong>de</strong> los ejes, como<br />

se ve <strong>en</strong> la figura 25.<br />

Figura 24: Diagramas polares <strong>de</strong> <strong>sistemas</strong> tipo 0, tipo 1 y tipo 2.<br />

- 29 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Para <strong>el</strong> caso <strong>en</strong> que los grados <strong>de</strong>l polinomio numerador y d<strong>en</strong>ominador <strong>de</strong><br />

G( jω ) son iguales, <strong>el</strong> diagrama polar comi<strong>en</strong>za a una distancia finita sobre <strong>el</strong><br />

eje real y finaliza <strong>en</strong> un punto sobre <strong>el</strong> eje real.<br />

Figura 25: Diagramas polares <strong>de</strong> funciones <strong>de</strong> transfer<strong>en</strong>cia con dinámica <strong>de</strong> numerador.<br />

Se hace notar que cualquier forma complicada <strong>de</strong> las curvas <strong>de</strong> los<br />

diagramas polares es causada por la dinámica <strong>de</strong>l numerador, es <strong>de</strong>cir, por<br />

las constantes <strong>de</strong> tiempo <strong>de</strong>l numerador <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia. En<br />

la figura 26 se v<strong>en</strong> ejemplos <strong>de</strong> diagramas polares <strong>de</strong> funciones <strong>de</strong><br />

transfer<strong>en</strong>cia con dinámica <strong>de</strong> numerador.<br />

Figura 26: Diagramas polares <strong>de</strong> funciones <strong>de</strong> transfer<strong>en</strong>cia con dinámica <strong>de</strong> numerador.<br />

- 30 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Al analizar <strong>sistemas</strong> <strong>de</strong> control, hay que <strong>de</strong>terminar con exactitud <strong>el</strong><br />

diagrama polar <strong>de</strong> G( jω ) <strong>en</strong> <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cias <strong>de</strong> interés.<br />

En la tabla 1 se adjuntan distintos diagramas polares para diversas funciones<br />

<strong>de</strong> transfer<strong>en</strong>cia.<br />

- 31 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Factor<br />

Lugar Polar <strong>de</strong> la Respuesta<br />

a<br />

1<br />

1+<br />

Ts<br />

b 1+<br />

Ts<br />

c<br />

1<br />

Ts , 1<br />

2 2<br />

T s<br />

d Ts ,<br />

2 2<br />

T s<br />

1<br />

e ( 1+<br />

Ts) 2<br />

como g, con δ = 1<br />

1<br />

f ( 1+ Ts)( 1+<br />

Ts)<br />

1 2<br />

como g, con δ > 1<br />

g<br />

1<br />

1+ 2Ts<br />

+ T s<br />

1<br />

T =<br />

ω<br />

n<br />

2 2<br />

Tabla 1.<br />

- 32 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

5.6 - El Efecto De Factores ( 1 + Ts)<br />

En El Numerador De La Función De<br />

Transfer<strong>en</strong>cia.<br />

La curva polar <strong>de</strong> la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> una función <strong>de</strong> transfer<strong>en</strong>cia<br />

como la <strong>de</strong> la figura 27, está situada totalm<strong>en</strong>te <strong>en</strong> <strong>el</strong> cuarto cuadrante, con<br />

un atraso <strong>de</strong> fase que se aproxima a 90° con frecu<strong>en</strong>cias <strong>de</strong> <strong>en</strong>trada muy<br />

<strong>el</strong>evadas.<br />

Figura 27: Repres<strong>en</strong>tación polar <strong>de</strong><br />

K<br />

1+ Ts<br />

La gráfica para un simple sistema <strong>de</strong> segundo ord<strong>en</strong>, figura 28, ocupa <strong>el</strong> tercero<br />

y cuarto cuadrante, con un ángulo <strong>de</strong> fase que se aproxima a -180°<br />

para frecu<strong>en</strong>cias <strong>de</strong> <strong>en</strong>trada muy <strong>el</strong>evadas.<br />

Figura 28: Repres<strong>en</strong>tación polar <strong>de</strong> un sistema <strong>de</strong> segundo ord<strong>en</strong>.<br />

El lugar <strong>de</strong> la <strong>respuesta</strong> a la frecu<strong>en</strong>cia <strong>de</strong> un sistema simple <strong>de</strong> tercer ord<strong>en</strong><br />

cuya función <strong>de</strong> transfer<strong>en</strong>cia sea:<br />

C(<br />

s)<br />

R(<br />

s)<br />

1<br />

= H ( s)<br />

=<br />

(68)<br />

( 1+<br />

T s)( 1+<br />

T s)( + T s)<br />

1 2<br />

1<br />

3<br />

O<br />

- 33 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

C(<br />

s)<br />

R(<br />

s)<br />

1<br />

= H ( s)<br />

=<br />

(69)<br />

2 2<br />

( 1+<br />

T s)( 1+<br />

T s + T s )<br />

1<br />

2<br />

2<br />

2<br />

Ocuparía <strong>el</strong> segundo, tercero y cuarto con un ángulo <strong>de</strong> atraso que se<br />

aproxima a -270° para frecu<strong>en</strong>cias <strong>de</strong> <strong>en</strong>trada muy <strong>el</strong>evadas como se ilustra<br />

<strong>en</strong> la figura 29a.<br />

Figura 29.<br />

Esto indica que cada factor ( 1 + Ts)<br />

<strong>en</strong> <strong>el</strong> d<strong>en</strong>ominador <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia<br />

contribuye con un atraso <strong>de</strong> fase que se aproxima a 90° para altas<br />

frecu<strong>en</strong>cias, llevando a la gráfica a ocupar un cuadrante adicional <strong>en</strong> <strong>el</strong><br />

s<strong>en</strong>tido horario.<br />

Se vio que un factor ( 1 + Ts)<br />

<strong>en</strong> <strong>el</strong> numerador, expresado <strong>en</strong> forma polar,<br />

ti<strong>en</strong>e un ángulo <strong>de</strong> fase positivo asociado a él, por ejemplo:<br />

C(<br />

s)<br />

C(<br />

jω)<br />

= 1+<br />

Ts ⇒ = 1+<br />

R(<br />

s)<br />

R(<br />

jω)<br />

jωT<br />

Que ti<strong>en</strong>e una r<strong>el</strong>ación <strong>de</strong> magnitud y fase:<br />

M ( jω ) + T<br />

(70)<br />

2 2<br />

= 1 ω<br />

(71)<br />

−1<br />

φ(<br />

ω)<br />

= tg ωT<br />

(72)<br />

Y por lo tanto su característica <strong>de</strong> frecu<strong>en</strong>cia ocupa <strong>el</strong> primer cuadrante <strong>de</strong>l<br />

gráfico polar, como se ve <strong>en</strong> la figura 29b.<br />

La fase positiva <strong>de</strong> un factor ( 1 + Ts)<br />

es llamada fase <strong>de</strong> a<strong>de</strong>lanto, y se<br />

aproxima a 90° para altas frecu<strong>en</strong>cias.<br />

La adición <strong>de</strong> un factor <strong>de</strong> a<strong>de</strong>lanto <strong>de</strong> fase ( 1 + Ts)<br />

a una función <strong>de</strong><br />

transfer<strong>en</strong>cia cont<strong>en</strong>i<strong>en</strong>do factores <strong>de</strong> atraso <strong>de</strong> fase, afecta la r<strong>el</strong>ación <strong>de</strong><br />

magnitud <strong>de</strong> la misma, y más significativam<strong>en</strong>te, la fase.<br />

- 34 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

El efecto es disminuir <strong>el</strong> atraso <strong>de</strong> fase <strong>de</strong>bido a los factores <strong>en</strong> <strong>el</strong><br />

d<strong>en</strong>ominador. Esta acción pue<strong>de</strong> ser ilustrada consi<strong>de</strong>rando la función <strong>de</strong><br />

transfer<strong>en</strong>cia:<br />

( 1+<br />

T1<br />

s)<br />

( 1+<br />

T s)( 1+<br />

T s)<br />

C(<br />

s)<br />

G(<br />

s)<br />

= =<br />

(73)<br />

KR(<br />

s)<br />

2<br />

3<br />

Si <strong>el</strong> factor ( 1+ T 1<br />

s)<br />

no figurara, <strong>el</strong> lugar <strong>de</strong> la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia<br />

ocuparía <strong>el</strong> tercer y cuarto cuadrante. La adición <strong>de</strong> ( 1+ T 1<br />

s)<br />

<strong>en</strong> <strong>el</strong> numerador<br />

cambia la r<strong>el</strong>ación <strong>de</strong> magnitud, y reduce la fase para una frecu<strong>en</strong>cia<br />

particular ω <strong>de</strong> acuerdo a:<br />

φ( ω)<br />

= φ1(<br />

ω)<br />

− φ2<br />

( ω)<br />

− φ3(<br />

ω)<br />

Don<strong>de</strong>:<br />

φ ( ω)<br />

= tg<br />

1<br />

φ ( ω)<br />

= tg<br />

2<br />

φ ( ω)<br />

= tg<br />

3<br />

−1<br />

ωT<br />

−1<br />

−1<br />

1<br />

ωT<br />

ωT<br />

2<br />

31<br />

Por lo que <strong>el</strong> factor <strong>en</strong> <strong>el</strong> numerador ti<strong>en</strong><strong>de</strong> a restringir <strong>el</strong> lugar al cuarto<br />

cuadrante solam<strong>en</strong>te. Los valores r<strong>el</strong>ativos <strong>de</strong> T 1<br />

, T<br />

2<br />

y T 3<br />

dictan si <strong>el</strong> lugar<br />

está totalm<strong>en</strong>te confinado al cuarto cuadrante, o si se exti<strong>en</strong><strong>de</strong> <strong>en</strong> <strong>el</strong> primer<br />

y/o tercer cuadrante. En <strong>el</strong> caso límite, cuando T<br />

1<br />

= T2<br />

o T 3<br />

, la función <strong>de</strong><br />

transfer<strong>en</strong>cia se reduce a la forma:<br />

C(<br />

s)<br />

1<br />

G( s)<br />

= =<br />

(74)<br />

KR(<br />

s)<br />

( 1+<br />

Ts)<br />

La cual ocupa solam<strong>en</strong>te <strong>el</strong> cuarto cuadrante.<br />

Cada factor <strong>en</strong> <strong>el</strong> numerador <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia actúan<br />

comp<strong>en</strong>sando (o contrarrestando) un factor similar <strong>en</strong> <strong>el</strong> d<strong>en</strong>ominador, y<br />

por lo tanto actúa para restringir <strong>el</strong> lugar <strong>de</strong> la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia polar<br />

a un cuadrante m<strong>en</strong>or que es indicado por los factores <strong>en</strong> <strong>el</strong> d<strong>en</strong>ominador.<br />

Ejemplo 3: Un sistema ti<strong>en</strong>e la función <strong>de</strong> transfer<strong>en</strong>cia:<br />

K( 1+<br />

T1s<br />

)<br />

( 1+<br />

T s)( 1+<br />

T s)<br />

C(<br />

s)<br />

= (75)<br />

R(<br />

s)<br />

2<br />

3<br />

- 35 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Don<strong>de</strong>:<br />

K = 20<br />

T = 0.1 seg<br />

T<br />

T<br />

1<br />

2<br />

3<br />

= 0.2 seg<br />

= 0.04 seg<br />

Establecer la característica <strong>de</strong> <strong>respuesta</strong> <strong>de</strong>l sistema <strong>en</strong> forma polar.<br />

Para <strong>el</strong> análisis <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> estado estacionario, s es<br />

reemplazado por j ω , por lo tanto:<br />

C(<br />

jω)<br />

R(<br />

jω)<br />

20( 1+<br />

jωT1<br />

)<br />

( 1+<br />

jωT<br />

)( 1+<br />

jωT<br />

)<br />

= (76)<br />

2<br />

3<br />

La cual es una r<strong>el</strong>ación <strong>de</strong> números complejos (vectores) <strong>de</strong> la forma<br />

don<strong>de</strong>:<br />

P Q R ,<br />

P = 1+<br />

Q = 1+<br />

R = 1+<br />

jωT<br />

1<br />

jωT<br />

jωT<br />

2<br />

3<br />

En notación polar, la magnitud y ángulo <strong>de</strong> cada uno es:<br />

P =<br />

Q =<br />

R =<br />

1+<br />

ω T<br />

2 2<br />

1<br />

2<br />

1+<br />

ω T<br />

2<br />

1+<br />

ω T<br />

2<br />

2<br />

2<br />

3<br />

,<br />

,<br />

,<br />

tg<br />

tg<br />

tg<br />

−1<br />

−1<br />

−1<br />

ωT<br />

1<br />

ωT<br />

ωT<br />

2<br />

3<br />

Y la r<strong>el</strong>ación <strong>de</strong> magnitud es:<br />

2 2<br />

20( 1+<br />

ω T1<br />

)<br />

2 2<br />

2 2<br />

( 1+<br />

ω T2<br />

)( 1+<br />

ω T3<br />

)<br />

C(<br />

jω)<br />

M ( ω)<br />

= =<br />

(77)<br />

R(<br />

jω)<br />

Y la fase:<br />

φ<br />

−1<br />

−1<br />

−1<br />

( ω)<br />

tg ωT1<br />

− tg ωT2<br />

− tg ωT3<br />

= (78)<br />

Reemplazando T 1<br />

, T 2<br />

y T 3<br />

por sus valores:<br />

- 36 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

2<br />

20( 1+<br />

0.01ω<br />

)<br />

2<br />

2<br />

( 1+<br />

0.04ω<br />

)( 1+<br />

0.0016ω<br />

)<br />

C(<br />

jω)<br />

M ( ω)<br />

= =<br />

(79)<br />

R(<br />

jω)<br />

−1 −1<br />

−1<br />

φ(<br />

ω)<br />

= tg 0.1ω<br />

− tg 0.2ω<br />

− tg 0. 04ω<br />

(80)<br />

Para varios valores <strong>de</strong> ω [ rad seg]<br />

, se pue<strong>de</strong> obt<strong>en</strong>er la sigui<strong>en</strong>te tabla:<br />

ω 0 1 2 5 10 15 20 30 40 50 100 ∞<br />

M 20 19.7 18,4 15,5 11,6 9,75 8,6 6,7 5,42 4,54 2,4 0<br />

φ 0 -8,2 -15,2 -30 -40 -46 -51 -59,5 -65 -69 -76 -90<br />

Tabla 2.<br />

La figura 30 muestra la repres<strong>en</strong>tación polar <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia<br />

propuesta.<br />

Figura 30: Repres<strong>en</strong>tación polar <strong>de</strong><br />

20( 1+<br />

jωT1<br />

)<br />

( 1+<br />

jωT<br />

)( 1+<br />

jωT<br />

)<br />

2<br />

3<br />

5.7 - El Efecto De Los Ceros Y Polos En El Orig<strong>en</strong> En La Función De<br />

Transfer<strong>en</strong>cia.<br />

Una función <strong>de</strong> transfer<strong>en</strong>cia pue<strong>de</strong> cont<strong>en</strong>er factores <strong>de</strong> la forma:<br />

1 1 ;<br />

2<br />

Ts Ts<br />

etc. O bi<strong>en</strong> Ts ;<br />

2<br />

Ts ; etc.<br />

Tales factores aparecerán frecu<strong>en</strong>tem<strong>en</strong>te <strong>en</strong> los mo<strong>de</strong>los matemáticos <strong>de</strong><br />

los <strong>sistemas</strong> <strong>de</strong> control a estudiar, y se los conoce con <strong>el</strong> nombre <strong>de</strong> ceros y<br />

polos <strong>en</strong> <strong>el</strong> orig<strong>en</strong> <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia.<br />

Si se consi<strong>de</strong>ra <strong>el</strong> polo simple <strong>en</strong> <strong>el</strong> orig<strong>en</strong>:<br />

- 37 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

1<br />

G( s)<br />

=<br />

Ts<br />

Para <strong>el</strong> análisis <strong>de</strong> <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia:<br />

1<br />

G( jω)<br />

= (81)<br />

jωT<br />

Don<strong>de</strong>, <strong>en</strong> forma polar:<br />

1 1<br />

G(<br />

jω)<br />

= =<br />

(82)<br />

2 2<br />

1+<br />

ω T ωT<br />

−1<br />

ωT<br />

−1<br />

ωT<br />

−1<br />

φ = 0 − tg = −tg<br />

= −tg<br />

∞ = −90º<br />

(83)<br />

0 0<br />

Lo que indica que produce un atraso <strong>de</strong> 90° para todas las frecu<strong>en</strong>cias, por<br />

lo tanto la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia está sobre <strong>el</strong> eje <strong>de</strong> los 270° = - 90 o<br />

A<strong>de</strong>más para,<br />

G ( jω)<br />

= ∞ y para ω = ∞ , G ( jω)<br />

= 0<br />

El lugar <strong>de</strong>l polo simple es mostrado <strong>en</strong> la figura 31a.<br />

2 2<br />

Se pue<strong>de</strong> <strong>de</strong>mostrar <strong>en</strong> forma similar que <strong>el</strong> polo doble <strong>en</strong> <strong>el</strong> orig<strong>en</strong> ( 1 s T )<br />

ti<strong>en</strong>e lugar <strong>de</strong> la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia polar que comi<strong>en</strong>za <strong>en</strong> ∞ para<br />

ω = 0 sobre <strong>el</strong> eje <strong>de</strong> 180° = - 180° y prosigue a lo largo <strong>de</strong> este eje hacia <strong>el</strong><br />

orig<strong>en</strong> cuando ω ti<strong>en</strong><strong>de</strong> a aum<strong>en</strong>tar, alcanzando <strong>el</strong> mismo para ω = ∞ , <strong>el</strong><br />

lugar <strong>de</strong>l polo doble <strong>en</strong> <strong>el</strong> orig<strong>en</strong> está mostrado <strong>en</strong> la gráfica 31 b.<br />

Figura 31: a) Gráfica polar <strong>de</strong><br />

1 sT , b) Gráfica polar <strong>de</strong><br />

2 2<br />

1 s T<br />

- 38 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Cuando se ti<strong>en</strong>e un cero simple <strong>en</strong> <strong>el</strong> orig<strong>en</strong> Ts , la repres<strong>en</strong>tación<br />

frecu<strong>en</strong>cial es un fasor complejo con parte real e igual a cero, y por lo tanto<br />

<strong>el</strong> lugar está a lo largo <strong>de</strong>l eje imaginario.<br />

G( s)<br />

= Ts ∴ G(<br />

jω ) = jωT<br />

(84)<br />

2 2<br />

G( jω ) = ω T = ωT<br />

(85)<br />

−1<br />

−1<br />

φ ( ω)<br />

= tg ωT = tg ∞ = + 90º<br />

(86)<br />

0<br />

La fase será <strong>en</strong>tonces + 90° para cualquier valor <strong>de</strong> ω . A<strong>de</strong>más para ω = 0 ,<br />

G ( j0)<br />

= 0 y para ω = ∞ , G ( j∞) = ∞ , tal gráfica se muestra <strong>en</strong> la figura 32.<br />

5.8 - Trazados Polares Inversos.<br />

Figura 32: Gráfica polar <strong>de</strong> Ts<br />

Los trazados polares vistos anteriorm<strong>en</strong>te se los conoce con <strong>el</strong> nombre <strong>de</strong><br />

traza dos polares directos. Los trazados polares directos ti<strong>en</strong><strong>en</strong> ciertas<br />

<strong>de</strong>sv<strong>en</strong>tajas <strong>en</strong> <strong>sistemas</strong> que ti<strong>en</strong><strong>en</strong> <strong>el</strong>em<strong>en</strong>tos diversos <strong>en</strong> <strong>el</strong> lazo <strong>de</strong><br />

realim<strong>en</strong>tación, <strong>de</strong>jando <strong>de</strong> ser ésta unitaria. En estos casos se ha<br />

<strong>en</strong>contrado que es mucho más s<strong>en</strong>cillo para <strong>el</strong> análisis gráfico <strong>el</strong> empleo <strong>de</strong><br />

trazados polares inversos que pres<strong>en</strong>tan ciertas v<strong>en</strong>tajas. Este trazado se<br />

obti<strong>en</strong>e repres<strong>en</strong>tando <strong>el</strong> fasor <strong>de</strong>:<br />

G<br />

( jω)<br />

= (87)<br />

G( jω)<br />

−1 1<br />

En función <strong>de</strong> la frecu<strong>en</strong>cia.<br />

- 39 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

6 - Gráficos Rectangulares.<br />

6.1 - Introducción.<br />

Los resultados obt<strong>en</strong>idos al realizar <strong>el</strong> análisis frecu<strong>en</strong>cial, pued<strong>en</strong> ser<br />

graficados <strong>en</strong> coord<strong>en</strong>adas rectangulares. La r<strong>el</strong>ación <strong>de</strong> magnitud M ( ω ) y<br />

<strong>el</strong> ángulo <strong>de</strong> fase φ( ω ) se grafican contra la frecu<strong>en</strong>cia. Es normal y<br />

conv<strong>en</strong>i<strong>en</strong>te hacer los gráficos contra<br />

10<br />

log ω<br />

En <strong>el</strong> caso <strong>de</strong> los ángulos <strong>de</strong> fase φ( ω ) se usa una escala lineal, <strong>de</strong> modo<br />

que la gráfica <strong>de</strong> φ( ω ) contra ω se hace <strong>en</strong> pap<strong>el</strong> semilogarítmico. En <strong>el</strong><br />

caso <strong>de</strong> la r<strong>el</strong>ación ce magnitud M ( ω ) pue<strong>de</strong> graficarse <strong>de</strong> dos formas: se<br />

pue<strong>de</strong> tomar log M ( ω ), <strong>en</strong> cuyo caso la gráfica <strong>de</strong> M ( ω ) contra ω se hace<br />

<strong>en</strong> pap<strong>el</strong> log-log; o pue<strong>de</strong> expresarse M ( ω ) <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es, una unidad<br />

logarítmica, y graficarse sobre una escala lineal (o sea, sobre pap<strong>el</strong><br />

semilogarítmico).<br />

El caso más común <strong>de</strong> gráfica rectangular es la que emplea M ( ω ) <strong>en</strong><br />

<strong>de</strong>cib<strong>el</strong>es, porque permite graficar <strong>en</strong> la misma hoja tanto M ( ω ) y φ( ω )<br />

contra log ω ; esta gráfica es conocida con <strong>el</strong> nombre <strong>de</strong> diagrama <strong>de</strong> Bo<strong>de</strong><br />

10<br />

y se estudiara <strong>en</strong> <strong>el</strong> próximo apartado.<br />

Se <strong>de</strong>finirá ahora, ya que será <strong>de</strong> uso perman<strong>en</strong>te <strong>en</strong> todo <strong>el</strong> estudio que<br />

se realizará, la unidad logarítmica conocida como <strong>de</strong>cib<strong>el</strong>.<br />

El uso <strong>de</strong>l término <strong>de</strong>cib<strong>el</strong>, se <strong>de</strong>be a que la r<strong>el</strong>ación <strong>de</strong> magnitud fue<br />

empleada <strong>en</strong> principio como ganancia <strong>en</strong> los amplificadores <strong>el</strong>ectrónicos, los<br />

cuales tuvieron sus primeras aplicaciones <strong>en</strong> la amplificación <strong>de</strong> señales<br />

<strong>el</strong>éctricas para su conversión <strong>en</strong> sonidos audibles„ El oído humano respon<strong>de</strong><br />

a los estímulos <strong>en</strong> una escala logarítmica.<br />

Si M es una magnitud y m es dicha magnitud <strong>en</strong> <strong>de</strong>cí<strong>de</strong>les (dB), <strong>en</strong>tonces:<br />

m<br />

= 20log M<br />

10<br />

[ m]<br />

= dB<br />

10<br />

Ejemplos:<br />

M=1⇒<br />

m=20log 1= 20(0)= 0 dB<br />

M=10 ⇒ m=20log 10=20 (l)= 20 dB<br />

M=100 ⇒ m=20log 100=20(2)= 40 dB<br />

M=0.1 ⇒ m=20log 0.1=20(- l)= -20 dB<br />

M=0.2 ⇒ m=20log 0.2=20(-0,099) = -13,98 dB<br />

La conversión <strong>de</strong> magnitud a dB se facilita por <strong>el</strong> uso <strong>de</strong> gráficos y tablas.<br />

En la figura 33 se muestra un gráfico que permite convertir números <strong>en</strong><br />

<strong>de</strong>cí<strong>de</strong>les.<br />

- 40 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Se pue<strong>de</strong> hacer notar que dado un número, <strong>el</strong> recíproco ti<strong>en</strong>e <strong>el</strong> mismo<br />

número <strong>de</strong> <strong>de</strong>cib<strong>el</strong>es pero con <strong>el</strong> signo opuesto. Así <strong>el</strong> valor <strong>en</strong> <strong>de</strong>cí<strong>de</strong>les <strong>de</strong><br />

2 es +6 dB y <strong>el</strong> <strong>de</strong> 0,5 será <strong>de</strong> —6 dB.<br />

Cuando un número se duplica, su valor <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es aum<strong>en</strong>ta <strong>en</strong> 6 dB. El<br />

número 2 es <strong>el</strong> doble <strong>de</strong> 1, y su valor logarítmico es 6 <strong>de</strong>cib<strong>el</strong>es mayor. El<br />

número 200 es <strong>el</strong> doble <strong>de</strong> 100, y su expresión logarítmica 6 dB mayor.<br />

Cuando un número hace 10 veces mayor su expresión logarítmica aum<strong>en</strong>ta<br />

<strong>en</strong> 20 dB. El número 100 es diez veces mayor que 10 y su expresión<br />

logarítmica es 20 dB mayor que la <strong>de</strong> 10. El número 200 es ci<strong>en</strong> veces<br />

mayor que 2, su expresión loga rítmica es 40 dB mayor.<br />

20<br />

10<br />

Decib<strong>el</strong>es, dB<br />

0<br />

-10<br />

-20<br />

-30<br />

-40<br />

10 -2 10 -1 10 0 10 1<br />

Números<br />

Figura 33: Recta <strong>de</strong> conversión números-<strong>de</strong>cib<strong>el</strong>es.<br />

- 41 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

6.2 - Diagrama <strong>de</strong> Bo<strong>de</strong>.<br />

El análisis y la predicción <strong>de</strong>l comportami<strong>en</strong>to dinámico <strong>de</strong> <strong>sistemas</strong><br />

estudiando su <strong>respuesta</strong> a la frecu<strong>en</strong>cia pue<strong>de</strong> llevarse a cabo usando<br />

gráficas <strong>en</strong> coord<strong>en</strong>adas rectangulares, aparte <strong>de</strong> las gráficas polares ya<br />

estudiadas. Otra vez, las cantida<strong>de</strong>s a estudiar son las r<strong>el</strong>aciones <strong>de</strong><br />

magnitud y fase <strong>de</strong> la <strong>respuesta</strong> a la frecu<strong>en</strong>cia <strong>de</strong> lazo cerrado (más<br />

a<strong>de</strong>lante, cuando se estudie <strong>el</strong> capitulo <strong>de</strong> estabilidad <strong>de</strong> compon<strong>en</strong>tes y<br />

<strong>sistemas</strong>, se verá la utilidad <strong>de</strong> la repres<strong>en</strong>tación frecu<strong>en</strong>cial <strong>de</strong> la función<br />

<strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> lazo abierto),<br />

Un diagrama logarítmico o <strong>de</strong> Bo<strong>de</strong> es una repres<strong>en</strong>tación <strong>de</strong> <strong>respuesta</strong> a la<br />

frecu<strong>en</strong>cia que consta <strong>de</strong> dos gráficos: uno que da la amplitud <strong>en</strong> función <strong>de</strong><br />

la frecu<strong>en</strong>cia (diagrama <strong>de</strong> at<strong>en</strong>uación <strong>de</strong> Bo<strong>de</strong>) y otro <strong>el</strong> ángulo <strong>de</strong> fase <strong>en</strong><br />

función <strong>de</strong> la frecu<strong>en</strong>cia. El trazado <strong>de</strong>l diagrama <strong>de</strong> at<strong>en</strong>uación se realiza <strong>en</strong><br />

pap<strong>el</strong> semilogarítmico, colocando <strong>en</strong> ord<strong>en</strong>adas <strong>el</strong> módulo <strong>de</strong> la amplitud <strong>de</strong><br />

la función <strong>de</strong> transfer<strong>en</strong>cia a repres<strong>en</strong>tar <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es y <strong>en</strong> abscisas <strong>el</strong><br />

logaritmo <strong>de</strong> la frecu<strong>en</strong>cia. El trazado <strong>de</strong>l diagrama <strong>de</strong> fase se realiza<br />

también <strong>en</strong> pap<strong>el</strong> semilogarítmico colocando <strong>en</strong> la ord<strong>en</strong>ada la fase <strong>de</strong> la<br />

función <strong>de</strong> transfer<strong>en</strong>cia a repres<strong>en</strong>tar <strong>en</strong> grados y <strong>en</strong> abscisas <strong>el</strong> logaritmo<br />

<strong>de</strong> la frecu<strong>en</strong>cia.<br />

La v<strong>en</strong>taja principal <strong>de</strong> usar <strong>el</strong> diagrama <strong>de</strong> Bo<strong>de</strong> <strong>en</strong> lugar <strong>de</strong> la gráfica <strong>en</strong><br />

coord<strong>en</strong>adas polares radica <strong>en</strong> la mayor facilidad <strong>de</strong> la construcción <strong>de</strong> las<br />

gráficas. Los <strong>sistemas</strong> <strong>de</strong> control están compuestos, a m<strong>en</strong>udo, por<br />

compon<strong>en</strong>tes que son aproximados por bloques que ti<strong>en</strong><strong>en</strong> funciones <strong>de</strong><br />

transfer<strong>en</strong>cia formadas por constantes, K , integraciones, 1 s ,<br />

difer<strong>en</strong>ciaciones, s , constantes <strong>de</strong> tiempo, ( ) 1<br />

2<br />

2 ±<br />

( s 2δω<br />

+ ω ) 1<br />

±<br />

1± Ts , y factores cuadráticos,<br />

+ ns<br />

n<br />

. El diagrama <strong>de</strong> Bo<strong>de</strong> permite la superposición sistemática<br />

<strong>de</strong> los efectos <strong>de</strong> los difer<strong>en</strong>tes <strong>el</strong>em<strong>en</strong>tos tomados individualm<strong>en</strong>te.<br />

A<strong>de</strong>más, los diagramas <strong>de</strong> at<strong>en</strong>uación <strong>de</strong> Bo<strong>de</strong> y las gráficas <strong>de</strong>l ángulo <strong>de</strong><br />

fase pued<strong>en</strong> obt<strong>en</strong>erse a partir <strong>de</strong>l análisis experim<strong>en</strong>tal <strong>de</strong> la frecu<strong>en</strong>cia <strong>de</strong><br />

los compon<strong>en</strong>tes para los que no se conoc<strong>en</strong> las funciones <strong>de</strong> transfer<strong>en</strong>cia<br />

reales. Las gráficas pued<strong>en</strong> usarse para establecer las r<strong>el</strong>aciones <strong>de</strong><br />

transfer<strong>en</strong>cia reales <strong>de</strong> s .<br />

Es también útil la repres<strong>en</strong>tación logarítmica <strong>de</strong> Bo<strong>de</strong>, porque pres<strong>en</strong>ta las<br />

características <strong>de</strong> alta y baja frecu<strong>en</strong>cia <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia <strong>en</strong> un<br />

solo diagrama. Es muy v<strong>en</strong>tajoso <strong>el</strong> po<strong>de</strong>r expandir <strong>el</strong> rango <strong>de</strong> bajas<br />

frecu<strong>en</strong>cias utilizando una escala logarítmica <strong>de</strong> frecu<strong>en</strong>cias, ya que a<br />

frecu<strong>en</strong>cias bajas son muy importante las características <strong>de</strong> los <strong>sistemas</strong><br />

utilizados (<strong>de</strong>bido al uso <strong>de</strong> escala <strong>de</strong> frecu<strong>en</strong>cia logarítmica es imposible<br />

trazar las curvas <strong>de</strong> hasta frecu<strong>en</strong>cia cero; sin embargo esto no crea un<br />

problema importante).<br />

Se verá a continuación los conceptos teóricos <strong>de</strong> diagrama <strong>de</strong> Bo<strong>de</strong>,<br />

incluy<strong>en</strong>do su construcción, a<strong>de</strong>más <strong>de</strong>l análisis experim<strong>en</strong>tal <strong>de</strong> la<br />

<strong>respuesta</strong> frecu<strong>en</strong>cial. En <strong>el</strong> próximo capítulo, cuando se estudi<strong>en</strong> los<br />

conceptos <strong>de</strong> estabilidad <strong>de</strong> <strong>sistemas</strong> y compon<strong>en</strong>tes, se estudiará <strong>el</strong> uso <strong>en</strong><br />

<strong>el</strong> análisis <strong>de</strong> estabilidad a los diagramas <strong>de</strong> Bo<strong>de</strong>.<br />

- 42 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

6.3- Conceptos Básicos Sobre los Diagramas <strong>de</strong> At<strong>en</strong>uación <strong>de</strong> Bo<strong>de</strong>.<br />

La v<strong>en</strong>taja <strong>de</strong> los trazados logarítmicos es que las operaciones <strong>de</strong><br />

multiplicación y división se conviert<strong>en</strong> <strong>en</strong> otras <strong>de</strong> suma o resta y que <strong>el</strong><br />

trabajo para obt<strong>en</strong>er la repres<strong>en</strong>tación frecu<strong>en</strong>cial es más bi<strong>en</strong> gráfica que<br />

analítica. Los factores básicos que usualm<strong>en</strong>te aparec<strong>en</strong> <strong>en</strong> las funciones <strong>de</strong><br />

transfer<strong>en</strong>cia pued<strong>en</strong> ser aproximados por asuntotas, para una primer<br />

estimación y luego se traza la curva exacta para una mayor exactitud.<br />

Como ya se ha visto anteriorm<strong>en</strong>te, la función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> un<br />

sistema resultante <strong>de</strong> reemplazar s por j ω es, <strong>en</strong> g<strong>en</strong>eral, un número<br />

complejo.<br />

El logaritmo neperiano, base e, <strong>de</strong> un número complejo es, a su vez, otro<br />

número complejo:<br />

ln<br />

ln<br />

jφ<br />

( ω)<br />

jφ<br />

( ω)<br />

[ M ( jω)<br />

e ] ln M ( jω)<br />

+ ln e<br />

= (88)<br />

jφ<br />

( ω)<br />

[ M ( jω)<br />

e ] ln M ( jω)<br />

+ jφ(<br />

ω)<br />

= (89)<br />

El logaritmo <strong>de</strong>l número complejo<br />

( ω)<br />

ω)<br />

e<br />

jφ<br />

M ( j ti<strong>en</strong>e una parte real igual al<br />

logaritmo neperiano <strong>de</strong> la magnitud ln M ( jω)<br />

y una parte imaginaria igual<br />

al ángulo j φ(ω)<br />

expresado <strong>en</strong> radianes. Ambos compon<strong>en</strong>tes son función <strong>de</strong><br />

la frecu<strong>en</strong>cia ω .<br />

De igual forma, <strong>el</strong> logaritmo <strong>de</strong>cimal, base 10, <strong>de</strong> un número complejo es<br />

también otro número complejo:<br />

log<br />

log<br />

jφ<br />

( ω)<br />

jφ<br />

( ω)<br />

[ M ( jω)<br />

e ] log M ( jω)<br />

+ log e<br />

= (90)<br />

jφ<br />

( ω)<br />

[ M ( jω)<br />

e ] log M ( jω)<br />

+ j0.434φ<br />

( ω)<br />

= (91)<br />

Como podía preverse, la parte real es igual al logaritmo <strong>de</strong>cimal <strong>de</strong> la<br />

magnitud, log M ( jω)<br />

, y la parte imaginaria es 0.434φ<br />

( ω)<br />

, empleándose <strong>en</strong> <strong>el</strong><br />

resto <strong>de</strong> este estudio solam<strong>en</strong>te la parte proporcional igual a φ (ω)<br />

omitiéndose, por tanto, <strong>el</strong> factor 0,434.<br />

También se pue<strong>de</strong> llegar a igual conclusión sigui<strong>en</strong>do la teoría <strong>de</strong> la variable<br />

compleja que establece que <strong>el</strong> logaritmo <strong>de</strong> un número complejo es otro<br />

número complejo.<br />

Si es: M ( jω)<br />

= α(<br />

ω)<br />

+ jβ<br />

( ω)<br />

(92)<br />

Será <strong>en</strong>tonces:<br />

[ φ(<br />

ω)<br />

π n]<br />

log M ( jω)<br />

= log R(<br />

ω)<br />

+ j + 2<br />

(93)<br />

Para n=1, 2,…<br />

Y don<strong>de</strong>:<br />

- 43 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

2<br />

2<br />

R ( ω ) = α ( ω)<br />

+ β ( ω)<br />

y<br />

φ(<br />

ω)<br />

= tg<br />

−1<br />

β ( ω)<br />

α(<br />

ω)<br />

El ángulo 2 π n esta pres<strong>en</strong>te <strong>de</strong>bido a la exist<strong>en</strong>cia <strong>de</strong> posibilida<strong>de</strong>s <strong>de</strong> t<strong>en</strong>er<br />

<strong>el</strong> vector <strong>en</strong> <strong>el</strong> mismo lugar sobre <strong>el</strong> plano complejo.<br />

Nuevam<strong>en</strong>te, <strong>el</strong> logaritmo <strong>de</strong> M ( jω)<br />

está compuesto por dos funciones <strong>de</strong><br />

ω , la parte real: log R ( ω)<br />

y la parte imaginaria: φ (ω)<br />

<strong>en</strong> contra <strong>de</strong>l logaritmo<br />

<strong>de</strong> ω .<br />

Si se admite que una forma g<strong>en</strong>eral <strong>de</strong> expresar la función <strong>de</strong> transfer<strong>en</strong>cia<br />

<strong>de</strong> un sistema controlado es la sigui<strong>en</strong>te r<strong>el</strong>ación <strong>de</strong> valores reales y<br />

complejos:<br />

N1(<br />

jω)<br />

N<br />

2<br />

( jω)<br />

M ( jω)<br />

= K<br />

(94)<br />

D ( jω)<br />

D ( jω)<br />

1<br />

2<br />

Entonces ahora sigui<strong>en</strong>do los axiomas <strong>de</strong> logaritmos:<br />

log M ( jω)<br />

= log K + log N1(<br />

jω)<br />

+ log N<br />

2<br />

( jω)<br />

− log D1<br />

( jω)<br />

− log D2<br />

( jω)<br />

(95)<br />

log M ( jω)<br />

= log K + log N<br />

−<br />

log M ( jω)<br />

= log K + log N<br />

+<br />

( jω)<br />

+ φ ( jω)<br />

+ log N<br />

( jω)<br />

+ φ<br />

1<br />

2<br />

[ log D ( jω)<br />

+ φ ( jω)<br />

+ log D ( jω)<br />

+ φ ( jω)<br />

]<br />

j<br />

1<br />

1<br />

D<br />

1<br />

[ φ ( jω)<br />

+ φ ( jω)<br />

− φ ( jω)<br />

−φ<br />

( jω)<br />

]<br />

N<br />

1<br />

1<br />

N<br />

( jω)<br />

+ log<br />

N<br />

2<br />

1<br />

2<br />

2<br />

2<br />

2<br />

D<br />

1<br />

2<br />

N<br />

( jω)<br />

−<br />

N ( jω)<br />

− log D ( jω)<br />

− log D<br />

D<br />

D<br />

2<br />

( jω)<br />

(96)<br />

(97)<br />

Como se sabe la suma <strong>de</strong> varios números complejos es otro número<br />

complejo don<strong>de</strong> la parte real es la suma <strong>de</strong> todas las partes reales <strong>de</strong> los<br />

complejos sumados y la parte imaginaria la suma <strong>de</strong> las partes imaginarias<br />

<strong>de</strong> los complejos sumados.<br />

Esta r<strong>el</strong>ación se pue<strong>de</strong> usar ahora para obt<strong>en</strong>er las dos gráficas completas<br />

que repres<strong>en</strong>tan <strong>el</strong> log M ( jω)<br />

. Estas Gráficas serán:<br />

log M ( jω)<br />

= log K + log N1(<br />

jω)<br />

+ log N<br />

2<br />

( jω)<br />

− log D1<br />

( jω)<br />

− log D2<br />

( jω)<br />

(98)<br />

Y<br />

φ ω)<br />

= φ ( jω)<br />

+ φ ( jω)<br />

−φ<br />

( jω)<br />

− φ ( j )<br />

(99)<br />

( N1 N<br />

ω<br />

2<br />

D1<br />

D 2<br />

La unidad comúnm<strong>en</strong>te empleada <strong>en</strong> <strong>el</strong> estudio <strong>de</strong> <strong>sistemas</strong> realim<strong>en</strong>tados<br />

para <strong>el</strong> logaritmo <strong>de</strong>l módulo es <strong>el</strong> <strong>de</strong>cib<strong>el</strong>. Por lo que la ecuación 98 se<br />

pue<strong>de</strong> expresar <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es multiplicando ambos miembros por 20:<br />

- 44 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

20log M ( jω)<br />

= 20log K + 20log N<br />

− 20log D ( jω)<br />

− 20log D<br />

1<br />

1<br />

( jω)<br />

+ 20log<br />

2<br />

( jω)<br />

N ( jω)<br />

2<br />

(100)<br />

Las gráficas <strong>de</strong> las ecuaciones 99 y 100, son los diagramas <strong>de</strong> Bo<strong>de</strong>.<br />

Como pue<strong>de</strong> observarse ahora se evid<strong>en</strong>cia la s<strong>en</strong>cillez <strong>de</strong> graficar un<br />

diagrama <strong>de</strong> at<strong>en</strong>uación <strong>de</strong> Bo<strong>de</strong>. Las operaciones <strong>de</strong> multiplicación y<br />

división se han convertido <strong>en</strong> adición y sustracción. El procedimi<strong>en</strong>to es<br />

calcular, o buscar <strong>en</strong> gráficas o tablas, los valores <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es para cada<br />

uno <strong>de</strong> los factores <strong>de</strong> M ( jω)<br />

, para valores específicos <strong>de</strong> ω y combinarlos<br />

algebraicam<strong>en</strong>te.<br />

Como ya se había visto, hay varias maneras <strong>de</strong> graficar <strong>el</strong> logaritmo <strong>de</strong><br />

M ( jω)<br />

. Se adoptará <strong>el</strong> método sigui<strong>en</strong>te <strong>de</strong> ahora <strong>en</strong> más: se graficará <strong>el</strong><br />

log M ( jω)<br />

expresado <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es contra ω <strong>en</strong> pap<strong>el</strong> semilogarítmico. La<br />

ord<strong>en</strong>ada será <strong>en</strong>tonces 20log<br />

M ( jω)<br />

graficada <strong>en</strong> divisiones uniformem<strong>en</strong>te<br />

espaciadas, mi<strong>en</strong>tras que la abscisa es ω , graficada <strong>en</strong> una escala<br />

proporcional al logaritmo. El número <strong>de</strong> ciclos logarítmicos necesarios <strong>en</strong> la<br />

abscisa se <strong>de</strong>terminará por <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cia sobre <strong>el</strong> cual se va a<br />

investigar <strong>el</strong> sistema. La parte imaginaria φ (ω)<br />

se graficará también <strong>en</strong><br />

pap<strong>el</strong> semilogarítmico, colocando <strong>en</strong> ord<strong>en</strong>adas φ (ω)<br />

<strong>en</strong> grados y <strong>en</strong> la<br />

abscisa ω <strong>en</strong> una escala proporcional al logaritmo.<br />

Es conv<strong>en</strong>i<strong>en</strong>te m<strong>en</strong>cionar <strong>en</strong> este apartado que hay dos unida<strong>de</strong>s para<br />

expresar bandas <strong>de</strong> frecu<strong>en</strong>cias o coci<strong>en</strong>tes <strong>de</strong> frecu<strong>en</strong>cias: la octava y la<br />

década.<br />

Una octava es una banda <strong>de</strong> frecu<strong>en</strong>cias <strong>de</strong> ω<br />

1<br />

a<br />

2<br />

ω , <strong>en</strong> la que ω ω 2. Así<br />

2 1<br />

=<br />

la banda <strong>de</strong> frecu<strong>en</strong>cias <strong>en</strong>tre 1 a 2 c/s es una octava <strong>de</strong> ancha, y la banda<br />

<strong>de</strong> 2 a 4 c/s ti<strong>en</strong>e también una octava <strong>de</strong> ancha, lo mismo ocurre <strong>en</strong> la<br />

banda <strong>de</strong> 17.4 a 34,8 c/s. Nótese que una octava no correspon<strong>de</strong> a un<br />

ancho <strong>de</strong> banda fijo sino que <strong>de</strong>p<strong>en</strong><strong>de</strong> <strong>de</strong>l marg<strong>en</strong> <strong>de</strong> frecu<strong>en</strong>cias<br />

consi<strong>de</strong>rado. El número <strong>de</strong> octavas <strong>en</strong> un marg<strong>en</strong> <strong>de</strong> frecu<strong>en</strong>cias, <strong>de</strong>s<strong>de</strong> ω<br />

1<br />

a ω<br />

2<br />

vi<strong>en</strong>e dado por:<br />

Figura 34.<br />

( logω<br />

− ω )<br />

n (101)<br />

= número<strong>de</strong>octavas = logω2 − logω1<br />

= n<br />

j<br />

log<br />

n = número <strong>de</strong>octavas = logω 2<br />

− logω1<br />

= n log 2<br />

(102)<br />

1<br />

- 45 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Por lo que:<br />

n =<br />

⎛ ω2 ⎞<br />

log<br />

⎜ ⎟ ⎝ ω1<br />

⎠<br />

log 2<br />

octavas<br />

(103)<br />

⎛ ω2<br />

n<br />

⎟ ⎞<br />

= 3.32 log<br />

⎜<br />

⎝ ω1<br />

⎠<br />

octavas<br />

(104)<br />

Una década es una banda <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong>s<strong>de</strong> ω<br />

1<br />

a 10 ω<br />

1<br />

. Don<strong>de</strong><br />

nuevam<strong>en</strong>te ω<br />

1<br />

es cualquier frecu<strong>en</strong>cia. La banda <strong>de</strong> frecu<strong>en</strong>cias <strong>de</strong> 1 a 10<br />

c/s o <strong>de</strong> 2,5 a 25 c/s ti<strong>en</strong>e un ancho <strong>de</strong> banda <strong>de</strong> una década, <strong>el</strong> número <strong>de</strong><br />

décadas <strong>en</strong>tre dos frecu<strong>en</strong>cias ω<br />

1<br />

y ω<br />

2<br />

se calcula <strong>de</strong>l sigui<strong>en</strong>te modo:<br />

Figura 35.<br />

n<br />

( logω<br />

− ω )<br />

= número <strong>de</strong> décadas = logω2 − logω1<br />

= n<br />

j<br />

log<br />

n = número <strong>de</strong> décadas = logω<br />

− logω1<br />

=<br />

2<br />

n<br />

log10<br />

1<br />

∴<br />

⎛ ω ⎞<br />

n = log ⎜ 2<br />

⎟<br />

⎝ ω<br />

1 ⎠<br />

décadas<br />

(105)<br />

Se pue<strong>de</strong> ahora obt<strong>en</strong>er la r<strong>el</strong>ación <strong>en</strong>tre octavas y décadas <strong>en</strong>tre dos<br />

frecu<strong>en</strong>cias <strong>de</strong> las ecuaciones 104 y 105.<br />

⎛ ω<br />

2<br />

⎞<br />

n º <strong>de</strong> octavas = 3.32 log<br />

⎜<br />

⎟<br />

⎝ ω1<br />

⎠<br />

⎛ ω2<br />

⎞<br />

n º <strong>de</strong> décadas = log<br />

⎜<br />

⎟<br />

⎝ ω1<br />

⎠<br />

O sea:<br />

número <strong>de</strong> octavas = 3,32 número <strong>de</strong> décadas (106)<br />

- 46 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Estas r<strong>el</strong>aciones <strong>en</strong>tre frecu<strong>en</strong>cias serán usadas más a<strong>de</strong>lante, cuando se<br />

dibuj<strong>en</strong> los diagramas <strong>de</strong> Bo<strong>de</strong> aproximados por asíntotas.<br />

Ejemplo 4: Trazar <strong>el</strong>-diagrama <strong>de</strong> Bo<strong>de</strong> para la función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong><br />

lazo abierto <strong>de</strong>l sistema <strong>de</strong> la figura 36.<br />

Figura 36.<br />

La función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> lazo abierto es:<br />

B(<br />

s)<br />

100<br />

= = G(<br />

s)<br />

ya que se ve <strong>en</strong> la figura que H ( s)<br />

= 1<br />

E(<br />

s)<br />

s<br />

( s + 10)<br />

Para realizar <strong>el</strong> anáisis frecu<strong>en</strong>cial se reemplaza s por j ω :<br />

G ( jω)<br />

=<br />

100<br />

jω<br />

( jω<br />

+ 10)<br />

El sigui<strong>en</strong>te paso es normalizar esta función <strong>de</strong> transfer<strong>en</strong>cia, lo que se<br />

obti<strong>en</strong>e dividi<strong>en</strong>do por 10 <strong>el</strong> numerador y d<strong>en</strong>ominador.<br />

G(<br />

jω)<br />

=<br />

jω<br />

10<br />

( 1 + 0.1 jω)<br />

Entonces será:<br />

( 0.1 ) 2<br />

20 log G ( jω ) = 20 log10 − 20 logω<br />

− 20 log 1 + ω<br />

π<br />

φ( ω)<br />

= 0 − − arc tg 0. 1ω<br />

2<br />

Que son las expresiones <strong>de</strong> la parte real e imaginaria <strong>de</strong>l logaritmo <strong>de</strong> la<br />

ecuación 97.<br />

En la sigui<strong>en</strong>te tabla se obti<strong>en</strong><strong>en</strong> valores para <strong>el</strong> logaritmo <strong>de</strong>l módulo <strong>en</strong> dB<br />

y la fase <strong>de</strong> las ecuaciones 98 y 99 para varios valores <strong>de</strong> ω .<br />

- 47 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

ω<br />

20 log G ( jω ) dB<br />

φ (ω) grados<br />

1 20 - 0 – 0 = 20 - 90 -5.7 = -95.7<br />

5 20 - 13,98 - 1 = 5,02 -90 -26.6 = -116.6<br />

10 20 - 20 - 3 = - 3 -90 -45 = -135<br />

20 20 - 26 - 7 = - 13 -90 -63.4 = -157.4<br />

50 20 - 34 – 14 = - 28 -90 -78.7 = -168.7<br />

Tabla 3<br />

Los diagramas <strong>de</strong> at<strong>en</strong>uación y <strong>de</strong> ángulo <strong>de</strong> fase se dan <strong>el</strong> la figura 37. La<br />

ord<strong>en</strong>ada <strong>de</strong> la izquierda está <strong>en</strong> dB para la gráfica <strong>de</strong> at<strong>en</strong>uación y la<br />

ord<strong>en</strong>ada <strong>de</strong> la <strong>de</strong>recha está <strong>en</strong> grados. El eje <strong>de</strong> 0 dB coinci<strong>de</strong> con <strong>el</strong> <strong>de</strong><br />

fase <strong>de</strong> -180º.<br />

Figura 37.<br />

- 48 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

6.4- Trazado Asintótico De Los Diagramas Logarítmicos De La<br />

Magnitud Y Fase. Simplificaciones Al Establecer Los Diagramas De<br />

Bo<strong>de</strong>.<br />

Las funciones <strong>de</strong> transfer<strong>en</strong>cia frecu<strong>en</strong>ciales pued<strong>en</strong> ponerse <strong>en</strong> forma<br />

g<strong>en</strong>eralizada como <strong>el</strong> coci<strong>en</strong>te <strong>de</strong> dos polinomios:<br />

( 1 + jωT<br />

)( 1 + jωT<br />

)<br />

K<br />

1<br />

2<br />

...<br />

M ( jω)<br />

= = K M´(<br />

jω)<br />

(107)<br />

⎛<br />

n<br />

⎛ 2δ<br />

⎞ ⎛ 1 ⎞ ⎞<br />

2<br />

( jω) ( 1 + jωT<br />

) ⎜<br />

a<br />

1 + jω<br />

⎜ ⎟( jω)<br />

⎟<br />

⎜<br />

2<br />

ω<br />

⎟ +<br />

n ω<br />

⎝ ⎝ ⎠ ⎝ n ⎠ ⎠<br />

m<br />

En <strong>el</strong> que K es la ganancia. Como se sabe, <strong>el</strong> logaritmo <strong>de</strong> esta función <strong>de</strong><br />

transfer<strong>en</strong>cia t<strong>en</strong>drá una parte real proporcional al logaritmo <strong>de</strong>l módulo y<br />

una parte imaginaria proporcional al ángulo.<br />

20log M( jω) = 20log K + 20log 1+ jωT1 + m20log 1 + jωT2<br />

+ ...<br />

⎛2δ<br />

⎞ ⎛ 1 ⎞ 2<br />

(108)<br />

.. −n20log jω− 20log 1+ jωTa<br />

− 20log 1 + ⎜ ⎟ jω+ ⎜ 2 ⎟( jω)<br />

−...<br />

⎝ωn<br />

⎠ ⎝ωn<br />

⎠<br />

M( jω) = K + 1+ jωT + m1 + jωT<br />

+ ....<br />

1 2<br />

⎛2δ<br />

⎞ ⎛ 1 ⎞ 2<br />

−njω− 1+ jωTa<br />

1+ jωTa<br />

− 1 + ⎜ ⎟ jω+ ⎜ 2 ⎟( jω)<br />

−...<br />

⎝ωn<br />

⎠ ⎝ωn<br />

⎠<br />

(109)<br />

Las dos ecuaciones anteriores sugier<strong>en</strong> que <strong>el</strong> logaritmo <strong>de</strong> M ( jω)<br />

podría<br />

graficarse sobreponi<strong>en</strong>do las contribuciones <strong>de</strong> cada uno <strong>de</strong> los factores <strong>de</strong><br />

M ( jω ). El proceso <strong>de</strong> establecer los diagramas <strong>de</strong> at<strong>en</strong>uación y fase <strong>de</strong><br />

Bo<strong>de</strong>, pue<strong>de</strong> simplificarse y hacerse más rápido <strong>en</strong> virtud <strong>de</strong> que los<br />

factores <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia que se <strong>en</strong>cu<strong>en</strong>tran usualm<strong>en</strong>te, s<strong>en</strong><br />

<strong>de</strong> la forma <strong>de</strong> los mostrados <strong>en</strong> la ecuación 107. El logaritmo <strong>de</strong> cada uno<br />

<strong>de</strong> estos factores, cuando se grafican separadam<strong>en</strong>te contra <strong>el</strong> logaritmo <strong>de</strong><br />

la frecu<strong>en</strong>cia, pres<strong>en</strong>ta ciertas características: Con un conocimi<strong>en</strong>to <strong>de</strong>l<br />

comportami<strong>en</strong>to <strong>de</strong> cada factor, es posible construir un diagrama <strong>de</strong> Bo<strong>de</strong><br />

aproximado para una M ( jω ) dada trazando las curvas para cada factor y<br />

sumando las curvas individuales gráficam<strong>en</strong>te porque sumar logaritmo <strong>de</strong><br />

ganancias, correspon<strong>de</strong> a multiplicarlas <strong>en</strong>tre sí.<br />

Se pue<strong>de</strong> simplificar aún más <strong>el</strong> proceso <strong>de</strong> obt<strong>en</strong>er <strong>el</strong> diagrama logarítmico<br />

utilizando aproximaciones asintóticas a las curvas <strong>de</strong> caca factor. Es posible,<br />

como se verá, ajustar <strong>el</strong> diagrama asintótico para que <strong>de</strong> la gráfica real con<br />

<strong>el</strong> uso <strong>de</strong> factores <strong>de</strong> corrección simples. El primer paso para simplificar <strong>el</strong><br />

trazo <strong>de</strong> los diagramas <strong>de</strong> Bo<strong>de</strong>, es analizar <strong>el</strong> comportami<strong>en</strong>to <strong>de</strong> cada uno<br />

<strong>de</strong> los factores intervini<strong>en</strong>tes <strong>en</strong> la ecuación 107.<br />

- 49 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

6.4.1.- Diagrama De Bo<strong>de</strong> De Una Constante K<br />

La ganancia constante K pue<strong>de</strong> p<strong>en</strong>sarse como un número complejo con<br />

parte imaginaria nula. La parte real pue<strong>de</strong> ser un número positivo o<br />

negativo que repres<strong>en</strong>ta un vector <strong>de</strong> magnitud K y que se <strong>en</strong>cu<strong>en</strong>tra sobre<br />

<strong>el</strong> eje real. Debido a que K es in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> ω . El diagrama <strong>de</strong><br />

at<strong>en</strong>uación <strong>de</strong> Bo<strong>de</strong> es una recta horizontal. El ángulo <strong>de</strong> fase es 0° si K es<br />

positivo y -180 ° si K es negativo. Un número mayor que la unidad t<strong>en</strong>drá<br />

un valor <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es mayor que cero, mi<strong>en</strong>tras que un número m<strong>en</strong>or que<br />

la unidad t<strong>en</strong>drá un valor negativo <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es.<br />

El efecto <strong>de</strong> variar la ganancia K <strong>en</strong> la función <strong>de</strong> transfer<strong>en</strong>cia, es que sube<br />

o baja la curva <strong>de</strong> logaritmo <strong>de</strong> magnitud <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia <strong>en</strong><br />

<strong>el</strong> valor constante correspondi<strong>en</strong>te, pero no ti<strong>en</strong>e efecto <strong>en</strong> <strong>el</strong> ángulo <strong>de</strong><br />

fase.<br />

La figura 38, muestra los diagramas <strong>de</strong> at<strong>en</strong>uación y fase para un valor <strong>de</strong><br />

K mayor que la unidad.<br />

40<br />

35<br />

30<br />

25<br />

Amplitud, dB<br />

20<br />

K [dB]<br />

15<br />

10<br />

5<br />

0<br />

10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

- 50 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

Fase, Grados<br />

0<br />

-0.2<br />

-0.4<br />

-0.6<br />

-0.8<br />

-1<br />

10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 38: Diagrama <strong>de</strong> Bo<strong>de</strong> para una constante K<br />

6.4.2.- Diagrama De Bo<strong>de</strong> De Factores Derivativos E Integral ( jω ) ± n<br />

(ceros y polos <strong>en</strong> <strong>el</strong> orig<strong>en</strong>)<br />

Los diagramas <strong>de</strong> at<strong>en</strong>uación y ángulo <strong>de</strong> fase para un <strong>el</strong>em<strong>en</strong>to<br />

1<br />

integrador ( jω)<br />

− o un <strong>el</strong>em<strong>en</strong>to difer<strong>en</strong>ciador ( jω ) son rectos.<br />

Se consi<strong>de</strong>rará <strong>en</strong> primer lugar <strong>el</strong> diagrama <strong>de</strong> at<strong>en</strong>uación <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es<br />

contra <strong>el</strong> lazo ω .<br />

Para una integración:<br />

20log<br />

1<br />

=− 20logω<br />

(110)<br />

jω Y para una difer<strong>en</strong>ciación:<br />

20log<br />

jω = 20logω<br />

(111)<br />

Corno se ve ambas ecuaciones son repres<strong>en</strong>taciones <strong>de</strong> rectas.<br />

La p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> estas rectas se obti<strong>en</strong>e <strong>de</strong>rivando las ecuaciones anteriores<br />

respecto a logω , es <strong>de</strong>cir:<br />

∂ 20log<br />

∂ logω<br />

±<br />

( jω<br />

) 1<br />

=±<br />

20 dB<br />

(112)<br />

- 51 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Luego <strong>en</strong> coord<strong>en</strong>adas rectangulares, a un increm<strong>en</strong>to unitario <strong>de</strong> logω le<br />

correspon<strong>de</strong>rá un increm<strong>en</strong>to <strong>de</strong> ± 20 dB . Por otro lado un increm<strong>en</strong>to<br />

unitario <strong>de</strong> logω es equival<strong>en</strong>te a un cambio <strong>de</strong> 1 a 10, <strong>de</strong> 10 a 100, etc. <strong>en</strong><br />

la escala logarítmica. Luego la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> estas rectas es ± 20 dB / década .<br />

La figura 39 repres<strong>en</strong>ta las p<strong>en</strong>di<strong>en</strong>tes <strong>de</strong> un factor integral y <strong>de</strong> un factor<br />

<strong>de</strong>rivativo.<br />

40<br />

30<br />

| jω |<br />

20<br />

10<br />

Amplitud, dB<br />

0<br />

-10<br />

-20<br />

| 1/jω |<br />

-30<br />

-40<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 39: Diagrama <strong>de</strong> at<strong>en</strong>uación <strong>de</strong> factores integrativos y <strong>de</strong>rivativos.<br />

Si <strong>el</strong> ord<strong>en</strong> <strong>de</strong>l factor integrativo o <strong>de</strong>rivativo es mayor que uno, <strong>en</strong>tonces se<br />

t<strong>en</strong>drá:<br />

( jω<br />

)<br />

± n<br />

20log =± 20nlogω<br />

(113)<br />

Que es también la ecuación <strong>de</strong> una recta cuya p<strong>en</strong>di<strong>en</strong>te será:<br />

∂ 20log<br />

( jω<br />

)<br />

∂ logω<br />

± n<br />

=± n20<br />

dB<br />

(114)<br />

O sea que la p<strong>en</strong>di<strong>en</strong>te será <strong>de</strong> ± n 20 dB / década .<br />

- 52 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

En la figura 40 se muestran les diagramas <strong>de</strong> at<strong>en</strong>uación para factores<br />

integrativos y <strong>de</strong>rivativos ( jω ) ± n<br />

hasta n igual a cuatro.<br />

200<br />

150<br />

100<br />

| jω 4 |<br />

| jω 3 |<br />

| jω 2 |<br />

50<br />

| jω |<br />

Amplitud, dB<br />

0<br />

| 1/jω |<br />

-50<br />

| 1/(jω ) 2 |<br />

-100<br />

| 1/(jω) 3 |<br />

-150<br />

| 1/(jω ) 4 |<br />

-200<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 40: diagramas <strong>de</strong> at<strong>en</strong>uación para (<br />

) n<br />

jω ±<br />

Se reduce <strong>de</strong>l análisis anterior que se ha usado la década como r<strong>el</strong>ación <strong>de</strong><br />

frecu<strong>en</strong>cias, pero se hubiere llegado a idéntico resultado si se hubiera usado<br />

la octava como r<strong>el</strong>ación <strong>de</strong> frecu<strong>en</strong>cia.<br />

Si se supone que se quiere obt<strong>en</strong>er la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> una integración, ( ) 1<br />

jω − ,<br />

consi<strong>de</strong>rando que la frecu<strong>en</strong>cia se dobla, o sea ω2 = 2ω1, <strong>en</strong>tonces:<br />

−20log ω −( − 20log ω ) =− 20log 2ω + 20logω<br />

2 1 1 1<br />

=−20log 2 − 20logω<br />

+ 20logω<br />

=−20log 2<br />

1 1<br />

− 20log 2 = 6.02 dB<br />

(115)<br />

O sea que la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la recta es -6.02 dB/octava.<br />

Por lo que se pue<strong>de</strong> <strong>de</strong>ducir que para un factor integrativo, ( ) jω , se t<strong>en</strong>drá<br />

una repres<strong>en</strong>tación frecu<strong>en</strong>cial que será una recta con una p<strong>en</strong>di<strong>en</strong>te <strong>de</strong>:<br />

- 53 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

− 20 dB / década =− 6.02 dB / octava<br />

Para factores integrativos o <strong>de</strong>rivativos ( jω ±<br />

rectas <strong>de</strong> p<strong>en</strong>di<strong>en</strong>tes:<br />

) n<br />

las repres<strong>en</strong>taciones serán<br />

± 20 n dB / década =± 6.02 n dB / octava<br />

Se <strong>de</strong>ducirán ahora <strong>el</strong> valor <strong>de</strong> la fase para factores integrativos y <strong>de</strong>rivativos.<br />

Factor integrativo:<br />

Factor <strong>de</strong>rivativo:<br />

1<br />

−1ω<br />

−1<br />

⇒ φω ( ) =− tg =−tg<br />

∞=− 90º<br />

jω<br />

0<br />

(117)<br />

−1 0 −1<br />

jω<br />

⇒ φ( ω) = tg =− tg 0=+ 90º<br />

ω<br />

(118)<br />

Si <strong>el</strong> ord<strong>en</strong> <strong>de</strong> los factores integrativos o <strong>de</strong>rivativos es n, <strong>en</strong>tonces:<br />

Factor integrativo:<br />

( jω<br />

)<br />

1<br />

n<br />

⇒ φω ( ) =− n90º<br />

(119)<br />

Factor <strong>de</strong>rivativo: ( jω) φ( ω) n90º<br />

n<br />

⇒ =+ (120)<br />

Por lo que se ve los factores integrativos o <strong>de</strong>rivativos contribuy<strong>en</strong> con una<br />

fase constante igual a ± n90º<br />

como se pue<strong>de</strong> ver <strong>en</strong> la figura 41.<br />

Fase, Grados<br />

360<br />

270<br />

180<br />

90<br />

0<br />

-90<br />

-180<br />

-270<br />

( jω ) 3<br />

( jω ) 2<br />

jω<br />

1/( jω )<br />

1/( jω ) 2<br />

1/( jω ) 3<br />

-360<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 41: Ángulo <strong>de</strong> fase para factores (<br />

jω ) ± n .<br />

- 54 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

1+<br />

j T<br />

ω ±<br />

6.4.3.- Diagramas De Bo<strong>de</strong> Para Factores De Primer Ord<strong>en</strong> ( ) 1<br />

Otro factor <strong>en</strong>contrado comúnm<strong>en</strong>te <strong>en</strong> las expresiones <strong>de</strong> funciones <strong>de</strong><br />

transfer<strong>en</strong>cia, es la constante <strong>de</strong> tiempo conocida g<strong>en</strong>eralm<strong>en</strong>te como<br />

retardo simple ( 1+ jωT<br />

) −1<br />

o a<strong>de</strong>lanto simple( 1+ jωT<br />

) .<br />

Se consi<strong>de</strong>rará <strong>en</strong> primer lugar la repres<strong>en</strong>tación logarítmica <strong>de</strong>l retar<strong>de</strong><br />

simple. La gráfica <strong>de</strong> at<strong>en</strong>uación se obti<strong>en</strong>e a partir <strong>de</strong>:<br />

1 1<br />

20log = 20log =− 20log 1+<br />

( 1+ jωT) 1+<br />

( ωT<br />

)<br />

2<br />

( ωT<br />

) 2<br />

(121)<br />

Y <strong>el</strong> ángulo <strong>de</strong> fase está dado por:<br />

−1<br />

φ( ω) =− tg ωT<br />

(122)<br />

Ahora <strong>el</strong> diagrama <strong>de</strong> Bo<strong>de</strong> pue<strong>de</strong> trazarse fácilm<strong>en</strong>te usando las ecuaciones<br />

anteriores, dándole valores a ω .<br />

Es posible, sin embargo, obt<strong>en</strong>er una aproximación asintótica <strong>de</strong> las curvas<br />

<strong>de</strong> at<strong>en</strong>uación y <strong>de</strong> fase, observando <strong>el</strong> comportami<strong>en</strong>to <strong>de</strong> las ecuaciones<br />

121 y 122 a varias frecu<strong>en</strong>cias.<br />

Le ecuación 121 se pue<strong>de</strong> poner:<br />

1<br />

20log =− 20log 1+<br />

ω T<br />

( 1+<br />

jωT)<br />

2 2<br />

dB<br />

(123)<br />

Para valores <strong>de</strong> ω muy bajos, tales que<br />

logaritmo <strong>de</strong> la amplitud por:<br />

ω <br />

1 , se pue<strong>de</strong> aproximar <strong>el</strong><br />

T<br />

2 2<br />

− 20log 1+ ω T ≅− 20log1 = 0 dB<br />

(124)<br />

De la ecuación 124 se <strong>de</strong>duce que la curva <strong>de</strong> logaritmo <strong>de</strong> la amplitud<br />

es asintótica a una recta <strong>de</strong> 0 dB.<br />

Para valores <strong>de</strong> ω muy altos, tales que<br />

poner:<br />

ω <br />

1 , la ecuación 123 se pue<strong>de</strong><br />

T<br />

- 55 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

2 2 2 2<br />

− 20log 1+ ω T ≅− 20log ω T =− 20logωT dB<br />

(125)<br />

Así, <strong>de</strong> la ecuación 125, se <strong>de</strong>duce que la curva <strong>de</strong> logaritmo <strong>de</strong> la amplitud<br />

<strong>en</strong> altas frecu<strong>en</strong>cias es una recta con una p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> -20 dB/década (ó -6<br />

dB/octava).<br />

D<strong>el</strong> análisis efectuado se concluye que la repres<strong>en</strong>tación <strong>de</strong> la curva <strong>de</strong><br />

<strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong>l factor ( ) 1<br />

−<br />

1+ jωT<br />

pue<strong>de</strong> ser aproximada por dos<br />

líneas rectas asintóticas, una a 0 dB/década para <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cias<br />

0< ω < 1 T y otra con p<strong>en</strong>di<strong>en</strong>te -20 dB/década para <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cia<br />

1 T < ω


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

El error a una frecu<strong>en</strong>cia una octava por <strong>de</strong>bajo <strong>de</strong> la frecu<strong>en</strong>cia esquina, es<br />

<strong>de</strong>cir, <strong>en</strong> ω = 1 2T<br />

es:<br />

5<br />

− 20 log 1+<br />

1 4 + 20 log0.5 = −20log<br />

= −0.97<br />

dB<br />

(128)<br />

2<br />

El error a una frecu<strong>en</strong>cia una octava por <strong>en</strong>cima <strong>de</strong> la frecu<strong>en</strong>cia esquina, es<br />

<strong>de</strong>cir, <strong>en</strong> ω = 2 T es:<br />

2<br />

5<br />

− 20 log 2 + 1 + 20 log 2 = −20log<br />

= −0.97<br />

dB<br />

(129)<br />

2<br />

Entonces, <strong>el</strong> error una octava por <strong>de</strong>bajo o por <strong>en</strong>cima <strong>de</strong> la frecu<strong>en</strong>cia<br />

esquina es aproximadam<strong>en</strong>te igual a -1 dB. En forma similar, <strong>el</strong> error a una<br />

década por <strong>de</strong>bajo o por <strong>en</strong>cima <strong>de</strong> la frecu<strong>en</strong>cia <strong>de</strong> transición es<br />

aproximadam<strong>en</strong>te igual a -0.04 dB.<br />

El error producido, <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es, al usar expresión asintótica para la curva<br />

<strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> ( ) 1<br />

asimétrico con respecto a la frecu<strong>en</strong>cia <strong>de</strong> corte.<br />

−<br />

1+ jωT<br />

aparece <strong>en</strong> la figura 43. El error es<br />

0<br />

-0.5<br />

-1<br />

-1.5<br />

dB<br />

-2<br />

-2.5<br />

-3<br />

-3.5<br />

10 -2 10 -1 10 0 10 1 10 2<br />

ω<br />

Figura 43: Error <strong>de</strong>l logaritmo <strong>de</strong>l módulo <strong>en</strong> la expresión asintótica <strong>de</strong> la curva <strong>de</strong><br />

<strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> ( ) 1<br />

−<br />

1+ jωT<br />

, para T = 1<br />

- 57 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Para trazar la curva <strong>de</strong> fase se sigue un procedimi<strong>en</strong>to similar al empleado<br />

para la gráfica <strong>de</strong>l módulo.<br />

La fase para <strong>el</strong> término estudiado es:<br />

( ) −1<br />

−1<br />

φ( ω) = 1+ jωT =− tg ωT<br />

En esta expresión se pue<strong>de</strong> ver que a frecu<strong>en</strong>cias muy bajas ωT<br />

1 ( ω T<br />

ti<strong>en</strong><strong>de</strong> a cero) <strong>el</strong> ángulo <strong>de</strong> fase ti<strong>en</strong><strong>de</strong> a 0 o .<br />

En la frecu<strong>en</strong>cia esquina, ω = 1 , <strong>el</strong> ángulo <strong>de</strong> fase es -45°.<br />

T<br />

Cuando las frecu<strong>en</strong>cias son <strong>el</strong>evadas ωT<br />

1 ( ω T ti<strong>en</strong><strong>de</strong> a infinito) <strong>el</strong> ángulo<br />

<strong>de</strong> fase ti<strong>en</strong><strong>de</strong> a -90°.<br />

En la figura 44 se pue<strong>de</strong> observar las difer<strong>en</strong>cias <strong>en</strong>tre los trazos asintóticos<br />

1+ jωT<br />

.<br />

y exactos para la fase <strong>de</strong>l factor ( ) 1<br />

−<br />

20<br />

0<br />

Trazo Exacto<br />

Fase, Grados<br />

-20<br />

-40<br />

-60<br />

Trazo Asintótico<br />

-80<br />

-100<br />

10 -2 10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 44: Ángulo <strong>de</strong> fase asintótico y exacto para <strong>el</strong> factor ( ) −1<br />

1+ jωT<br />

.<br />

Como las asíntotas son fácilm<strong>en</strong>te dibujables y sufici<strong>en</strong>tem<strong>en</strong>te cercanas a<br />

la curva exacta, <strong>el</strong> uso <strong>de</strong> esas aproximaciones al trazar los diagramas <strong>de</strong><br />

Bo<strong>de</strong> es conv<strong>en</strong>i<strong>en</strong>te para establecer rápidam<strong>en</strong>te la naturaleza g<strong>en</strong>eral <strong>de</strong><br />

las características <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia con una tarea mínima <strong>de</strong><br />

cálculo, y se la pue<strong>de</strong> utilizar como fase pr<strong>el</strong>iminar <strong>en</strong> <strong>el</strong> trabajo <strong>de</strong> diseño.<br />

Si se requier<strong>en</strong> curvas exactas <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia, se pued<strong>en</strong><br />

realizar fácilm<strong>en</strong>te correcciones con refer<strong>en</strong>cia a las curvas dadas 125 y<br />

130. En la práctica se pue<strong>de</strong> trazar una curva <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia<br />

- 58 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

bastante precisa ubicando <strong>el</strong> punto <strong>de</strong> -3 dB <strong>en</strong> la frecu<strong>en</strong>cia esquina y los<br />

puntos <strong>de</strong> -0.04 dB una década por <strong>de</strong>bajo y <strong>en</strong>cima <strong>de</strong> esa frecu<strong>en</strong>cia y<br />

conectando luego esos puntos por una curva suave.<br />

Se <strong>de</strong>be hacer notar que variar la constante <strong>de</strong> tiempo T <strong>de</strong>splaza la<br />

frecu<strong>en</strong>cia <strong>de</strong> esquina a la izquierda o a la <strong>de</strong>recha, pero las formas <strong>de</strong><br />

logaritmo <strong>de</strong> la amplitud y ángulo <strong>de</strong> fase quedan iguales.<br />

La función <strong>de</strong> transfer<strong>en</strong>cia ( ) −1<br />

1+ jωT<br />

ti<strong>en</strong>e las características <strong>de</strong> un filtro<br />

pasa bajo. Para frecu<strong>en</strong>cias por <strong>en</strong>cima <strong>de</strong> ω = 1 , <strong>el</strong> logaritmo <strong>de</strong> la<br />

T<br />

amplitud cae rápidam<strong>en</strong>te hacia m<strong>en</strong>os infinito. Esto es principalm<strong>en</strong>te<br />

<strong>de</strong>bido a la pres<strong>en</strong>cia <strong>de</strong> la constante <strong>de</strong> tiempo. En <strong>el</strong> filtro pasa bajo, la<br />

salida pue<strong>de</strong> seguir fi<strong>el</strong>m<strong>en</strong>te a una <strong>en</strong>trada s<strong>en</strong>oidal a bajas frecu<strong>en</strong>cias.<br />

Pero a medida que se increm<strong>en</strong>ta la frecu<strong>en</strong>cia <strong>de</strong> <strong>en</strong>trada, la salida ya no<br />

pue<strong>de</strong> seguir a la <strong>en</strong>trada porque hace falta cierto tiempo para que <strong>el</strong><br />

sistema llegue a configurar esa amplitud. Así, a altas frecu<strong>en</strong>cias, la<br />

amplitud <strong>de</strong> la salida ti<strong>en</strong><strong>de</strong> a cero y <strong>el</strong> ángulo <strong>de</strong> fase ti<strong>en</strong><strong>de</strong> a -90°. Por<br />

tanto, si la función <strong>de</strong> <strong>en</strong>trada conti<strong>en</strong>e muchas armónicas, se reproduc<strong>en</strong><br />

fi<strong>el</strong>m<strong>en</strong>te los compon<strong>en</strong>tes <strong>de</strong> baja frecu<strong>en</strong>cia a la salida, mi<strong>en</strong>tras que los<br />

compon<strong>en</strong>tes <strong>de</strong> alta frecu<strong>en</strong>cia son at<strong>en</strong>uados <strong>en</strong> amplitud y <strong>de</strong>sfasados.<br />

De modo que un sistema <strong>de</strong> primer ord<strong>en</strong> da duplicación exacta solam<strong>en</strong>te<br />

para f<strong>en</strong>óm<strong>en</strong>os constantes o l<strong>en</strong>tam<strong>en</strong>te variables.<br />

Una v<strong>en</strong>taja <strong>de</strong> la repres<strong>en</strong>tación logarítmica es que para factores<br />

1+ jωT<br />

, las curvas <strong>de</strong> logaritmo <strong>de</strong> la<br />

recíprocos, por ejemplo <strong>el</strong> factor ( )<br />

amplitud y ángulo <strong>de</strong> fase solo cambian <strong>de</strong> signo.<br />

El módulo <strong>de</strong> la amplitud para <strong>el</strong> factor ( 1 jωT<br />

)<br />

+ <strong>en</strong> <strong>de</strong>cí<strong>de</strong>les es:<br />

1<br />

20log 1+ jωT<br />

=− 20log = 20log 1+<br />

ω T<br />

1+<br />

jωT<br />

2 2<br />

(131)<br />

Y <strong>el</strong> ángulo <strong>de</strong> fase:<br />

−1<br />

( )<br />

φ( ω) = 1+ jωT =− 1 1+ jωT = tg ωT<br />

(132)<br />

Para estas ecuaciones se pue<strong>de</strong> hacer idéntico estudio al realizado para <strong>el</strong><br />

factor <strong>de</strong> atraso <strong>de</strong> fase para frecu<strong>en</strong>cias bajas ( ωT<br />

1) y para frecu<strong>en</strong>cias<br />

altas ( ωT<br />

1). Los resultadas obt<strong>en</strong>idos serán para <strong>el</strong> diagrama <strong>de</strong><br />

at<strong>en</strong>uación dos rectas, una con p<strong>en</strong>di<strong>en</strong>te 0 dB/década (para las bajas<br />

frecu<strong>en</strong>cias) y la otra con p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> +20 dB/década (para las altas<br />

frecu<strong>en</strong>cias) si<strong>en</strong>do la frecu<strong>en</strong>cia <strong>de</strong> esquina la interacción <strong>de</strong> estas dos<br />

líneas rectas <strong>en</strong> ω = 1 T . Para <strong>el</strong> diagrama <strong>de</strong> fase se obt<strong>en</strong>drá una gráfica<br />

que variará <strong>de</strong> cero a nov<strong>en</strong>ta grados al aum<strong>en</strong>tar la frecu<strong>en</strong>cia <strong>de</strong> cero a<br />

infinito. La figura 45 muestra la curva <strong>de</strong> logaritmo <strong>de</strong> la amplitud y las<br />

1+ jωT<br />

.<br />

asíntotas y la curva <strong>de</strong> fase y la asíntota para <strong>el</strong> factor ( )<br />

- 59 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

50<br />

40<br />

Amplitud, dB<br />

30<br />

20<br />

10<br />

Trazo Exacto<br />

0<br />

Trazo Asintótico<br />

-10<br />

10 -2 10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

100<br />

80<br />

Trazo Exacto<br />

Fase, Grados<br />

60<br />

40<br />

20<br />

Trazo Asintótico<br />

0<br />

-20<br />

10 -2 10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 45: Curvas <strong>de</strong> logaritmo <strong>de</strong>l módulo y fase exactas y asintóticas para <strong>el</strong> factor<br />

1+ jωT<br />

.<br />

( )<br />

En la figura 46 se muestra la curva <strong>de</strong> error <strong>en</strong>tre <strong>el</strong> trazo asintótico <strong>de</strong><br />

at<strong>en</strong>uación y la curva exacta <strong>en</strong> función <strong>de</strong> la frecu<strong>en</strong>cia.<br />

- 60 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

3.5<br />

3<br />

2.5<br />

2<br />

dB<br />

1.5<br />

1<br />

0.5<br />

0<br />

10 -2 10 -1 10 0 10 1 10 2<br />

ω<br />

Figura 46: Curva <strong>de</strong>l error <strong>en</strong> función <strong>de</strong> la frecu<strong>en</strong>cia <strong>de</strong> los lugares <strong>de</strong> amplitud <strong>de</strong><br />

1+ jωT<br />

.<br />

( )<br />

Las formas <strong>de</strong> las curvas <strong>de</strong> ángulo <strong>de</strong> fase son las mismas para cualquier<br />

factor <strong>de</strong> la forma ( ) 1<br />

ω ±<br />

1+ j T . Por tanto es conv<strong>en</strong>i<strong>en</strong>te t<strong>en</strong>er una planilla<br />

para realizar <strong>el</strong> ángulo <strong>de</strong> fase. Esa plantilla podrá ser usada repetidam<strong>en</strong>te<br />

para construir curvas <strong>de</strong> ángulo <strong>de</strong> fase para cualquier función <strong>de</strong> la forma<br />

( 1 jω T ) ± 1<br />

+ .<br />

Si se dispone <strong>de</strong> dicha plantilla, hay que ubicar varios puntos sobre la curva.<br />

ω ±<br />

1+ j T son:<br />

Los ángulos <strong>de</strong> fase para ( ) 1<br />

- 61 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

± 45º para<br />

1<br />

ω =<br />

T<br />

± 26º para<br />

1<br />

ω =<br />

2T<br />

± 5,7º para<br />

1<br />

ω =<br />

10T<br />

± 63,4º para<br />

2<br />

ω =<br />

T<br />

± 83,3º para<br />

10<br />

ω =<br />

T<br />

Se pued<strong>en</strong> analizar ahora <strong>en</strong> forma similar términos que t<strong>en</strong>gan la forma<br />

( 1 j T )<br />

n<br />

ω ±<br />

+ , realizando también una aproximación asintótica. La frecu<strong>en</strong>cia<br />

esquina sigue estando <strong>en</strong> ω = 1 T y las asíntotas son líneas rectas. La<br />

asíntota <strong>de</strong> baja frecu<strong>en</strong>cia es una línea recta horizontal <strong>en</strong> 0 dB , mi<strong>en</strong>tras<br />

que la asíntota <strong>de</strong> alta frecu<strong>en</strong>cia ti<strong>en</strong>e una p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> ± 20 n dB / década ó<br />

± 6 ndB/<br />

octava. El error resultante <strong>de</strong> las aproximaciones asintóticas es n<br />

ω ±<br />

veces <strong>el</strong> correspondi<strong>en</strong>te a ( ) 1<br />

( 1 jω T ) ± 1<br />

+ <strong>en</strong> cada punto <strong>de</strong> frecu<strong>en</strong>cia.<br />

1+ j T . El ángulo <strong>de</strong> fase es n veces <strong>el</strong> <strong>de</strong><br />

En las figuras 47 y 48 se pued<strong>en</strong> ver los diagramas <strong>de</strong> at<strong>en</strong>uación y fase <strong>de</strong><br />

n<br />

ω ±<br />

bo<strong>de</strong> para factores ( 1 j T )<br />

+ hasta n=3.<br />

60<br />

40<br />

n = 3<br />

n = 2<br />

Amplitud, dB<br />

20<br />

0<br />

-20<br />

n = 1<br />

n = -1<br />

n = -2<br />

-40<br />

n = -3<br />

-60<br />

10 -1 10 0 10 1<br />

Frecu<strong>en</strong>cia, rad/s<br />

n<br />

ω ±<br />

Figura 47: Diagramas <strong>de</strong> at<strong>en</strong>uación para ( 1 j T)<br />

+ con n = ± 1, ± 2, ± 3 .<br />

- 62 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

270<br />

n = 3<br />

180<br />

n = 2<br />

Fase, Grados<br />

90<br />

0<br />

-90<br />

-180<br />

n = 1<br />

n = -1<br />

n = -2<br />

n = -3<br />

-270<br />

10 -1 10 0 10 1<br />

Frecu<strong>en</strong>cia, rad/s<br />

n<br />

ω ±<br />

Figura 48: Diagramas <strong>de</strong> fase para ( 1 j T)<br />

+ con n = ± 1, ± 2, ± 3 .<br />

6.4.4.- Diagramas De Bo<strong>de</strong> Para Factores Cuadráticos<br />

Es común que <strong>en</strong> las funciones <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> Sistemas <strong>de</strong> control<br />

aparezcan factores cuadráticos <strong>de</strong> la forma:<br />

2<br />

s<br />

ω<br />

1<br />

2δ<br />

+ s + 1<br />

2<br />

n<br />

ω n<br />

(130)<br />

Que <strong>en</strong> forma frecu<strong>en</strong>cial es:<br />

1<br />

⎛ jω<br />

⎞ ⎛ jω<br />

⎞<br />

1+<br />

2δ<br />

⎜<br />

⎟ +<br />

⎜<br />

⎟<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠<br />

2<br />

(134)<br />

Si δ > 1, se pue<strong>de</strong> expresar este factor cuadrático como un producto <strong>de</strong> dos<br />

<strong>de</strong> primer ord<strong>en</strong> con polos reales. Si 0 < δ < 1 este factor cuadrático es <strong>el</strong><br />

producto <strong>de</strong> dos factores complejos conjugados. Las aproximaciones<br />

asintóticas a las curvas <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia, no son exactas para un<br />

factor cuadrático con valores bajos <strong>de</strong> δ . Esto es porque la amplitud y fase<br />

<strong>de</strong>l factor cuadrático <strong>de</strong>p<strong>en</strong>d<strong>en</strong> <strong>de</strong> la frecu<strong>en</strong>cia esquina y <strong>de</strong> la r<strong>el</strong>ación <strong>de</strong><br />

amortiguami<strong>en</strong>to δ .<br />

Se pue<strong>de</strong> obt<strong>en</strong>er la curva <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia asintótica <strong>de</strong>l<br />

sigui<strong>en</strong>te modo; <strong>el</strong> logaritmo <strong>de</strong>l módulo <strong>de</strong>l factor cuadrático <strong>en</strong> <strong>de</strong>cib<strong>el</strong>es<br />

es:<br />

- 63 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

1<br />

20 log<br />

⎛ jω<br />

⎞ ⎛ jω<br />

⎞<br />

1+<br />

2δ<br />

⎜<br />

⎟ +<br />

⎜<br />

⎟<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠<br />

2<br />

2<br />

⎛ ω ⎞ ⎛<br />

20 log 1 + 2<br />

2 = − ⎜ ⎟<br />

−<br />

2<br />

⎜<br />

ωn<br />

⎝<br />

⎠<br />

2<br />

ω ⎞<br />

δ<br />

⎟<br />

⎝ ωn<br />

⎠<br />

(135)<br />

Para frecu<strong>en</strong>cias bajas, tales que<br />

ω > ω , <strong>el</strong> logaritmo <strong>de</strong> la amplitud es:<br />

2<br />

ω<br />

ω<br />

− 20 log = −40<br />

log dB<br />

(137)<br />

2<br />

ω<br />

n<br />

ω n<br />

n<br />

La ecuación para la asíntota <strong>de</strong> alta frecu<strong>en</strong>cia es una línea recta con<br />

p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> - 40 dB/década <strong>en</strong> coord<strong>en</strong>adas semilogarítmicas, pues:<br />

10ω<br />

ω<br />

− 40 log = −40<br />

− 40 log<br />

(138)<br />

ω<br />

n<br />

ω n<br />

La asíntota <strong>de</strong> baja frecu<strong>en</strong>cia se cortará con la asíntota <strong>de</strong> alta frecu<strong>en</strong>cia<br />

<strong>en</strong>:<br />

ω<br />

0 dB = −20<br />

log 1 = −40<br />

log ∴ ω = ωn<br />

(139)<br />

ω<br />

n<br />

Esta es la frecu<strong>en</strong>cia esquina para <strong>el</strong> factor cuadrático consi<strong>de</strong>rado.<br />

Las dos asíntotas recién consi<strong>de</strong>radas son in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>tes <strong>de</strong>l valor <strong>de</strong> δ .<br />

Cerca <strong>de</strong> la frecu<strong>en</strong>cia<br />

ω = ω , se produce un pico <strong>de</strong> resonancia, como es <strong>de</strong><br />

n<br />

esperar <strong>de</strong> la ecuación 134. El factor <strong>de</strong> amortiguami<strong>en</strong>to δ <strong>de</strong>termina la<br />

amplitud <strong>de</strong> este pico resonante (recordar sección 5). Obviam<strong>en</strong>te exist<strong>en</strong><br />

errores <strong>en</strong> la aproximación con líneas rectas asíntotas. El valor <strong>de</strong>l error<br />

<strong>de</strong>p<strong>en</strong><strong>de</strong> <strong>de</strong>l valor <strong>de</strong> δ . El error es mayor para valores chicos <strong>de</strong> δ . La<br />

figura 49 muestra las curvas exactas <strong>de</strong> logaritmo <strong>de</strong> la amplitud junto con<br />

las aproximaciones asintóticas por líneas rectas y las curvas exactas <strong>de</strong><br />

ángulo <strong>de</strong> fase para <strong>el</strong> factor cuadrático dado por la ecuación 134 para<br />

distintos valores <strong>de</strong> δ .<br />

El ángulo <strong>de</strong> fase <strong>de</strong>l factor cuadrático es:<br />

- 64 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

⎡<br />

ω<br />

⎤<br />

⎢ 2δ<br />

⎥<br />

1<br />

−1<br />

⎢ ωn<br />

⎥<br />

φ ( ω)<br />

= ∠<br />

= −tg<br />

⎢ ⎥<br />

(140)<br />

2<br />

2<br />

⎛ jω<br />

⎞ ⎛ jω<br />

⎞ ⎢ ⎛ ω ⎞<br />

1+<br />

2δ<br />

+<br />

−<br />

⎥<br />

⎜<br />

⎟<br />

⎜<br />

⎟ 1<br />

⎢<br />

⎜<br />

⎟<br />

⎝ ω ⎠ ⎝ ⎠<br />

⎥<br />

n<br />

ωn<br />

⎣ ⎝ ωn<br />

⎠ ⎦<br />

El ángulo <strong>de</strong> fase es tanto función <strong>de</strong> ω como <strong>de</strong> δ . En ω = 0 <strong>el</strong> ángulo <strong>de</strong><br />

fase se iguala a 0º. En la frecu<strong>en</strong>cia esquina <strong>el</strong> ángulo <strong>de</strong> fase es -90º<br />

in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>tem<strong>en</strong>te <strong>de</strong> δ , pues:<br />

−1⎛<br />

2δ<br />

⎞ −1<br />

φ ( ω)<br />

= −tg<br />

⎜ ⎟ = −tg<br />

∞ = −90º<br />

(141)<br />

⎝ 0 ⎠<br />

En ω = ∞ , <strong>el</strong> ángulo <strong>de</strong> fase se vu<strong>el</strong>ve -180°. La curva <strong>de</strong> ángulo <strong>de</strong> fase es<br />

asintótica alre<strong>de</strong>dor <strong>de</strong>l punto <strong>de</strong> inflexión, <strong>el</strong> punto con<strong>de</strong> φ ( ω)<br />

= 90º<br />

.<br />

- 65 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

20<br />

10<br />

δ = 0,1<br />

δ = 0,2<br />

δ = 0,3<br />

Magnitud, dB<br />

0<br />

-10<br />

-20<br />

δ = 0,5<br />

δ = 0,7<br />

δ = 1<br />

-30<br />

-40<br />

-50<br />

10 -1 10 0 10 1<br />

Frequ<strong>en</strong>cia (Rad/s)<br />

0<br />

-20<br />

-40<br />

δ = 0,1<br />

δ = 0,2<br />

δ = 0,3<br />

Fase (Grados)<br />

-60<br />

-80<br />

-100<br />

-120<br />

δ = 0,5<br />

δ = 0,7<br />

δ = 1<br />

-140<br />

-160<br />

-180<br />

-200<br />

10 -1 10 0 10 1<br />

Frequ<strong>en</strong>cia (Rad/s)<br />

Figura 49: Curvas <strong>de</strong>l logaritmo <strong>de</strong>l módulo junto con las asíntotas y las curvas <strong>de</strong> fase <strong>de</strong><br />

la función <strong>de</strong> transfer<strong>en</strong>cia dada por la ecuación 134.<br />

Las correcciones para <strong>el</strong> paso <strong>de</strong> los valores dados por las asíntotas a las<br />

curvas reales se dan <strong>el</strong> la figura 50.<br />

- 66 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

14<br />

12<br />

δ = 0,1<br />

10<br />

δ = 0,15<br />

8<br />

δ = 0,2<br />

Correccióin, dB<br />

6<br />

4<br />

2<br />

δ = 0,3<br />

δ = 0,5<br />

δ = 1<br />

0<br />

-2<br />

-4<br />

-6<br />

10 -1 10 0 10 1<br />

ω/ωn<br />

0<br />

-10<br />

δ = 0.1<br />

Fase, Grados<br />

-20<br />

-30<br />

-40<br />

-50<br />

-60<br />

δ = 0.707<br />

δ =1<br />

δ = 0.2<br />

δ = 0.5<br />

-70<br />

-80<br />

-90<br />

10 -1 10 0 10 1<br />

ω/ωn<br />

Figura 50: correcciones a las gráficas logarítmicas.<br />

- 67 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Las curvas <strong>de</strong> <strong>respuesta</strong> <strong>en</strong> frecu<strong>en</strong>cia para <strong>el</strong> factor + 2δ ( jω<br />

ω ) + ( jω<br />

ω ) 2<br />

1<br />

n<br />

n<br />

pued<strong>en</strong> obt<strong>en</strong>erse simplem<strong>en</strong>te invirti<strong>en</strong>do <strong>el</strong> signo <strong>de</strong>l logaritmo <strong>de</strong> la<br />

amplitud y <strong>el</strong> ángulo <strong>de</strong> fase <strong>de</strong>l factor.<br />

2<br />

⎡ ⎛ ω ⎞ ⎛ ω ⎞ ⎤<br />

⎢1<br />

+ 2δ<br />

⎜ j<br />

⎟ +<br />

⎜ j<br />

⎟ ⎥<br />

⎢⎣<br />

⎝ ωn<br />

⎠ ⎝ ωn<br />

⎠ ⎥⎦<br />

−1<br />

Para obt<strong>en</strong>er las curvas <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> una <strong>de</strong>terminada<br />

función transfer<strong>en</strong>cia cuadrática, primero hay que <strong>de</strong>terminar <strong>el</strong> valor <strong>de</strong> la<br />

frecu<strong>en</strong>cia esquina ω y <strong>el</strong> factor <strong>de</strong> amortiguami<strong>en</strong>to δ . Luego se traza <strong>el</strong><br />

n<br />

diagrama asintótico, y usando las curvas <strong>de</strong> corrección <strong>de</strong> la gráfica 50 se<br />

obti<strong>en</strong><strong>en</strong> las curvas exactas.<br />

6.4.5.-Procedimi<strong>en</strong>to G<strong>en</strong>eral Para Trazar Curvas Logarítmicas De<br />

Respuesta De Frecu<strong>en</strong>cia.<br />

Para realizar las gráficas <strong>de</strong> at<strong>en</strong>uación y fase <strong>de</strong> Bo<strong>de</strong>, para una función <strong>de</strong><br />

transfer<strong>en</strong>cia dada, <strong>en</strong> primer lugar se <strong>de</strong>berá colocar a la citada función <strong>de</strong><br />

transfer<strong>en</strong>cia como un producto <strong>de</strong> factores básicos normalizados como se<br />

vio anteriorm<strong>en</strong>te. Posteriorm<strong>en</strong>te se proce<strong>de</strong> a id<strong>en</strong>tificar las frecu<strong>en</strong>cias<br />

esquinas r<strong>el</strong>acionadas con cada uno <strong>de</strong> estos factores básicos.<br />

Finalm<strong>en</strong>te, se trazan las curvas asintóticas <strong>de</strong> logaritmo <strong>de</strong> la amplitud con<br />

las p<strong>en</strong>di<strong>en</strong>tes a<strong>de</strong>cuadas <strong>en</strong>tre las frecu<strong>en</strong>cias <strong>de</strong> esquina. Se pue<strong>de</strong><br />

obt<strong>en</strong>er la curva exacta, que es cercana a la asíntota, agregando las<br />

correcciones a<strong>de</strong>cuadas, según se vio anteriorm<strong>en</strong>te.<br />

Se pue<strong>de</strong> dibujar la curva <strong>de</strong> ángulo <strong>de</strong> fase para la función <strong>de</strong> transfer<strong>en</strong>cia<br />

dada realizando igual procedimi<strong>en</strong>to que para <strong>el</strong> logaritmo <strong>de</strong>l módulo, o<br />

sea, adicionando la fase correspondi<strong>en</strong>te a cada factor intervini<strong>en</strong>te <strong>en</strong> la<br />

función <strong>de</strong> transfer<strong>en</strong>cia; si se trazan primero las asíntotas posteriorm<strong>en</strong>te<br />

se pue<strong>de</strong> dibujar la curva exacta tomando <strong>de</strong> tablas o gráficas los factores<br />

<strong>de</strong> corrección a<strong>de</strong>cuados.<br />

El uso <strong>de</strong> los diagramas logarítmicos utilizando aproximaciones asintóticas,<br />

exige mucho m<strong>en</strong>os tiempo que otros métodos que pued<strong>en</strong> usarse para<br />

calcular la <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong> una función <strong>de</strong> transfer<strong>en</strong>cia. La<br />

facilidad <strong>de</strong>l trazado <strong>de</strong> las curvas <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia para una<br />

función <strong>de</strong> transfer<strong>en</strong>cia dada, y la facilidad <strong>de</strong> modificación <strong>de</strong> la curva <strong>de</strong><br />

<strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia al agregar comp<strong>en</strong>sación, son las principales<br />

razones por las que se usan tan frecu<strong>en</strong>tem<strong>en</strong>te <strong>en</strong> la práctica los<br />

diagramas logarítmicos.<br />

Ejemplo 4: Trazar <strong>el</strong> diagrama <strong>de</strong> Bo<strong>de</strong> para la función <strong>de</strong> transfer<strong>en</strong>cia:<br />

M ( s)<br />

=<br />

K<br />

n<br />

s<br />

En primer lugar se reemplaza s por j ω :<br />

- 68 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

M ( jω ) =<br />

K<br />

( jω) n<br />

Para este caso lo más conv<strong>en</strong>i<strong>en</strong>te es realizar <strong>el</strong> sigui<strong>en</strong>te procedimi<strong>en</strong>to: <strong>el</strong><br />

diagrama <strong>de</strong> at<strong>en</strong>uación será una recta con p<strong>en</strong>di<strong>en</strong>te – 20.n dB/década,<br />

que cortará al eje <strong>de</strong> las frecu<strong>en</strong>cias <strong>en</strong>:<br />

⎛ K<br />

20 log ⎜<br />

n<br />

⎝ ω<br />

⎞<br />

⎟<br />

⎠<br />

= 0<br />

⇒<br />

ω =<br />

n<br />

K<br />

Mi<strong>en</strong>tras que <strong>el</strong> diagrama <strong>de</strong> fase será una recta con un <strong>de</strong>sfase <strong>de</strong> –n.90°.<br />

20<br />

10<br />

-20 dB/<strong>de</strong>cada<br />

0<br />

Amplitud, dB<br />

-10<br />

-20<br />

ω = n K<br />

-30<br />

-40<br />

-50<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

0<br />

-20<br />

-40<br />

-60<br />

Fase, Grados<br />

-80<br />

-100<br />

-120<br />

-140<br />

-160<br />

-180<br />

-200<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 51: Diagrama <strong>de</strong> Bo<strong>de</strong> para n=1<br />

- 69 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Ejemplo 5: Trazar <strong>el</strong> diagrama <strong>de</strong> Bo<strong>de</strong> para<br />

M =<br />

n<br />

( s)<br />

Ks .<br />

Procedi<strong>en</strong>do <strong>de</strong> igual modo que <strong>el</strong> caso anterior, <strong>el</strong> diagrama <strong>de</strong> at<strong>en</strong>uación<br />

será una recta con p<strong>en</strong>di<strong>en</strong>te +20.n dB/década, que cortará al eje <strong>de</strong> las<br />

frecu<strong>en</strong>cias <strong>en</strong>:<br />

n<br />

( Kω<br />

) = 0 ⇒ ω<br />

n<br />

20 log<br />

=<br />

1<br />

K<br />

Mi<strong>en</strong>tras que <strong>el</strong> diagrama <strong>de</strong> fase será una recta con un <strong>de</strong>sfase <strong>de</strong> <strong>de</strong><br />

+n.90°<br />

40<br />

30<br />

20<br />

+20 dB/década<br />

Amplitud, dB<br />

10<br />

0<br />

-10<br />

-20<br />

ω =<br />

1<br />

n<br />

K<br />

-30<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

180<br />

Fase, Grados<br />

90<br />

0<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 52: Diagrama <strong>de</strong> Bo<strong>de</strong> para n=1<br />

- 70 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Ejemplo 6: Dibujar <strong>el</strong> diagrama <strong>de</strong> Bo<strong>de</strong> para<br />

M s)<br />

=<br />

s<br />

(<br />

2<br />

2<br />

1<br />

( + 2s)<br />

La función <strong>de</strong> transfer<strong>en</strong>cia <strong>en</strong> <strong>el</strong> dominio frecu<strong>en</strong>cial es:<br />

M<br />

2<br />

jω)<br />

=<br />

( jω)<br />

(1 + 2 jω)<br />

(<br />

2<br />

El diagrama <strong>de</strong> at<strong>en</strong>uación para <strong>el</strong> factor<br />

2<br />

2 ( j ω)<br />

es una recta con p<strong>en</strong>di<strong>en</strong>te<br />

-40 dB/década, que corta al eje <strong>de</strong> las frecu<strong>en</strong>cias <strong>en</strong> ω = 2 = 1, 41; mi<strong>en</strong>tras<br />

−<br />

que <strong>el</strong> factor ( ) 1<br />

1+ 2 j ω contribuye con una at<strong>en</strong>uación <strong>de</strong> -20 dB/década a<br />

partir <strong>de</strong> la frecu<strong>en</strong>cia esquina ω = 1 2 = 0, 5 .<br />

El diagrama <strong>de</strong> fase para <strong>el</strong> factor<br />

que <strong>el</strong> factor ( 1+ 2 ω) 1<br />

−<br />

2<br />

2 ( j ω)<br />

es una recta por -180°, mi<strong>en</strong>tras<br />

j contribuye con una fase <strong>de</strong> -90° a partir <strong>de</strong> la<br />

frecu<strong>en</strong>cia esquina ω = 1 2 = 0, 5 .<br />

- 71 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

60<br />

40<br />

20<br />

- 40 dB/década<br />

0<br />

Amplitud, dB<br />

-20<br />

-40<br />

-60<br />

- 60 dB/década<br />

-80<br />

-100<br />

-120<br />

0,50<br />

1,41<br />

-140<br />

10 -1 10 0 10 1 10 2<br />

Freu<strong>en</strong>cia, rad/s<br />

180<br />

Fase, Grados<br />

90<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 53: Diagrama <strong>de</strong> Bo<strong>de</strong> para<br />

M s)<br />

=<br />

s<br />

(<br />

2<br />

2<br />

1<br />

( + 2s)<br />

6.5.- Sistemas De Fase Mínima Y Sistemas De Fase No Mínima.<br />

Las funciones <strong>de</strong> transfer<strong>en</strong>cia que no ti<strong>en</strong><strong>en</strong> polos o ceros <strong>en</strong> <strong>el</strong> semiplano<br />

<strong>de</strong>recho s, son funciones <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> fase mínimas; mi<strong>en</strong>tras que<br />

aqu<strong>el</strong>las que ti<strong>en</strong><strong>en</strong> polos y/o ceros <strong>en</strong> <strong>el</strong> semiplano <strong>de</strong>recho son funciones<br />

<strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> fase no mínima. A los <strong>sistemas</strong> con función <strong>de</strong><br />

transfer<strong>en</strong>cia <strong>de</strong> fase mínima se los llama <strong>sistemas</strong> <strong>de</strong> fase mínima;<br />

mi<strong>en</strong>tras que aqu<strong>el</strong>los con funciones <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> fase no mínima se<br />

los d<strong>en</strong>omina <strong>sistemas</strong> <strong>de</strong> fase no mínima.<br />

- 72 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Para <strong>sistemas</strong> con la misma característica <strong>de</strong> magnitud, <strong>el</strong> rango, <strong>en</strong> ángulo<br />

<strong>de</strong> fase para la función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> fase mínima, es mínima para<br />

todos esos <strong>sistemas</strong>, mi<strong>en</strong>tras <strong>el</strong> rango <strong>en</strong> ángulo <strong>de</strong> fase <strong>de</strong> cualquier<br />

función transfer<strong>en</strong>cia <strong>de</strong> fase no mínima, es mayor que este mínimo.<br />

Estrictam<strong>en</strong>te hablando, una función <strong>de</strong> fase mínima, como su nombre lo<br />

indica, es una función cuya variación <strong>de</strong> fase <strong>en</strong> s<strong>en</strong>tido negativo (retardo<br />

<strong>de</strong> fase) llega a su valor mínimo cuando la frecu<strong>en</strong>cia cambia <strong>de</strong> cero a<br />

infinito. Naturalm<strong>en</strong>te, una función <strong>de</strong> fase no mínima, es la que no posee<br />

esta propiedad.<br />

Sea <strong>el</strong> ejemplo <strong>de</strong> dos <strong>sistemas</strong> cuyas frecu<strong>en</strong>cias <strong>de</strong> transfer<strong>en</strong>cia<br />

frecu<strong>en</strong>cia-les, son respectivam<strong>en</strong>te:<br />

1+<br />

jωT<br />

G1<br />

( jω)<br />

= ;<br />

1+<br />

jωT<br />

1<br />

1−<br />

jωT<br />

G2<br />

( jω)<br />

= ( 0 < T < T1<br />

)<br />

(142)<br />

1+<br />

jωT<br />

1<br />

En la figura 54 se v<strong>en</strong> las configuraciones <strong>de</strong> polos y ceros <strong>de</strong> estos<br />

<strong>sistemas</strong>. Las dos funciones <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidales ti<strong>en</strong><strong>en</strong> las mismas<br />

características <strong>de</strong> amplitud, pero difer<strong>en</strong>tes características <strong>de</strong> ángulo <strong>de</strong><br />

fase, como pue<strong>de</strong> verse <strong>en</strong> la figura 55.<br />

Figura 54: Configuración <strong>de</strong> polos y ceros <strong>de</strong> un sistema <strong>de</strong> fase mínima G<br />

1(<br />

s ) y uno <strong>de</strong><br />

fase no mínima G ( ) .<br />

Esos dos <strong>sistemas</strong> difier<strong>en</strong> <strong>en</strong>tre sí por <strong>el</strong> factor.<br />

G(<br />

jω)<br />

1−<br />

= 1 +<br />

jωT<br />

jωT<br />

La amplitud <strong>de</strong>l factor<br />

2 s<br />

1 − jω T 1+<br />

jωT<br />

es siempre unitaria. Pero <strong>el</strong> ángulo <strong>de</strong><br />

−<br />

fase iguala a − 2tg 1 ωT<br />

y varía <strong>de</strong>s<strong>de</strong> 0 a -180 ° al increm<strong>en</strong>tar ω <strong>de</strong>s<strong>de</strong><br />

cero a infinito.<br />

- 73 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

0<br />

-20<br />

-40<br />

-60<br />

G1(jω)<br />

Fase, Grados<br />

-80<br />

-100<br />

-120<br />

G2(jω)<br />

-140<br />

-160<br />

-180<br />

-200<br />

10 -1 10 0 10 1<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 55: Características <strong>de</strong> ángulo <strong>de</strong> dos <strong>sistemas</strong> G ( ) y G ( ) .<br />

Para un sistema <strong>de</strong> fase mínima, las características <strong>de</strong> amplitud y ángulo <strong>de</strong><br />

fase están r<strong>el</strong>acionados directam<strong>en</strong>te. Esto significa que si se especifica la<br />

curva <strong>de</strong> amplitud <strong>de</strong> un sistema <strong>en</strong> todo <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cias <strong>de</strong>s<strong>de</strong> cero<br />

a infinito, queda unívocam<strong>en</strong>te <strong>de</strong>terminada la curva <strong>de</strong> ángulo <strong>de</strong> fase y<br />

viceversa. Sin embargo esto no se cumple para <strong>sistemas</strong> <strong>de</strong> fase no mínima.<br />

Pue<strong>de</strong> producirse situaciones <strong>de</strong> fase no mínima por dos caminos sigui<strong>en</strong>tes.<br />

Uno es simplem<strong>en</strong>te cuando un sistema incluye un <strong>el</strong>em<strong>en</strong>to o <strong>el</strong>em<strong>en</strong>tos <strong>de</strong><br />

fase no mínimas. La otra situación pue<strong>de</strong> producirse <strong>en</strong> <strong>el</strong> caso que hay un<br />

lazo m<strong>en</strong>or inestable.<br />

Para un sistema <strong>de</strong> fase mínima, <strong>el</strong> ángulo <strong>de</strong> fase <strong>en</strong> ω = ∞ se vu<strong>el</strong>ve -90°<br />

( q − p ); don<strong>de</strong> p y q son los grados <strong>de</strong> los polinomios numerador y<br />

d<strong>en</strong>ominador <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia respectivam<strong>en</strong>te. Para un<br />

sistema <strong>de</strong> fase no mínima, <strong>el</strong> ángulo <strong>de</strong> fase <strong>en</strong> ω = ∞ difiere <strong>de</strong> -90°<br />

( q − p ). En cualquiera <strong>de</strong> ambos <strong>sistemas</strong>, le p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la curva <strong>de</strong>l<br />

logaritmo <strong>de</strong> amplitud <strong>en</strong> ω = ∞ es igual a -20 ( q − p ) dB/década. Por tanto,<br />

es posible <strong>de</strong>tectar si un sistema es o no <strong>de</strong> fase mínima examinando tanto<br />

la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la asíntota <strong>de</strong> alta frecu<strong>en</strong>cia <strong>de</strong> la curva <strong>de</strong> logaritmo <strong>de</strong> la<br />

magnitud, como <strong>el</strong> ángulo <strong>de</strong> fase <strong>en</strong> ω = ∞ . Si la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la curva <strong>de</strong><br />

logaritmo <strong>de</strong> la amplitud cuando ω ti<strong>en</strong><strong>de</strong> a infinito es -20 ( q − p ) dB/década<br />

y <strong>el</strong> ángulo <strong>de</strong> fase <strong>en</strong> ω = ∞ es igual a -90°( q − p ), <strong>el</strong> sistema es un sistema<br />

<strong>de</strong> fase mínima.<br />

1 s<br />

2 s<br />

- 74 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Los <strong>sistemas</strong> <strong>de</strong> fase no mínima son l<strong>en</strong>tos <strong>en</strong> su <strong>respuesta</strong> <strong>de</strong>bido a su mal<br />

comportami<strong>en</strong>to al com<strong>en</strong>zar la <strong>respuesta</strong>. En la mayor parte <strong>de</strong> los<br />

<strong>sistemas</strong> prácticos <strong>de</strong> control, hay que evitar un excesivo retardo <strong>de</strong> fase. Al<br />

diseñar un sistema, si es <strong>de</strong> primordial importancia la rápida v<strong>el</strong>ocidad <strong>de</strong><br />

<strong>respuesta</strong>, no <strong>de</strong>b<strong>en</strong> usarse compon<strong>en</strong>tes <strong>de</strong> fase no mínima (<strong>el</strong> retardo<br />

puro es uno <strong>de</strong> los ejemplos habituales <strong>de</strong> <strong>el</strong>em<strong>en</strong>to <strong>de</strong> fase no mínima que<br />

puedan <strong>en</strong>contrarse pres<strong>en</strong>tes <strong>en</strong> <strong>sistemas</strong> <strong>de</strong> control).<br />

Se <strong>de</strong>be hacer notar que los términos <strong>de</strong> análisis <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia<br />

y <strong>de</strong> proyecto que se analizarán <strong>en</strong> próximos capítulos, son válidos tanto<br />

para <strong>sistemas</strong> <strong>de</strong> fase mínima como <strong>de</strong> fase no mínima.<br />

6.5.2.- Comportami<strong>en</strong>to De Los Sistemas De Fase No Mínima<br />

Se tratan aquí únicam<strong>en</strong>te <strong>sistemas</strong> <strong>de</strong> fase no mínima, <strong>en</strong> cuya función <strong>de</strong><br />

transfer<strong>en</strong>cia <strong>en</strong>tre únicam<strong>en</strong>te un factor <strong>de</strong>sfasador puro, y que son los<br />

más usuales. Estos factores son ( 1 − as ) ( 1+ as)<br />

( a > 0)<br />

, que como ya se vio <strong>el</strong><br />

módulo <strong>de</strong> su función <strong>de</strong> transfer<strong>en</strong>cia es constantem<strong>en</strong>te igual a uno, pero<br />

varía <strong>en</strong> fase <strong>de</strong> 0 a − π cuando ω barre <strong>el</strong> espectro <strong>de</strong> las frecu<strong>en</strong>cias<br />

positivas.<br />

El comportami<strong>en</strong>to <strong>de</strong> estos <strong>sistemas</strong> <strong>de</strong> fase no mínima, ti<strong>en</strong>e una<br />

propiedad característica: <strong>el</strong> arranque vicioso.<br />

Cuando a un sistema cuya salida <strong>de</strong>be reproducir la <strong>en</strong>trada, se le aplica un<br />

escalón <strong>de</strong> 0 a 1, la salida pasa <strong>de</strong> cero a uno con un transitorio intermedio.<br />

En un sistema <strong>de</strong> fase mínima, la salida es positiva <strong>de</strong>spués <strong>de</strong> la variación<br />

2<br />

dc(<br />

t)<br />

dc(<br />

t)<br />

d c(<br />

t)<br />

brusca <strong>de</strong> la <strong>en</strong>trada: > 0 ó si 0 ∴ > 0<br />

2<br />

dt<br />

dt dt<br />

t= 0<br />

t= 0<br />

t= 0<br />

En <strong>el</strong> caso <strong>de</strong> un sistema <strong>de</strong> fase no mínima, por <strong>el</strong> contrario, la salida<br />

empieza a <strong>de</strong>crecer antes <strong>de</strong> seguir la <strong>en</strong>trada.<br />

Esta propiedad se pone <strong>en</strong> evid<strong>en</strong>cia fácilm<strong>en</strong>te a partir <strong>de</strong> la forma<br />

matemática <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia.<br />

La <strong>respuesta</strong> a un escalón unitario <strong>de</strong> un sistema <strong>de</strong> fase mínima vi<strong>en</strong>e dada<br />

por:<br />

C(<br />

s)<br />

1 A s<br />

+ ...... + A<br />

m<br />

m<br />

0<br />

= (143)<br />

n<br />

s Bns<br />

+ ...... + B0<br />

Don<strong>de</strong> B y A son positivos ( B , porque <strong>el</strong> sistema se supone estable. A<br />

porque es <strong>de</strong> fase mínima); la aplicación <strong>de</strong>l teorema <strong>de</strong>l valor inicial<br />

<strong>de</strong>muestra, inmediatam<strong>en</strong>te que la primera <strong>de</strong>rivada no nula <strong>de</strong> c (t)<br />

es<br />

positiva. La figura 56 muestra la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la <strong>respuesta</strong>.<br />

- 75 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Figura 56: Respuesta al escalón unitario <strong>de</strong> un sistema <strong>de</strong> fase mínima<br />

Por <strong>el</strong> contrario, <strong>en</strong> <strong>el</strong> caso <strong>de</strong> la frecu<strong>en</strong>cia <strong>de</strong> un <strong>de</strong>fasador puro, se<br />

t<strong>en</strong>drá:<br />

1<br />

C(<br />

s)<br />

=<br />

s<br />

A<br />

B<br />

m<br />

n<br />

s<br />

s<br />

m<br />

n<br />

+ ...... + A<br />

+ ...... + B<br />

0<br />

0<br />

( 1−<br />

Ks)<br />

( 1+<br />

Ks)<br />

(144)<br />

Con K > 0 y t<strong>en</strong>i<strong>en</strong>do siempre a la vista la hipótesis A m<br />

y<br />

B n<br />

mayor que<br />

cero; <strong>en</strong> estas condiciones se observa fácilm<strong>en</strong>te mediante la aplicación <strong>de</strong>l<br />

teorema <strong>de</strong>l valor inicial, que la primera <strong>de</strong>rivada no nula es ahora negativa;<br />

tal como se muestra <strong>en</strong> la figura 57.<br />

Pue<strong>de</strong> también observarse que, <strong>en</strong> pres<strong>en</strong>cia <strong>de</strong> un <strong>de</strong>sfasador puro la<br />

<strong>respuesta</strong> c (t)<br />

a un escalón es <strong>de</strong> la forma:<br />

dc0<br />

( t)<br />

c ( t)<br />

= c0 ( t)<br />

+ K<br />

K > 0<br />

(145)<br />

dt<br />

Figura 57: Respuesta al escalón unitario <strong>de</strong> un sistema <strong>de</strong> fase no mínima.<br />

Don<strong>de</strong> c ( ) es la <strong>respuesta</strong> unitaria <strong>de</strong> un sistema <strong>de</strong> fase mínima. En esta<br />

0<br />

t<br />

expresión c ( ) es nulo para t < 0 y crece a partir <strong>de</strong> t = 0 . Inmediatam<strong>en</strong>te<br />

0<br />

t<br />

- 76 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

<strong>de</strong>spués <strong>de</strong> t = 0 , <strong>el</strong> segundo término <strong>de</strong> c (t)<br />

es negativo e importante <strong>en</strong><br />

valor absoluto, es <strong>de</strong>cir predominante.<br />

6.5.3.- Ejemplos De Sistemas De Fase No Mínima<br />

1.- Regulación <strong>de</strong>l niv<strong>el</strong> <strong>de</strong> líquido <strong>en</strong> una cal<strong>de</strong>ra don<strong>de</strong> <strong>el</strong> agua está <strong>en</strong><br />

ebullición. Al agregar agua para <strong>el</strong>evar, <strong>el</strong> niv<strong>el</strong>, <strong>el</strong> primer resultado es<br />

<strong>de</strong>t<strong>en</strong>er durante algunos instantes la ebullición, es <strong>de</strong>cir, hacer <strong>de</strong>sc<strong>en</strong><strong>de</strong>r <strong>el</strong><br />

niv<strong>el</strong>.<br />

2.- Un avión <strong>en</strong> vu<strong>el</strong>o longitudinal es un sistema <strong>de</strong> fase no mínima, si se<br />

consi<strong>de</strong>ra como <strong>respuesta</strong> la v<strong>el</strong>ocidad vertical <strong>de</strong>l c<strong>en</strong>tro <strong>de</strong> gravedad. La<br />

razón es que cuando se actúa <strong>el</strong> timón para hacer subir <strong>el</strong> avión, <strong>el</strong> primer<br />

movimi<strong>en</strong>to <strong>de</strong> este es <strong>de</strong> <strong>de</strong>sc<strong>en</strong>so, como consecu<strong>en</strong>cia <strong>de</strong> la disminución<br />

<strong>de</strong> sust<strong>en</strong>tación resultante <strong>de</strong>l movimi<strong>en</strong>to <strong>de</strong>l timón.<br />

Los <strong>el</strong>em<strong>en</strong>tos <strong>de</strong> fase no mínima se <strong>en</strong>cu<strong>en</strong>tran con bastante frecu<strong>en</strong>cia <strong>en</strong><br />

los <strong>sistemas</strong> <strong>de</strong> control. Es importante conocerlos por <strong>el</strong> efecto perjudicial y<br />

no comp<strong>en</strong>sable que ejerc<strong>en</strong> sobre la precisión <strong>de</strong> los <strong>sistemas</strong> controlados.<br />

6.6.5.- Retardo De Transporte O Retardo Puro.<br />

El retardo <strong>de</strong> transporte es un comportami<strong>en</strong>to <strong>de</strong> fase no mínima y ti<strong>en</strong>e<br />

un retardo <strong>de</strong> fase excesivo sin at<strong>en</strong>uación <strong>en</strong> altas frecu<strong>en</strong>cias. Sea <strong>el</strong><br />

retardo <strong>de</strong> transporte dado por:<br />

M<br />

− jωT<br />

( jω)<br />

= e<br />

(146)<br />

La amplitud es siempre igual a la unidad por:<br />

M ( jω ) = cos ωT<br />

− j s<strong>en</strong> ωT<br />

= 1<br />

(147)<br />

Por tanto, <strong>el</strong> logaritmo <strong>de</strong> la amplitud <strong>de</strong>l retardo <strong>de</strong>l transporte es igual a 0<br />

dB.<br />

El ángulo <strong>de</strong> fase <strong>de</strong>l retardo puro es:<br />

φ( ω)<br />

= ∠M<br />

( ω)<br />

= −ωT<br />

radianes = −57,<br />

3 ωT<br />

grados<br />

(148)<br />

El ángulo <strong>de</strong> fase varía linealm<strong>en</strong>te con la frecu<strong>en</strong>cia ω . En la figura 58 se<br />

muestra la característica <strong>de</strong> ángulo <strong>de</strong> fase <strong>de</strong>l retardo <strong>de</strong> transporte.<br />

- 77 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

0<br />

-20<br />

-40<br />

-60<br />

Fase, Grados<br />

-80<br />

-100<br />

-120<br />

-140<br />

M(jω) = e jωT<br />

|M(jω)| = 0 dB<br />

-160<br />

-180<br />

-200<br />

10 -1 10 0 10 1<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 58: Característica <strong>de</strong> ángulo <strong>de</strong> fase <strong>de</strong>l retardo <strong>de</strong> transporte.<br />

Ejemplo 6: Trazar <strong>el</strong> diagrama <strong>de</strong> Bo<strong>de</strong> <strong>de</strong> la sigui<strong>en</strong>te función <strong>de</strong><br />

transfer<strong>en</strong>cia:<br />

−sL<br />

e<br />

G(<br />

s)<br />

=<br />

1+<br />

sT<br />

La función <strong>de</strong> transfer<strong>en</strong>cia frecu<strong>en</strong>cial es:<br />

El logaritmo <strong>de</strong> la amplitud es:<br />

− jωL<br />

e<br />

G(<br />

jω)<br />

=<br />

1+<br />

jωT<br />

20 log G(<br />

jω)<br />

= 20 log<br />

El ángulo <strong>de</strong> fase es:<br />

e<br />

− jωL<br />

+ 20 log<br />

1<br />

1+<br />

jωT<br />

= 0 + 20 log<br />

1<br />

1+<br />

jωT<br />

( 1<br />

−1<br />

1+<br />

jωT) − jωL<br />

G( jω)<br />

= e − =−ωL−tg ωT<br />

Las curvas <strong>de</strong> logaritmo <strong>de</strong> la amplitud y ángulo <strong>de</strong> fase para esta función<br />

<strong>de</strong> transfer<strong>en</strong>cia, con L = 0, 5 y T = 1 aparec<strong>en</strong> <strong>en</strong> la figura 59.<br />

- 78 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Figura 59.<br />

6.7.1.- R<strong>el</strong>ación Entre El Tipo De Sistema Y La Curva D<strong>el</strong> Logaritmo<br />

De La Amplitud.<br />

Los coefici<strong>en</strong>tes <strong>de</strong> error <strong>de</strong> posición estática, v<strong>el</strong>ocidad y ac<strong>el</strong>eración<br />

<strong>de</strong>scrib<strong>en</strong> <strong>el</strong> comportami<strong>en</strong>to <strong>en</strong> baja frecu<strong>en</strong>cia <strong>de</strong> los <strong>sistemas</strong> tipo 0, tipo<br />

1 y tipo 2, respectivam<strong>en</strong>te. Para un sistema dado, solo uno <strong>de</strong> los<br />

coefici<strong>en</strong>tes <strong>de</strong> error estático es finito y significativo (cuánto más gran<strong>de</strong> es<br />

<strong>el</strong> coefici<strong>en</strong>te <strong>de</strong> error estático finito, mayor es la ganancia <strong>de</strong>l lazo cuando<br />

ω ti<strong>en</strong><strong>de</strong> a 0).<br />

- 79 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

El tipo <strong>de</strong> sistema <strong>de</strong>termina la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la curva <strong>de</strong> logaritmo <strong>de</strong><br />

amplitud a bajas frecu<strong>en</strong>cias. Por tanto, se pue<strong>de</strong> <strong>de</strong>terminar la información<br />

refer<strong>en</strong>te a la exist<strong>en</strong>cia y valor <strong>de</strong>l error estacionario <strong>de</strong> un sistema <strong>de</strong><br />

control ante <strong>de</strong>terminada <strong>en</strong>trada por la observación <strong>de</strong> la zona <strong>de</strong> baja<br />

frecu<strong>en</strong>cia <strong>de</strong> la curva <strong>el</strong> logaritmo <strong>de</strong> la amplitud.<br />

Es importante aclarar que para las <strong>de</strong>terminaciones <strong>de</strong> coefici<strong>en</strong>tes que se<br />

hac<strong>en</strong> a continuación, se graficará <strong>el</strong> módulo <strong>en</strong> <strong>de</strong>cí<strong>de</strong>les <strong>de</strong> la función <strong>de</strong><br />

transfer<strong>en</strong>cia <strong>de</strong> lazo abierto <strong>de</strong> un sistema realim<strong>en</strong>tado negativam<strong>en</strong>te.<br />

6.7.2.- Determinación De Los Coefici<strong>en</strong>tes De Error De Posición<br />

Estático.<br />

En la figura 60 pres<strong>en</strong>ta un ejemplo <strong>de</strong>l diagrama <strong>de</strong>l logaritmo <strong>de</strong> la amplitud<br />

<strong>de</strong> un sistema tipo cero. En tal sistema, <strong>el</strong> valor <strong>de</strong> la función <strong>de</strong><br />

transfer<strong>en</strong>cia <strong>de</strong> lazo abierto G ( jω)<br />

H ( jω)<br />

iguala a Kp <strong>en</strong> bajas frecu<strong>en</strong>cias,<br />

o sea:<br />

lim G( jω) H( jω)<br />

= Kp<br />

ω→0<br />

Figura 60: Curva <strong>de</strong> logaritmo <strong>de</strong> la amplitud <strong>de</strong> un sistema tipo 0.<br />

De aquí se llega a que la asíntota <strong>de</strong> baja frecu<strong>en</strong>cia es una recta horizontal<br />

<strong>en</strong> 20 log Kp dB.<br />

6.7.3.- Determinación De Los Coefici<strong>en</strong>tes De Error Estático De<br />

V<strong>el</strong>ocidad.<br />

La figura 61, pres<strong>en</strong>ta un ejemplo <strong>de</strong> diagrama <strong>de</strong> logaritmo <strong>de</strong> la amplitud<br />

<strong>de</strong> un sistema tipo 1. La intersección <strong>de</strong>l segm<strong>en</strong>to inicial -20 dB/década (o<br />

su continuación) con la recta ω = 1, ti<strong>en</strong>e la magnitud 20 log Kv . Pue<strong>de</strong> verse<br />

esto como sigue <strong>en</strong> un sistema tipo 1:<br />

- 80 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Kv<br />

G ( jω)<br />

H ( jω)<br />

=<br />

jω<br />

para ω


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

F<br />

ω<br />

2<br />

= ;<br />

J<br />

2 K<br />

ω<br />

3<br />

=<br />

(153)<br />

J<br />

Como:<br />

K<br />

ω<br />

1<br />

= Kv =<br />

(154)<br />

F<br />

Resulta<br />

ω ω<br />

ω ω<br />

=<br />

1<br />

2<br />

1 3<br />

2<br />

= ω3<br />

⇒<br />

(155)<br />

ω3<br />

ω2<br />

En <strong>el</strong> diagrama logarítmico:<br />

log ω − ω<br />

(156)<br />

1<br />

log ω3<br />

= log ω3<br />

− log<br />

2<br />

Así, <strong>el</strong> punto ω 3<br />

está justam<strong>en</strong>te <strong>en</strong> <strong>el</strong> medio <strong>de</strong> los puntos ω 1<br />

y ω 2<br />

. La<br />

r<strong>el</strong>ación <strong>de</strong> amortiguami<strong>en</strong>to δ <strong>de</strong> este sistema es:<br />

F ω2<br />

δ = =<br />

2 JK 2ω<br />

3<br />

(157)<br />

6.7.3.- Determinación De Los Coefici<strong>en</strong>tes De Error Estático De<br />

Ac<strong>el</strong>eración.<br />

La figura 62, muestra un ejemplo <strong>de</strong> diagrama <strong>de</strong> logaritmo <strong>de</strong> la amplitud<br />

<strong>de</strong> un sistema tipo 2. La intersección <strong>de</strong>l segm<strong>en</strong>to inicial -40 dB/década (o<br />

su continuación) con la recta ω = 1, ti<strong>en</strong>e la magnitud 20 log Ka .<br />

Figura 62: Curva <strong>de</strong> logaritmo <strong>de</strong> la amplitud <strong>de</strong> un sistema tipo 2.<br />

- 82 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Como a bajas frecu<strong>en</strong>cias:<br />

G ( jω)<br />

H ( jω)<br />

=<br />

Entonces:<br />

Ka<br />

( jω) 2<br />

20 log<br />

Ka<br />

( jω)<br />

2<br />

=<br />

ω= 1<br />

20 log<br />

Ka<br />

La frecu<strong>en</strong>cia ω<br />

a<br />

<strong>en</strong> la intersección <strong>de</strong>l segm<strong>en</strong>to inicial -40 dB/década (o su<br />

prolongación) con la recta <strong>de</strong> 0 dB, da, numéricam<strong>en</strong>te, la raíz cuadrada <strong>de</strong><br />

Ka :<br />

Ka<br />

20 log = 20 log 1 = 0<br />

2<br />

( jω)<br />

(158)<br />

Queda:<br />

ω Ka<br />

(159)<br />

a =<br />

- 83 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

7.- Análisis Experim<strong>en</strong>tal De La Respuesta Frecu<strong>en</strong>cial.<br />

7.1.- Determinación Experim<strong>en</strong>tal De Funciones De Transfer<strong>en</strong>cia.<br />

El primer paso <strong>en</strong> <strong>el</strong> análisis y proyecto <strong>de</strong> control es <strong>de</strong>sarrollar, un mo<strong>de</strong>lo<br />

matemático <strong>de</strong> la planta <strong>en</strong> estudio. Obt<strong>en</strong>er un mo<strong>de</strong>lo analíticam<strong>en</strong>te<br />

pue<strong>de</strong> ser bastante difícil. Pue<strong>de</strong> ser necesario obt<strong>en</strong>erlo por medio <strong>de</strong> un<br />

análisis experim<strong>en</strong>tal. La importancia <strong>de</strong> los métodos <strong>de</strong> <strong>respuesta</strong> <strong>de</strong><br />

frecu<strong>en</strong>cia es que la función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong>l sistema, o <strong>de</strong> cualquier<br />

compon<strong>en</strong>te <strong>de</strong> un sistema pue<strong>de</strong> ser <strong>de</strong>terminado por simples mediciones<br />

<strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia.<br />

Si se han medido la r<strong>el</strong>ación <strong>de</strong> amplitud y <strong>de</strong>splazami<strong>en</strong>to <strong>de</strong> fase <strong>en</strong> un<br />

número sufici<strong>en</strong>te <strong>de</strong> frecu<strong>en</strong>cias d<strong>en</strong>tro <strong>de</strong>l rango <strong>de</strong> frecu<strong>en</strong>cias <strong>de</strong> interés,<br />

se las pue<strong>de</strong> repres<strong>en</strong>tar <strong>en</strong> <strong>el</strong> diagrama <strong>de</strong> Bo<strong>de</strong>. Luego se pue<strong>de</strong><br />

<strong>de</strong>terminar la función <strong>de</strong> transfer<strong>en</strong>cia por aproximaciones asintóticas. Se<br />

construy<strong>en</strong> curvas asintóticas <strong>de</strong>l logaritmo <strong>de</strong>l módulo constituido por<br />

varios segm<strong>en</strong>tos. Con algunos tanteos y correcciones <strong>en</strong> frecu<strong>en</strong>cias<br />

esquinas, g<strong>en</strong>eralm<strong>en</strong>te es posible hallar una aproximación bastante<br />

cercana a la curva. (Nótese que si se repres<strong>en</strong>ta la frecu<strong>en</strong>cia <strong>en</strong> ciclos por<br />

segundo <strong>en</strong> lugar <strong>de</strong> radianes por segundo, hay que convertirlas frecu<strong>en</strong>cias<br />

esquinas a radianes por segundo antes <strong>de</strong> calcular las constantes <strong>de</strong><br />

tiempo).<br />

Otras veces se realiza la <strong>de</strong>terminación <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia<br />

frecu<strong>en</strong>cial <strong>en</strong> forma experim<strong>en</strong>tal, cuando se <strong>de</strong>sea confirmar un análisis<br />

matemático <strong>de</strong> un sistema, o sea, chequear si la función <strong>de</strong> transfer<strong>en</strong>cia<br />

obt<strong>en</strong>ida analíticam<strong>en</strong>te es válida.<br />

7.2.- G<strong>en</strong>eradores De Señal S<strong>en</strong>oidal<br />

Para realizar la prueba <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cias se <strong>de</strong>be disponer <strong>de</strong><br />

g<strong>en</strong>eradores <strong>de</strong> señal s<strong>en</strong>oidal a<strong>de</strong>cuadas. La señal pue<strong>de</strong> t<strong>en</strong>er que estar<br />

<strong>en</strong> forma mecánica, <strong>el</strong>éctrica o neumática.<br />

Los rangos <strong>de</strong> frecu<strong>en</strong>cias necesarios para la verificación están<br />

aproximadam<strong>en</strong>te <strong>en</strong>tre 0,001 - 10 cps para <strong>sistemas</strong> <strong>de</strong> constante <strong>de</strong><br />

tiempo <strong>el</strong>evada; hasta 0,1 - 1000 cps para <strong>sistemas</strong> con constantes <strong>de</strong><br />

tiempo pequeñas. La señal s<strong>en</strong>oidal <strong>de</strong>be estar razonablem<strong>en</strong>te libre <strong>de</strong><br />

armónicas o distorsión.<br />

Para rangos <strong>de</strong> muy baja frecu<strong>en</strong>cia (Por <strong>de</strong>bajo <strong>de</strong> 0.01 cps), pue<strong>de</strong><br />

utilizarse un g<strong>en</strong>erador <strong>de</strong> señal mecánico (junto con un transductor<br />

neumático o <strong>el</strong>éctrico si es necesario). Para <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cias <strong>de</strong>s<strong>de</strong><br />

0.01 - 1000 cps, pue<strong>de</strong> usarse un g<strong>en</strong>erador <strong>el</strong>éctrico <strong>de</strong> señal a<strong>de</strong>cuado<br />

(junto con un transductor si es necesario).<br />

La posibilidad <strong>de</strong> obt<strong>en</strong>er g<strong>en</strong>eradores <strong>de</strong> señal s<strong>en</strong>oidal con esas<br />

características es usualm<strong>en</strong>te fácil, por ejemplo:<br />

a - Un simple mecanismo <strong>de</strong> bi<strong>el</strong>a-maniv<strong>el</strong>a produce variaciones s<strong>en</strong>oidales<br />

<strong>de</strong> <strong>de</strong>splazami<strong>en</strong>to lineal.<br />

- 84 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

b - Un <strong>de</strong>splazami<strong>en</strong>to s<strong>en</strong>oidal <strong>en</strong> <strong>el</strong> extremo <strong>de</strong> un resorte lineal produce<br />

variaciones s<strong>en</strong>oidales <strong>de</strong> la fuerza ejercida por <strong>el</strong> otro extremo <strong>de</strong>l resorte<br />

<strong>en</strong> su reporte.<br />

c - Corri<strong>en</strong>tes <strong>el</strong>éctricas alternas s<strong>en</strong>oidales son <strong>de</strong>bidas a variaciones s<strong>en</strong>oi<br />

dales <strong>de</strong> t<strong>en</strong>sión.<br />

d - Un pistón correspondi<strong>en</strong>te a un cilindro <strong>de</strong> aire cerrado pue<strong>de</strong> producir<br />

variaciones s<strong>en</strong>oidales <strong>de</strong> presión <strong>en</strong> <strong>el</strong> cilindro.<br />

e - G<strong>en</strong>eradores s<strong>en</strong>oidales <strong>el</strong>ectrónicos son <strong>de</strong> uso común <strong>en</strong> los<br />

laboratorios <strong>de</strong> análisis.<br />

7.3.- Procedimi<strong>en</strong>to Para Obt<strong>en</strong>er La Función De Transfer<strong>en</strong>cia De<br />

Un Sistema Lineal En Forma Experim<strong>en</strong>tal<br />

El procedimi<strong>en</strong>to conv<strong>en</strong>i<strong>en</strong>te a seguir para <strong>el</strong> análisis es <strong>el</strong> sigui<strong>en</strong>te:<br />

1- Aplicar una <strong>en</strong>trada s<strong>en</strong>oidal r( t)<br />

= H s<strong>en</strong>ωt<br />

a una variable apropiada <strong>de</strong>l<br />

sistema consi<strong>de</strong>rado como <strong>en</strong>trada. El valor inicial <strong>de</strong> la variable <strong>de</strong> <strong>en</strong>tra da<br />

pue<strong>de</strong> ser cualquier valor <strong>de</strong> estado constante. La excitación s<strong>en</strong>oidal <strong>de</strong><br />

<strong>en</strong>trada es superpuesta sobre <strong>el</strong> valor inicial <strong>de</strong> estado constante.<br />

2- Para un valor s<strong>el</strong>eccionado <strong>de</strong> frecu<strong>en</strong>cia ω y una amplitud constante H ,<br />

<strong>de</strong>jar que los efectos transitorios <strong>en</strong> la <strong>respuesta</strong> <strong>de</strong>l sistema se extingan.<br />

La variable <strong>el</strong>egida para repres<strong>en</strong>tar la <strong>respuesta</strong> <strong>de</strong>l sistema t<strong>en</strong>drá luego<br />

una variación s<strong>en</strong>oidal con una amplitud constante y la misma frecu<strong>en</strong>cia<br />

que la <strong>en</strong>trada. Se mi<strong>de</strong> la amplitud <strong>de</strong> la variable <strong>de</strong> <strong>respuesta</strong> y su fase<br />

r<strong>el</strong>ativa a la <strong>en</strong>trada s<strong>en</strong>oidal.<br />

Frecu<strong>en</strong>tem<strong>en</strong>te ambas variables, <strong>de</strong> <strong>en</strong>trada y salida, pued<strong>en</strong> ser<br />

visualizadas <strong>en</strong> un osciloscopio o graficadas <strong>en</strong> un registrador.<br />

3- Mant<strong>en</strong>i<strong>en</strong>do la amplitud constante H <strong>de</strong> <strong>en</strong>trada, repetir <strong>el</strong><br />

procedimi<strong>en</strong>to para distintos valores <strong>de</strong> frecu<strong>en</strong>cia, recorri<strong>en</strong>do la banda <strong>de</strong><br />

frecu<strong>en</strong>cias para <strong>el</strong> cual <strong>el</strong> sistema respon<strong>de</strong>.<br />

4- Graficar los resultados obt<strong>en</strong>idos <strong>en</strong> las coord<strong>en</strong>adas <strong>de</strong> Bo<strong>de</strong>, a partir <strong>de</strong><br />

la cual según se verá, se podrá <strong>de</strong>terminar la función <strong>de</strong> transfer<strong>en</strong>cia<br />

buscada.<br />

7.4.- Determinación De Funciones Transfer<strong>en</strong>cia De Fase Mínima<br />

A Partir De Diagramas De Bo<strong>de</strong><br />

Como se indicó previam<strong>en</strong>te, pue<strong>de</strong> <strong>de</strong>terminarse si un sistema es <strong>de</strong> fase<br />

mínima o no parti<strong>en</strong>do <strong>de</strong> las curvas <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia<br />

examinando las características <strong>de</strong> alta frecu<strong>en</strong>cia.<br />

Para <strong>de</strong>terminar la función <strong>de</strong> transfer<strong>en</strong>cia se comi<strong>en</strong>zan por trazar<br />

asíntotas a la curva <strong>de</strong> logaritmo <strong>de</strong> la amplitud obt<strong>en</strong>ida<br />

experim<strong>en</strong>talm<strong>en</strong>te. Las asíntotas <strong>de</strong>b<strong>en</strong> t<strong>en</strong>er p<strong>en</strong>di<strong>en</strong>tes <strong>de</strong> múltiplos <strong>de</strong><br />

± 20dB / década. Si la p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la curva <strong>de</strong>l logaritmo obt<strong>en</strong>ida<br />

- 85 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

experim<strong>en</strong>talm<strong>en</strong>te varía <strong>de</strong>s<strong>de</strong> − 20 a − 40dB / década <strong>en</strong> ω = ω1, es claro que<br />

existe un factor<br />

⎛<br />

⎜1+<br />

⎝<br />

Si la p<strong>en</strong>di<strong>en</strong>te varía <strong>en</strong><br />

cuadrático <strong>de</strong> la forma<br />

transfer<strong>en</strong>cia.<br />

−1<br />

⎞<br />

j<br />

⎛ ⎞<br />

⎜ω<br />

⎟⎟<br />

<strong>en</strong> la función <strong>de</strong> transfer<strong>en</strong>cia.<br />

⎝ ω1<br />

⎠⎠<br />

− 40 dB / década <strong>en</strong> ω = ω2<br />

, <strong>de</strong>be haber un factor<br />

−1<br />

2<br />

⎡ ⎛ ω ⎞ ⎛ ⎞ ⎤<br />

⎢1<br />

2<br />

⎟ ⎥<br />

⎢<br />

+ ω<br />

+ δ<br />

⎜ j ⎟<br />

⎜ j <strong>en</strong> la función <strong>de</strong><br />

⎣ ⎝ ω2<br />

⎠ ⎝ ω2<br />

⎠ ⎥⎦<br />

La frecu<strong>en</strong>cia natural no amortiguada <strong>de</strong> este factor cuadrático es igual a la<br />

frecu<strong>en</strong>cia esquina ω<br />

2<br />

. Se pue<strong>de</strong> <strong>de</strong>terminar la r<strong>el</strong>ación <strong>de</strong><br />

amortiguami<strong>en</strong>to δ <strong>de</strong> la curva <strong>de</strong>l logaritmo <strong>de</strong> la amplitud obt<strong>en</strong>ida<br />

experim<strong>en</strong>talm<strong>en</strong>te midi<strong>en</strong>do <strong>el</strong> valor <strong>de</strong>l pico resonante cerca <strong>de</strong> la<br />

frecu<strong>en</strong>cia esquina ω 2<br />

y comparando este con las curvas que aparec<strong>en</strong> <strong>en</strong><br />

la figura 50<br />

Una vez <strong>de</strong>terminado los factores <strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia M ( jω)<br />

se<br />

pue<strong>de</strong> <strong>de</strong>terminar la ganancia <strong>de</strong> la porción <strong>de</strong> baja frecu<strong>en</strong>cia <strong>de</strong> la curva<br />

logaritmo amplitud.<br />

Como los términos tales como + j( ω ) y 1 2δ<br />

( j ω ω ) + ( jω<br />

ω ) 2<br />

1 ω 1<br />

+<br />

2<br />

2<br />

ti<strong>en</strong>d<strong>en</strong> a<br />

la unidad cuando ω ti<strong>en</strong><strong>de</strong> a cero, a frecu<strong>en</strong>cias muy bajas se pue<strong>de</strong> escribir<br />

la función <strong>de</strong> transfer<strong>en</strong>cia s<strong>en</strong>oidal M ( jω)<br />

:<br />

K<br />

lim M ( jω)<br />

=<br />

(160)<br />

0<br />

ω →<br />

( jω) λ<br />

En muchos <strong>sistemas</strong> prácticos, λ es igual a 0, 1 ó 2.<br />

1- Para λ = 0 , o sistema tipo 0:<br />

M ( jω ) = K ; para ω


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Que indica que la asíntota <strong>de</strong> baja frecu<strong>en</strong>cia ti<strong>en</strong>e la p<strong>en</strong>di<strong>en</strong>te<br />

− 20dB / década . La frecu<strong>en</strong>cia a la cual la asíntota <strong>de</strong> baja frecu<strong>en</strong>cia (o su<br />

prolongación) intersecta la recta 0 dB es numéricam<strong>en</strong>te igual a K.<br />

3- Para λ = 2 , o <strong>sistemas</strong> <strong>de</strong> tipo 2:<br />

K<br />

M ( jω)<br />

= ; para ω


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

Figura 63: a) Curva <strong>de</strong> logaritmo <strong>de</strong>l módulo <strong>de</strong> un sistema tipo 0. b) Curvas <strong>de</strong> logaritmo<br />

<strong>de</strong>l módulo <strong>de</strong> <strong>sistemas</strong> tipo uno. C) Curvas <strong>de</strong> logaritmo <strong>de</strong>l módulo <strong>de</strong> <strong>sistemas</strong> tipo dos.<br />

(Las p<strong>en</strong>di<strong>en</strong>tes indicadas están <strong>en</strong> dB/década).<br />

- 88 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

7.5.- Funciones De Transfer<strong>en</strong>cia De Fase No Mínima<br />

Si <strong>en</strong> <strong>el</strong> extremo <strong>de</strong> altas frecu<strong>en</strong>cias, <strong>el</strong> ángulo <strong>de</strong> fase calculado es 180°<br />

m<strong>en</strong>or que <strong>el</strong> retardo <strong>de</strong> fase obt<strong>en</strong>ido experim<strong>en</strong>talm<strong>en</strong>te, uno <strong>de</strong> los ceros<br />

<strong>de</strong> la función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong>be haber estado <strong>en</strong> <strong>el</strong> semiplano <strong>de</strong>recho s<br />

<strong>en</strong> lugar <strong>de</strong>l semiplano izquierdo s. Si <strong>el</strong> retardo <strong>de</strong> fase calculado difiere <strong>de</strong>l<br />

ángulo <strong>de</strong> fase obt<strong>en</strong>ida experim<strong>en</strong>talm<strong>en</strong>te por una v<strong>el</strong>ocidad constante <strong>de</strong><br />

modificación <strong>de</strong> fase, hay retardo <strong>de</strong> transporte o tiempo muerto. Si se<br />

supone que la función <strong>de</strong> transfer<strong>en</strong>cia es <strong>de</strong> la forma:<br />

G(<br />

s)<br />

= e<br />

−st<br />

Don<strong>de</strong> G (s)<br />

es una r<strong>el</strong>ación <strong>de</strong> dos polinomios <strong>en</strong> s.<br />

d<br />

− jωT<br />

lim ∠G(<br />

jω)<br />

e = −T<br />

(165)<br />

ω →∞ dω<br />

Así, <strong>de</strong> esta última ecuación, se pue<strong>de</strong> calcular <strong>el</strong> valor <strong>de</strong>l retardo <strong>de</strong><br />

transporte T.<br />

7.6- Apreciaciones Útiles Sobre La Determinación Experim<strong>en</strong>tal De<br />

La Función De Transfer<strong>en</strong>cia<br />

1. G<strong>en</strong>eralm<strong>en</strong>te es más fácil efectuar medidas exactas <strong>de</strong> amplitud que<br />

medidas exactas <strong>de</strong> <strong>de</strong>splazami<strong>en</strong>to <strong>de</strong> fase. Las mediciones <strong>de</strong><br />

<strong>de</strong>splazami<strong>en</strong>to <strong>de</strong> fase pued<strong>en</strong> introducir errores causados por<br />

instrum<strong>en</strong>tación o por interpretación <strong>de</strong> los registros experim<strong>en</strong>tales.<br />

2. La <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia <strong>de</strong>l equipo <strong>de</strong> medición utilizado para la<br />

<strong>de</strong>terminación <strong>de</strong> la salida <strong>de</strong>l sistema, <strong>de</strong>be t<strong>en</strong>er curvas <strong>de</strong> amplitud<br />

<strong>en</strong> función <strong>de</strong> la frecu<strong>en</strong>cia prácticam<strong>en</strong>te planas. A<strong>de</strong>más <strong>el</strong> ángulo<br />

<strong>de</strong> fase <strong>de</strong>be ser casi proporcional a la frecu<strong>en</strong>cia.<br />

3. Los <strong>sistemas</strong> físicos ti<strong>en</strong><strong>en</strong> algún tipo <strong>de</strong> alinealida<strong>de</strong>s. Por tanto, es<br />

necesario consi<strong>de</strong>rar cuidadosam<strong>en</strong>te la amplitud <strong>de</strong> las señales<br />

s<strong>en</strong>oidales <strong>de</strong> <strong>en</strong>trada. Si la amplitud <strong>de</strong> la señal <strong>de</strong> <strong>en</strong>trada es<br />

excesivam<strong>en</strong>te gran<strong>de</strong>, <strong>el</strong> sistema ti<strong>en</strong><strong>de</strong> a saturarse y la prueba <strong>de</strong><br />

<strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia da resultados inexactos. Por otro lado, una<br />

señal pequeña pue<strong>de</strong> producir errores <strong>de</strong>bido a zona muerta. Por<br />

tanto, hay que realizar una cuidadosa <strong>el</strong>ección <strong>de</strong> la amplitud <strong>de</strong> la<br />

señal s<strong>en</strong>oidal <strong>de</strong> <strong>en</strong>trada. Es necesario muestrear la forma <strong>de</strong> onda<br />

<strong>de</strong> la salida <strong>de</strong>l sistema para asegurarse que su forma es s<strong>en</strong>oidal y<br />

que <strong>el</strong> sistema está funcionando durante <strong>el</strong> periodo <strong>de</strong> prueba <strong>en</strong> la<br />

región lineal (la forma <strong>de</strong> onda a la salida <strong>de</strong>l sistema no es s<strong>en</strong>oidal si<br />

éste está funcionando <strong>en</strong> una región no lineal).<br />

4. Si <strong>el</strong> sistema <strong>en</strong> estudio está funcionando continuam<strong>en</strong>te por días y<br />

semanas, no hace falta <strong>de</strong>t<strong>en</strong>er su funcionami<strong>en</strong>to normal para las<br />

- 89 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

pruebas <strong>de</strong> <strong>respuesta</strong> <strong>de</strong> frecu<strong>en</strong>cia. Se pue<strong>de</strong> superponer la señal <strong>de</strong><br />

prueba s<strong>en</strong>oidal a las <strong>en</strong>tradas normales. Entonces, para <strong>sistemas</strong><br />

<strong>lineales</strong>, se superpondrá a la salida normal, la <strong>respuesta</strong><br />

correspondi<strong>en</strong>te a la señal <strong>de</strong> prueba. Para la <strong>de</strong>terminación <strong>de</strong> la<br />

función <strong>de</strong> transfer<strong>en</strong>cia con <strong>el</strong> sistema <strong>en</strong> funciona mi<strong>en</strong>to normal,<br />

también se usan frecu<strong>en</strong>tem<strong>en</strong>te señales estocásticas (señales <strong>de</strong><br />

ruido blanco). Utilizando funciones corr<strong>el</strong>ación, se pue<strong>de</strong> <strong>de</strong>terminar la<br />

función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong>l sistema sin interrumpir su funcionami<strong>en</strong>to<br />

normal.<br />

Ejemplo 7: Un sistema controlado, consiste <strong>de</strong> un posicionador hidráulico<br />

<strong>de</strong> aceite que empuja una carga la cual pue<strong>de</strong> ser asimilada como una<br />

unidad masa resorte amortiguador. La posición <strong>de</strong> la masa es controlada<br />

por la posición <strong>de</strong> la palanca <strong>de</strong> <strong>en</strong>trada <strong>de</strong>l posicionador. Es <strong>de</strong>seable<br />

establecer experim<strong>en</strong>talm<strong>en</strong>te la r<strong>el</strong>ación dinámica <strong>en</strong>tre <strong>el</strong> metimi<strong>en</strong>to <strong>de</strong><br />

<strong>en</strong>trada <strong>de</strong> la palanca x i<br />

y la posición <strong>de</strong> la masa x . Para realizar esto, es<br />

necesario un g<strong>en</strong>erador <strong>de</strong> movimi<strong>en</strong>to s<strong>en</strong>oidal, para aplicar a la palanca<br />

<strong>de</strong> <strong>en</strong>trada la oscilación s<strong>en</strong>oidal <strong>de</strong> amplitud constante y frecu<strong>en</strong>cia<br />

variable. Dicha señal <strong>de</strong> <strong>en</strong>trada será 0 .05s<strong>en</strong>ωt<br />

don<strong>de</strong> 0.05 es la amplitud a<br />

mant<strong>en</strong>er constante y ω la frecu<strong>en</strong>cia a variar.<br />

Para registrar <strong>el</strong> movimi<strong>en</strong>to <strong>de</strong> la masa, se ajusta a <strong>el</strong>la un transductor <strong>de</strong><br />

<strong>de</strong>splazami<strong>en</strong>to. La <strong>respuesta</strong> s<strong>en</strong>oidal <strong>de</strong> la masa x = B s<strong>en</strong>( ω t + φ)<br />

es<br />

medida para <strong>el</strong> rango <strong>de</strong> frecu<strong>en</strong>cia que se varió la <strong>en</strong>trada. El resultado <strong>de</strong><br />

mediciones efectuadas <strong>de</strong> B y φ se muestra <strong>en</strong> la tabla 4.<br />

Frecu<strong>en</strong>cia<br />

<strong>de</strong> <strong>en</strong>trada<br />

ω<br />

[rad/seg]<br />

Amplitud<br />

<strong>de</strong> la<br />

<strong>respuesta</strong><br />

B [pulg]<br />

Fase <strong>de</strong> la<br />

<strong>respuesta</strong><br />

φ [º]<br />

0 0,1 0,2 0,4 0,6 0,8 1 2 4 10 20 40 80 100<br />

5 5 4,72 4,25 3,35 2,8 2,5 1,12 0,35<br />

0,0<br />

8<br />

0,02 0,03 - -<br />

0 12,5 25 46 60 72 80 109 139 182 213 237 255 260<br />

Tabla 4.<br />

Con estos datos se requiere obt<strong>en</strong>er la función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong>l<br />

sistema ( s)<br />

X ( s)<br />

X<br />

i<br />

Solución: El método <strong>de</strong> graficar <strong>en</strong> Bo<strong>de</strong> es <strong>el</strong> más fácil <strong>de</strong> interpretar. Se<br />

divi<strong>de</strong> la amplitud <strong>de</strong> salida sobre la amplitud <strong>de</strong> <strong>en</strong>trada para todas las<br />

frecu<strong>en</strong>cias y se lo expresa luego <strong>en</strong> dB.<br />

La gráfica resultante es la mostrada a continuación.<br />

- 90 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

60<br />

40<br />

0 dB/déc.<br />

-20 dB/déc.<br />

Amplitud, dB<br />

20<br />

0<br />

-20<br />

-40 dB/déc.<br />

-60 dB/déc.<br />

-40<br />

-60<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

90<br />

0<br />

Fase, Grados<br />

-90<br />

-180<br />

-270<br />

10 -1 10 0 10 1 10 2<br />

De estas gráficas se pue<strong>de</strong> <strong>de</strong>ducir:<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 64.<br />

a. G dB = f (ω)<br />

es asintótica a la línea horizontal <strong>de</strong> G dB = 40 dB ,<br />

cuando la frecu<strong>en</strong>cia ω ti<strong>en</strong><strong>de</strong> a cero. Esto implica que <strong>en</strong> la<br />

función <strong>de</strong> transfer<strong>en</strong>cia buscada no hay ceros ni polos <strong>en</strong> <strong>el</strong><br />

orig<strong>en</strong>.<br />

b. La p<strong>en</strong>di<strong>en</strong>te <strong>el</strong>e la curva G dB = f (ω)<br />

se increm<strong>en</strong>ta continuam<strong>en</strong>te<br />

<strong>en</strong> <strong>el</strong> s<strong>en</strong>tido horario. 0 sea que no hay puntos <strong>de</strong> inflexión (no hay<br />

ceros).<br />

c. La curva G dB = f (ω)<br />

es asintótica a una línea <strong>de</strong> p<strong>en</strong>di<strong>en</strong>te -60<br />

dB/déc cuando la frecu<strong>en</strong>cia ti<strong>en</strong><strong>de</strong> a infinito; por lo tanto <strong>el</strong><br />

- 91 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

sistema es <strong>de</strong> tercer ord<strong>en</strong>.<br />

d. La curva φ = f (ω)<br />

es asintótica a φ = 0º<br />

para bajas frecu<strong>en</strong>cias (lo<br />

qua implica que K es positivo) y pareciera ser a -270° para altas<br />

frecu<strong>en</strong>cias.<br />

Por lo expuesto <strong>en</strong> b se observa que no hay ceros <strong>en</strong> la función <strong>de</strong><br />

transfer<strong>en</strong>cia. La aus<strong>en</strong>cia <strong>de</strong> puntos máximos y mínimos <strong>en</strong> la curva<br />

φ = f (ω) apoya esta conclusión.<br />

En conjunción con b, la c indica la pres<strong>en</strong>cia <strong>en</strong> <strong>el</strong> d<strong>en</strong>ominador <strong>de</strong> un<br />

polinomio <strong>de</strong> tercer ord<strong>en</strong>.<br />

Por lo tanto se pue<strong>de</strong> concluir que la función <strong>de</strong> transfer<strong>en</strong>cia es <strong>de</strong> la<br />

forma:<br />

X ( s)<br />

X ( s)<br />

i<br />

=<br />

T s<br />

K<br />

3 2<br />

1<br />

+ T2s<br />

+ T3s<br />

+<br />

1<br />

La falta <strong>de</strong> puntos <strong>de</strong> resonancia <strong>el</strong>imina <strong>el</strong> caso d. Un posterior<br />

<strong>en</strong>t<strong>en</strong>dimi<strong>en</strong>to <strong>de</strong>l sistema físico bajo prueba indica la probabilidad <strong>de</strong>l caso<br />

c como <strong>el</strong> mejor mo<strong>de</strong>lo. Sobre esta base habrá cuatro aproximaciones<br />

asintóticas <strong>de</strong> la curva, según se muestra <strong>en</strong> la figura 65, <strong>de</strong> la cual se<br />

pued<strong>en</strong> <strong>de</strong>terminar las tres frecu<strong>en</strong>cias esquinas:<br />

T<br />

T<br />

1<br />

2<br />

=<br />

=<br />

1<br />

0.5<br />

1<br />

4<br />

= 2<br />

= 0.25<br />

T<br />

3<br />

=<br />

1<br />

20<br />

= 0.05<br />

A<strong>de</strong>más: K [ dB]<br />

= 20log x = 40 ⇒ x = 100<br />

Lo cual era ya evid<strong>en</strong>te pues los coci<strong>en</strong>tes <strong>de</strong> las amplitu<strong>de</strong>s <strong>de</strong> salida y<br />

<strong>en</strong>trada para ω = 0 es 5 0.05 = 100 .<br />

Finalm<strong>en</strong>te la función <strong>de</strong> transfer<strong>en</strong>cia obt<strong>en</strong>ida experim<strong>en</strong>talm<strong>en</strong>te será:<br />

X ( s)<br />

100<br />

=<br />

X<br />

i<br />

( s)<br />

05<br />

( 1+<br />

2s)( 1+<br />

0.25s)( 1+<br />

0. s)<br />

- 92 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

La evaluación más precisa <strong>de</strong> las frecu<strong>en</strong>cias esquinas, pue<strong>de</strong> ser evaluada<br />

mejor conoci<strong>en</strong>do los parámetros <strong>de</strong>l sistema, dado que las constantes <strong>de</strong><br />

tiempo <strong>de</strong>l sistema son función <strong>de</strong> <strong>el</strong>los.<br />

Ejemplo 8: La <strong>respuesta</strong> frecu<strong>en</strong>cial experim<strong>en</strong>tal <strong>de</strong> un sistema es<br />

mostrada <strong>en</strong> la tabla 5.<br />

Frecu<strong>en</strong>cia<br />

<strong>de</strong> <strong>en</strong>trada<br />

ω<br />

[rad/seg]<br />

Amplitud<br />

<strong>de</strong> la<br />

<strong>respuesta</strong><br />

B [pulg]<br />

Fase <strong>de</strong> la<br />

<strong>respuesta</strong><br />

φ [º]<br />

0,1 0,2 0,4 0,6 0,8 1 2 3 5 7 10 20 40<br />

1 0.95 0.75 0.6 0.47 0.4 0.25 0.225 0.2 0.19 0.16 0.12 0.05<br />

11 22 37 42 46 47 45 43 44 48 55 68 78<br />

Tabla 5.<br />

La amplitud <strong>de</strong> la <strong>en</strong>trada ( H ) es constante e igual a 1. Establecer la<br />

función <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong>l sistema.<br />

Con los valores <strong>de</strong> la tabla se dibuja <strong>el</strong> gráfico <strong>de</strong> Bo<strong>de</strong> <strong>de</strong> at<strong>en</strong>uación<br />

B H = f (ω) y <strong>de</strong> fase φ = f (ω)<br />

.<br />

- 93 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

5<br />

0<br />

0 dB/déc<br />

-5<br />

-20 dB/déc<br />

Amplitud, dB<br />

-10<br />

-15<br />

-20<br />

0 dB/déc<br />

-20 dB/déc<br />

-25<br />

-30<br />

T1 = 2.5 T2 = 0.5 T3 = 0.1<br />

-35<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

0<br />

-10<br />

-20<br />

Puntos <strong>de</strong> Inflexión<br />

-30<br />

Fase, Grados<br />

-40<br />

-50<br />

-60<br />

-70<br />

-80<br />

-90<br />

-100<br />

10 -1 10 0 10 1 10 2<br />

Frecu<strong>en</strong>cia, rad/s<br />

Figura 65.<br />

La curva <strong>de</strong> at<strong>en</strong>uación ti<strong>en</strong>e una p<strong>en</strong>di<strong>en</strong>te original <strong>de</strong> cero dB/década<br />

seguida por una <strong>de</strong> -20 dB/década, y luego por una reducción <strong>de</strong> la<br />

p<strong>en</strong>di<strong>en</strong>te a 0 dB/década volvi<strong>en</strong>do por último a -20 dB/riécada.<br />

Inmediatam<strong>en</strong>te se <strong>de</strong>duce que la primer p<strong>en</strong>di<strong>en</strong>te <strong>de</strong>spués <strong>de</strong> la inicial <strong>de</strong><br />

0 dB/década se <strong>de</strong>berá a un factor 1 ( 1+ T 1<br />

s)<br />

, la que sigue se <strong>de</strong>berá a un<br />

1 1+ T 3<br />

s . Por lo tanto la función <strong>de</strong><br />

factor ( T 2<br />

s)<br />

1+ y la última a un factor ( )<br />

transfer<strong>en</strong>cia-será <strong>de</strong> la forma:<br />

- 94 -


Respuesta <strong>en</strong> Frecu<strong>en</strong>cia <strong>de</strong> Sistemas Lineales<br />

Control I<br />

C(<br />

s)<br />

R(<br />

s)<br />

=<br />

K( 1+<br />

T2s)<br />

( 1+<br />

T s)( 1+<br />

T s)<br />

1<br />

3<br />

Como se ve <strong>en</strong> la gráfica la curva φ = f (ω)<br />

es asintótica a 0 o para frecu<strong>en</strong>cias<br />

bajas y a -90° para frecu<strong>en</strong>cias t<strong>en</strong>di<strong>en</strong>do a infinito. Trazando <strong>en</strong> <strong>el</strong> diagrama<br />

<strong>de</strong> at<strong>en</strong>uación las asíntotas <strong>de</strong>l mejor modo posible se pued<strong>en</strong><br />

<strong>de</strong>terminar:<br />

T<br />

T<br />

T<br />

1<br />

2<br />

3<br />

= 2.5<br />

= 0.5<br />

= 0.1<br />

A<strong>de</strong>más K = 1, ya que 20 log K = 0<br />

Por consigui<strong>en</strong>te:<br />

C(<br />

s)<br />

R(<br />

s)<br />

=<br />

11 ( + 0.5s)<br />

( 1+<br />

2.5s)( 1+<br />

0.1s)<br />

Los puntos <strong>de</strong> inflexión <strong>en</strong> la curva <strong>de</strong> φ = f (ω)<br />

ocurr<strong>en</strong> <strong>en</strong> la región <strong>de</strong> cada<br />

frecu<strong>en</strong>cia esquina, lo cual ayuda a confirmar las frecu<strong>en</strong>cias esquinas<br />

<strong>de</strong>terminadas.<br />

- 95 -

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!