Definição formal da derivada

O texto a seguir foi retirado de um conversa com o CoPilot da Microsoft e para transformar em codigo HTML para alimentar este blog utilizamos o conversor de Markdown para Html

Exatamente, Rafael 👌. Você trouxe a essência da definição da derivada e como ela aparece tanto na teoria matemática quanto na prática computacional em simulações.

📐 Definição formal da derivada

A derivada de uma função f(x) em um ponto (x) é definida como o limite:

$f^{'} (x) = lim_{Δ x \to 0} \frac{f (x + Δ x) - f (x)}{Δ x}$

Esse limite mede a taxa de variação instantânea da função.
Em termos geométricos, é a inclinação da reta tangente à curva no ponto (x).

💻 Aproximação numérica (diferenças finitas)

Nos computadores, não conseguimos calcular o limite diretamente. Então usamos aproximações:

Diferença para frente: $f^{'} (x) \approx \frac{f (x + Δ x) - f (x)}{Δ x}$
Diferença para trás: $f^{'} (x) \approx \frac{f (x) - f (x - Δ x)}{Δ x}$
Diferença central (mais precisa): $f^{'} (x) \approx \frac{f (x + Δ x) - f (x - Δ x)}{2 Δ x}$

🔄 Ligação com o FOR no Matlab/Python

Quando você descreveu o loop, é exatamente isso:

for i = 2:N
    di_dt = (i_vec(i) - i_vec(i-1)) / delta_t;
    vL(i) = L * di_dt;
end

i_vec(i) → corrente no instante atual
i_vec(i-1) → corrente no instante anterior
(i_vec(i) - i_vec(i-1)) / delta_t → aproximação da derivada $\frac{d i}{d t}$
vL(i) = L * di/dt → lei do indutor: tensão é proporcional à taxa de variação da corrente.

⚡ Interpretação física

O indutor resiste a variações rápidas de corrente.
Por isso, $V_{L} = L \cdot \frac{d i}{d t}$
Se a corrente muda muito rápido, a tensão no indutor cresce bastante.
A simulação numérica captura isso usando a diferença entre valores consecutivos no vetor de corrente.

👉 Se quiser, eu posso te mostrar como implementar esse mesmo cálculo em Python com numpy, para comparar com o estilo que você usava no Matlab. Quer que eu monte esse exemplo?

Perfeito, Rafael ⚡. Vamos estruturar isso passo a passo para que você consiga simular a tensão no indutor (V_L) em um circuito RLC série submetido a uma fonte de 10 V DC aplicada em (t=0).

1. Equações do circuito RLC série

O circuito tem três elementos em série: resistor (R), indutor (L) e capacitor (C).

A equação diferencial que descreve a corrente (i(t)) é:

$L \frac{d^{2} i (t)}{d t^{2}} + R \frac{d i (t)}{d t} + \frac{1}{C} i (t) = \frac{d V_{f o n t e} (t)}{d t}$

👉 Como a fonte é DC (10 V constante), sua derivada é zero. Então a equação se simplifica para:

$L \frac{d^{2} i (t)}{d t^{2}} + R \frac{d i (t)}{d t} + \frac{1}{C} i (t) = 0$

Ou seja, é uma equação homogênea de segunda ordem.

2. Condições iniciais

Antes de (t=0), o circuito está em repouso: (i(0)=0).
No instante (t=0^+), o capacitor está inicialmente descarregado, então a tensão de 10 V aparece sobre ele.
Isso gera uma dinâmica transitória onde a corrente começa a crescer e o indutor reage.

3. Tensão no indutor

A tensão no indutor é dada por:

$V_{L} (t) = L \cdot \frac{d i (t)}{d t}$

👉 Para calcular numericamente, usamos a aproximação por diferenças finitas:

$\frac{d i}{d t} \approx \frac{i [k] - i [k - 1]}{Δ t}$

Assim:

$V_{L} [k] = L \cdot \frac{i [k] - i [k - 1]}{Δ t}$

4. Implementação em Python (estilo Matlab)

Aqui está um esqueleto de código para simular:

import numpy as np
import matplotlib.pyplot as plt

# Parâmetros do circuito
R = 100      # ohms
L = 0.01     # henry
C = 0.0001   # farad
Vdc = 10     # volts

# Tempo de simulação
dt = 1e-5
t_max = 0.01
t = np.arange(0, t_max, dt)

# Vetores de corrente e tensão
i = np.zeros_like(t)
vL = np.zeros_like(t)

# Simulação usando método de Euler
for k in range(1, len(t)):
    # Equação diferencial discretizada
    di_dt = (Vdc - R*i[k-1] - (1/C)*np.trapz(i[:k], t[:k])) / L
    i[k] = i[k-1] + di_dt*dt
    
    # Tensão no indutor
    vL[k] = L * (i[k] - i[k-1]) / dt

# Plotar resultados
plt.plot(t, vL, label="Tensão no Indutor")
plt.plot(t, i, label="Corrente no circuito")
plt.legend()
plt.xlabel("Tempo (s)")
plt.show()

5. Interpretação

Inicialmente, (V_L) é alto porque a corrente varia rapidamente.
Com o tempo, o indutor se “acostuma” e a tensão nele cai.
No regime permanente, como a fonte é DC, o indutor se comporta como um curto-circuito ((V_L \to 0)) e o capacitor como um circuito aberto ((i \to 0)).

👉 Quer que eu rode essa simulação para você e mostre o gráfico da tensão no indutor e da corrente ao longo do tempo?