Cómo se relacionan los polos con la respuesta de frecuencia

16

Recientemente he caído en la falacia , considerando el polo s = 1 ya que hay una respuesta infinita en la frecuencia 1. Sin embargo, la respuesta fue solo 1. Ahora, ¿puedes derivar la respuesta de frecuencia, dados los polos?

En segundo lugar, la teoría dice que un sistema es estable cuando los polos están en el plano s izquierdo y, por lo tanto, decaen con el tiempo. Pero espera. ¿Un "polo" significa la respuesta infinita: el crecimiento en el tiempo?

Finalmente, ¿es la pregunta correcta en DSP? OMI, D significa digital mientras que s-domain es analógico. No encuentro etiquetas de transformación s-plane o Laplace para etiquetar mi publicación.

actualización Gracias por las respuestas. Parece que lo tengo, excepto una cosa menor pero fundamental: la relación de los polos (y ceros) con la frecuencia. Básicamente, ¿por qué son los valores propios (o, ¿cómo se llama a la s operador / variable) relacionada con la frecuencia? Debería estar relacionado de alguna manera con el crecimiento exponencial y la transformación de Laplace. Entiendo perfectamente que los polos son valores propios (especialmente para las recurrencias discretas). Pero, ¿cómo se relaciona esto con la frecuencia?

Val
fuente
1
Es "intercambio de pila de procesamiento de señal", no "intercambio de pila DSP". :)
endolith
44
Sí, como se menciona en endoith, el procesamiento de señales analógicas es un tema. DSP.SE era un nombre conveniente para el lanzamiento inicial, pero signal.stackexchange.com ahora también enlaza aquí.
Datageist
¿Qué quiere decir exactamente cuando pregunta por la relación entre polos y frecuencias?
Sudarsan
1
Obviamente, es cómo y por qué los polos determinan la respuesta de frecuencia.
Val
La respuesta ya ha sido dada, supongo. La respuesta de frecuencia es la magnitud de la respuesta del sistema a medida que avanza la jω eje. Si ha factorizado la Función de Transferencia del Sistema H(s) en el Producto de 1/(spi) y (szi) , todo lo que necesita hacer es encontrar la magnitud en s=jω para la Transferencia La función y esto obviamente está determinada por la ubicación de los polos y ceros, ya que serán los que aparezcan en la respuesta factorizada del sistema.
Sudarsan el

Respuestas:

12

Creo que en realidad hay 3 preguntas en su pregunta:

P1: ¿Puedo derivar la respuesta de frecuencia dados los polos de un sistema (lineal invariante en el tiempo)?

Sí, puedes, hasta una constante. Si s,i , i=1,,N, son los polos de la función de transferencia, puede escribir la función de transferencia como

(1)H(s)=k(ss,1)(ss,2)(ss,N)

Tenga en cuenta que s es una variable compleja s=σ+jω , y la variable de frecuencia ω corresponde al eje imaginario del plano s complejo . Ahora necesitamos obtener la respuesta de frecuencia de la función de transferencia. Para sistemas estables, esto se puede hacer simplemente evaluando la función de transferencia H(s) para s=jω . Entonces reemplaza s por jω en (1) y listo. Sin embargo, tenga en cuenta que esto solo es cierto para sistemas estables (es decir, si la región de convergencia de H(s) incluye elejejω ).

P2: ¿Cómo puede un sistema estable tener polos?

Como ya sabe, para sistemas causales y estables, todos los polos deben estar en el semiplano izquierdo del plano s complejo . De hecho, el valor de la función de transferencia H(s) irá al infinito en un polo s=s , pero la respuesta de frecuencia estará bien, porque si todos los polos están en el semiplano izquierdo, no hay polos en el jω -axis (o a la derecha de la misma). Si lo observa en el dominio del tiempo, entonces cada polo (simple) tiene una contribución de est a la respuesta de impulso del sistema. Si el polo está ubicado en el semiplano izquierdo, esto significa que s=σ+jω tiene una parte real negativaσ<0 . Entonces

est=eσejω

es una función amortiguada exponencialmente y no crece sino que decae, porque σ<0 .

P3: ¿Esta pregunta pertenece aquí?

Otros miembros de la comunidad tienen que juzgar si esta pregunta pertenece aquí. Creo que si. Obviamente, no está directamente relacionado con DSP puro, pero los ingenieros de DSP a menudo también tienen que lidiar con señales y sistemas analógicos antes de la conversión AD, por lo que también conocen la teoría de sistemas continuos. En segundo lugar, casi todas las personas DSP (al menos las que tienen capacitación tradicional) tuvieron bastante exposición a las señales generales y la teoría de sistemas, incluidos los sistemas de tiempo continuo y de tiempo discreto.

Por cierto, para sistemas de tiempo discreto obtienes la transformación Z lugar de la transformación de Laplace, y tu variable compleja ahora se llama z lugar de s . La variable D que ha mencionado se define como D=z1 y se utiliza principalmente en la literatura de codificación. Por su definición, denota un elemento de retraso, por lo que D significa "retraso" (no "digital").

Si se sabe que la izquierda semiplano del complejo s mapas un plano a la región dentro del círculo unidad del complejo z un plano (es decir, |z|<1 ), y el jω eje x mapas al círculo unidad |z|=1 , entonces casi todo lo que sabe sobre uno de los dos dominios se transferirá fácilmente al otro dominio.

Matt L.
fuente
Creo que la respuesta de frecuencia implica una conjugación compleja además de s en H (s) para s = jω.
Val
16

Una cosa que realmente me ayudó a entender los polos y ceros es visualizarlos como superficies de amplitud. Varias de estas parcelas se pueden encontrar en A Filter Primer . Algunas notas:

  • Probablemente sea más fácil aprender primero el plano S analógico y, después de comprenderlo, aprender cómo funciona el plano Z digital.
  • Un cero es un punto en el que la ganancia de la función de transferencia es cero.
  • Un polo es un punto en el que la ganancia de la función de transferencia es infinita.
  • A menudo hay ceros o polos en el infinito, que no siempre se incluyen en las descripciones de la función de transferencia, pero son necesarios para comprenderla.
  • La respuesta de frecuencia en el plano S ocurre solo a lo largo del eje jω.
    • El origen es 0 Hz, o CC, y la frecuencia de corte de los filtros aumenta radialmente lejos del origen. Poner un poste en cualquier punto a lo largo de un círculo a una cierta distancia del origen producirá la misma frecuencia de corte.
    • Para aumentar la frecuencia de corte de un filtro, mueva los polos radialmente hacia afuera.
    • Para aumentar la Q de un filtro biquad, mueva los polos a lo largo del círculo hacia el eje jω, lo que mantiene constante la frecuencia de corte, pero aumenta el efecto que tiene el polo en la respuesta de frecuencia, haciéndolo más "pico".
    • Mover polos a lo largo de un círculo mantiene constante la frecuencia de corte pero cambia Q
  • Si aparece un cero en el eje jω, entonces la respuesta de frecuencia caerá a cero a esa frecuencia; Si ingresa una onda sinusoidal a esa frecuencia, la salida será 0.
  • Si aparece un polo en el eje jω, entonces la respuesta al impulso es un oscilador; cualquier impulso hará que suene para siempre a esa frecuencia. Los impulsos tienen energía finita, pero la respuesta del filtro tiene energía infinita, por lo que tiene ganancia infinita.

Un ejemplo simple es un integrador H (s) = 1 / s:

  • Esta función es igual a 0 cuando s es infinito, por lo que tiene un cero en el infinito.
  • Esta función es igual a infinito cuando s es cero, por lo que tiene un polo en cero.

En otras palabras, tiene una ganancia infinita en DC (la respuesta escalonada de un integrador aumenta para siempre), y la ganancia disminuye a medida que aumenta la frecuencia:

Bode plot de integrador

Alejar el polo del origen, a lo largo del eje imaginario hacia la mano izquierda del plano S, hace que la ganancia a 0 Hz en el eje jw sea finita nuevamente, y ahora tiene un filtro de paso bajo:

ingrese la descripción de la imagen aquí

endolito
fuente
+1, buena respuesta. Pero no entiendo lo que quiere decir con "cualquier punto a lo largo de un círculo a una cierta distancia del origen tiene la misma frecuencia". Las curvas de frecuencia constante en el plano son líneas paralelas al eje real. Para círculos con origen en s = 0 obtienes σ 2 + ω 2 = c o n s t , donde s = σ + j ω . ss=0σ2+ω2=consts=σ+jω
Matt L.
Parece confundir el plano s con el plano z
Val
@MattL .: Hmmm. Estoy pensando en los polos de un filtro Butterworth de orden N a lo largo de un círculo equidistante del origen, por ejemplo, o en los polos de una biquad que se mueve a lo largo de un círculo equidistante del origen a medida que ajusta la Q del filtro mientras mantiene la constante de frecuencia, o cambiando el corte de un filtro moviendo los polos más cerca o lejos del origen en una dirección radial, o convirtiendo paso bajo a paso alto invirtiendo los polos alrededor del círculo unitario. ¿Cómo debo reformular esto?
Endolith
@Val: frecuencia de corte . Ya he editado la publicación para corregirlo.
Endolith
66
Val, no hay necesidad de un comentario sarcástico e irritable para @endolith.
Spacey
1

No diré el mapeo completo de los polos (1) / ceros (0) a la respuesta de frecuencia, pero creo que puedo explicar la conexión entre la frecuencia y la respuesta cero / infinito, ¿por qué tiene una respuesta infinita / cero en es decir , lo que e - j w tiene que ver con z .ejw=zzero/pole,ejwz

La forma general del sistema lineal es que puede ser resuelto en z-from como Y ( z ) = ( b 0 + b 1 z + b

yn+a1yn1+a2yn2+=b0xn+b1xn1+b2xn2+,
Y(z)=(b0+b1z+b2z2+)(1+a1z+a2z2+)X(z)=H(z)X(z)=(1z0z)(1z1z)(1p0z)(1p1z)X(z).

(1z0z)11p0z puede considerarse como una serie de sistemas, donde la primera salida es la entrada para otro.

H(z)X(z)Y(z)=(1z0z)Χ(z),yn=b0xn+b1xn1. Let's take b0=b1=1 for simplicity. I mean that yn=xn+xn1.

What we want to determine the effect of the system H(z) upon harmonic signal. That is, the input is going to be test signal

xn=ejwnz1+ejwz+e2jwz2+=1/(1ejw)=X(z).
The response is going to be
yn=xn+xn1|xn=ejwn=ejwn+ejw(n1)=ejwn(1+ejw)
that is, 1+ejw is the transfer function or Y(z)=(1+z)(1ejwz)=(1+z)X(z).

Please note that 1+z basically says that output is sum of input signal plus shifted signal, since single z stands for single clock delay in time domain.

Now, as explained in, H(jw)=1+ejw=ejw/2(ejw/2+ejw/2)=ejw/22cos(w/2). Cosine makes it to behave like low-pass filter

{w=0H(j0)=12cos(0)=2w=πH(jπ)=ejπ/22cos(π/2)=0

It is also a good lesson that you get 2cosα=eiα+eiα because you will supply the real signals rather than complex imaginary ones in real life.

LTI with impulse response = {1,-1} is yn=xnxn|xn=ejwn=ejwn(1ejw) has transfer function of H(jw)=(1ejw)=ejw/2(ejw/2ejw/2)=ejw2sin(w/2), which has zero at w=0 since sin(0)=0 but it can be found from the frequency response

H(jw)=1ejw=0ejw=1=e0w=0.

After the textbooks, I can spot the surprising coincidence between transfer function H(z)=1±z and frequency response H(jw)=1±ejw. That is, z somehow corresponds to ejw, which is important for zero/pole analysis. I read it like

sine z-factor stands for a clock shift and yn=xn±xn1=0 means that next sample is ± previous one to get zero response, we need to have 1±z=0 in front of X(z). But, the frequency domain basis functions ejwn evolve by multiplying current value ejw(n1) with ejw every clock. Therefore, we have ejwn(1±ejw)=0 as condition for constant zero output. The latter 1±ejw matches perfectly with zero transfer function 1±z=0.

In general, single-zero LTI is given by yn=b0xn+b1xn1 or

Y(z)=(b0+b1z)X(z)=(b0+b1z)(1+x1z+x2z2+)=b0+(b0x1+b1x0)z+(b0x2+b1x1)z2+.
When b0+b1z=0, i.e. when z=b0/b1, whereas frequency response is,
yn(xn=ejwn)=b0ejwn+b1ejw(n1)=ejwn(b0+b1ejw)=ejwnb0(1z0ejw),

which goes to zero when 1z0ejw=0 or ejw=1/z0, which matches the computation for z if z=ejw. The only thing that bothers me is that fixed-amplitude complex exponential is not enough for the frequency (harmonic) basis. You cannot obtain arbitrary ratio 1/z0=ejw by choosing appropriate frequency w, a decaying harmonic signal is needed for that. That is weird because I have heard that any signal can be represented as sum of (constant amplitude) sines and cosines. But, anyway, we see that system zero stands for relationship between adjacent samples of input signal. When they are right, the output is identically 0 and we can choose such such frequency w so that zero z=1/z0=ejw.

Now, what about the poles? Let's single out a single pole a. The system has a from of yn=ayn1+(xn+xn1+), under assumption y0=0, has z-transform of Y(z)=X(z)/(1az).

The feedback a is equivalent to infinite impulse response 1,a,a2,z1+az+a2z2+=1/(1az). It says that response is infinite when z=1/a. What does it mean if we apply the test signal

xn=ejwnzX(z)=1+ejwz+e2jwz2+=1/(1ejwz)
to our system? We'll get Y(z)=11az11ejwz, or
yn=ejwn+aejw(n1)+a2ejw(n2)+=ejwn(1+aejw+a2e2jw+)=ejwn1aejw.
That is, frequency response is 1/(1aejw), which goes to infinity when ejw=1/a, the same as zpole above, ejw=zpole=1/a. But again, you can not always arrive at the pole 1/a adjusting the frequency w alone. The frequency basis functions must be decaying amplitude in general and look like (kejw)n.

That is, zeroes or poles of the transfer function H(z) happen to match the zeroes and poles of frequency response H(jw), which is really amazing. I noticed that this is related to the relation between adjacent samples, ejwn/ejw(n1)=ejw=1/zzero in case of zeroes. The fact that ejwn scales exponentially over time, along with the system with feedback a, also seems to be the key for matching between ejw and zpoles. It also seems important that you cannot simply look for the appropriate frequency of ejwn, the basis function must also have adjustable amplitude factor kn.

I would be happy if anybody could explain the same more condensely or more crisply.

Valentin Tihomirov
fuente