Compreender o enredo de Nyquist

O gráfico Nyquist é uma ferramenta fundamental na engenharia de controle, usada para analisar a estabilidade dos sistemas de feedback. Nomeado após Harry Nyquist, esta representação gráfica retrata a resposta de frequência de um sistema, plotando a função de transferência complexa G(jω) como a frequência angular ω varia de 0 a . O gráfico é desenhado no plano complexo, com a parte real de G(jω)[]] no eixo horizontal e na parte imaginária no eixo vertical. O seu utilitário primário consiste em avaliar a estabilidade de circuito fechado a partir da resposta de frequência de circuito aberto, sem computar explicitamente os pólos de circuito fechado. O gráfico Nyquist também fornece a percepção sobre o ganho e as margens de fase, tornando- o indispensável para o design robusto do controlador.

Os fundamentos matemáticos das parcelas de Nyquist requerem uma compreensão completa da análise complexa, particularmente a avaliação das funções de transferência no eixo imaginário e a interpretação geométrica do mapeamento de contornos. Os engenheiros devem apreender o conceito do ] contorno de Nyquist, um caminho fechado no plano complexo que encerra todo o meio plano direito (RHP), e como este traçado de contornos através da função de transferência de circuito aberto G(s)H(s)[] para produzir o Caminho de Nyquist[. O gráfico resultante está a girar em torno do ponto crítico −1 + 0j[[] indica a estabilidade do sistema, seguindo o critério de estabilidade de Nyquist. Este artigo expande-se nestas bases matemáticas, percorrendo a lógica de construção, o princípio de argumentação e as pitfalls práticas.

Avaliação da função complexa de transferência

A função de transferência de circuito aberto G(s)H(s) é uma função racional da variável complexa s[ = σ + jω. Para análise de resposta de frequência, definimos σ = 0, então s[] = []jω. A expressão resultante ]G(jω)H(jω)] é um número complexo que depende de ω:

G(jω)H(jω) = Re(ω) + j Im(ω) = M(ω) e^{jω(ω)}[

Aqui, M(ω)] é a magnitude (também chamada ganho) e

Para calcular matematicamente G(jω)H(jω), considere uma função de transferência geral fatorada em pólos e zeros:

[[FLT: 0]]]G( s)H( s) = K \frac{\prod {i=1}^{m} (s - z i)}{\prod {k=1}^{n} (s - p k)}

Quando s = j[ω, cada termo torna-se um vetor complexo do ponto zero ou do pólo até o ponto j[ω no eixo imaginário. A magnitude total é o produto dos comprimentos dos vetores zero divididos pelo produto dos comprimentos dos vetores pólos, multiplicado por . A fase total é a soma dos ângulos dos vetores zero menos a soma dos ângulos dos vetores pólos, mais a fase de K (que é 0 se K é real e positiva). Esta interpretação geométrica é central para entender como a curva Nyquist se comporta, especialmente perto das frequências correspondentes às frequências naturais do sistema.

Passos de Construção Matemática

A construção de um gráfico Nyquist à mão (ou usando software) envolve a avaliação sistemática de G(jω)H(jω) sobre uma faixa de frequência. Os passos principais são:

  1. Identidade da função de transferência e escreva-a em Bode (pólos e zeros).Para sistemas com atrasos de tempo, incluir o fator ]e^{-jωT}].
  2. Computar a magnitude e a fase em um conjunto de frequências. Tipicamente, são escolhidas frequências baixas (perto de 0) e altas (perto de ..) com pontos intermediários perto de frequências de quebra (onde as contribuições de pólo ou zero mudam sinal).
  3. [[FLT: 0]]Iniciar em ω = 0[[FLT: 1]]: O ponto depende do tipo de sistema. Para o Tipo 0 (sem integradores), o gráfico começa em [[FLT: 2]]K[[[FLT: 3]] (o ganho DC) no eixo real positivo. Para o Tipo 1 (um integrador), o gráfico começa em uma magnitude infinita com uma fase de −90°. Para o Tipo 2 (dois integradores), a fase é de −180° a ω = 0.
  4. Plota pontos em frequências crescentes, conectando-os suavemente. À medida que a frequência aumenta, a curva geralmente se move em direção à origem (desde que a magnitude diminui). A fase pode cruzar -180° na frequência de crossover de fase onde a magnitude é crítica para análise de estabilidade.
  5. [[FLT: 0]]Mirror para frequências negativas (se usar o contorno completo de Nyquist). O gráfico para ω de − . . para 0 é o complexo conjugado do gráfico para ω positivo, refletido no eixo real. Combinado com o ramo de frequência positiva (e o semicírculo infinito se o sistema tiver um excesso de pólo- zero líquido), obtém- se o caminho completo de Nyquist.
  6. Identifique o número de circunscrições de −1 + 0j para aplicar o critério de estabilidade.

Matematicamente, a curva de Nyquist para frequências positivas pode ser parametrizada como ω . Para sistemas com grau relativo (n − m) ≥ 1, a curva termina na origem, aproximando-se dela com uma fase de −(n − m)× 90°. Para grau relativo 0 (mesmo número de pólos e zeros), o ponto final é um número finito não- zero complexo em ω → ♦.

O critério de estabilidade de Nyquist em detalhe

O critério de estabilidade de Nyquist é um método poderoso para determinar a estabilidade de circuito fechado a partir da resposta de frequência de circuito aberto. Ele alavanca Cauchy's article principle-plane é mapeado através de uma função F(s)] para um contorno fechado F F no ] F-plane (s) F(s) F(s) F) F[FLT] F[FLT]] é mapeado através de uma função [FLT:] F[FLT] (s) f(s) f), o número de vezes F) F[FLT] [F] (s[FLT]] [F] [F] é [FLT] [F] [F]] F] F [F] é [FLT(s] F] F] [F] [F] F] [F

Porque ]F(s)] está relacionado com a função de transferência de circuito aberto por F(s) = 1 + G(s)H(s), o mapa de F(s)[] através do contorno de Nyquist é simplesmente o plano Nyquist de [G(s)H(s)[] deslocado para a direita por 1. Equivalentemente, os encirclementos da origem no F-plane correspondem exatamente aos encirclementos do ponto −1 + 0j[][[G(s)H(s)-plane.

N = Z − P

em que:

  • N = número de circunscrições de −1 + 0j pelo gráfico de Nyquist (contado no sentido horário como positivo, anti-horário como negativo).
  • Z = número de pólos de circuito fechado no RHP (zeros de 1 + G(s)H(s)] no RHP).
  • P = número de pólos de circuito aberto no RHP.

Para um sistema de circuito fechado estável, nós exigimos ]Z = 0. Assim, o sistema de circuito fechado é estável se e somente se o gráfico de Nyquist de G(s)H(s) circundar o ponto −1 + 0j[[ exatamente −P[[] vezes (i.e., encirclementos no sentido anti-horário iguais P[j[[[][[[][] = 0 (estável em circuito aberto), então o gráfico de Nyquist não deve [][[[(FLT:13]]] encircle −1 em toda estabilidade de loop fechado.

Aplicando o Critério: Contagem de Encirculamento

Contar circunstancias rigorosamente envolve analisar o caminho de Nyquist sobre o contorno completo, incluindo o mapeamento do semicírculo infinito. Para sistemas com grau relativo ≥ 1, o semicírculo infinito mapeia para um pequeno círculo em torno da origem no plano G(s)H(s) (ou até um ponto fixo se grau relativo zero). A abordagem prática para engenheiros é:

  1. Desenhar o gráfico de Nyquist para ω de 0 para .
  2. Refletir através do eixo real para ω de 0 a − .
  3. Se o sistema tiver pólos no eixo imaginário (por exemplo, integradores), indentar o contorno de Nyquist com pequenos semicírculos à direita para evitar esses pólos, e mapear essas indentações para grandes arcos no plano G(s)H(s)].
  4. Contar o número de cercos no sentido horário de −1 traçando a curva fechada completa (incluindo o mapeamento infinito do semicírculo).

Matematicamente, a contagem de cercos N pode ser calculada através da integral do número de enrolamento:

N = \frac{1}{2πj} \oint {Nyquist \, contorno} \frac{G'(s) H(s)}{1+G(s) H(s)} \, ds

mas isso raramente é avaliado analiticamente; engenheiros dependem de ferramentas de inspeção geométrica ou software. O principal insight é que os cercos ocorrem quando a fase de G(jω)H(jω)[ cruza –180° (i.e., o ponto –1 no eixo real) enquanto a magnitude é ≥1. Esta condição define a margem de ganho [] e ]margem de fase.

Exemplos práticos de construção

Para solidificar a compreensão, consideramos um sistema de primeira ordem simples e um sistema de segunda ordem com amortecimento.

Sistema de primeira ordem: G(s) = 1/(s+1)

Aqui G(jω) = 1/(jω+1). A magnitude é [M(ω) = 1/ğ(1+ω2)[ e a fase é

Sistema de segunda ordem com adaptação: G(s) = ω n2/(s2 + 2ω n s + ω n2)[

Para um sistema de segunda ordem com frequência natural ωn e a relação de amortecimento

Para ver um crossover de fase adequado, considere um sistema com um pólo adicional em baixa frequência: G(s) = 1/(s+1)(s2 + 29,5%s + 1)].Aqui a fase pode cruzar-se -180° em algum finito ω, e o gráfico pode circundar -1 se o ganho for suficientemente alto.

Conectando Nyquist, Bode e Nichols Plots

O gráfico de Nyquist está intimamente relacionado com o gráfico de Bode (magnitude e fase) e o gráfico de Nichols. O gráfico de Bode fornece magnitude em dB e fase em graus como funções de frequência de log. O gráfico de Nyquist é uma combinação paramétrica da magnitude e fase: cada ponto na curva de Nyquist corresponde a uma magnitude e fase em um dado ω. A margem de ganho [] é o fator pelo qual o ganho pode ser aumentado antes da magnitude na frequência de crossover de fase se torna 1 (i.e., a distância de −1 para a curva de Nyquist ao longo do eixo real negativo). A margem de fase [ é o defasamento de fase adicional necessário na frequência de ganho crossover (onde magnitude = 1) para fazer a fase −180°. Estas margens são diretamente legíveis do gráfico de Nyquist: a margem de ganho é 1/ .

O gráfico de Nichols transforma o gráfico de Nyquist num gráfico cartesiano de magnitude (dB) vs. fase (deg). Esta transformação é útil para desenhar controladores, porque as margens de ganho e fase são facilmente vistas, e a resposta de circuito fechado pode ser lida usando os círculos M. Todas as três representações são equivalentes matematicamente; escolher uma depende do contexto e da preferência do engenheiro.

Limitações e Considerações Avançadas

Embora poderoso, o gráfico de Nyquist tem limitações que os engenheiros devem reconhecer. Primeiro, para sistemas com atrasos de tempo, a curva de Nyquist torna-se infinita e altamente curvada, exigindo um cuidadoso manuseio do mapeamento semicírculo infinito. O atraso adiciona uma mudança de fase −ωT que se acumula linearmente, fazendo com que o gráfico de Nyquist se torne espiral para dentro ou para fora, dependendo do atraso. A análise de estabilidade requer então o critério modificado de Nyquist] que responde pelo atraso usando o conceito do contorno de Nyquist que exclui o pólo de atraso no infinito.

Em segundo lugar, o gráfico de Nyquist, como descrito, é para sistemas de tempo contínuo. Para sistemas de tempo discreto (domínio z), o contorno de Nyquist torna-se o círculo unitário, e o critério de estabilidade examina circunscrições de −1 no plano w ou o uso da transformada bilinear para mapear o plano s. As fundações matemáticas permanecem semelhantes, mas as mudanças de contorno.

Terceiro, o critério Nyquist assume que o sistema de circuito fechado é estável se e somente se o gráfico de ciclo aberto Nyquist circundar −1 adequadamente. Contudo, isto requer que o carácter dos pólos de circuito aberto no RHP seja conhecido ] a priori. Se a função de transferência de circuito aberto tiver pólos no eixo imaginário (incluindo integradores puros), o contorno de Nyquist deve ser indentado para evitar estes pólos, e o mapeamento dessas indentações produz arcos infinitos que devem ser incluídos na contagem de encirculamento. Isto pode ser complicado para engenheiros novos ao método.

Finalmente, o gráfico de Nyquist é uma ferramenta de domínio de frequência; fornece condições necessárias e suficientes para estabilidade apenas para sistemas lineares invariantes do tempo (LTI). Para sistemas não lineares ou variantes do tempo, o critério de Nyquist não se aplica diretamente, embora extensões como o critério ] do círculo existem para certas classes de não linearidades.

Conclusão

A construção de um gráfico Nyquist está profundamente enraizada em análises complexas, particularmente no mapeamento do contorno Nyquist através da função de transferência de malha aberta e do princípio do argumento. Ao dominar as etapas matemáticas – avaliar ]G(jω)H(jω), plotar os pontos complexos e contar os círculos de −1 – os engenheiros ganham um método poderoso para avaliar a estabilidade de alça fechada e projetar controladores robustos.O critério Nyquist fornece intuitivo insight em margens de ganho e fase, e sua conexão com os gráficos Bode e Nichols torna-o uma ferramenta versátil no kit de ferramentas do engenheiro de controle.Compreender essas fundações permite aos engenheiros não só aplicar o gráfico Nyquist corretamente, mas também interpretar seus resultados com confiança em aplicações do mundo real, desde sistemas aeroespaciais até controle de processos industriais.

Para leitura posterior, consulte livros didáticos padrão como O artigo de Wikipedia sobre o critério de estabilidade de Nyquist, Tutoriais de Controle para MATLAB e Simulink[, ou o texto clássico Engenharia de Controle Moderno[] por Ogata. Esses recursos fornecem exemplos adicionais e derivações matemáticas mais profundas.