Jak simulovat odporový dělič napětí pomocí PySpice

Následující kód simuluje odporový dělič napětí 10kΩ / 1kΩ pomocí PySpice a může sloužit jako dobrý výchozí bod pro simulaci jednoduchých obvodů.

Tento příspěvek ukazuje, jak simulovat dělič napětí pomocí transientní analýzy. Viz také alternativní variantu tohoto příspěvku s využitím DC sweep analýzy: Jak provést DC sweep odporového děliče napětí pomocí PySpice

pyspice_voltage_divider.py

import PySpice.Logging.Logging as Logging
logger = Logging.setup_logging()

from PySpice.Probe.Plot import plot
from PySpice.Spice.Netlist import Circuit
from PySpice.Unit import *

circuit = Circuit("MyCircuit")
# Vytvoření zdroje napětí: 5V DC
source = circuit.VoltageSource('V1', 'in', circuit.gnd, dc_value=5@u_V)
# Vytvoření odporového děliče
r1 = circuit.R('R1', 'in', 'n1', 10@u_kΩ)
r2 = circuit.R('R2', 'n1', circuit.gnd, 1@u_kΩ)
# Simulace po dobu 1 sekundy v krocích 1 milisekundy
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=1@u_ms, end_time=1@u_s)

K poli výstupních napětí děliče (tj. uzel n1) můžete přistupovat pomocí analysis['n1']:

PySpice transientní simulace odporového děliče napětí zobrazující vstupní a výstupní napětí v čase

Toto je kód, který jsme použili pro vykreslení grafu:

plot_voltage_divider.py

import matplotlib.ticker as mtick
import matplotlib.pyplot as plt
from UliEngineering.EngineerIO import format_value

def format_volts(value, pos=None):
    return format_value(value, 'V')

plt.style.use("ggplot")
plt.xlabel("Čas [ms]")
plt.ylim([0.0, 5.5])
plt.gca().yaxis.set_major_formatter(mtick.FuncFormatter(format_volts))
plt.gcf().set_size_inches(8,5)
plt.plot(analysis["in"], label="Vstupní napětí")
plt.plot(analysis["n1"], label="Výstupní napětí děliče")
plt.gca().legend()
plt.savefig("/ram/PySpice-Voltage-Divider.svg")