Como os pólos estão relacionados à resposta em frequência

16

Recentemente, caí na falácia , considerando o polo s = 1, pois há uma resposta infinita na frequência 1. No entanto, a resposta foi de apenas 1. Agora, você pode derivar a resposta de frequência, dados os polos?

Segundo, a teoria diz que um sistema é estável quando os polos estão no plano s esquerdo e, portanto, decaem no tempo. Mas espere. Não significa "pólo" a resposta infinita - o crescimento no tempo?

Finalmente, é a pergunta certa no DSP? IMO, D significa digital, enquanto o domínio s é analógico. Não encontro tags s-plane ou Laplace transform para rotular minha postagem.

atualizar Obrigado pelas respostas. Parece que eu entendi, exceto a coisa menor, mas fundamental - a relação dos polos (e zeros) com a frequência. Basicamente, por que os valores próprios (ou, como você chama o operador / variável ) estão relacionados à frequência? De alguma forma, deve estar relacionado ao crescimento exponencial e à transformação de Laplace. Entendo perfeitamente que os pólos são autovalores (especialmente para recorrências discretas). Mas, como isso está relacionado à frequência?s

Val
fonte
1
É "Troca de pilha de processamento de sinal", não "Troca de pilha DSP". :)
endolith 12/06
4
Sim, como endoith mencionou, o processamento de sinal analógico está no tópico. DSP.SE era um nome conveniente para o lançamento inicial, mas o signs.stackexchange.com agora também está vinculado aqui.
datageist
O que exatamente você quer dizer quando pede a relação entre poloneses e frequências?
Sudarsan 04/08/13
1
Obviamente, é como e por que os polos determinam a resposta em frequência.
Val
A resposta já foi dada, eu acho. A resposta em frequência é a magnitude da resposta do sistema à medida que você se move ao longo do eixo . Se você fatorou a Função de transferência do sistema no produto de e , tudo o que você precisa fazer é encontrar a magnitude em para a função de transferência e isso obviamente é determinado pelo local dos pólos e zeros, pois eles serão os que aparecerão na resposta do sistema fatorado. jωH(s)1/(spi)(szi)s=jω
Sudarsan 04/08/13

Respostas:

12

Eu acho que existem 3 perguntas na sua pergunta:

Q1: Posso derivar a resposta de frequência, considerando os polos de um sistema (linear invariante no tempo)?

Sim, você pode, até uma constante. Se , são os pólos da função de transferência, você pode escrever a função de transferência comos,ii=1,,N,

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

Observe que é uma variável complexa , e a variável de frequência corresponde ao eixo imaginário do plano complexo . Agora precisamos obter a resposta de frequência da função de transferência. Para sistemas estáveis, isso pode ser feito simplesmente avaliando a função de transferência para . Então você substitui por em (1) e pronto. Observe, no entanto, que isso só é verdade para sistemas estáveis ​​(ou seja, se a região de convergência de incluir o eixo ).ss=σ+jωωsH(s)s=jωsjωH(s)jω

P2: Como um sistema estável pode ter postes?

Como você já sabe, para sistemas causais e estáveis, todos os pólos devem estar no semiplano esquerdo do plano complexo . De fato, o valor da função de transferência irá para o infinito no polo , mas a resposta de frequência será boa, porque se todos os polos estiverem no semiplano esquerdo, não haverá pólos no eixo (ou à direita). Se você observar no domínio do tempo, cada pólo (simples) terá uma contribuição de para a resposta de impulso do sistema. Se o polo estiver localizado no semi-plano esquerdo, isso significa que tem uma parte real negativa . entãosH(s)s=sjωests=σ+jωσ<0 0

est=eσejω

é uma função exponencialmente amortecida e não cresce, mas decai, porque .σ<0 0

Q3: Esta pergunta pertence aqui?

Outros membros da comunidade devem julgar se esta pergunta pertence a este lugar. Eu acho que sim. Obviamente, isso não está diretamente relacionado ao DSP puro, mas os engenheiros do DSP muitas vezes também precisam lidar com sinais e sistemas analógicos antes da conversão do AD, para que também saibam sobre a teoria contínua dos sistemas. Segundo, quase todas as pessoas com DSP (pelo menos aquelas com treinamento tradicional) tiveram bastante exposição a sinais gerais e à teoria de sistemas, incluindo sistemas de tempo contínuo e de tempo discreto.

A propósito, para sistemas de tempo discreto, você obtém a transformação vez da transformação de Laplace, e sua variável complexa agora é chamada vez de . A variável que você mencionou é definida como e é usada principalmente na literatura de codificação. Por sua definição, denota um elemento de atraso, então significa "atraso" (não "digital").ZzsDD=z-1D

Se você souber que o meio plano esquerdo do plano complexo mapeado para a região dentro do círculo unitário do plano complexo (ou seja ), e o eixo mapeia o círculo unitário , quase tudo o que você sabe sobre um dos dois domínios será facilmente transferido para o outro domínio.sz|z|<1jω|z|=1

Matt L.
fonte
Penso que a resposta em frequência envolve conjugação complexa além de s em H (s) para s = jω.
Val
16

Uma coisa que realmente me ajudou a entender pólos e zeros é visualizá-los como superfícies de amplitude. Várias dessas plotagens podem ser encontradas em A Filter Primer . Algumas notas:

  • Provavelmente é mais fácil aprender primeiro o plano S analógico e, depois de entendê-lo, aprenda como o plano Z digital funciona.
  • Um zero é um ponto no qual o ganho da função de transferência é zero.
  • Um polo é um ponto no qual o ganho da função de transferência é infinito.
  • Frequentemente, existem zeros ou pólos no infinito, que nem sempre são incluídos nas descrições da função de transferência, mas são necessários para entendê-la.
  • A resposta de frequência no plano S ocorre apenas ao longo do eixo jω.
    • A origem é 0 Hz, ou DC, e a frequência de corte dos filtros aumenta radialmente para longe da origem. Colocar um polo em qualquer ponto ao longo de um círculo a uma certa distância da origem produzirá a mesma frequência de corte.
    • Para aumentar a frequência de corte de um filtro, mova os pólos radialmente para fora.
    • Para aumentar o Q de um filtro de biquad, mova os pólos ao longo do círculo em direção ao eixo jω, o que mantém a frequência de corte constante, mas aumenta o efeito que o pólo tem na resposta de frequência, tornando-o mais "pico".
    • Mover polos ao longo de um círculo mantém a frequência de corte constante, mas muda Q
  • Se um zero aparecer no eixo jω, a resposta da frequência cairá para zero nessa frequência; se você inserir uma onda senoidal nessa frequência, a saída será 0.
  • Se um polo aparecer no eixo jω, a resposta ao impulso será um oscilador; qualquer impulso fará com que ele toque para sempre nessa frequência. Os impulsos têm energia finita, mas a resposta do filtro tem energia infinita, por isso possui ganho infinito.

Um exemplo simples é um integrador H (s) = 1 / s:

  • Esta função é igual a 0 quando s é infinito, portanto, ele tem um zero no infinito.
  • Essa função é igual ao infinito quando s é zero, então ele tem um polo em zero.

Em outras palavras, ele possui um ganho infinito no DC (a resposta de etapa de um integrador aumenta para sempre) e o ganho diminui à medida que a frequência aumenta:

Gráfico Bode do integrador

Afastar o polo da origem, ao longo do eixo imaginário para a mão esquerda do plano S, torna o ganho a 0 Hz no eixo jw finito novamente e agora você tem um filtro passa-baixo:

insira a descrição da imagem aqui

endólito
fonte
+1, boa resposta. Mas não entendo o que você quer dizer com "Qualquer ponto ao longo de um círculo a uma certa distância da origem tem a mesma frequência". Curvas de frequência constante no plano são linhas paralelas ao eixo real. Para círculos com origem em s = 0 Você começa σ 2 + ω 2 = c O n s t , onde s = σ + j ω . ss=0σ2+ω2=consts=σ+jω
Matt L.
Ele parece confundir plano s com z-plano
Val
@MattL .: Hmmm. Estou pensando nos pólos de um filtro Butterworth de ordem n de estar em um círculo equidistante da origem, por exemplo, ou nos pólos de um biquad se movendo ao longo de um círculo equidistante da origem, conforme você ajusta o Q do filtro enquanto mantém a constante de frequência, ou alterando o corte de um filtro, movendo os pólos para mais perto ou para longe da origem em uma direção radial, ou convertendo a passagem baixa em passagem alta invertendo os pólos em torno do círculo unitário. Como devo reformular isso?
Endolith 13/06/2013
@ Val: frequência de corte . Já editei a postagem para corrigi-la.
endolith 13/06
6
Val, não há necessidade de um comentário sarcástico para @endolith.
Spacey
1

Não vou contar o mapeamento completo dos pólos (1) / zeros (0) para a resposta de frequência, mas acho que posso explicar a conexão entre frequência e resposta zero / infinita, por que você tem resposta infinita / zero em ou seja , o que e - j w tem a ver com z .ejw=zzero/pole,ejwz

A forma geral do sistema linear é que pode resolva em z-de 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).

No final, a série de produtos binomiais pode ser considerado como uma série de sistemas, onde a primeira saída é a entrada para outra.(1-z0 0z)11-p0 0z

Eu gostaria de analisar o efeito de pólo único e zero. Vamos destacar o primeiro zero, considerando a função de transferência para que o restante de seja o sinal de entrada, Y ( z ) = ( 1 - z 0 z ) Χ ( z ) , que corresponde a alguns y n = b 0 x n + b 1 x n - 1 . Vamos dar b 0H(z)X(z)Y(z)=(1-z0 0z)Χ(z),yn=b0xn+b1xn1. para simplificar. Quero dizer que y n = x n + x n - 1 .b0=b1=1yn=xn+xn1

O que queremos determinar o efeito do sistema H (z) sobre o sinal harmônico. Ou seja, a entrada será o sinal de teste A resposta será y n = x n + x n

xn=ejwnz1+ejwz+e2jwz2+=1/(1ejw)=X(z).
ou seja,1+ e - j w é a função de transferência ouY(z)= ( 1 + z )
yn=xn+xn1|xn=ejwn=ejwn+ejw(n1)=ejwn(1+ejw)
1+ejw.Y(z)=(1+z)(1ejwz)=(1+z)X(z)

Observe que basicamente diz que a saída é a soma do sinal de entrada mais o sinal deslocado, uma vez que z único representa atraso de relógio único no domínio do tempo.1+zz

Agora, como explicado em , . O cosseno faz com que ele se comporte como um filtro passa-baixo { w = 0 H(jw)=1+ejw=ejw/2(ejw/2+ejw/2)=ejw/22cos(w/2)

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

2cosα=eiα+eiα

yn=xnxn|xn=ejwn=ejwn(1ejw)H(jw)=(1ejw)=ejw/2(ejw/2ejw/2)=ejw2sin(w/2)w=0sin(0)=0

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

H(z)=1±zH(jw)=1±ejwejw

yn=xn±xn1=0±1±z=0ejwnejw(n1)ejwejwn(1±ejw)=01±ejw1±z=0

yn=b0xn+b1xn1

Y(z)=(b0+b1z)X(z)=(b0+b1z)(1+x1z+x2z2+)=b0+(b0x1+b1x0)z+(b0x2+b1x1)z2+.
b0+b1z=0z=b0/b1,
yn(xn=ejwn)=b0ejwn+b1ejw(n1)=ejwn(b0+b1ejw)=ejwnb0(1z0ejw),

1z0ejw=0ejw=1/z0zz=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)/(1-umaz).

O retorno uma é equivalente a resposta de impulso infinito 1,uma,uma2,...z1+umaz+uma2z2+=1/(1-umaz). Diz que a resposta é infinita quandoz=1/uma. O que significa se aplicarmos o sinal de teste

xn=ejWnzX(z)=1+ejWz+e2jWz2+=1/(1-ejWz)
ao nosso sistema? Nós conseguiremosY(z)=11-umaz11-ejWz, ou
yn=ejWn+umaejW(n-1)+uma2ejW(n-2)+=ejWn(1+umae-jW+uma2e-2jW+)=ejWn1-umae-jW.
Ou seja, a resposta em frequência é 1/(1-umae-jW), que vai para o infinito quando e-jW=1/uma, o mesmo que zpoeue acima, e-jW=zpoeue=1/uma. Mas, novamente, você nem sempre pode chegar ao poste1/uma ajustando a frequência Wsozinho. As funções de base de frequência devem estar em geral diminuindo a amplitude e parecer(kejW)n.

Ou seja, zeros ou pólos da função de transferência 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.

Ficaria feliz se alguém pudesse explicar o mesmo de forma mais condensada ou mais clara.

Valentin Tihomirov
fonte