Respuesta escalonada de HPF de segundo orden

En el análisis de un filtro de paso alto RC pasivo de segundo orden, encuentro que la relación de amortiguamiento zeta = 3/2 y el factor de calidad Q = 1/3. Estos valores no dependen de R y C. Tenía la impresión de que si la relación de amortiguamiento es mayor que uno, entonces la respuesta de paso del sistema no sobrepasa su valor de estado estable. Sin embargo, cuando realizo la respuesta escalonada en MATLAB para mi sistema, encuentro que el circuito se sobrepasa ligeramente. Estoy analizando el circuito a continuación para R y C iguales (10k y 100nF respectivamente).

Código MATLAB

Código MATLAB:

RC HPF de segundo orden

SalidaMATLAB:Respuesta escalonada de RC HPF de segundo orden

Creo que los cálculos de zeta y Q son correctos, algo debe estar mal en mi comprensión de su papel en la respuesta al escalón. ¿Por qué el circuito se sobrepasa si está "sobreamortiguado"?

He repetido este procedimiento con un RC LPF pasivo de segundo orden (mismos valores R y C) y no veo ningún exceso. Esto está de acuerdo con mi intuición. ¿Por qué el HPF es diferente?

Gracias por cualquier ayuda que se pueda proporcionar.

La Q de un filtro de dos o más órdenes ciertamente está relacionada con los valores de los componentes, y la zeta = 1/2Q, por lo que también está relacionada con los valores de los componentes.
Pansy: gracias por esta pregunta básica pero interesante.
@EJP. Para esta configuración particular, no estoy de acuerdo. La forma estándar para el denominador es s^2 + (wn/Q)*s + wn^2. Se puede ver en el TF que la frecuencia natural es sqrt(1/(R^2*C^2))=1/RC. Entonces, desde la teoría, se supone que el coeficiente de s en el denominador es igual a wn/Q, y en este circuito, el coeficiente de s es 3/RC. Entonces, si deja wn = 1/RC, establezca 3/RC = wn/Q = 1/QRC. Entonces, 3/RC = 1/QRC, por lo tanto, los valores de RC se cancelan y Q siempre es 1/3 para este circuito, independientemente de los valores de R y C. Zeta será 1/2Q = 3/2, independientemente de R y C.

Respuestas (3)

El doble cero en s = 0 da lugar al rebrote.

Para facilitar el análisis, es mejor normalizar el TF a ω norte = 1 , así dividir ω norte por 1000 dar:

GRAMO ( s ) = s 2 s 2 + 3 s + 1

Multiplicar por 1 s para obtener la respuesta al escalón: R ( s ) = s ( 1 s 2 + 3 s + 1 )

Ahora encuentre la transformada inversa de Laplace del término entre paréntesis y luego diferencie (multiplique por s = diferenciación) para determinar r ( t ) , de este modo:

1 s 2 + 3 s + 1 = 0,45 s + 0.38 0,45 s + 2.62 0,45 mi 0.38 t 0,45 mi 2.62 t
diferencie para dar la respuesta al escalón:

r ( t ) = d d t ( 0,45 mi 0.38 t 0,45 mi 2.62 t ) = 0.17 mi 0.38 t + 1.17 mi 2.62 t

esto comienza en r ( 0 ) = 1.0 , luego cae a un mínimo (undershoot), r ( 1.72 ) = 0.0755 , y se establece en r ( ) = 0

Finalmente, escala el eje de tiempo por el factor de normalización, 1 / 1000 , dando el rebrote de 0.0755 en t = 1.72 metro s

Chu, buen trabajo!
@Chu. Gracias, es una forma interesante de resolverlo. Estoy recibiendo una respuesta diferente. Obtengo que las raíces del denominador sean -0.38 y -2.62. Después de derivar, r(t) = -0.17e^(-0.38t) + 1.17e^(-2.62t). La función de residuos de MATLAB proporciona los mismos números para el sistema. Esto proporciona una expresión para la respuesta en el dominio del tiempo, pero ¿qué nos dice zeta sobre este sistema?
2.62 es correcto, lo editaré. Cuando hay un cero (o más de uno), a menudo da lugar a un sobreimpulso o un subimpulso en un sistema sobreamortiguado, ya que tiene el efecto de agregar el diferencial a la respuesta natural. Si hay un cero en el plano s positivo, la respuesta transitoria puede comenzar a ir en la dirección "opuesta" antes de seguir a la respuesta natural. En tales casos, el modelo clásico de segundo orden no es aplicable y zeta no es un buen indicador de la forma transitoria.
... con un cero o dos en el TF, la forma de la respuesta transitoria está determinada por las posiciones de los polos y ceros. Los polos determinan la estabilidad relativa y definen la respuesta 'natural', pero los ceros pueden aumentar significativamente la respuesta natural.

He repetido este procedimiento con un RC LPF pasivo de segundo orden (mismos valores R y C) y no veo ningún exceso. Esto está de acuerdo con mi intuición. ¿Por qué el HPF es diferente?

No es lo mismo que el sobreimpulso que obtendría de un filtro de paso bajo de segundo orden con zeta por debajo de la unidad.

Ha creado un filtro de paso alto y no puede pasar CC, por lo que encontrará que el área de la forma de onda por encima de cero es exactamente igual al área de la forma de onda por debajo de cero. Esto es lo que hace un filtro de paso alto y no tiene nada que ver con la amortiguación.

¿Diría que un filtro de paso alto simple de primer orden tiene un sobreimpulso debido a problemas de relación de amortiguación: -

ingrese la descripción de la imagen aquí

No, no diría que un HPF de primer orden tiene problemas de sobreimpulso. La salida decae exponencialmente a cero, sin sobrepasarlo. Con un LPF de primer orden, aumenta exponencialmente hasta el valor de estado estable (sin sobreimpulso). Ambos circuitos de primer orden parecen ser pseudoinversos entre sí. ¿Qué cambia en los sistemas de segundo orden? ¿Está diciendo que zeta no tiene sentido cuando se consideran HPF de segundo orden?
@pansy, los dos filtros de paso alto de primer orden que ha conectado en cascada podrían separarse/almacenarse entre sí y seguir dando el mismo tipo de respuesta. La salida del primero dará la respuesta de caída a cero sin negatividad, pero esta salida luego forma la entrada al segundo filtro de paso alto (primer orden) y esto, naturalmente, genera la parte negativa de la forma de onda con la que pensó que estaba relacionado. zeta. Es solo una coincidencia.

En mi opinión, la respuesta no es tan sencilla como parece a primera vista. Las siguientes consideraciones son más o menos generales y no se adaptan únicamente al circuito dado.

Creo que podríamos (¿debemos?) usar la transformación de paso bajo a paso alto para ver qué sucede. Esta transformación consiste en una inversión simple de la variable de frecuencia compleja "s". Eso significa: en la función de transferencia de paso bajo correspondiente reemplazamos "s" bei "1/s".

Como siguiente paso, consideramos la respuesta al impulso h(t) que no es más que la transformada inversa de LAPLACE de la función de transferencia de paso alto H(1/s). Más que eso, sabemos que la respuesta escalonada g(t) es igual a la integral de tiempo sobre h(t) que, en el dominio de la frecuencia, es equivalente a una multiplicación con (1/s). Por lo tanto, tenemos que encontrar la transformada inversa de LAPLACE para una función " H(1/s)/s ".

Ahora, puede demostrarse que la transformada inversa de LAPLACE para tal expresión contiene el producto de (a) la respuesta de impulso de paso bajo h(t) y (b) la función de Bessel de orden cero Io, que es la causa del subimpulso observado ( la función de Bessel Io exhibe un comportamiento oscilatorio). La derivación matemática exacta es bastante complicada (y se puede encontrar, por ejemplo, en Claude S. Lindquist: Active Network Design).

EDITAR: para una comprensión mejor y más "intuitiva", es útil recordar que cualquier paso alto tiene, en principio, un comportamiento "diferenciador" que no sea un paso bajo con propiedades de integración. Eso significa que es la PENDIENTE de las variables de estado dentro del circuito lo que importa y determina la forma de la respuesta al escalón. Por lo tanto, dependiendo de las constantes de tiempo, podemos tener rebasamientos insuficientes, pequeños o incluso grandes.

No sigo bien su explicación / notación. Entiendo que la respuesta de impulso del sistema es la transformada inversa de Laplace del sistema sin multiplicar el sistema por la fuente del dominio de frecuencia (dado que la función de impulso se transforma a 1 en el dominio de frecuencia, entonces 1*sistema = sistema). Cuando dices "H(1/s)", ¿significa "H de 1/s"? Si mi sistema es H(s) "H de s", entonces h(t) = Laplace inversa (H(s)), ¿verdad? Entonces, ¿está diciendo que reemplace cada "s" en la función con "1/s"?
Sí, escribir H(1/s) significa reemplazar en la función de paso bajo la variable de s por 1/s. Este es el contenido de la transformación de paso bajo a paso alto.