3a

Prévia do material em texto

import numpy as np
import matplotlib.pyplot as plt
# Parâmetros do circuito
R1 = 1e3 # Resistência R1 em ohms (1kΩ)
R2 = 1e3 # Resistência R2 em ohms (1kΩ)
V_m = 10.0 # Amplitude da tensão senoidal (Volts)
f = 1e3 # Frequência da tensão senoidal (1 kHz)
# Impedância equivalente do circuito (resistências em série)
Z_eq = R1 + R2 # Soma direta pois resistências estão em série
# Gerar vetor de tempo
t = np.linspace(0, 2/f, 1000) # 2 períodos completos da onda
# Frequência angular
omega = 2 * np.pi * f
# Tensão aplicada v(t) = V_m * sin(ωt)
v_t = V_m * np.sin(omega * t)
# Corrente no circuito i(t) = v(t) / Z_eq (Lei de Ohm)
i_t = v_t / Z_eq
# Tensão no resistor R2 (V_R = i(t) * R2)
v_R2 = i_t * R2
# Plotar a tensão e corrente no circuito
plt.figure(figsize=(10, 6))
# Tensão aplicada
plt.subplot(2, 1, 1)
plt.plot(t, v_t, label="Tensão aplicada v(t)")
plt.plot(t, v_R2, label="Tensão no resistor R (v_R2)")
plt.title("Tensão no Circuito")
plt.xlabel("Tempo (s)")
plt.ylabel("Tensão (V)")
plt.grid(True)
plt.legend()
# Corrente no circuito
plt.subplot(2, 1, 2)
plt.plot(t, i_t, label="Corrente i(t)", color='r')
plt.title("Corrente no Circuito")
plt.xlabel("Tempo (s)")
plt.ylabel("Corrente (A)")
plt.grid(True)
plt.legend()
plt.tight_layout()
plt.show()
# Imprimir resultados
print(f"Impedância equivalente Z_eq = {Z_eq} ohms")
Impedância equivalente Z_eq = 2000.0 ohms