Significado del "polo" de una función de transferencia con retardo de tiempo

Aprendí en teoría básica de control que podemos determinar la estabilidad de un sistema LTI por los signos de la parte real de los polos de su función de transferencia H ( s ) = norte ( s ) D ( s ) . Para un sistema racional, donde ambos norte ( s ) y D ( s ) son polinomios es s, es fácilmente comprensible porque una parte real positiva del polo provocará un aumento exponencial en el dominio del tiempo, siendo así inestable. Mi pregunta es, si el sistema tiene un retraso de T segundos en alguna parte, entonces H ( s ) contendrá un factor Exp ( s T ) , a veces en el denominador. En este caso, ¿por qué todavía podemos determinar la estabilidad mirando la parte real de los "polos" de H ( s ) (puntos en el plano s donde D ( s ) = 0 )?

Un ejemplo:

norte ( s ) = s

D ( s ) = a s 2 + b s + C + d Exp ( s T )

donde a, b, c, d y T son constantes.

En este caso, no podemos hacer expansión en fracciones parciales como lo hacemos para funciones de transferencia racionales. Entonces, ¿cómo vincular el comportamiento en el dominio del tiempo con las raíces de D ( s ) ?

Respuestas (2)

Creo que la propiedad a la que te refieres es

L { F ( t τ ) } = mi s τ L { F ( t ) }

Esto significa que su función de transferencia realmente no sigue esa regla y realmente no tiene un retraso de tiempo:

H ( s ) = s a s 2 + b s + C + d mi s T

Resolver la transformada inversa de Laplace para esta función probablemente implicará métodos numéricos. Lo que es posible, es tener algo como esto:

H ( s ) = s + d mi s T a s 2 + b s + C

Porque se puede separar en:

H ( s ) = s a s 2 + b s + C + mi s T d a s 2 + b s + C

Puedes ver que el segundo término también se puede separar en fracciones parciales, como el primer término. Los polos del denominador dirán algo acerca de la estabilidad de la función de transferencia.

Estabilidad en general

La fórmula de la transformada inversa de Laplace es

F ( t ) = 1 2 π j límite T σ j T σ + j T F ( s ) mi s t d s

dónde σ es elegido para ser mayor que todas las singularidades de F ( s ) en el plano complejo (en nuestro caso tal que incluye todos los polos).

Esta integral se resuelve de manera equivalente mediante el teorema del residuo de Cauchy :

F ( t ) = L 1 { F ( s ) } = a yo yo   pag o yo mi s   o F   F ( s ) R mi s [ F ( s ) mi s t ]
Recuerde, esto es general , es decir. ¡siempre funciona para una función de transferencia con cualquier singularidad!

Entonces no importa si los polos provienen o no de un polinomio o de una función trascendental. Siempre que las singularidades estén todas en el semiplano izquierdo, su residuo siempre contendrá un exponencial que decae a 0 en lugar de infinito. Cualquier singularidad en el RHP siempre conducirá a un exponencial que explota.


Apéndice

Se puede notar que mi s T se puede ampliar en su serie Taylor:

mi z = norte = 0 + z norte norte !

Entonces esto significa que su función de transferencia de ejemplo se puede escribir como

H ( s ) = s a s 2 + b s + C + d ( norte = 0 + ( s T ) norte norte ! ) = s norte = 0 + A norte s norte

Dónde A 0 = C + d , A 1 = b d T , A 2 = a + d T 2 2 , A norte = d ( T ) norte norte ! , norte > 2 .

Entonces, este tipo de función de transferencia tiene "polos" por todas partes. Esto es para ilustrar que definitivamente no puede manejar esto como una función de transferencia normal de segundo orden.


Apéndice

El residuo de un polo (simple) es

R mi s s = a [ H ( s ) ] = límite s a ( s a ) F ( s )

Para polos con multiplicidad norte :

R mi s s = a [ H ( s ) ] = 1 ( norte 1 ) ! límite s a d norte 1 d s norte 1 ( ( s a ) norte H ( s ) )

Creo que la expansión de Taylor en su Apéndice es un buen punto, ya que revela que habrá muchos polos sobre el espacio como sugirió. En mi caso, la exponencial está en efecto en el denominador debido a que el retraso se produce en un ciclo de retroalimentación. Generalmente se requieren métodos numéricos, estoy de acuerdo, pero me gustaría saber por qué conceptualmente podemos determinar la estabilidad mirando las raíces del denominador.
@GeorgeC Edité mi respuesta a lo que creo que quisiste decir.
George, las raíces del denominador son los polos de la función de transferencia. Y ninguno de estos polos no debe aparecer en la mitad derecha del plano s (RHP). Por lo tanto, observamos las ubicaciones de todos los polos mientras variamos la ganancia. Usando este método, podemos ver a qué ganancia los polos se mueven hacia el RHP (inestabilidad). Esto indica dónde está el límite de estabilidad.

La expresión para un retardo de tiempo (tiempo muerto) en el dominio de la frecuencia exp(-sT) es una función trascendental. No se puede modelar con elementos agrupados y, por lo tanto, no se puede combinar directamente con una función de transferencia que consta de polinomos en "s". Sin embargo, para este propósito, se puede aproximar mediante la siguiente expresión (aproximación de Pade de segundo orden):

H d ( s ) = 1 s T 2 + s 2 T 2 12 1 + s T 2 + s 2 T 2 12

Esta es una función racional (allpass) y aparece como un factor en la ganancia de bucle y aparece, por supuesto, en la función de transferencia de bucle cerrado. Aquí se puede tratar como cualquier otro bloque en el circuito cerrado.

Sí, he aprendido que el Pade es una técnica común para aproximar retrasos. Sin embargo, incluso si no hacemos Pade, todavía determinamos la estabilidad encontrando los "polos" de la función de transferencia original. Como ejemplo, la función de transferencia de bucle abierto es H(s)=KG(s), donde G(s)=exp(-sT)/(s(s+1)). El libro de texto simplemente dibuja el lugar geométrico de las raíces del polo en bucle cerrado frente a la ganancia K (sin aproximación de Pade). ¿Por qué podemos hacer esto?
Está bien, sé lo que quieres decir. En mi respuesta he dado la aproximación de Pade de segundo orden. Por supuesto, el orden tercero o incluso el orden 15 es más exacto. Un simulador puede manejar órdenes más altas que n = 2, pero no sé exactamente qué hace el simulador. Supongo que el orden de aproximación está limitado a cualquier valor grande y, finalmente, el simulador conecta todos los polos calculados y, por lo tanto, es capaz de producir un lugar de raíces.