Introdução ao Critério de Estabilidade Routh-Hurwitz

O critério de estabilidade de Routh-Hurwitz é uma pedra angular da teoria de controle e análise polinomial. Ele fornece um método sistemático para determinar o número de raízes de um polinômio real que se encontram na metade direita do plano complexo - sem exigir computação de raiz explícita. Para engenheiros e matemáticos que trabalham com polinômios complexos, este critério transforma um problema algébrico de outra forma difícil em um procedimento tabular simples.

Um “polinômio complexo” aqui se refere a um polinômio com coeficientes reais (como é padrão em aplicações de controle), mas cujas raízes podem ser complexas. A estabilidade de muitos sistemas físicos – de circuitos eletrônicos a ligações mecânicas – depende do sinal das partes reais das raízes características do polinômio. O critério Routh-Hurwitz responde a essa pergunta de forma eficiente, mesmo para polinômios de alta ordem que são impraticáveis para fatorar.

Este artigo amplia o resumo original, fornecendo uma análise detalhada do método, exemplos, casos especiais e aplicações práticas. No final, os leitores poderão aplicar a técnica de Routh-Hurwitz para simplificar e analisar polinômios complexos com confiança.

O que são polinômios complexos?

Em engenharia de controle e matemática, um polinômio complexo é tipicamente uma expressão da forma:

P(s) = ansn[ + an−1[s[n−1[ + ... + a1[]s + a0][

onde os coeficientes ai são números reais e s[] é uma variável complexa. O termo “complexo” não se refere aos coeficientes (que permanecem reais), mas ao fato de que as raízes da equação P(s)=0[] são geralmente números complexos. Essas raízes determinam o comportamento dinâmico de um sistema descrito pelo polinômio.

Por exemplo, a equação característica de um sistema de segunda ordem como um mass-spring-damper é ]ms[2[ + bs + k = 0. Dependendo da relação amortecimento, as raízes podem ser reais e negativas (sobrepujadas) ou um par conjugado complexo com parte real negativa (submersa). O critério Routh-Hurwitz fornece uma maneira de verificar se todas as raízes têm partes reais negativas sem resolver a equação.

Por que o cálculo direto da raiz é impraticável

Para polinômios da ordem cinco ou mais, raramente existem soluções de raiz analítica exatas (teorema de Abel- Ruffini). Os métodos numéricos podem aproximar as raízes, mas são sensíveis à incerteza do coeficiente e podem perder sinais sutis de instabilidade. O critério de Routh-Hurwitz, por contraste, dá uma resposta definitiva sim/não sobre o sinal de partes reais, usando apenas os coeficientes polinomiais. Isto torna indispensável para análise de estabilidade preliminar e para verificação de margens de projeto.

O Critério de Estabilidade Routh–Hurwitz

Desenvolvido independentemente por Edward John Routh (1876) e Adolf Hurwitz (1895), o critério afirma que para um polinômio com coeficientes reais, todas as raízes têm partes reais negativas se e somente se todos os principais menores da matriz Hurwitz (ou, equivalentemente, a primeira coluna do array Routh) forem positivos. O método é geralmente ensinado usando o array Routh porque é mais simples construir à mão.

Construindo o Array de Routh

Dado o polinômio:

P(s) = ansn[ + an−1s[n−1[ + ... + a1[s + a[]0[]]] (]an[ > 0]),

O conjunto Routh é construído da seguinte forma:

  1. ]Label the first two rows.] Linha s[n contém coeficientes de poderes equilibrados (começando com a[n[, ]an−2[an−1n[, ]n−2[FLT[FT:18]a]n[FLT[Flt][Flt]]]]]n:2]]]s.
  2. Computar linhas subsequentes. Para a linha sn−2[, cada elemento é derivado das duas linhas acima usando uma fórmula determinante:
    ]b[n] = (a]n−1 /n]] - a[FLT:n:(FLT:20] =
    n ] [FLT [FLT[FT[:]] [FT:3]] T[FT:2T[FT:2T[FT:2T[:2T:2]]]]]]]]]
  3. Continua até que apenas a linha constante permaneça. A última linha corresponde a s[0 e contém um único elemento, que é o termo constante (normalmente a0]).
  4. Examine a primeira coluna. Se todas as entradas na primeira coluna do array são positivas (e nenhuma é zero), então todas as raízes estão no semiplano esquerdo (o sistema é estável). Qualquer mudança de sinal indica o número de raízes com partes reais positivas – cada mudança de sinal conta uma dessas raízes.

Casos Especiais na Linha de Routh

Surgem duas situações especiais que requerem tratamento cuidadoso:

  • Zero na primeira coluna. Se o primeiro elemento de uma linha for zero, mas o resto da linha não for zero, substitua o zero por um pequeno número positivo ε e prossiga. O sinal de ε determina o resultado da estabilidade. Alternativamente, pode-se usar o método auxiliar [] (ver abaixo).
  • Fila inteira de zeros. Isto indica que o polinômio contém raízes simétricas localizadas (pares de raízes que são negativas umas das outras ou conjugados complexos com parte real zero). Neste caso, formam um polinômio auxiliar da linha acima dos zeros, pegam a sua derivada e usam os coeficientes dessa derivada para substituir a linha zero. Em seguida, continuem o array.

Exemplo: Um Polinômio de Quarta Ordem

Considerar P(s) = s4 + 2s[3[ + 3s2[ + 4s + 5[].

  • Linha s4: 1, 3, 5
  • Linha s3: 2, 4, 0
  • Linha s2: b1 = (2·3 - 1·4)/2 = (6-4)/2 = 1, ]b[2[] = (2·5 - 1·0)/2 = 10/2 = 5]]
  • Linha s1: c[1 = (1·4 - 2·5)/1 = (4-10)/1 = -6, então zero
  • Linha s0: d[1] = (-6·5 - 1·0)/(-6) = (-30)/(-6) = 5[

Primeira coluna: 1, 2, 1, -6, 5 → o sinal muda uma vez (de 1 para -6, e de -6 para trás para 5), indicando uma raiz com parte real positiva. De fato, o polinômio tem raízes: -1,65 ± 0,35i, 0,65 ± 1,72i — o último par tem parte real positiva.

Simplificar Polinomiais Complexos Usando o Critério

O método de Routh-Hurwitz simplifica a análise de estabilidade de polinômios complexos de várias maneiras:

  1. Avoids root computation. Para polinômios de alta ordem, fatorar ou procurar numericamente por raízes é caro e às vezes não confiável.O array requer apenas adição, subtração, multiplicação e divisão - operações facilmente executadas à mão ou planilha.
  2. Certeza dos parâmetros de mãos. Ao expressar coeficientes em termos de parâmetros do sistema (ganhos, constantes de tempo), pode-se derivar condições de estabilidade em termos de desigualdades. Por exemplo, o intervalo de um ganho de controlador K[ para o qual o sistema de circuito fechado permanece estável pode ser encontrado aplicando o array Routh e exigindo que todos os elementos de primeira coluna sejam positivos.
  3. Identifica a estabilidade marginal. Um zero na primeira coluna (o ε caso) ou uma linha de zeros indica raízes no eixo imaginário ou simétricomente colocadas sobre a origem. Estas condições definem frequentemente o limite entre comportamento estável e instável.

Trabalhar com Coeficientes Complexos?

Embora o critério clássico de Routh-Hurwitz exija coeficientes reais, existem muitas extensões. Para polinômios com coeficientes complexos (ai , modificações tais como o ] modificado Routh array ou o [Critério de Michailov[]] podem ser empregadas. No entanto, na maioria dos problemas práticos de engenharia de controle, a equação característica surge de parâmetros físicos reais, então os coeficientes reais são a norma.

Aplicações em Sistemas de Controle e Além

A aplicação primária do critério de estabilidade Routh-Hurwitz é determinar a estabilidade dos sistemas linear-invariantes (LTI) descritos por funções de transferência. O denominador da função de transferência é o polinômio característico, e a estabilidade requer que todas as raízes tenham partes reais negativas.

Concepção de sistemas de controle

Durante o desenho do controlador, os engenheiros frequentemente ajustam ganhos (por exemplo, os ganhos proporcionais, integrais e derivados de um controlador PID) para alcançar o desempenho desejado. O critério Routh fornece uma maneira rápida de encontrar ] margens de ganho[—a gama de ganhos que mantêm o sistema estável. Por exemplo, dado um sistema de feedback de unidade com função de transferência de frente G(s) = K / [(s+1)(s+2)(s+3)],, o polinomial característico é s[[3[[ + 6s[[2[ + 11s + (6+K) = 0. Aplicando o arranjo de Routh produz a condição K [)]] [[F:11]] para estabilidade.

Análise da Rede

Na engenharia elétrica, o polinômio característico de um circuito RLC determina se as oscilações transientes decaem. O critério Routh-Hurwitz ajuda os engenheiros a escolherem valores de componentes para evitar oscilações sustentadas (instabilidade). Da mesma forma, em sistemas mecânicos, a estabilidade de robôs ou aeronaves controlados por feedback é rotineiramente avaliada usando este método.

Controle de Processo Industrial

Nas indústrias química e de petróleo, os controladores de processo regulam a temperatura, pressão e vazão. A matriz Routh é usada durante a afinação de controladores proporcional-integral-derivativos (PID) para garantir que a planta de malha fechada permaneça estável sob variações em condições operacionais.

Benefícios e Limitações

Como qualquer ferramenta, o critério Routh-Hurwitz tem pontos fortes e fracos.

Pontos fortes

  • Não é necessário cálculo de raiz. O array opera exclusivamente em coeficientes.
  • Scalable. O método funciona para polinômios de qualquer ordem, embora ordens mais altas se tornem tediosas à mão (fácilmente automatizadas).
  • Capacidade simbólica. Os coeficientes podem ser expressões algébricas, permitindo análise de estabilidade paramétrica.
  • Fornece o número de raízes instáveis. Cada mudança de sinal na primeira coluna corresponde a uma raiz com parte real positiva.

Limitações

  • Coeficientes reais apenas. A forma clássica não lida diretamente com polinômios com coeficientes complexos, embora existam extensões.
  • Nenhuma informação sobre as localizações das raízes. O critério apenas diz quantas raízes instáveis existem, não os seus valores exatos. Para análise transitória detalhada, é necessário um lócus raiz ou um radiculamento numérico.
  • Sensibilidade a erros numéricos. Quando os coeficientes são quase zero ou contêm erros de arredondamento, a construção do array pode produzir ambiguidades de sinais. O uso do método ε é uma solução alternativa, mas requer interpretação cuidadosa.
  • Não lida com sistemas de atraso. Os sistemas com atrasos de tempo têm equações características transcendentais, não polinômios. O critério Routh–Hurwitz não é diretamente aplicável; em vez disso, métodos como aproximações de Padé ou o preditor de Smith são usados.

Contexto Histórico e Leitura Adicional

O critério Routh-Hurwitz é nomeado em homenagem a Edward John Routh, um matemático inglês que publicou o método de matriz em 1876, e Adolf Hurwitz, um matemático alemão que independentemente derivava uma formulação de matriz semelhante em 1895. Seu trabalho construído sobre contribuições anteriores por Cauchy e Hermite. O critério permanece um básico nos livros didáticos de controle modernos - veja o artigo da Wikipédia para uma visão geral. Para um tratamento matemático mais profundo, A página de tópico da ScienceDirect fornece referências às fontes originais.

Engenheiros que estudam controle avançado também encontrarão a matriz de Hurwitz no contexto da estabilidade de Lyapunov. A conexão entre o array Routh e o determinante da matriz de Hurwitz é um resultado clássico na álgebra linear e teoria de controle.

Exemplo ilustrado passo a passo

Vamos percorrer um exemplo completo com um polinômio complexo (de alta ordem) para demonstrar o poder do método. Considere o polinômio característico de um sistema de controle:

P(s) = 2 s5 + 4 s[4[ + 6 s3 + 8 s[]2[[ + 10 s + 12]

Passo 1: Organize as duas primeiras linhas

  • Linha s5 (potência 5): coeficientes de s5[, ]s[3[][, ]]]s[1[ → 2, 6, 10
  • Linha s4 (potência 4): coeficientes de s4, ]s[2[][, ]]s[0[] → 4, 8, 12

Passo 2: Calcular a terceira linha (]s3)

Utilizando a fórmula determinante:

  • Primeiro elemento: (4·6 - 2·8)/4 = (24 - 16)/4 = 8/4 = 2
  • Segundo elemento: (4·10 - 2·12)/4 = (40 - 24)/4 = 16/4 = 4
  • Terceiro elemento: zero (desde que não haja mais coeficientes na primeira linha)

Linha s3: 2, 4, 0

Passo 3: Calcular a quarta linha (]s2)

Usando a linha s4[(4,8,12] e a linha ]s3[ (2,4,0):

  • Primeiro elemento: (2·8 - 4·4)/2 = (16 - 16)/2 = 0/2 = 0
  • Segundo elemento: (2·12 - 4·0)/2 = (24 - 0)/2 = 12

Linha s2[: 0, 12

Aqui temos um zero na primeira coluna. O resto da linha é não-zero (12). Substitua o zero por ε (um pequeno número positivo).

Passo 4: Calcular a quinta linha (]]1)

Usando a linha s3[2,4,0]] e a linha s2[](ε,12):

  • Primeiro elemento: (ε·4 - 2·12)/ε = (4ε - 24)/ε = 4 - 24/ε
  • Segundo elemento: (ε·0 - 2·0)/ε = 0

Linha s1: (4 - 24/ε), 0

Como ε → 0[+, o termo -24/ε domina, tornando este elemento negativo (− .).

Passo 5: Calcular a sexta linha (]s0)

Usando a linha s2 (ε,12] e a linha ]s1[[ (4 - 24/ε, 0):

  • Primeiro elemento: [(4 - 24/ε)·12 - ε·0] / (4 - 24/ε) = 12

Linha s0: 12

Passo 6: Analise a primeira coluna

As entradas da primeira coluna são: linha s5: 2 (positivo), linha s4[: 4 (positivo), linha s3: 2 (positivo), linha s[2[[: ε (positivo, como ε>0), linha s[[1: 4 - 24/ε (negativo para pequeno ε), linha s[0[: 12 (positivo).

Existem duas alterações de sinais: de positivo (linha s2] para negativo (linha s1), e de negativo para positivo (linha s0]). Portanto, o polinômio tem duas raízes com partes reais positivas[].

Verificação: Usando ferramentas numéricas, as raízes de 2 s5[ + ... + 12 = 0 são aproximadamente: -1,09, -0,50 ± 0,87i, 0,79 ± 1,04i. De fato, as duas raízes conjugadas complexas têm parte real positiva (0,79), confirmando a predição de Routh-Hurwitz.

Abordagem alternativa: o polinômio auxiliar para o caso de linha zero

Se tivéssemos encontrado uma linha completa de zeros em vez de apenas um zero na primeira coluna, formaríamos um polinômio auxiliar da linha acima. Essa técnica é essencial para casos como s[4[ + 5s[2[ + 4[[] (que tem raízes puramente imaginárias ±j, ±2j). A linha zero indica raízes simétricas sobre a origem, e o polinômio auxiliar (derivado da linha acima) dá as raízes reais.

Conclusão

O critério de estabilidade Routh-Hurwitz transforma a análise complexa de polinômios de alta ordem em um algoritmo simples. Ao construir a matriz Routh e inspecionar os sinais da primeira coluna, engenheiros e matemáticos podem determinar a estabilidade do sistema sem resolver as raízes. O método é indispensável no projeto do sistema de controle, análise de circuito e qualquer campo onde as raízes de um comportamento polinômio real ditam.

Dominar esta técnica permite lidar com polinômios complexos com confiança, seja para estudos acadêmicos, design industrial ou pesquisa. Combinado com ferramentas computacionais modernas (por exemplo, MATLAB, Python), o critério Routh-Hurwitz continua a ser um bloco fundamental da análise dinâmica do sistema. Para aqueles que desejam explorar mais, o Control Systems Principles whitepaper[ fornece exemplos adicionais trabalhados e extensões avançadas.