Diagrama de Bode de polos que están cerca uno del otro

Tengo una función de transferencia de la siguiente forma:

  T ( j w ) = a 0 / ( 1 + s / w 1 ) ( 1 + s / w 2 ) ( 1 + s / w 3 )

donde a0 = 3600 w1= 1MHz w2= 4MHz w3= 40MHz

Dibujé el diagrama de Bode y verifiqué dos veces con matlab y encontré algunas discrepancias. Descubrí que el problema ocurre entre 1 MHz y 4 MHz ya que los polos están cerca uno del otro (menos de 1 década de diferencia).

Entonces, dado que a 0 Hz, la ganancia es de alrededor de 71 db, esperaba que a 1 MHz, la gráfica comenzara a disminuir con 20 db/dec. Entre 4Mhz y 1Mhz, hay (   yo o gramo ( 4 METRO / 1 METRO ) = 0.6 década (no una década completa) y, por lo tanto, a 4MHz la ganancia es 71 - log(4M/1M)*20 = 71 - 12 = 59 dB.

Sin embargo, en MATLAB, a 4 MHz la ganancia es de 69 db. Lo que significa que la ganancia se redujo en 8 dB, no en 12 dB. ¿Puedes decirme dónde está el flujo en mi comprensión?

Sé que este problema ocurre porque los polos están cerca uno del otro con menos de 1 década de diferencia, y no estoy seguro de cómo se manejan estos casos.

Entonces, el propósito de Mi punto principal es ¿cuál es la ganancia a 4M (2pi) rad/s? O en otras palabras, ¿cuánto caerá la ganancia entre 2pi*1M rad/sy 2pi*4M rad/sy por qué?

Mi trama (El eje x se multiplica por 2pi y está en rad/s):

ingrese la descripción de la imagen aquí

diagrama de MATLAB:

ingrese la descripción de la imagen aquí

Entonces, ¿qué dice MATLAB y por qué no usar Excel y trazarlo con más detalle? ¿Por qué está preguntando el intercambio de pila?
Verifique la escala de dB (su dibujo no es lineal) así como la escala de frecuencia (tampoco es correcta).
@Andyaka En MATLAB, a 4MHz la ganancia es 71 - 8 = 63dB... así que dice básicamente de 1MHz a 4MHz, son solo 0,4 décadas, no 0,6 décadas como en mis cálculos
Descargue papel cuadriculado log-Lin de papersnake
Debe mostrar lo que matlab dice en forma de diagrama. Sé lo que está sucediendo, pero debe trabajar un poco para que esta sea una mejor pregunta que pueda ser útil para otros.
@Chu perdón por la mala trama. Es solo un boceto. El eje x se multiplica por 2pi y está en rad/s. Mi punto principal es ¿cuál es la ganancia a 4M (2pi) rad/s? O en otras palabras, ¿cuánto caerá la ganancia entre 2pi*1M rad/s y 2pi*4M rad/s?
@Andyaka Espero que sea mejor ahora y que sea útil para otros más adelante

Respuestas (2)

¿Ha considerado que la trama de MatLab ha sido mal interpretada? Por ejemplo, todas las ubicaciones de los polos en la pregunta están etiquetadas en MHz, pero el eje del gráfico de MatLab muestra rad/seg. Si multiplicas las frecuencias por 2 π , las ubicaciones de los polos pasan a ser 6,28 Mrad/seg, 25,1 Mrad/seg y 251 Mrad/seg.

Ahora, está realizando un análisis asintótico, por lo que los números que obtenga estarán un poco alejados del resultado de MatLab, pero debe realizarse correctamente después del cambio en la escala de frecuencia. Por ejemplo, después de corregir el error asintótico, los dos primeros polos deberían tener magnitudes de 68,1 dB a 6,28 Mrad/seg y 55,8 dB a 25,1 Mrad/seg.

Tenga en cuenta que no tendrá un diagrama de Bode completo hasta que agregue también un diagrama de fase.

@Andy, gracias. No eres el único que puede necesitar una prueba de la vista, tuve que hacer zoom en el diagrama de MatLab en un 200% o más para leer las unidades del eje.

Esta respuesta perdió el punto (imagen muy pequeña en la pregunta del operador) de que la respuesta de frecuencia estaba en radianes por segundo en lugar de Hz. Sin embargo, lo dejo aquí porque sigue siendo un punto bastante importante sobre lo que sucede cuando dos puntos de ruptura de frecuencia están cerca: tienden a fusionarse en un solo punto de ruptura de segundo orden a mitad de camino logarítmicamente entre los dos puntos de ruptura originales.

Respuesta original: -


Bien, ahora que tiene la imagen de MATLAB, puedo confirmar que potencialmente está obteniendo un efecto de "pico" de segundo orden. Olvídese del término a 40 MHz; esto es irrelevante para el rango de 1 MHz a 4MHz. Así que ahora tu TF se convierte en: -

A 0 ( 1 + s / ω 1 ) ( 1 + s / ω 2 )

Si multiplicas esto obtienes: -

ω 1 ω 2 A 0 s 2 + s ( ω 1 + ω 2 ) + ω 1 ω 2

En otras palabras, tiene la forma de un filtro de paso bajo clásico de segundo orden que tiene una resonancia neta en ω 1 ω 2 o 2 MHz. Cuánto emerge ese pico por encima de la línea de base de pendiente descendente depende del término medio ω 1 + ω 2 : -

ingrese la descripción de la imagen aquí

El término medio ( ω 1 + ω 2 ) incorpora zeta ( ζ ) con la ecuación entera de la forma: -

H ( s ) = ω 0 2 s 2 + s ( 2 ζ ω 0 ) + ω 0 2

Dónde ω 0 es ω 1 ω 2 o 2 MHz

Básicamente, si las dos frecuencias están cerca, la zeta se reduce y se notará un efecto de combinación. En su caso, zeta es 1.25 y esto no producirá un pico notable para el ojo desinformado, pero tenderá a juntar los dos puntos de 3dB para convertirse en un punto de 3dB a más de 2MHz, es decir, lo suficiente como para confundirlo en lo predicho (pero demasiado simple ) diagrama de Bode.

Tenga en cuenta que aunque he mostrado la ecuación general para este tipo de filtro de paso bajo de segundo orden, debido a la naturaleza de la fórmula original, zeta nunca puede ser inferior a 1, por lo tanto, los dos puntos separados de -3dB se combinan para dar un nuevo punto de ruptura. logarítmicamente entre 1 MHz y 4 MHz, es decir, 2 MHz.

Básicamente, para las frecuencias 1 y 4, el efecto es pequeño. Numéricamente, esas dos frecuencias dan una zeta de 5/(2*2) que es 1,25 y, a medida que la zeta crece, el efecto es cada vez menor. Zeta solo puede crecer a medida que las frecuencias se amplían, así que supongo que en 1 y 4 podría argumentar que este es un punto de umbral, pero para otras aplicaciones podría estar más lejos. Es lo que es y diferente dependiendo de su punto de vista.
Oye, calculaste zeta LOL
Sí, pero lo eliminé, pensé que estaba mal. ¡Gracias por la confirmación!
Por supuesto, se pueden escribir libros completos sobre este tema y las implicaciones son profundas en todos los ámbitos de la vida: la resonancia mecánica y el desastre del puente de Tacoma ( youtube.com/watch?v=3mclp9QmCGs ) salta a 1 minuto.
Pero este TF tiene tres polos reales de primer orden, por lo que no puede haber un pico de resonancia.
@Chu No digo que tenga un pico resonante. Al final llegué a la conclusión (sin sesgo) de que zeta es mayor que 1, por lo tanto, los puntos de ruptura en 1 MHz y 4 MHz forman un punto de ruptura combinado en 2 MHz. Tener un tercer punto de interrupción a 40 MHz es en gran medida irrelevante, pero estoy preparado para dejarme influenciar por la sabiduría.