Fluxos em uma dimensão

Fase 1 — Capítulo 1 (Strogatz, caps. 1–2)

Semana1 de 24
FaseStrogatz
Tempo60–90 min

Joana, segunda-feira, 8h.

Joana chegou ao polo da UAB-UNITINS em Augustinópolis com o caderno aberto na contagem de matrículas. A turma de Pedagogia que entrou em 2022 tinha 38 alunos. No quarto semestre, 22. No oitavo, projetava-se algo entre 12 e 15. A turma de 2023, com 41 ingressantes, parecia caminhar para 25 no quarto semestre — mais devagar.

Ela não procurava ainda uma fórmula. Procurava um padrão. Existe um número de equilíbrio para cada polo? Uma turma que estabiliza vira “estável”? Uma que despenca colapsa para zero?

À tarde, abriu o RStudio e em seis linhas de geobr e censobr puxou o shape do município de Augustinópolis e a faixa etária 18–24 anos do Censo 2022. Antes de qualquer cálculo, tinha o território em frente — pixel por pixel — do sistema que ia diagnosticar.

Sentou-se com o mapa aberto e começou a desenhar setas: matrícula → trancamento → reativação → conclusão → evasão. Foi o primeiro retrato de fase da vida dela.

Pense numa caixa d’água com torneira aberta na entrada e ralo aberto na saída. Se entra mais água do que sai, o nível sobe. Se sai mais do que entra, desce. Em algum momento, se as duas torneiras ficam constantes, o nível para de mexer — não está vazio, não está cheio, está num ponto de equilíbrio. Esse ponto é o que vamos estudar.

A matemática do capítulo só formaliza essa intuição. Quem entender a caixa d’água já entendeu metade.

A pergunta operacional

A pergunta central da Fase 1 é geometricamente simples: dado um sistema autônomo de primeira ordem na reta, \(\dot{x} = f(x)\), com \(f\) suficientemente regular, o que é possível dizer sobre o comportamento de longo prazo de qualquer trajetória, sabendo apenas onde \(f\) se anula e qual é o sinal de \(f\) entre seus zeros? A resposta de Strogatz (2015) é desproporcional à modéstia da pergunta: para a dinâmica unidimensional autônoma, tudo o que importa está nos zeros de \(f\) (os pontos fixos) e na monotonicidade local de \(f\) em torno deles.

A razão é topológica. Em um espaço de estados unidimensional, uma trajetória é uma curva monótona até atingir um ponto fixo ou divergir; não há laços fechados, não há oscilações sustentadas, não há caos. Esse fato — que parece pobre — funda toda a Fase 1: é por ser pobre que é tratável, e tratabilidade é a condição para que possamos depois reintroduzir, ordenadamente, complexidade.

Ponto fixo

Estado x* em que f(x*) = 0. Sistema iniciado ali não se move.

Estabilidade linear

Se f'(x*) < 0: atrai. Se f'(x*) > 0: repele. Se zero: degenerado.

Bacia de atração

Conjunto de condições iniciais cujas trajetórias convergem para um x* estável.

Potencial V(x)

V = -∫f. Trajetórias descem o gradiente. Pontos fixos estáveis = mínimos locais.

Pontos fixos e estabilidade linear

Um ponto fixo \(x^*\) de \(\dot{x}=f(x)\) é um zero de \(f\), isto é, \(f(x^*)=0\). Iniciado em \(x^*\), o sistema permanece em \(x^*\) para todo \(t \geq 0\). A questão delicada é: iniciado próximo a \(x^*\), o sistema retorna ou se afasta?

A resposta é dada pela análise de estabilidade linear. Expandindo \(f\) em série de Taylor em torno de \(x^*\),

\[ \dot{x} \;=\; f(x^*) + f'(x^*)(x-x^*) + O\big((x-x^*)^2\big) \;=\; f'(x^*)\,\eta + O(\eta^2), \]

com \(\eta := x - x^*\), e desprezando os termos de ordem superior, obtém-se a equação linear \(\dot{\eta} = f'(x^*)\,\eta\), cuja solução é exponencial: \(\eta(t) = \eta(0)\,e^{f'(x^*)t}\). Daí o critério:

\[ \boxed{\;f'(x^*) < 0 \;\Rightarrow\; x^*\text{ estável};\qquad f'(x^*) > 0 \;\Rightarrow\; x^*\text{ instável}.\;} \tag{1}\]

O caso \(f'(x^*)=0\) é degenerado e exige análise não-linear. Geometricamente: se \(f\) cruza o eixo \(x\) descendo no zero \(x^*\), então \(f\) é positiva à esquerda (sistema vai para a direita) e negativa à direita (sistema vai para a esquerda) — atrai. Se \(f\) cruza subindo, repele.

Existência e unicidade (Picard)

Antes de tratar comportamento de longo prazo, é prudente garantir que as soluções existem e são únicas. O teorema de Picard-Lindelöf estabelece que, se \(f\) é Lipschitz contínua em uma vizinhança de \(x_0\), o problema de Cauchy

\[ \dot{x} = f(x), \qquad x(t_0) = x_0, \]

admite solução única em algum intervalo de tempo aberto contendo \(t_0\). A consequência operacional é forte: trajetórias não se cruzam. Em 1D, isso já força o comportamento monótono — duas trajetórias que se cruzassem violariam a unicidade de Picard.

O potencial \(V(x)\) e a interpretação energética

A função \(V(x) := -\int^x f(s)\,ds\) define um potencial fictício sobre o espaço de estados. Pela definição, \(\dot{x} = f(x) = -V'(x)\), ou seja, o sistema desce o gradiente de \(V\). Pontos fixos estáveis correspondem a mínimos locais; instáveis, a máximos. A “bacia de atração” do dinamicismo é literalmente uma bacia do potencial.

Gráfico do potencial V(N) com bola descendo até o mínimo local em N igual a K
Figura 1: Bola descendo o gradiente de \(V(N)\) até o mínimo local em \(N=K\). A trajetória é monótona; o ponto fixo estável é literalmente o fundo da bacia.

Voz dos personagens

Exemplos canônicos

Modelo logístico (Verhulst)

\(\dot{N} = r\,N\!\left(1 - N/K\right)\), com \(r,K > 0\). Pontos fixos: \(N^* = 0\) (instável, \(f'=r>0\)) e \(N^* = K\) (estável, \(f'=-r<0\)). A bacia de atração de \(K\) é \((0, \infty)\). Tempo característico: \(\sim 1/r\). Modelo canônico de capacidade de carga — cresce exponencialmente longe do limite, satura próximo a ele.

Carregamento RC

\(\dot{V} = (V_0 - V)/\tau\), com \(\tau = RC\). Único ponto fixo \(V^* = V_0\), globalmente estável. Solução explícita: \(V(t) = V_0 + (V(0) - V_0)\,e^{-t/\tau}\). O tempo característico \(\tau = |f'(x^*)|^{-1}\) é o tempo de meia-vida local da perturbação.

Curva exponencial crescente saturando em V0 com marcação de tau
Figura 2: Carregamento RC: \(V(t)\) aproxima-se exponencialmente de \(V_0\) com tempo característico \(\tau = RC\). Em \(t = \tau\), o sistema cobriu 63% do caminho.

Polo UAB-UNITINS de Augustinópolis (aplicação institucional)

Considere \(S(t)\) = número de estudantes ativos no polo num dado instante. Em modelagem agregada (capítulo F3-01), o balanço aproxima-se por

\[ \dot{S} = i(t) - \frac{S}{\tau_{\text{permanência}}}, \]

onde \(i(t)\) é o fluxo de novas matrículas e \(\tau_{\text{permanência}}\) é o tempo médio que um estudante permanece ativo. Em regime estacionário com \(i\) constante: \(S^* = i\,\tau_{\text{permanência}}\), estável com tempo característico \(\tau_{\text{permanência}}\). Forma RC aplicada à demografia institucional. O sistema filtra distúrbios com banda passante \(1/\tau\). Para Joana: comparar polos UAB-UNITINS de Augustinópolis e Dianópolis pelo \(\tau\) estimado é o primeiro passo de diagnóstico — polos com \(\tau\) pequeno (alta evasão) e polos com \(\tau\) grande (baixa rotatividade) operam em regimes dinâmicos qualitativamente diferentes.

Conexão com Ashby/Beer

Em linguagem cibernética, um ponto fixo estável é a definição rigorosa de homeostase: estado para o qual o sistema retorna após perturbações pequenas. A bacia de atração delimita o conjunto de distúrbios que o sistema absorve sem mudar de regime — uma instância concreta da entropia \(H(D)\) que aparecerá em F4-01. O critério \(|f'(x^*)|\) mede a taxa de absorção; a recíproca, \(1/|f'(x^*)|\), mede o tempo de absorção. Em VSM (capítulo F4-01), esse tempo característico é o que limita a banda passante de S2 (antioscilação) e S3 (alocação aqui-agora): se um distúrbio chega mais rápido do que \(\tau = 1/|f'|\), o regulador simplesmente não tem tempo de responder.

Importante

A diferença entre “regulador com variedade requisita” e “regulador rápido o suficiente” é o que separa o desenho institucional bom no papel do desenho institucional que funciona na prática.

Pergunta de verificação

Eixo horizontal com três pontos fixos marcados: x igual a menos um, zero e mais um, com setas convergindo aos extremos
Figura 3: Retrato de fase de \(\dot{x} = x - x^3\). Pontos fixos estáveis em \(x = \pm 1\) (círculos cheios), instável em \(x = 0\) (círculo vazio). As setas seguem o sinal de \(f\) e convergem aos estáveis; cada um tem sua bacia de atração separada por \(x = 0\).

Para o sistema \(\dot{x} = x - x^3\):

  1. Encontre todos os pontos fixos. (Três; Figura 3 antecipa.)
  2. Classifique cada um via Equação 1.
  3. Verifique o retrato de fase contra Figura 3.
  4. Identifique a bacia de atração de cada ponto fixo estável.
  5. Calcule o potencial \(V(x)\) e confirme que pontos fixos estáveis são mínimos locais.
  6. Em linguagem cibernética: qual o “regulador implícito” e qual sua “variedade requisita”?
Dica

Importe o deck Anki (vide revisão) — os cartões strogatz fase1 cobrem precisamente os conceitos deste capítulo. Quinze minutos diários por seis dias na semana, e a Fase 1 inteira fica fixada.

Próxima parada

Semana 2 — Bifurcações em 1D
O que acontece quando um parâmetro muda e a topologia colapsa?
Continuar →

Referências

STROGATZ, Steven H. Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering. 2. ed. [S.l.]: Westview Press, 2015.