Usando Routh-hurwitz para Prever Resposta do Sistema e Comportamento Transiente

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

O critério Routh-Hurwitz é uma pedra angular da teoria clássica do controle, fornecendo um método sistemático para avaliar a estabilidade de sistemas lineares invariantes do tempo (LTI) sem computar explicitamente as raízes da equação característica.Desenvolvido independentemente por Edward John Routh em 1874 e Adolf Hurwitz em 1895, esta técnica algébrica permanece amplamente utilizada na engenharia de controle, pois oferece uma maneira rápida e confiável de determinar se um sistema irá responder de forma estável, previsível ou exibir oscilações e instabilidade.A estabilidade é fundamental para o design do sistema – um sistema instável pode levar a um super-respondência perigosa, oscilações persistentes ou mesmo falha catastrófica, especialmente em aplicações como aeroespacial, controles automotivos, processos químicos e robótica.

Ao analisar os coeficientes do polinômio característico de um sistema, o critério Routh-Hurwitz identifica o número de pólos (raiz da equação característica) que se encontram na metade direita do plano complexo. Pólos com partes reais positivas causam comportamento instável, enquanto pólos com partes reais negativas contribuem para respostas convergentes estáveis. O método não requer a resolução de polinômios de alto grau, tornando-o prático para análise manual e projeto de estágios iniciais. Além disso, suas percepções sobre comportamento transitório – como tempo de ajuste, sobreposição e frequência de oscilação – ajudam os engenheiros a ajustar controladores para atender as especificações de desempenho antes de se comprometerem com simulações detalhadas ou protótipos de hardware.

Este artigo irá primeiro explicar a base matemática do critério Routh-Hurwitz, e depois guiá-lo através da construção do array Routh com um exemplo passo a passo. Depois disso, vamos explorar como o critério prevê a resposta do sistema e o comportamento transitório, incluindo características-chave como sobreposição e tempo de ajuste. Um exemplo numérico prático ilustrará o processo, e discutiremos as limitações do método e as ferramentas complementares. Ao final, você terá uma compreensão sólida de como aplicar o critério Routh-Hurwitz como uma ferramenta prática e autoritária na análise do sistema de controle.

O que é o Critério Routh-Hurwitz?

O critério Routh-Hurwitz é uma condição necessária e suficiente para a estabilidade nos sistemas de LTI. Opera na equação característica do sistema, que é derivada do denominador da função de transferência de circuito fechado após a sua fixação em zero. Para um sistema descrito pela função de transferência , a equação característica é , onde é um polinômio em s de grau n:

O critério indica que para que o sistema seja estável (isto é, todos os pólos têm partes reais negativas), todos os coeficientes a_i devem ser não zero e ter o mesmo sinal (tipicamente positivo), e o array Routh construído a partir desses coeficientes não deve ter alterações de sinal na primeira coluna. Se ocorrerem alterações de sinal, o número de mudanças de sinal é igual ao número de pólos com partes reais positivas, indicando instabilidade.

É importante notar que o critério Routh-Hurwitz é puramente algébrico, não dá as localizações reais do pólo. Em vez disso, fornece uma contagem de pólos de meio plano direito (RHP). Isto torna inestimável quando o polinômio característico é de alta ordem (por exemplo, quarto, quinto ou superior) onde o achado de raiz é computacionalmente intensivo ou analiticamente inviável.

O critério está diretamente ligado à matriz Hurwitz, uma matriz quadrada formada a partir dos coeficientes do polinômio característico. Todos os principais menores da matriz Hurwitz devem ser positivos para estabilidade – esta é a condição Hurwitz. O método tabular de Routh simplifica este processo de verificação determinante para o array Routh mais simples.

Construindo o Array de Routh: Passo a passo

A construção do array Routh é um procedimento sistemático. Vamos ilustra-lo usando um exemplo concreto de polinomial:

Este polinômio tem raízes conhecidas em s = -1, -2, -3, -4 (todos negativos), então o sistema é estável. Nós vamos usá-lo para verificar que o array Routh não dá nenhum sinal de mudanças.

Passo 1: Escreva o Polinomial em Formulário Padrão

Garantir que o polinômio está organizado em poderes descendentes de s com todos os coeficientes escritos. Por exemplo, os coeficientes são:[
a₄ = 1 (s⁴, a_3[ = 10 (s³, a_2[ = 35 (s^]2[, a[_1[] = 50 (s^]1, a[_0[[[ = 24 (constant term).

Passo 2: Inicializar as primeiras duas linhas

A primeira linha do array Routh contém os coeficientes de termos indexados (começando a partir de s^]n): linha s[^4[: [1, 35, 24]. A segunda linha contém os coeficientes de termos indexados ímpares: linha s^3[]: [10, 50, 0]. (Preencher termos em falta com zeros.)

Assim:
Linha s⁴: 1, 35, 24
Linha s³: 10, 50, 0

Passo 3: Calcular linhas subsequentes usando a regra Rhombus

Cada elemento da linha seguinte é calculado usando as duas linhas acima dela. Para um elemento na linha s², primeira coluna, a fórmula é:

Mas é mais fácil aplicar a fórmula padrão: para a linha k, coluna j (começando da coluna 1), o elemento é:

Vamos fazer isso manualmente para clareza.

Row s²:
Primeiro elemento: Usando linhas s ^-4 e s³][
R(s²,1) = [(1×50) - (10×35)]][10 - 0(50]/10 = -[FLT][-300]10 = 30.[FLT: 12]]] Segundo elemento: próxima coluna (j=2). R(s 4,] [35, R(s ^{4[FLT] 4[FLT] 4[FT: 16] 10 = 2 [FLT: 24 = R(s[FT: 17] 3] [F[FT[FT: 18] 24[F[FT[F: 20] = FT}

Assim, a linha s²: 30, 24, 0.

Row s¹:
Usar linhas s³ e s[^2[.
]R(s^1[,1) = -[ (10×24) - (30×50)] / 30 = -[240 - 1500]/30 = -[-1260]/30 = 42.
] R(s[[¹2,2) = -[ (10×0) - (30×0)]] / 30 = 0.

Assim, a linha s¹: 42, 0.

Row s⁰:
Use linhas s² e s[^1[.
] R(s[^0[,1) = [(30×0) - (42×0)] / 42 = 0. Mas a última linha é sempre o termo constante do polinômio, que deve ser 24. Na verdade, existe uma regra conhecida: o último elemento da última linha é igual a[(FLT:12]0[ (o termo constante). Podemos calcular como: R(s⁰],1) = -[30×0) - (42×0) - (42×10) = 24.

Então, o conjunto completo de Routh é:

s4	1	35	24
s3	10	50	0
s2	30	24	0
s1	42	0	0
s0	24	0	0

Passo 4: Examine a primeira coluna

A primeira coluna contém: 1, 10, 30, 42, 24. Todos são positivos, e não há mudanças de sinal. Portanto, o sistema é estável, que corresponde às nossas raízes conhecidas.

Casos Especiais

Dois casos especiais podem surgir durante a construção do array Routh:

• Zero na primeira coluna: Se o primeiro elemento de uma linha for zero, mas a linha não for todos os zeros, substitua o zero por um pequeno épsilon (ε), então prossiga. O sinal de ε determina as alterações de sinais. Isto indica estabilidade marginal ou instabilidade, dependendo do resto do array.
• [[FLT: 0]] Fila inteira de zeros: Isto ocorre quando existem raízes simétricamente colocadas (por exemplo, pares de raízes imaginárias ou pares simétricos sobre o eixo real). Neste caso, forme um polinômio auxiliar da linha acima da linha zero e diferencie-a para obter coeficientes para a linha em falta. Isto indica a presença de um par de raízes puramente imaginárias ou de um par de raízes reais com sinais opostos (isto é, instabilidade).

Esses casos especiais estendem a aplicabilidade do critério a sistemas com pólos no eixo imaginário ou padrões simétricos.

Prevendo a Resposta do Sistema e Comportamento Transiente

O critério Routh-Hurwitz não se limita a um veredicto binário estável/instável. Ao revelar o número de pólos com partes reais positivas, ele informa diretamente o comportamento transitório do sistema. Cada pólo no meio plano direito contribui para um modo exponencialmente crescente, fazendo com que a resposta transitória diverja. Por outro lado, todos os pólos no meio plano esquerdo levam a uma deterioração exponencial, atingindo eventualmente um estado estável.

No entanto, o critério não dá as localizações exatas do pólo, portanto não pode quantificar precisamente o tempo de superação, de fixação ou de elevação. No entanto, fornece pistas essenciais:

• Se a primeira coluna não tiver alterações de sinal: Todos os pólos têm partes reais negativas. A resposta transitória será estável, eventualmente se fixando em um estado estacionário. Mas pode haver decadência lenta ou oscilatória (pólos quase imaginários) que o critério não detecta – apenas o sinal da parte real é garantido.
• Se houver alterações de sinais: O número de alterações de sinais é igual ao número de pólos RHP. O sistema é instável e a resposta transitória irá crescer sem limites. O número de pólos RHP indica a natureza da instabilidade: um único pólo positivo leva a um crescimento exponencial monotónico, enquanto um par de pólos RHP complexos produz oscilações crescentes.
• Estabilidade marginal: Se a equação característica tiver raízes no eixo imaginário (puramente imaginário ou múltiplos pólos no eixo), o array Routh produzirá uma linha de zeros. Isto indica que o sistema está no limite da estabilidade, e a resposta transitória pode apresentar oscilações sustentadas (se simples pólos imaginários) ou oscilações rampa-up (se pólos imaginários repetidos).

Ao combinar o resultado de Routh-Hurwitz com o conhecimento da ordem do sistema, os engenheiros podem fazer previsões informadas. Para um sistema de segunda ordem, o critério reduz-se a exigir que todos os coeficientes sejam positivos (que é necessário e suficiente para a estabilidade). Para sistemas de ordem superior, o array fornece um indicador de margem de estabilidade – quanto mais próximo os elementos da primeira coluna estão de zero, mais próximo o sistema é da instabilidade. Isto é especialmente útil quando ajustar parâmetros de ganho em um loop de controle.

Relação com Características Transientes

A resposta transitória de um sistema estável é determinada principalmente pelos pólos dominantes, os mais próximos do eixo imaginário. O critério Routh- Hurwitz não localiza diretamente estes pólos, mas pode ajudar na afinação. Por exemplo, se o critério mostrar que aumentar um ganho de controlador irá eventualmente levar a uma mudança de sinal na primeira coluna, ele marca a margem de ganho. Os engenheiros usam isto frequentemente para definir ganhos proporcionais em controladores P, PI ou PID sem precisar simular o sistema inteiro a cada vez.

As principais métricas transitórias que são influenciadas pela estabilidade incluem:

• Sobreposição – O valor máximo de pico da resposta em relação ao valor final. Superação é tipicamente causada por pólos conjugados complexos com uma baixa relação de amortecimento. Embora o critério Routh-Hurwitz não compute amortecimento, um sistema com uma linha de zeros (estabilidade marginal) muitas vezes corresponde a uma relação de amortecimento perto de zero, implicando alto sobreposição.
• Tempo de ajuste – O tempo necessário para que a resposta permaneça dentro de uma determinada percentagem (por exemplo, 2%) do valor final. O tempo de ajuste é inversamente proporcional à parte real do pólo dominante. O array Routh não pode dar esta parte real, mas se o array mostrar que todos os pólos estão no meio plano esquerdo, o tempo de fixação é finito; se algum pólo estiver muito próximo do eixo imaginário, o tempo de fixação será longo.
• Rise Time – O tempo para ir de 10% para 90% do valor final. Isto está relacionado com a frequência natural dos pólos complexos dominantes. Novamente, o critério Routh-Hurwitz não fornece nenhum valor numérico direto, mas indica se o sistema é estável o suficiente para ter um tempo de ascensão bem definido.
• Oscilações – A frequência de oscilações transitórias é determinada pelas partes imaginárias dos pólos complexos.O critério não revela partes imaginárias, mas a presença de uma linha de zeros indica pólos puramente imaginários, que produzem oscilações sustentadas na resposta transitória.

Assim, o critério Routh-Hurwitz atua como um gatekeeper: garante estabilidade como um pré-requisito para análise significativa transiente usando outras ferramentas como locus de raiz, gráficos Bode, ou simulações de resposta de passo.

Exemplo prático: Sistema de Terceira Ordem

Vamos aplicar o critério Routh-Hurwitz a um exemplo mais complexo com parâmetros desconhecidos. Considere a função de transferência de circuito fechado com o feedback de unidade. A equação característica é:

Queremos determinar o intervalo de K para o qual o sistema é estável. Construir o array Routh para polinomial .

Linhas:[
s³: 1, 10[
s[²: 7, K
s[^1[: ^[
s[0[]: Na verdade, mais simples: o elemento da última linha é o termo constante K (desdevido que a linha acima não é todos os zeros). Mais sistematicamente: a linha s^]0[[][[] é derivada das linhas s[^2]2[2][]][[[[FLT16]]]]]][[[[[[[[[FLT16

s3	1	10
s2	7	K
s1	(K-70)/7	0
s0	K	0

Para a estabilidade, todos os elementos da primeira coluna devem ser positivos: 1 > 0, 7 > 0, (K- 70)/7 > 0, K > 0. Isto dá K > 70 e K > 0, por isso K > 70. O limite inferior é K > 70. Se K = 70, a linha s[[ FLT: 0]] 1[[ FLT: 1]] torna- se zero, o que é um caso especial que indica estabilidade marginal (raízes puramente imaginárias). Assim, o sistema é estável para K > 70. Isto diz ao engenheiro que ganhos acima de 70 produzem um transiente estável, mas à medida que K se aproxima de 70, as oscilações transitórias tornam- se inalteradas. A resposta transitória para K = 200, por exemplo, terá dois pólos conjugados complexos no meio plano esquerdo e um pólo real, levando a uma resposta de sobreposição, mas eventualmente de liquidação. O critério Routh- Hurwitzz prediz explicitamente a fronteira de estabilidade.

Implicações de resposta transientes para o exemplo

Para K pouco acima de 70, o sistema é estável, mas os pólos dominantes estão muito próximos do eixo imaginário. A resposta transitória mostrará uma deterioração muito lenta e uma grande superação. O tempo de fixação será longo. À medida que o K aumenta mais, os pólos movem- se para a esquerda (geralmente), melhorando o tempo de amortecimento e de fixação. Contudo, se o K for demasiado elevado, a matriz Routh para este sistema em particular não mostrará quaisquer alterações adicionais de sinais — a estabilidade persiste. Na prática, o K muito elevado poderá causar outros problemas, como saturação do atuador ou amplificação de ruído, mas o critério Routh- Hurwitz só verifica a estabilidade do ciclo.

Assim, o critério fornece informações críticas: o intervalo de ganho para estabilidade e a proximidade com instabilidade.Ao avaliar os elementos da primeira coluna, os engenheiros podem calcular margens de ganho e controladores de projeto que evitam estabilidade marginal.

Limitações do Critério Routh-Hurwitz

Apesar do seu poder, o critério Routh-Hurwitz tem limitações importantes que cada engenheiro de controlo deve compreender:

• Aplicável apenas aos sistemas LTI: Não pode ser utilizado diretamente para sistemas de dados variáveis, não lineares ou amostrados (embora existam extensões para sistemas de tempo discreto através de transformadas bilineares).
• Não quantifica o desempenho transitório: O critério não fornece informações sobre o excesso de tempo, o tempo de pico, o tempo de subida ou o tempo exato de fixação. Para isso, são necessárias locus de raiz, resposta de frequência ou simulações de domínio temporal.
• Requer um polinômio característico conhecido: Os coeficientes devem ser conhecidos exatamente. Na prática, incertezas em parâmetros podem exigir uma análise robusta da estabilidade (por exemplo, usando o teorema de Kharitonov).
• Torna-se complicado para sistemas de alta ordem: Cálculos manuais para polinômios de décima ordem são propensas a erros, embora software (MATLAB, Python) automatize a construção de array.
• Não é possível lidar com atrasos de tempo diretamente: Os sistemas com atrasos de tempo têm equações características transcendentais – o array Routh não pode ser diretamente aplicado a menos que use aproximações de Padé.

Para superar essas limitações, os engenheiros usam comumente o critério Routh-Hurwitz em conjunto com outros métodos. Por exemplo, gráficos de locus de raiz mostram como pólos de loop fechado se movem com ganho, fornecendo informações transitórias exatas. O critério Nyquist lida com atrasos de tempo e não linearidades. No entanto, o array Routh muitas vezes serve como uma primeira verificação rápida, especialmente durante o projeto preliminar ou solução de problemas.

Ferramentas Complementares

O critério de Routh-Hurwitz é “necessário e suficiente” para a estabilidade dos polinômios de LTI, mas para a análise de desempenho, os engenheiros dependem:

• Locus de Root:] A trajetória dos pólos como parâmetro (por exemplo, ganho) varia. Mostra diretamente a relação amortecimento e a frequência natural.
• Bode Gráficos e Diagramas Nyquist: Fornecer margem de ganho e margem de fase, que são indicadores de estabilidade relativa e características de resposta transientes.
• Simulação de Resposta de Passo: Dá as métricas transientes reais (superação, tempo de ajuste) para um determinado sistema.
• Métodos de Estado-Espaço:] Os autovalores da matriz de estado (que são os pólos) podem ser calculados numericamente, e o método de Lyapunov estende a análise de estabilidade a sistemas não lineares.

O critério Routh-Hurwitz ocupa assim um ponto doce histórico e prático: é matemático o suficiente para ser rigoroso, mas simples o suficiente para ser feito em um quadro branco durante uma revisão de design. Continua a ser uma habilidade essencial para qualquer engenheiro de controle.

Conclusão

O critério Routh-Hurwitz oferece uma forma clara e algébrica de avaliar a estabilidade do sistema, analisando o polinômio característico sem resolver as raízes. Através da construção do arranjo Routh, os engenheiros podem determinar o número de pólos instáveis, entender a sensibilidade do sistema às mudanças de parâmetros e identificar limites de estabilidade que afetam diretamente o comportamento transitório. Embora não forneça métricas transientes exatas como overshoot ou o tempo de fixação, ele serve como uma ferramenta de triagem fundamental que informa uma análise mais profunda usando locus raiz, resposta de frequência ou simulação.

O domínio do critério Routh-Hurwitz equipa-o com a capacidade de avaliar rapidamente se um sistema de controlo se irá comportar de forma fiável em condições transitórias. Quer esteja a desenhar um controlador proporcional simples ou a ajustar um sistema multi- loop complexo, o array Routh é um ponto de partida fiável. Para mais leitura, considere explorar o critério de estabilidade Routh–Hurwitz na Wikipedia[ para profundidade matemática adicional, ou consulte a Universidade de Michigan Control Tutoriais em Routh Array para exemplos mais trabalhados. Adicionalmente, MathWorks documentação[[] fornece insights em implementações de software. Ao integrar o critério Routh-Hurwitz no seu fluxo de design de controlo, você garante que os sistemas que constrói não são apenas estáveis, mas também exibem respostas transitórias previsíveis e de alto desempenho.