O critério de estabilidade de Routh-Hurwitz é uma pedra angular da teoria clássica de controle, fornecendo um método algébrico direto para determinar se um sistema linear invariante de tempo é estável sem resolver seus pólos. Para sistemas de ordem inteira, onde o polinômio característico tem apenas potências inteiras da variável Laplace s- o teste é tanto eficiente quanto elegante. No entanto, um número crescente de aplicações do mundo real, desde o amortecimento viscoelástico até a espectroscopia de impedância eletroquímica, são descritos por equações diferenciais de ordem fracionária. Estes sistemas introduzem derivados de ordem não-inteira, e suas equações características contêm poderes fracionais de s, tornando o tradicional conjunto Routh impossível de construir diretamente. Engenheiros e pesquisadores devem, portanto, adaptar a filosofia subjacente da análise Routh-Hurwitz para lidar com dinâmicas de ordem fracionária. Este artigo explica como estender essas ideias, apresentando métodos práticos que preservam o critério original enquanto que a simplicidade dos sistemas fracionais mais ricos.

Compreendendo o clássico Routh–Hurwitz Criterion

Antes de explorar extensões, é essencial recapitular o método clássico. Para um sistema com um polinômio característico

o critério Routh-Hurwitz constrói uma tabela (o array Routh) a partir dos coeficientes. A condição de estabilidade é que cada elemento da primeira coluna do array deve ter o mesmo sinal (geralmente positivo). Se ocorrerem alterações de sinais, o número de mudanças de sinais é igual ao número de pólos de meio plano à direita. A beleza do teste reside na sua natureza puramente algébrica: sem descoberta de raiz, sem análise complexa. Para sistemas de ordem inteira, isto funciona perfeitamente porque os polinômios têm um número finito de raízes, e a condição de meio plano à esquerda é inequívoca.

O critério também se estende aos sistemas com atrasos de tempo (através da aproximação Padé) e aos sistemas de tempo discreto (usando a transformada bilinear), mas todas essas extensões dependem do fracionamento polinomial. Os sistemas de ordem fracionária quebram esse framework, exigindo uma nova perspectiva.

Sistemas de ordem fraccional: uma breve visão geral

Um sistema de ordem fracionária é aquele cuja dinâmica é descrita por uma equação diferencial fracionária da forma

onde as ordens de diferenciação αi e βi não são necessariamente inteiros. A transformada de Laplace de tal equação (em condições iniciais zero) produz uma função de transferência contendo termos como s[γ[] com fraccional γ. Por exemplo, um simples integrador de ordem fracionária tem a função de transferência 1/s[[α] com 0 < α < 2. Estes sistemas exibem propriedades não vistas em sistemas de ordem inteira: representação de estado-espaço de dimensão infinita, relaxamento não-exponencial e memória de longo alcance.

Cálculo fracionário encontrou aplicações no controle de estruturas flexíveis, modelagem de tecidos biológicos, caracterização de impedância de bateria e até mesmo séries temporais financeiras. A necessidade de verificar a estabilidade para tais sistemas é aguda, mas a tabela clássica Routh-Hurwitz não pode ser construída porque a equação característica não é um polinômio.

Por que o teste clássico Routh-Hurwitz falha para sistemas fracionais

O obstáculo é fundamental: uma equação característica de ordem fracionária, como

não é um polinômio. Os poderes de s não são múltiplos inteiros um do outro. O conjunto Routh requer poderes inteiros para que uma eliminação sistemática possa ser realizada. Além disso, os poderes fracionários levam a cortes de ramificações no plano complexo; o conceito de “pole de meio plano direito” torna-se mais sutil porque um poder fracionário pode ter múltiplas folhas de Riemann. Um sistema fracionário pode ter um número infinito de pólos (se a ordem for irracional) ou um número finito se as ordens forem compatíveis.

Como o critério Routh-Hurwitz é essencialmente um teste algébrico para a propriedade Hurwitz de um polinômio, não pode ser aplicado diretamente a uma expressão não polinomial. Ao invés disso, são necessárias modificações que preservem o espírito do teste — determinando estabilidade examinando alguma condição algébrica transformada — enquanto se trata da natureza fracionária.

O critério de estabilidade Matignon: uma substituição algébrica direta

A extensão mais conhecida é o critério de estabilidade Matignon (também chamado de teorema de Matignon). Fornece uma condição necessária e suficiente para uma grande classe de sistemas de ordem fracionária com ordens proporcionais. Um sistema fracionário é dito ter ordens proporcionais se todos os poderes fracionários na equação característica são múltiplos inteiros de uma ordem de base α (por exemplo, α, 2α, 3α, ...). Muitos sistemas práticos, especialmente aqueles derivados de controladores PID fracionários ou de modelos viscoelásticos fracionários, são compatíveis.

Critério de Mattignon:Para um sistema de ordem fraccionada com equação característica

onde α é um número real positivo, o sistema é estável se e somente se

..arg(si]] .. > απ/2[] para cada raiz si] de P[.

Em outras palavras, as raízes devem estar em uma região do plano complexo fora de um setor de ângulo απ/2 centrado no eixo real negativo. Esta condição reduz para a condição clássica de meio plano esquerdo quando α = 1 (desde então απ/2 = π/2), mas para α < 1 a região de estabilidade é uma cunha.

O critério Matignon é poderoso porque transforma o problema de estabilidade de volta em um problema de localização-raiz para o polinomial transformado Pw = sα. Pode-se aplicar o teste de Routh-Hurwitz ao polinomial P(w)[[] após a substituição ]w = sα[[, mas a condição do argumento de s]] não é capturada apenas pelo conjunto clássico de Routh. No entanto, para sistemas de commensuração, a condição de estabilidade pode ser verificada pela primeira construção do arranjo Routh para o polinomial ][FTP][F][F] e então verificar a condição.

Isto conduz a um procedimento prático:

  1. Identifique a ordem fraccionada base α de modo que a equação característica se torne um polinômio em sα].
  2. Definir w = sα e reescrever a equação característica como P(w) = 0], um polinômio convencional.
  3. Aplicar o teste clássico de Routh–Hurwitz para P(w). Se ocorrerem alterações de sinais na primeira coluna, o sistema é instável sob a definição polinomial clássica em w[.
  4. Verifique a condição do argumento:] para cada raiz wi[, calcule o correspondente s[i]i[w[i[1/α[][] e verifique se .arg(s[i]) . > απ/2. Este passo pode exigir computação numérica porque as raízes wi][[[[]]] .

O critério Matignon é frequentemente utilizado como a principal ferramenta para a análise de estabilidade de sistemas fraccionados proporcionais. Para sistemas não-compensados (onde as ordens fraccionadas não são múltiplas de uma base comum), devem ser utilizados métodos alternativos, como a abordagem de domínio de frequência ou a inversão numérica.

Abordagens analíticas alternativas

Função Mittag–Leffler e Representação Estatal-Espaço

Os sistemas de ordem fraccional podem ser representados de forma pseudo-estatal-espacial utilizando o derivado fraccional da ordem α. A condição de estabilidade torna-se uma restrição nos autovalores da matriz do sistema. Especificamente, para um sistema fraccional linear descrito por

o sistema é estável se e somente se

.arg(λ(A)) . . .π/2

onde λ(A) são os autovalores da matriz A. Este é o análogo multivariável direto do critério Matignon. Para sistemas compatíveis, esta condição reduz-se para a mesma região de estabilidade em forma de cunha. Os engenheiros podem calcular os autovalores e verificar seus argumentos usando álgebra linear numérica padrão.

Transformação para um sistema de ordem integral através da substituição variável

Outro caminho analítico é introduzir uma nova variável ]z = sα[ e, em seguida, aplicar o teste clássico de Routh–Hurwitz ao polinômio transformado. Contudo, esta abordagem só verifica a “parte polinomial” da estabilidade. Uma condição necessária é que o polinômio em z[ tem todas as raízes no semiplano esquerdo, mas a condição suficiente requer o teste de argumento adicional. Na prática, os engenheiros frequentemente executam ambos: primeiro execute o arranjo Routh em P(z)[[] para detectar instabilidade polinomial óbvia; então, para os candidatos estáveis restantes, aplique a condição de argumento.

Técnicas de Domínio de Frequência

A estabilidade pode também ser avaliada pelos critérios Nyquist ou Bode adaptados às funções de transferência de ordem fraccionada. A função de transferência de circuito aberto G(s)[] de um sistema de ordem fraccionada é uma função não-racional, mas a sua resposta de frequência pode ser avaliada numericamente. O critério de estabilidade de Nyquist ainda contém: o número de circunscrições do ponto −1 determina a estabilidade de circuito fechado. No entanto, a aplicação direta do método Routh-Hurwitz é frequentemente preferida porque evita a plotagem e proporciona uma condição algébrica.

Ferramentas numéricas e computacionais

Dada a complexidade dos sistemas de ordem fraccionada, a prática de engenharia depende fortemente de ferramentas de software. A caixa de ferramentas FOMCON (Modelagem e Controle de Ordem Fraccional) para MATLAB é uma escolha popular. Fornece funções para calcular verificações de estabilidade usando o critério Matignon, construir aproximações (por exemplo, aproximação recursiva Oustaloup) de operadores fraccionais e realizar simulações de domínio temporal.

Outras ferramentas incluem:

  • CRONE toolbox (Comandante Robuste d’Ordre Non Entier): desenvolvido na Universidade de Bordeaux, suporta o design de controlador de ordem fraccionada e a análise de estabilidade.
  • Bibliotecas de python: fraccional (para cálculo fraccional) e ]control[ (para análise do sistema de controlo geral) podem ser utilizadas em conjunto para aplicar o critério Matignon.
  • Computação simbólica no Mathematica ou Maple pode lidar com manipulações algébricas de poderes fracionários e de busca de raízes no plano complexo.

Para uma verificação rápida, muitos engenheiros usam o seguinte algoritmo numérico:

  1. Defina a equação característica como uma expressão simbólica em s com poderes fracionários.
  2. Converte para um polinômio em w = sα se o sistema for proporcional.
  3. Calcular as raízes do polinômio em w utilizando um solucionador numérico.
  4. Para cada raiz, computar o correspondente s e seu argumento.
  5. Verificar se .arg(s) . > απ/2] para todas as raízes.

Este algoritmo pode ser implementado em algumas linhas de código e é robusto para a maioria dos sistemas proporcionais. Para sistemas não-compensadores, é preciso contar com métodos de domínio de frequência ou com avaliação numérica baseada em grade da equação característica no plano complexo.

Procedimento prático passo a passo para engenheiros

Abaixo está um fluxo de trabalho pragmático que combina o espírito clássico de Routh-Hurwitz com as necessárias modificações fracionárias.

  1. Obtenha a equação característica da função de transferência do sistema. Para um sistema de ordem fracionária, isto envolverá termos como s]γ][] com γ não-integer.
  2. [[FLT: 0]] Verifique se existe uma commensurabilidade. Determinar se todas as ordens fracionárias são inteiros múltiplos de uma base comum α. Por exemplo, as ordens 0.5, 1.5, 2.0 são múltiplos de 0.5; as ordens 0.7 e 1.2 não são múltiplos de uma única base (a menos que a base seja 0,1, mas raramente racional). Se o sistema for proporcional, prossiga para a transformação polinomial. Se não, salte para o passo 6 (método de domínio de frequência).
  3. ]Definir o polinômio em w = sα. Escrever a equação característica como P(w) = 0. Por exemplo, s[2.5[ + 3 s1.21.5] + 3 s0.5 + 2 = 0[FLT: 19][FLT][0.5[FLT] + 0[FLT(0.5] [FLT] + 0 [FLT(16) 0.5[FT=0] + 0.5[FLT(0.)[FLTT:2[FT:2] + + + 0] + 0.2[FT:
  4. Aplicar o teste clássico de Routh–Hurwitz para P(w). Construir o array de Routh. Se ocorrer qualquer alteração de sinal na primeira coluna, o sistema é instável (sob a interpretação polinomial). Se todas as entradas de primeira coluna tiverem o mesmo sinal, prossiga para o próximo passo.
  5. Verifique a condição do argumento.Computar as raízes de P(w]. Para cada raiz w[i, computar s[i] = w[i[[]i[[. O poder complexo é multivalorizado; use o valor principal. Então, verifique ?arg(s[[ii[α > απ/2. Se isto se mantiver para todas as raízes, o sistema fracional é estável.
  6. Para sistemas não-commensuráveis, use uma abordagem de domínio de frequência. Avaliar a equação característica ao longo do limite da região de estabilidade, que para sistemas fracionados é um setor. O gráfico Nyquist da função de transferência de loop aberto pode indicar estabilidade de loop fechado. Alternativamente, use o root-finding numérico no s-plane discretizando o operador fracionário, por exemplo usando a aproximação Oustaloup. Em seguida, aplique o teste clássico de Routh-Hurwitz ao sistema aproximado de ordem inteira. Este método dá um resultado aproximado e deve ser validado com simulações.

Muitos engenheiros descobrem que o caso proporcional abrange uma grande parte dos projetos práticos de controle fracionário, especialmente controladores fracionários de PID (onde as ordens derivadas e integrais são múltiplas de uma α comum, muitas vezes 0,5 ou 0,1).

Estudo de caso ilustrativo: Controlador de PID Fractional-Order

Considere um controlador PID de ordem fraccionada (PI[αDβ) que controla uma planta de primeira ordem.

Este sistema não é compatível porque as ordens fraccionadas 2.5, 1.2, 0. 1 não são múltiplos inteiros de uma base comum. Não podemos aplicar directamente o critério Matignon. Em vez disso, aproximamos os operadores fraccionados usando um filtro de ordem Oustaloup, digamos, 5 sobre uma faixa de frequências [0,01, 100] rad/s. Isto produz uma função de transferência de ordem inteira de alta ordem. O polinômio característico aproximado será de ordem 5×5 = 25 para cada termo fraccionado, mas podemos reduzir a ordem usando truncamento equilibrado. Então, o teste clássico de Routh-Hurwitz pode ser aplicado ao polinômio reduzido. Uma verificação de estabilidade via simulação confirma o resultado. Este método aproximado é o padrão de facto na indústria quando as ferramentas analíticas não estiverem disponíveis.

Em contraste, se o controlador for desenhado com ordens proporcionais (por exemplo, α = 0,5, por isso as ordens são 0.5, 1.0, 1.5, 2.0), o sistema torna-se proporcional. Então, podemos usar o critério Matignon diretamente, obtendo uma condição exata necessária e suficiente (módulo a interpretação de corte de ramificação). Isto ilustra a vantagem de projetar controladores fracionários com ordens proporcionais quando possível.

Limitações e Pesquisa em andamento

A extensão da análise de Routh-Hurwitz aos sistemas de ordem fraccionada não é uma única fórmula, mas uma coleção de técnicas, cada uma com o seu próprio domínio de aplicabilidade. O critério Matignon é exato apenas para sistemas proporcionais. Para sistemas não-commensuráveis, a abordagem aproximada de ordem inteira pode introduzir erros, especialmente se o sistema tiver modos de alta frequência levemente amortecidos. Pesquisadores estão explorando critérios algébricos diretos para sistemas não-commensurados usando o conceito de matrizes de ordem fraccionada (por exemplo, usando a transformação variável de ordem fraccionada ]s = (z)1/α[ e aplicando um array modificado). Outra direção promissora é o uso de desigualdades de matriz linear (LMIs) adaptadas a sistemas de ordem fraccionada.

Além disso, a região de estabilidade para sistemas fracionários com ordens irracionais não é uma simples cunha; pode se tornar uma região mais complicada na superfície de Riemann. Os engenheiros práticos geralmente dependem da avaliação numérica da magnitude e da fase da equação característica ao longo do limite de estabilidade, que é computacionalmente intensiva, mas produz resultados confiáveis para valores de parâmetros específicos.

Conclusão

A extensão da análise de Routh-Hurwitz aos sistemas de ordem fraccionada requer que se mova para além da estrutura polinomial clássica. O critério de estabilidade de Matignon para sistemas proporcionais re-lança o problema numa condição de argumento raiz que pode ser verificada com a ajuda da matriz clássica de Routh e do enquadramento numérico. Para sistemas não-compensadores, os métodos de domínio de frequência e as aproximações numéricas continuam a ser as ferramentas primárias. Ao compreender tanto a base teórica como o fluxo de trabalho numérico prático, os engenheiros podem realizar com confiança a análise de estabilidade para uma ampla gama de sistemas de ordem fraccionada — aproximando o fosso entre as aplicações de controlo clássico e as aplicações de ordem fraccionada emergentes.

Para leitura posterior, os leitores podem consultar a entrada Wikipédia sobre o critério Routh–Hurwitz e a visão geral do cálculo fracionário. O trabalho original de Matignon é detalhado em este artigo[ (Link DOI). Além disso, a caixa de ferramentas FOMCON para MATLAB[] fornece a implementação prática dos métodos descritos aqui.