¿Cómo estimar el tiempo de asentamiento de un sistema sobreamortiguado?

He estado tratando de averiguar cómo estimar el tiempo de establecimiento de un sistema de segundo orden en respuesta a un paso de entrada de magnitud 5. La función de transferencia del sistema es

GRAMO ( s ) = 1 ( s + 2 ) ( s + 4 )
y ya he determinado el tiempo de respuesta con la entrada de paso R(s):
C ( s ) = R ( s ) GRAMO ( s ) C ( t ) = 5 8 + 5 8 mi 4 t 5 4 mi 2 t
Ahora necesito estimar el 2 % del tiempo de establecimiento de la respuesta utilizando esta información, pero no estoy seguro de cómo hacerlo. Sé que el sistema está sobreamortiguado como ζ>1, por lo que no puedo usar la ecuación de tiempo de establecimiento normal
T s = 4 ζ ω norte
Investigué esta publicación: ( tiempo de establecimiento de sistemas con amortiguación excesiva y crítica ) pero las respuestas solo explican formas largas de obtener un resultado preciso. Ya he usado MATLAB para obtener un resultado exacto de 2,3 segundos, pero necesito poder estimarlo sin MATLAB.

Estaba pensando que podría probar prueba y error con diferentes valores de t hasta que c (t) esté dentro del 2% del valor de estado estable (que es 0.625) pero aunque esto funcionaría, dudo que sea la forma correcta de hacerlo, también lo es ¿Alguien puede ayudarme con un mejor método?

Respuestas (3)

Para los sistemas con polos reales del semiplano izquierdo, generalmente puede estimarlo considerando solo el polo dominante (el polo con la frecuencia más baja). En tu caso esto seria pag d = 2 . El resultado se vuelve más preciso a medida que el polo no dominante ( pag norte d ) se aleja más del polo dominante.

Al considerar solo el polo dominante, obtienes una ecuación bastante simple:

5 4 mi 2 t = 0.02 5 8 t = 1 2 en ( 0.02 5 8 4 5 ) 2.30258509299 s

La idea es que el polo no dominante en pag norte d = 4 conduce a un término mi 4 t que se amortiguará tan rápidamente que no afectará el tiempo total de asentamiento. La ventaja es la simplicidad de la ecuación y el hecho de que en realidad es bastante común tener un polo muy dominante y polos no dominantes lejanos en los circuitos electrónicos.

En su caso específico, es posible calcular analíticamente el tiempo de asentamiento. El tiempo que tardan los términos dependientes del tiempo en amortiguarse al 2 % del valor final se puede calcular usando (similar a la respuesta de Andy, pero usando el valor absoluto):

| mi 4 t 2 mi 2 t | = 0.02 ( y = mi 2 t ) y 2 2 y = ± 0.02 ( Hay 4 soluciones distintas, pero solo tomo la relevante ) y = mi 2 t = 1 7 5 2 t = 1 2 en ( 1 7 5 2 ) 2.30006613189 s

Entonces un factor de 2 para pag norte d / pag d conduce a un error de alrededor del 0,1 % en el tiempo de establecimiento calculado cuando se utiliza la aproximación del polo dominante. Si esto es suficiente o no, se lo dejo a usted.

Eso es perfecto gracias!

Sí, su cálculo de Laplace inverso es correcto.

El valor de estado estable final será 5/8: este es el valor de CC después de un largo período de tiempo. Entonces, realmente está buscando que el resto de la ecuación caiga en magnitud al 2% de 5/8: -

5 8 mi 4 t 5 4 mi 2 t = 5 8 0.02

= 8 8 mi 4 t 8 4 mi 2 t = 8 8 0.02

= mi 4 t 2 mi 2 t = 0.02

¿Eso ayuda?

Estaba bastante cerca entonces ya que estaba haciendo algo similar a esto, pero estaba usando 0.98 en lugar de 0.02. ¿Por qué usas el 2% y no el 98%? Seguramente si queremos un 2% menos que el valor final, ¿sería 0,98 x 5/8?
@Sam, tu enfoque y el mío son igualmente válidos, es solo que pensé de una manera y tú pensaste de la otra.
Tendríamos que usar 0.02 en lugar de + 0.02 para encontrar el tiempo de asentamiento correcto (o invertir los coeficientes de las exponenciales, o usar el 0.98 @Sam propuesto e incluir el término DC). Probablemente un descuido, pero ¿tal vez la fuente de la confusión?

Bueno, resolvamos esto en un caso más general. Tenemos la siguiente función de transferencia (suponiendo un valor real positivo para ϵ ):

(1) H ( s ) := Y ( s ) X ( s ) = 1 ( s + ϵ ) ( s + 2 ϵ )

Cuando miramos la respuesta escalonada que estamos usando X ( s ) = L t [ θ ( t ) ] ( s ) = 1 s , por lo que la salida viene dada por:

(2) Y ( s ) = 1 s 1 ( s + ϵ ) ( s + 2 ϵ )

Usando la transformada inversa de Laplace, encontramos:

(3) y ( t ) = L s 1 [ 1 s 1 ( s + ϵ ) ( s + 2 ϵ ) ] ( t ) = Exp ( 2 ϵ t ) ( Exp ( ϵ t ) 1 ) 2 2 ϵ 2

No es difícil demostrar que cuando t (suponiendo un valor real positivo para ϵ ), obtenemos:

(4) límite t y ( t ) = 1 2 ϵ 2

Ahora, para el tiempo de establecimiento, queremos encontrar el tiempo t cuando norte % del valor final se alcanza:

(5) y ( t norte ) = norte % 100 1 2 ϵ 2   t norte =

Resolviendo eso da:

(6) t norte = 1 ϵ en ( 10 10 norte % )