Respuesta a Excitaciones Exponenciales
Respuesta a Excitaciones Exponenciales
Exponenciales
Dentro de la categoría de excitaciones exponenciales involucramos las dadas por:
= ∙ ∙ cos " +
(1)
∙ ∙
= ∙ ∙ + = ∙ + ∙ = ′ + ′′ (2)
donde
∙ =Ė
∙ =Ễ
son números complejos (fasores) que se caracterizan por la amplitud y por la fase al origen de la excitación.
Las constantes de atenuación de las exponenciales temporales son: = + !" y # $ = − !" que, como
quedará evidenciado más adelante, son los polos de la función excitación (modo forzado). De esta manera, la
función queda expresada como suma de dos funciones:
Ė& Ễ&
= ∙ = ∙
' ∙ )∙
'( (3)
′ (a) y ′′ (b)
Ė& Ễ&
= *
+ **
= ∙+ + ∙+ )
' '(
(4)
que si la escribiéramos con común denominador, quedaría en evidencia que $ son los polos de la excitación.
y #
Conocidos y # $ , junto a los fasores Ė y Ễ, queda definida la función excitación propuesta de (1), pero puede
concluirse que toda la información necesaria está contenida en ′ dada por (3)(a), donde figuran Ė y , pues en
(3)(b)
tenemos iguales parámetros pero en su forma conjugada. Luego se verá que, calculada la respuesta en régimen
= ℜ- * . y si es /0*
permanente a ′ , la correspondiente a ′′ será la conjugada de la primera y, por lo tanto, puede escribirse sólo
por inspección. Otra manera de ver lo mismo sería considerando que es la respuesta
*
permanente (forzada) a , también será:
1
Ahora podemos detenernos a analizar (1) para ver que encierra toda una gama de funciones, considerando que "
puede ser nulo o tomar cualquier valor positivo (frecuencias " negativas no tienen significado físico, puede haber
valores negativos de " , pero nunca frecuencias negativas), y que puede asumir cualquier valor real. Asi con;
a) 78 ≠ :
2
b) 78 = : → ?8@ AB = ?CB. = EA → B = EA ∙ F8 ∙ BG8 C
b2) = 0, Escalón
Naturalmente debe entenderse que es la inversa de la constante de tiempo de las exponenciales y caracteriza la
H
velocidad de extinción o crecimiento de las mismas, y " es la frecuencia de las oscilaciones " = siendo J el
I
,
período de ellas.
Habiendo relacionado (1) con los polos que la identifican, y en virtud de la posibilidad de expresarla como
superposición de las (3), que son “funciones exponenciales” es que al tema lo denominamos de búsqueda de
“respuesta a excitaciones exponenciales”, aunque debemos tener presente que cada una de las (3) no son funciones
reales del tiempo.
La linealidad de los sistemas que estudiamos permite que se busque la respuesta por separado a ′ y a ′′ , y
luego superponerlas para encontrar la respuesta a . Sin embargo, para mayor claridad aplicaremos ambas
simultáneamente.
Para ello supondremos que el sistema viene dado por su función Transferencia H(s) que en términos generales es:
3
donde R = 0 es el polinomio característico y sus raíces ST , ST , . . , STV son los polos del sistema (modos
naturales) pudiendo escribirse:
K =
M '
' 0vQ ' 0vw .. ' 0vO (7)
Impondremos la condición de que los polos ST tengan todos parte real negativa para evitar que arrojen respuesta
libre divergente o contante: con esto nos aseguramos que la respuesta libre se extinga en el tiempo, constituyéndose
la “respuesta forzada” en el “régimen permanente”.
$ = − !"
resulta ser un número real, y además, resuelta en
valores conjugados de resultan números complejos conjugados: si = + !" y #
.
b ,
con ] = ] ," = ^/_ ` ac ,
d = ^/`-K (10)
V
e* e** gh
t = + +u
s− $
s− # − STh
hk (12)
e′ = i s − ∙ R s |+k' = i lĖm +
Ễ& ∙ + '
+ ')(
n∙H s p = ∙ Ėm ∙ H s = (13)
+k'
e′′ = iqs − #
$r∙R s s = i l &+ ' + Ễm n ∙ H s p $r
Ė q+ '(r )
+k')
(
= ∙ Ễm ∙ Hqs = # (14)
+k')
(
e′ = ∙ Ėm ∙ |H |∙ [
y $ rs ∙
e′′ = ∙ Ễm ∙ sHq # [ (15)
4
de igual manera resultarán los
Que serán reales si el polo STh es real, y serán complejos conjugados en polos complejos conjugados.
V
1 Ėm ∙ [
Ễm ∙ [
gh
t = ∙ |H |… + †+u
2 s− $
s− # − STh
hk
(17)
Evidentemente la parte forzada puede desglosarse en dos términos, comparables con (4);
1 Ėm 1 Ėm ∙ |H |∙ [
1 t•
t *y = ∙ ∙H = E* s ∙ H = ∙ = ∙
2 s− 2 s− 2 s−
(18)
1 Ễm $ rs ∙
1 Ễ ∙ sHq # [
1 t‚ƒ
t **y = ∙ $ r = E ** s ∙ Hq #
∙ Hq # $r = ∙ m = ∙
$
2 s− # 2 s− # $ 2 s− #$
(19)
Esto indica que los respectivos fasores respuesta y excitación se relacionan por:
t•
t• = Ėm ∙ H → =H
Ėm
tƒ‚
tƒ‚ = Ễm ∙ Hq #
$r → $r
= Hq #
Ễm
(20) y (21)
V
1 1
/ = /y* z{P|P + /y** z{P|P + /}h~z = ∙ t• ∙ '
+ ∙ t‚ƒ ∙ ')
(
+ u gh ∙ 0vx
2 2
hk (22)
considerando que la respuesta libre se extingue con el tiempo, la respuesta en régimen permanente coincide con la
forzada;
luego, sustituyendo en (22) las (20) y (21), y teniendo en cuenta las (18) y (19), salen respuestas también exponenciales;
5
/0 = ∙ Ėm ∙ H ∙ ' $r∙
+ ∙ Ễm ∙ Hq #
)
'( (24)
1
/0 = ∙ Ė ∙ |H |∙ ∙ [
∙ + ∙ [
∙
2 m
[
+ [
= Ėm ∙ |H |∙ … †
2
= Ėm ∙ |H |∙ ∙ cos " + +]
(25)
y por (20);
/0 =t ∙ ∙ cos " + z
(27)
Observando la (27) y comparándola con (1), en cualquiera de sus formas (a1, a2, a3, b1, b2 y b3) vemos que la
respuesta forzada es del mismo tipo de la excitación; es la misma función armónica, de igual frecuencia " , y su
amplitud está modulada por la misma exponencial . Las únicas diferencias están entre las amplitudes
= |H | ) y en la fase en el origen
L
(verificándose que z (siendo, z = + ], o sea el desfasaje entre /0
y es z − = ] = ^/`-K . ).
excitación es una de las seis descriptas, identificamos = + !" , y por (27) escribimos directamente la
Como conclusión, conocida la función de transferencia entre una respuesta y una excitación “H ”, si la
Electrotecnia de Corriente Alterna es sólo un caso particular del aquí presentado, donde = !" =0 .
De lo visto con el tratamiento fasorial surge que el método operativo simbólico o polar empleado en la
Interpretando las impedancias operacionales de elementos simples (R, L y C) como relaciones de transferencia
entre corriente y tensión en bornes, resultarán:
a1 a2 a3 b1 b2 b3
6
Dados en régimen permanente senoidal, ‹ = ℜ X• yŒ = ℜ •• , la impedancia en régimen
c• •
permanente senoidal es el complejo Ž0 = ••
≠ h
. O sea para régimen permanente es el valor de la impedancia de
entrada Z(s) resuelta en = !" .
En la columna a2 (excitación armónica pura) se nota que aparecen los valores !" ‡ y
‘
de tratamiento usual en
Corriente Alterna. Si se tratara de corriente continua aparecen los valores de columna b2 (recordar que se dedujeron
conclusiones para régimen permanente).
Si las variables respuesta y excitación son alejados entre sí y pertenecen a un circuito con cualquier tipo
interconexión de elementos R, L y C, se puede llegar a interrelación entre respuesta y excitación (o sea, la función
columna correspondiente en el cuadro anterior, con lo cual resulta directamente el módulo |H | y la fase ].
de transferencia) por agrupamientos serie, paralelo, transfiguraciones, Thevenin, etc. de las impedancias dadas en la
/0 = t• ∙ '
= Ėm ∙ H ∙ '
= Ėm ∙ |H |∙ [ ' (29)
t• = Ėm ∙ H = Ėm ∙ X + !Y = W + !R = st• s ’
estando contenida en t• toda la información necesaria para definir la respuesta temporal en régimen permanente
(módulo y fase).
H
Frecuencia fundamental: " =
Período: a
P
7
˜
e
– = + u“eV ∙ _™ ”" + —V ∙ ” ”" •
2
Vk
(30)
donde,
œ P
2
e = š – ›
^
œ
œ P
2
eV = š – _™ ”" ›
^
œ
œ P
2
—V = š – ” ”" ›
^
œ (31)
Si se conoce la función de trasferencia K , aplicando como excitación tal f(t), la respuesta resulta / ′ ;
˜
e
/0 = K 0 + u“eV |K !”" |_™ ”" + ]V + —V |K !”" | ” ”" + ]V •
2
Vk (32)
donde;
]V = ^/`“K = !”" •
Las series de Fourier también pueden expresarse en su forma compleja, aplicando Euler (o escribiendo el seno y
el coseno en su forma compleja) en (30):
˜
e V
+ V V
− V
– = + u …eV ∙ + —V ∙ †
2 2 2!
Vk
Definimos:
8
eV − !—V
•V =
2
eV + !—V
• =
V
2
e
• =
2
(33)
entonces:
˜ ˜
– = • + u •V ∙ V
+• V ∙ V
= u •V ∙ V
Vk ˜
(34)
œ P
1
•V = š – V
›
^
œ
(35)
dada la f(t) según (34) y la función de Transferencia H(s), su respuesta ahora se calcula según (28) y (29) empleando
relaciones fasoriales:
/0 = u •V ∙ K = !”" ∙ V
Vk ˜ (36)
tanto la (32) como la (36) convergen a una función /0 que será periódica, de igual período que la f(t), esto es:
9
Respuesta en régimen permanente a excitaciones periódicas generales
aplicando T. de Laplace
La transformada de Laplace permite, a través de un número finito de operaciones, encontrar la respuesta en
régimen permanente rp(t) de forma exacta. Mientras que aplicando Fourier, se requería un número infinito de
operaciones (o bien obtener una expresión aproximada si se utilizan los primeros términos de la serie).
Sea tal excitación aplicada a través de una función de transferencia H(s) “incondicionalmente estable” (todos sus
polos tienen parte real negativa).
Sea e(t) tal excitación periódica (con un número finito de discontinuidades finitas, sin impulsos y absolutamente
integrable).
e(t) si 0 ≤ < ^
e1(t) =
0 si ≥ ^
Gráficamente:
H
Frecuencia fundamental: " =
Período: a
P
De las ecuaciones (34), (35) y (36), podemos escribir de acuerdo a la serie de Fourier en forma compleja excitación y
respuesta en régimen permanente;
= u •V ∙ V
˜ (37)
P P
1 1
•V = š ∙ V
› = š ∙ V
›
^ ^
ž ž (38)
10
˜
/0 = u •V ∙ K !”" ∙ V
Vk ˜ (39)
Recordamos que rp(t) es periódica: rp(t+k·a)= rp(t), por lo que es suficiente hallarla en el primer período, para lo
cual es indistinto excitar con e(t) o con e1(t). Sin embargo el procedimiento nos exigirá trabajar con ambas
excitaciones.
Aplicando el segundo Teorema de Traslación en T. de Laplace, se puede expresar e(t) en función de e1(t):
=u − ª^ ∙ ¡ − ª^
¥kž
¡ : _^£™” › K ^‹ Œ›
(40)
˜ ˜
1
= ℒ“ •=u ∙ ¥P'
= u ¥P'
= ∙
1− P'
¥kž ¥kž
Luego,
=
1− P'
(41)
La función se obtiene tras realizar una cantidad finita de operaciones, mientras que contiene infinitos
términos.
/ =ℒ “ K •=ℒ … K †, ∀
1− P'
(42)
/ =ℒ “ K •, 0≤ ≤^
(43)
(1er período)
/ =/ 0≤ ≤^ (44)
11
Con esto queremos decir que en el primer período, ambas tienen igual componente de régimen permanente (la
buscada) e igual componente transitoria.
/ se obtiene anti-transformando a través de un número finito de operaciones, pues H(s) tiene “n” polos y E1(s)
colabora con algunos términos más.
Comparando (42) con (43), vemos que son similares salvo por el denominador 1 − P' , por lo que sus raíces son
polos de r(t) pero no de r1(t) (no olvidar que se están estudiando las funciones en el intervalo del primer período).
a) los polos de H(s): que dan el transitorio, por tener una parte real negativa
Caso a
Caso b
˜ P
=š '
› =š '
›
ž ž
(45)
≤ gš '
›
ž
que tiene valor finito para toda “s”, con lo que deducimos que “no tiene polos”, vale la regla de
diferenciación bajo el signo integral:
P P
› «
=š “ ' •› = − š“ · · ' •›
› «
ž ž
la cual existe para todo valor finito de “s”, por lo que es regular en todo el plano “s”: es una función analítica
en todo el plano complejo, es entera.
Caso c
12
H
1− P'
= 0, esto da como resultado ¥ = !2ª = !ª" ;
P
para todo “k” entero, pues ez es una función periódica, por lo tanto me va a aportar infinitos polos en R(s).
Aquí queda al descubierto que E(s) (y por lo tanto R(s)) tienen infinitos polos, que No son polos de R1(s) (la
transformada de r1(t)).
˜ V
/ = u e¥ ¥
+ u gh · 0vx
¥k ˜ hk
(48)
Se deberían agregar a (48) los términos correspondientes a los polos de , pero como se demostró, no existen,
por lo que (48) es completa, y la podemos sub-dividir en dos respuestas:
/ = /¥ + /T (49)
con, /¥ = ∑¥̃k ˜ e¥
¥ (50)
K
e¥ = it − !ª" |'k =i − !ª" ®
¥
1− P'
'k ¥
ž
ž
el que resulta una indeterminación por los dos términos entre paréntesis. La indeterminación se puede resolver
por L’Hopital haciendo el cociente de las respectivas derivadas:
›
− !ª" 1
i› ¯ =i
› ^ P' p
1 − P' 'k ¥
› 'k ¥
Volviendo a e¥ ;
1 1
e¥ = i K P' p
= !ª" K !ª"
^ 'k ¥ ^
(52)
13
P P
1 1
•V = š ∙ V
› = š ∙ V
›
^ ^
ž ž (38)
˜ P
=š '
› =š '
›
ž ž
(45)
!”" =š V
›
ž
(53)
tanto “k” como “n” son enteros, y varían entre −∞ y +∞, por ello la (52) es:
e¥ = •¥ K !ª"
(54)
luego en (50):
˜ ˜
/¥ = u e¥ ¥
= u •¥ K !ª" ¥
= /0
¥k ˜ ¥k ˜
O sea en (49), la parte de la respuesta r(t) correspondiente a la excitación e(t), calculada en los ceros de 1 − P'
es
la respuesta en régimen permanente, pues coincide con la que nos da la serie de Fourier. Precisamente es la
respuesta que pretendemos calcular (por otro método).
/ = /0 + /T (55)
/0 =/ − /T (56)
14
Determinación de -° C
K
gh = i − STh ®
1− P'
'k0vx
(57)
Sin embargo, en (56) subsiste el inconveniente para conocer / , para obtener luego /0 .
Es aquí que nos apoyaremos en la (43), que nos da, según (44), una / =/ , pero sólo coincidentes en el primer
período.
/0 =/ − /T (58)
Pero como es sabido, la respuesta en régimen permanente es periódica, de período “a”. Entonces de (58)
calcularemos la /0 para 0 ≤ ≤ ^ y luego la extendemos como réplica al resto de los períodos en virtud de
˜
/0 = u /0 − ª^ ∙ ¡ − ª^
¥kž
Debemos de aclarar que / calculada según (43) se obtiene por cantidad finita de operaciones, por lo que
concluimos que la (58), que nos da la solución exacta para el régimen permanente en el primer período, se obtiene
por un número finito de pasos.
15