def test_transformer(self):
import os
from PySpice.Spice.Netlist import Circuit
circuit = Circuit('Diode Characteristic Curve')
circuit.L('primary', 'Vlp', 'Vdrain', '{l_trf}')
circuit.C('resonance', 'Vlv', 'Vdrain', '{cap_r}')
circuit.L('secondary', 'Vls', 'ghv', '{Ls}')
circuit.R('secondary', 'Vls', 1, 5.15)
circuit.K('flyback', 'Lprimary', 'Lsecondary', 1)
def test(self):
self._test_spice_declaration(
Resistor(Circuit(''), '1', 'n1', 'n2', 100),
'R1 n1 n2 100'.lower())
self._test_spice_declaration(
Resistor(Circuit(''), '1', 'n1', 'n2', kilo(1)),
'R1 n1 n2 1k'.lower())
self._test_spice_declaration(
Resistor(Circuit(''),
'1',
'n1',
'n2',
kilo(1),
ac=kilo(2),
multiplier=2,
scale=1.5,
temperature=25,
device_temperature=26,
noisy=True),
'R1 n1 n2 1k ac=2k dtemp=26 m=2 noisy=1 scale=1.5 temp=25'.lower())
self._test_spice_declaration(
Resistor(Circuit(''), '1', 'n1', 'n2', kilo(1), noisy=False),
'R1 n1 n2 1k'.lower())
self._test_spice_declaration(
XSpiceElement(Circuit(''), '1', 1, 0, model='cap'),
'A1 1 0 cap'.lower())
def test_ground_node(self):
circuit = Circuit('')
circuit.V('input', 'in', circuit.gnd, '10V')
circuit.R(1, 'in', 'out', 9 @ u_kΩ)
circuit.R(2, 'out', circuit.gnd, 1 @ u_kΩ)
self.assertTrue(circuit.has_ground_node())
circuit = Circuit('')
circuit.V('input', 'in', 'fake_ground', '10V')
circuit.R(1, 'in', 'out', 9 @ u_kΩ)
circuit.R(2, 'out', 'fake_ground', 1 @ u_kΩ)
self.assertFalse(circuit.has_ground_node())
def nmos_characteristics(Vd, Vg):
circuit = Circuit('NMOS Input Characteristic')
circuit.include(
'./app/circuits/libraries/transistor/ptm_65nm_nmos_bulk.mod')
# Define the DC supply voltage value
Vdd = float(Vd)
# Instanciate circuit elements
Vgate = circuit.V('gate', 'gatenode', circuit.gnd, u_V(float(Vg)))
Vdrain = circuit.V('drain', 'vdd', circuit.gnd, u_V(Vdd))
# M <name> <drain node> <gate node> <source node> <bulk/substrate node>
circuit.MOSFET(1,
'vdd',
'gatenode',
circuit.gnd,
circuit.gnd,
model='ptm65nm_nmos')
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.dc(Vdrain=slice(0, Vdd, .01))
figure, ax = plt.subplots(figsize=(20, 10))
ax.plot(analysis['vdd'], u_A(-analysis.Vdrain))
ax.legend('NMOS characteristic')
ax.grid()
ax.set_xlabel('Vds [V]')
ax.set_ylabel('Id [A]')
return circuit, analysis, plt
def negative_clamper(v, r, c, f):
print(v, r, c, f)
circuit = Circuit('Negative Clamper')
circuit.include('./app/circuits/libraries/diode/switching/1N4148.lib')
source = circuit.SinusoidalVoltageSource('input',
'in',
circuit.gnd,
amplitude=u_V(float(v)),
frequency=u_Hz(float(f)))
circuit.C('C1', 'in', 'output', u_mF(float(c)))
circuit.X('D1', '1N4148', 'output', circuit.gnd)
circuit.R('load', 'output', circuit.gnd, u_Ohm(float(r)))
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=source.period / 200,
end_time=source.period * 2)
plt.close('all')
plt.title('Negative Clamper')
plt.xlabel('Time [s]')
plt.ylabel('Voltage [V]')
plt.grid()
plt.plot(analysis['in'])
plt.plot(analysis.output)
plt.legend(('input', 'output'), loc=(.05, .1))
plt.tight_layout()
return circuit, analysis, plt
def high_pass_rc_filter(v, r, c):
circuit = Circuit('Low-Pass RC Filter')
circuit.SinusoidalVoltageSource('input',
'in',
circuit.gnd,
amplitude=u_V(float(v)))
C1 = circuit.C(1, 'in', 'out', u_uF(float(c)))
R1 = circuit.R(1, 'out', circuit.gnd, u_kΩ(float(r)))
break_frequency = 1 / (2 * math.pi * float(R1.resistance * C1.capacitance))
print("Break frequency = {:.1f} Hz".format(break_frequency))
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.ac(start_frequency=1 @ u_Hz,
stop_frequency=1 @ u_MHz,
number_of_points=10,
variation='dec')
# print(analysis.out)
plt.close('all')
figure, axes = plt.subplots(2, figsize=(20, 10))
plt.title("Bode Diagram of a Low-Pass RC Filter")
bode_diagram(
axes=axes,
frequency=analysis.frequency,
gain=20 * np.log10(np.absolute(analysis.out)),
phase=np.angle(analysis.out, deg=False),
marker='.',
color='blue',
linestyle='-',
)
for ax in axes:
ax.axvline(x=break_frequency, color='red')
plt.tight_layout()
return circuit, analysis, plt
def simulate_RLC(w=2000,V=2,R=10,L=3,C=320,**kwargs):
circuit = Circuit('Name me please')
circuit.R('1', 1, 2,R@u_kΩ)
circuit.L('1', 2, 3,L@u_H)
circuit.C('1', 3, 0,C@u_nF)
#circuit.V('1', circuit.gnd,1 ,f'SIN(0 {V} {w})')
circuit.V('1',0,1, f'DC 0 AC {V} SIN(0 {V} {w})')
dt = 2*np.pi/w*4
tf = 5/w
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
#analysis = simulator.transient(step_time=0.1@u_ms, end_time=tf@u_s)
analysis = simulator.transient(step_time=dt@u_ms, end_time=tf@u_s)
if kwargs.get('view',False):
#print(str(circuit))
fig = plt.figure(figsize=(20,4)) # create a figure object
ax = fig.add_subplot(1, 1, 1)
plt.plot(analysis['1']-analysis['2'],label='R')
plt.plot(analysis['2']-analysis['3'],label='L')
ax.plot(analysis['3'],label='C')
#ax.set_ylim(-int(V*1.1)-10,int(V*1.1)+10)
ax.legend()
ax.set_title(f'freq : {w}, voltage: {V}')
ax.set_xlabel('time (ms)')
ax.set_ylabel('current (ma)')
plt.savefig('gallery/RLC_example.png')
print(f'resonance was at {np.sqrt(1/(L*C))}')
return max(abs(analysis['1']-analysis['2']))/R
def create_circuit(genome, config):
libraries_path = '/home/alan/ngspice/libraries' # os.path.join(os.path.dirname(os.path.dirname(__file__)), 'libraries')
spice_library = SpiceLibrary(libraries_path)
circuit = Circuit('NEAT')
circuit.include(spice_library['1N4148'])
Vbase = circuit.V('base', 'input1', circuit.gnd, 2)
Vcc = circuit.V('cc', 'input2', circuit.gnd, 5)
Vgnd = circuit.V('gnd', 'input3', circuit.gnd, 0)
#circuit.R('test1', 'node0', circuit.gnd, 1e6)
#circuit.R('test2', 'node0', 'input1', 1e6)
ridx = 1
xidx = 1
for key, c in iteritems(genome.connections):
if c.component == 'resistor':
pin0, pin1 = get_pins(key)
R = 10**c.value
circuit.R(ridx, pin1, pin0, R)
ridx += 1
elif c.component == 'diode':
pin0, pin1 = get_pins(key)
circuit.X(xidx, '1N4148', pin1, pin0)
xidx += 1
return circuit
def test_spinit(self):
from PySpice.Spice.Netlist import Circuit
import PySpice.Unit as U
circuit = Circuit('Astable Multivibrator')
source = circuit.V('cc', 'vcc', circuit.gnd, 15 @ U.u_V)
circuit.R(1, 'output', 'comparator', 1 @ U.u_kΩ)
circuit.C(1, 'comparator', circuit.gnd, 100 @ U.u_nF)
circuit.R(2, 'output', 'reference', 100 @ U.u_kΩ)
circuit.R(3, 'vcc', 'reference', 100 @ U.u_kΩ)
circuit.R(4, 'reference', circuit.gnd, 100 @ U.u_kΩ)
circuit.NonLinearVoltageSource(1,
'output',
circuit.gnd,
expression='V(reference, comparator)',
table=((-U.micro(1), 0),
(U.micro(1), source.dc_value)))
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
simulator.initial_condition(
comparator=0) # Fixme: simulator.nodes.comparator == 0
analysis = simulator.transient(step_time=1 @ U.u_us,
end_time=500 @ U.u_us)
if (len(analysis.output)) < 500:
raise NameError('Simualtion failed')
def freq_resp():
shutil.copyfile('PA13.LIB', '/tmp/PA13.LIB')
shutil.copyfile('op27.cir', '/tmp/op27.cir')
circuit = Circuit('Freq Response')
circuit.include('/tmp/PA13.LIB')
circuit.include('/tmp/op27.cir')
circuit.subcircuit(PA13Amplifier(Z_1, Z_2))
circuit.subcircuit(OP27Amplifier(R_4, C_4))
# V_Ir input
circuit.SinusoidalVoltageSource('input', 'vr', circuit.gnd, amplitude=0.05)
circuit.R('r3', 'vir', 'vir2', Z_3)
# Current controller stage
circuit.X('curr', 'op27_amplifier', 'vir2', 'vr')
# Voltage stage
circuit.X('volt', 'pa13_amplifier', 'vr', 'vo')
# Motor stage
circuit.R('rm', 'vo', 'vm1', R_m)
circuit.L('lm', 'vm1', 'vio', L_m)
circuit.R('rs', 'vio', circuit.gnd, R_s)
circuit.R('r5', 'vio', 'vir2', Z_5)
simulator = circuit.simulator()
# Force ngspice to have shorter time steps near sharp transitions
simulator.options(trtol=0.0001)
import pdb
pdb.set_trace()
analysis = simulator.ac(
start_frequency=50e0,
stop_frequency=5 * f_h,
number_of_points=1000, # Lab manual suggests 20, might as well do more
variation='dec')
return analysis
def create_test_circuit(fet_type, iparam, fet_L, fet_W, coner_path):
c = Circuit('gm_id')
c.include(
'/home/tclarke/skywater-pdk/libraries/sky130_fd_pr/latest/models/corners/tt.spice'
)
fet_L = 0.15
fet_W = 1
# create the circuit
c.V('gg', 1, c.gnd, 0 @ u_V)
c.V('dd', 2, c.gnd, 1.8 @ u_V)
c.X('M1',
fet_type,
2,
1,
c.gnd,
c.gnd,
L=fet_L,
W=fet_W,
ad="'W*0.29'",
pd="'2*(W+0.29)'",
as_="'W*0.29'",
ps="'2*(W+0.29)'",
nrd="'0.29/W'",
nrs="'0.29/W'",
sa=0,
sb=0,
sd=0,
nf=1,
mult=1)
return c
def test_spinit(self):
from PySpice.Spice.Netlist import Circuit
import PySpice.Unit as U
circuit = Circuit('Test')
# Fixme: On Windows
# Supplies reduced to 2.5749% Supplies reduced to 1.7100% Warning: source stepping failed
# doAnalyses: Too many iterations without convergence
source = circuit.V('cc', 'vcc', circuit.gnd, 15 @ U.u_V)
circuit.R(1, 'output', 'comparator', 1 @ U.u_kΩ)
circuit.C(1, 'comparator', circuit.gnd, 100 @ U.u_nF)
circuit.R(2, 'output', 'reference', 100 @ U.u_kΩ)
circuit.R(3, 'vcc', 'reference', 100 @ U.u_kΩ)
circuit.R(4, 'reference', circuit.gnd, 100 @ U.u_kΩ)
# circuit.NonLinearVoltageSource(1, 'output', circuit.gnd,
# expression='V(reference, comparator)',
# table=((-U.micro(1), 0),
# (U.micro(1), source.dc_value))
# )
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
simulator.initial_condition(
comparator=0) # Fixme: simulator.nodes.comparator == 0
analysis = simulator.transient(step_time=1 @ U.u_us,
end_time=500 @ U.u_us)
if (len(analysis.output)) < 500:
raise NameError('Simualtion failed')
def simulate_attenuation_factor_core(w, R_1, R_2, C_1, C_2, optboolean=False):
circuit = Circuit('Name me please')
circuit.include(spice_library['D1N4148'])
# TLV3201 is a 0-5 T.I. OpAmp
circuit.include(spice_library['LM741'])
V = 7
#w = 20E3
#circuit.V('1', circuit.gnd,1,f'DC 0 AC {V} SIN(0 {V} {w})')
circuit.V('1', 1, circuit.gnd, f'DC 0 AC {V} SIN(0 {V} {w})')
#R_1 = 11.2
#R_2 = 11.2
R_A = 100000
R_B = 0.001
#C_1 = 2000
#C_2 = 1000
links = [
(1, 2),
(2, 3),
(circuit.gnd, 4),
(4, 5),
]
R_vector = [
R_1,
R_2,
R_A,
R_B,
]
for x in range(len(R_vector)):
circuit.R(str(x + 1), links[x][0], links[x][1], R_vector[x] @ u_kΩ)
circuit.C('1', 2, 5, C_1 @ u_pF)
circuit.C('2', circuit.gnd, 3, C_2 @ u_pF)
circuit.X('opamp', 'LM741', 3, 4, 'Vcc', 'Vee', 5)
circuit.V('2', 'Vcc', circuit.gnd, 'DC +15')
circuit.V('3', 'Vee', circuit.gnd, 'DC -15')
if optboolean:
circuit.X('diodus', 'D1N4148', 2, 1)
f = w
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=((1 / f) / 10) @ u_s,
end_time=(4 / f) @ u_s) # Printer?
if False:
#print(str(circuit))
fig = plt.figure(figsize=(20, 4)) # create a figure object
ax = fig.add_subplot(1, 1, 1)
plt.plot(analysis['1'] - analysis['2'], label='R')
plt.plot(analysis['2'] - analysis['3'], label='L')
ax.plot(analysis['3'], label='C')
#ax.set_ylim(-int(V*1.1)-10,int(V*1.1)+10)
ax.legend()
ax.set_title(f'freq : {w}')
print(f'resonance was at {np.sqrt(1/(L*C))}')
if True:
outputter = max(analysis['5'])
del (circuit)
return outputter
return
def test_subcircuit(self):
circuit = Circuit('')
circuit.include('.../mosdriver.lib')
circuit.X('test', 'mosdriver', '0', '1', '2', '3', '4', '5')
circuit.BehavioralSource('test',
'1',
'0',
voltage_expression='if(0, 0, 1)',
smoothbsrc=1)
print(circuit)
def instantiate_circuit(self): ''' class method that instantiates a circuit model from an inherited netlist ''' self.circuit = Circuit(self.circuitname) for line in self.netlist[1:]: try: self.add_circuit_element_from_netlist_line(line) except: print(traceback.format_exc())
def test(self):
self._test_spice_declaration(
PieceWiseLinearVoltageSource(
Circuit(''),
'pwl1', '1', '0',
values=[(0, 0), (10@u_ms, 0), (11@u_ms, 5@u_V), (20@u_ms, 5@u_V)],
),
'Vpwl1 1 0 PWL(0s 0V 10ms 0V 11ms 5V 20ms 5V r=0s td=0.0s)',
)
def test(self):
self._test_spice_declaration(
PieceWiseLinearVoltageSource(
Circuit(''),
'pwl1',
'1',
'0',
values=[(0, 0), (10 @ u_ms, 0), (11 @ u_ms, 5 @ u_V),
(20 @ u_ms, 5 @ u_V)],
),
'Vpwl1 1 0 PWL(0s 0V 10ms 0V 11ms 5V 20ms 5V)',
)
self._test_spice_declaration(
PieceWiseLinearVoltageSource(
Circuit(''),
'pwl1',
'1',
'0',
values=[(0, 0), (10 @ u_ms, 0), (11 @ u_ms, 5 @ u_V),
(20 @ u_ms, 5 @ u_V)],
repeat_time=12 @ u_ms,
delay_time=34 @ u_ms,
),
'Vpwl1 1 0 PWL(0s 0V 10ms 0V 11ms 5V 20ms 5V r=12ms td=34ms)',
)
self._test_spice_declaration(
PieceWiseLinearVoltageSource(
Circuit(''),
'pwl1',
'1',
'0',
values=[(0, 0), (10 @ u_ms, 0), (11 @ u_ms, 5 @ u_V),
(20 @ u_ms, 5 @ u_V)],
repeat_time=12 @ u_ms,
delay_time=34 @ u_ms,
dc=50 @ u_V,
),
'Vpwl1 1 0 DC 50V PWL(0s 0V 10ms 0V 11ms 5V 20ms 5V r=12ms td=34ms)',
)
def init():
logger = Logging.setup_logging()
# libraries_path = find_libraries()
libraries_path = '/home/anurag/Workspace/DTU/CourseProj/AE/website/DTU-VLAB/dtuvlab/lab/simulations/libraries'
print(libraries_path)
spice_library = SpiceLibrary(libraries_path)
print(spice_library)
circuit = Circuit('Diode Characteristic Curve')
circuit.include(spice_library['1N4148'])
print(spice_library['1N4148'])
return circuit
def voltage_divider(v, r1, r2):
circuit = Circuit('Voltage Divider')
circuit.V('input', 1, circuit.gnd, u_V(float(v)))
circuit.R(1, 1, 2, u_kOhm(float(r1)))
circuit.R(2, 2, circuit.gnd, u_kOhm(float(r2)))
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.operating_point()
output = {}
for node in analysis.nodes.values():
output[str(node)] = str(round(float(node), 2)) + "V"
return circuit, analysis, output
def simulate_attenuation_factor_core(w, V, R, C):
circuit = Circuit('RC low-pass filter')
fc = 1 / 2 / 3.14 / (R * 1000) / (C * 1E-9)
circuit.V('1', 1, circuit.gnd, f'DC 0 AC {V} SIN(0 {V} {w})')
circuit.R('1', 1, 2, R @ u_kΩ)
circuit.C('1', 2, circuit.gnd, C @ u_nF)
f = w / 2 / np.pi
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=((1 / f) / 10) @ u_s,
end_time=(4 / f) @ u_s)
if True:
outputter = max(analysis['2'])
del (circuit)
return outputter
return
def setup_circuit(kind="biased"):
if kind=="biased":
circuit = Circuit('Biased Envelope Circuit')
circuit.include(spice_library['hsms'])
circuit.V('in', 'input', circuit.gnd, 'dc 0 external')
# bias portion
circuit.C(2, 'input', 1, 10@u_nF)
circuit.R(2, 1, 2, 1@u_kOhm)
circuit.X('D2', 'hsms', 2, circuit.gnd)
circuit.R(3, 2, 'bias', 1@u_kOhm)
circuit.V('bias', 'bias', circuit.gnd, 5@u_V)
# envelope portion
circuit.X('D1', 'hsms', 1, 'output')
circuit.C(1, 'output', circuit.gnd, 220@u_pF)
circuit.R(1, 'output', circuit.gnd, 50@u_Ohm)
return circuit
elif kind=="basic":
circuit = Circuit('Biased Envelope Circuit')
circuit.include(spice_library['hsms'])
circuit.V('in', 'input', circuit.gnd, 'dc 0 external')
# envelope portion
circuit.X('D1', 'hsms', 'input', 'output')
circuit.C(1, 'output', circuit.gnd, 220@u_pF)
circuit.R(1, 'output', circuit.gnd, 50@u_Ohm)
return circuit
elif kind=="diode":
circuit = Circuit('Diode Output')
circuit.include(spice_library['hsms'])
circuit.V('in', 'input', circuit.gnd, 'dc 0 external')
circuit.X('D1', 'hsms', 'input', 'output')
return circuit
def simple_bjt_amp():
circuit = Circuit('test circuit')
model_npn = set_model_qbc847b(circuit)
n = NodeNames('n1', 'n2', 'n3', 'n4', 'n5')
gnd = 0
circuit.R('1', n.n3, n.n2, 2.2e6)
circuit.Q('1', n.n3, n.n2, gnd, model=model_npn)
circuit.R('3', n.n1, n.n3, 1e3)
circuit.C('1', n.n4, n.n2, 100e-9)
circuit.C('2', n.n3, n.n5, 1e-6)
circuit.C('4', n.n1, n.n3, 10e-9)
circuit.R('5', n.n5, gnd, 10e3)
circuit.V('pwr', n.n1, gnd, 6)
circuit.V('in', n.n4, gnd, 'dc 0 ac 1 distof1 1 distof2 0.1')
return circuit, n
def make_database_core(w=2000,V=12,R=10,L=3,C=320):
threshold = 4
circuit = Circuit('RLC-series')
circuit.V('1',0,1, f'DC 0 AC {V} SIN(0 {V} {w})')
circuit.R('1', 1, 2,R@u_kΩ)
circuit.L('1', 2, 3,L@u_H)
circuit.X('diodus','D1N4148', 2, 3)
circuit.include(spice_library['D1N4148'])
circuit.C('1', 3, 0,C@u_nF)
dt = 2*np.pi/w*4
tau = 2*np.pi/w*4
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=dt@u_ms, end_time=tau@u_s)
auxiliar = analysis['2']-analysis['3']
return [w,V,R,L,C,int(float(max(auxiliar))>threshold)]
def test_raw_spice(self):
spice_declaration = """
.title Voltage Divider
R2 out 0 1kOhm
vinput in 0 10v
r1 in out 9kohm
"""
# .end
circuit = Circuit('Voltage Divider')
circuit.V('input', 'in', circuit.gnd, '10V')
circuit.R(1, 'in', 'out', raw_spice='9kOhm')
circuit.raw_spice += 'R2 out 0 1kOhm'
self._test_spice_declaration(circuit, spice_declaration)
def createCircuit(filter1o, filter2o):
circuit = Circuit('Filter')
circuit.include('Models/BasicOpamp.cir')
circuit.include('Models/AD8619.cir')
circuit.include('Models/TL084.cir')
circuit.subcircuit(filter1o)
circuit.subcircuit(filter2o)
circuit.V('1', '5V', circuit.gnd, '5')
circuit.V('2', 'VRef', circuit.gnd, '2.5')
circuit.SinusoidalVoltageSource('In', 'In', 'VRef', amplitude=1)
circuit.X('1', filter1o.name, 'In', 'out1o', 'VRef', '5V', circuit.gnd)
circuit.X('2', filter2o.name, 'In', 'out2o', 'VRef', '5V', circuit.gnd)
print(circuit)
return circuit
def current_divider(i, r1, r2):
circuit = Circuit('Current Divider')
circuit.I('input', 1, circuit.gnd, u_A(float(i))) # Fixme: current value
circuit.R(1, 1, circuit.gnd, u_kOhm(float(r1)))
circuit.R(2, 1, circuit.gnd, u_kOhm(float(r2)))
for resistance in (circuit.R1, circuit.R2):
resistance.minus.add_current_probe(circuit) # to get positive value
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.operating_point()
output = {}
for node in analysis.branches.values():
output[str(node)] = str(round(float(node), 2)) + "A"
return circuit, analysis, output
def full_wave_rectifier(v, r, c, f):
circuit = Circuit('Full-wave rectification')
circuit.include('./app/circuits/libraries/diode/switching/1N4148.lib')
source = circuit.SinusoidalVoltageSource('input',
'in',
circuit.gnd,
amplitude=u_V(float(v)),
frequency=u_Hz(float(f)))
circuit.X('D1', '1N4148', 'in', 'output_plus')
circuit.R('load', 'output_plus', 'output_minus', u_Ω(float(r)))
circuit.X('D2', '1N4148', 'output_minus', circuit.gnd)
circuit.X('D3', '1N4148', circuit.gnd, 'output_plus')
circuit.X('D4', '1N4148', 'output_minus', 'in')
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=source.period / 200,
end_time=source.period * 2)
plt.close('all')
figure, (ax3, ax4) = plt.subplots(2, figsize=(20, 10))
ax3.set_title('Full-Wave Rectification')
ax3.set_xlabel('Time [s]')
ax3.set_ylabel('Voltage [V]')
ax3.grid()
ax3.plot(analysis['in'])
ax3.plot(analysis.output_plus - analysis.output_minus)
ax3.legend(('input', 'output'), loc=(.05, .1))
ax3.set_ylim(float(-source.amplitude * 1.1), float(source.amplitude * 1.1))
circuit.C('1', 'output_plus', 'output_minus', u_mF(float(c)))
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=source.period / 200,
end_time=source.period * 2)
ax4.set_title('Full-Wave Rectification with filtering')
ax4.set_xlabel('Time [s]')
ax4.set_ylabel('Voltage [V]')
ax4.grid()
ax4.plot(analysis['in'])
ax4.plot(analysis.output_plus - analysis.output_minus)
ax4.legend(('input', 'output'), loc=(.05, .1))
ax4.set_ylim(float(-source.amplitude * 1.1), float(source.amplitude * 1.1))
plt.tight_layout()
return circuit, analysis, plt
def make_database_core(w, V, R_1, R_2, C_1, C_2, ax=False):
circuit = Circuit('Sallen-key low-pass-filter')
circuit.include(spice_library['D1N4148'])
circuit.include(spice_library['LM741'])
circuit.V('1', 1, circuit.gnd, f'DC 0 AC {V} SIN(0 {V} {w})')
R_A = 100000
R_B = 0.001
links = [
(1, 2),
(2, 3),
(circuit.gnd, 4),
(4, 5),
]
R_vector = [
R_1,
R_2,
R_A,
R_B,
]
for x in range(len(R_vector)):
circuit.R(str(x + 1), links[x][0], links[x][1], R_vector[x] @ u_kΩ)
circuit.C('1', 2, 5, C_1 @ u_pF)
circuit.C('2', circuit.gnd, 3, C_2 @ u_pF)
circuit.X('opamp', 'LM741', 3, 4, 'Vcc', 'Vee', 5)
circuit.V('2', 'Vcc', circuit.gnd, 'DC +15')
circuit.V('3', 'Vee', circuit.gnd, 'DC -15')
circuit.X('diodus', 'D1N4148', 2, 1)
f = w
tf = 3 / w
t0 = 1
treshold = 4
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=((1 / f) / 10) @ u_s,
end_time=3 / w @ u_s)
if False:
ax.plot(analysis['5'], label='output')
ax.plot(analysis['1'], label='input')
ax.plot(analysis['1'] - analysis['2'], label='diodo')
ax.hlines(4, 0, 5 / w)
ax.legend()
ax.set_title(
f'R1: {R_1}, R2: {R_2}, C1: {C_1}, C2: {C_2}, V: {V}, W: {w}')
ax.set_ylim(-15, 15)
auxiliar = analysis['2'] - analysis['1']
del (circuit)
return [w, V, R_1, R_2, C_1, C_2, int(float(max(auxiliar)) > treshold)]
def simulate_amplification_factor_core(R1,R2,RC,RE,RL,C1,C2,V,w):
circuit = Circuit('Amplifier')
circuit.R(1, 5, 2, R1@u_kΩ) #kOhm
circuit.R(2, 2, 0, R2@u_kΩ) #kOhm
circuit.R('C', 5, 4, RC@u_kΩ) #kOhm
circuit.R('E', 3, 0, RE@u_kΩ) #kOhm
circuit.R('Load', 'out', 0, RL@u_MΩ) #MOhm
circuit.C(1, 'inp', 2, C1@u_uF) #uF
circuit.C(2, 4, 'out', C2@u_uF) #uF
circuit.BJT(1, 4, 2, 3, model='bjt') # Q is mapped to BJT !
circuit.model('bjt', 'npn', bf=80, cjc=pico(5), rb=100)
circuit.V('power', 5, circuit.gnd, 15@u_V)
circuit.V('var','inp',circuit.gnd, f'DC 0 AC {V} SIN(0 {V}V {w})')
T = 25
simulator = circuit.simulator(temperature=T, nominal_temperature=T)
analysis = simulator.transient(step_time = (2*np.pi/w/20)@u_s,end_time=(2*np.pi/w*4)@u_s)
return [float(min(analysis['out']))]
def simulate_amplifier():
circuit = Circuit('Amplifier')
R_1 = [100]
R_2 = [20]
R_C = [10]
R_E = [2]
R_L = [1]
C_1 = [10]
C_2 = [10]
R1 = random.choice(R_1)
R2 = random.choice(R_2)
RC = random.choice(R_C)
RE = random.choice(R_E)
RL = random.choice(R_L)
C1 = random.choice(C_1)
C2 = random.choice(C_2)
circuit.R(1, 5, 2, R1@u_kΩ) #kOhm
circuit.R(2, 2, 0, R2@u_kΩ) #kOhm
circuit.R('C', 5, 4, RC@u_kΩ) #kOhm
circuit.R('E', 3, 0, RE@u_kΩ) #kOhm
circuit.R('Load', 'out', 0, RL@u_MΩ) #MOhm
circuit.C(1, 'in', 2, C1@u_uF) #uF
circuit.C(2, 4, 'out', C2@u_uF) #uF
circuit.BJT(1, 4, 2, 3, model='bjt') # Q is mapped to BJT !
circuit.model('bjt', 'npn', bf=80, cjc=pico(5), rb=100)
V_vector = [0.5]
w_vector = [1E3]
V = random.choice(V_vector)
w = random.choice(w_vector)
circuit.V('power', 5, circuit.gnd, 15@u_V)
circuit.V('var','in',circuit.gnd, f'DC 0 AC {V} SIN(0 {V}V {w})')
simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time = (1/w/25)@u_s,end_time=(4/w)@u_s)
figure = plt.figure(1, (20, 10))
axe = plt.subplot(111)
plt.title('')
plt.xlabel('Time [s]')
plt.ylabel('Voltage [V]')
plt.grid()
plot(analysis['in'], axis=axe)
plot(analysis.out, axis=axe)
plt.legend(('input', 'output'), loc=(.05,.1))
plt.tight_layout()
plt.savefig('gallery/amplifier.png')
return
