示例#1
0
def test_pulse():
    """Test time_functions.pulse"""
    # as we cannot define a periodic function in sympy, we will only check one
    # period
    v1, v2, td, tr, pw, tf, per, time = sympy.symbols('v1, v2, td, tr, pw, ' +
                                                      'tf, per, time')
    F = v1*Heaviside(td - time) + \
        ((v2 - v1)/tr*time + v1 +(v1 - v2)*td/tr)* \
        Heaviside(time - td)*Heaviside(td + tr - time) + \
        v2*Heaviside(time - td - tr)*Heaviside(td + tr + pw - time) + \
        ((v1 - v2)/tf*time + v2 + (-v1*(pw + td + tr) + v2*(pw + td + tr))/tf)* \
        Heaviside(time - td - tr - pw)*Heaviside(td + tr + pw + tf - time) + \
        v1*Heaviside(time - td - tr - pw - tf)

    v1n, v2n, tdn, trn, pwn, tfn, pern = (-2, 3, 0.1e-3, 0.05e-3, .5e-3,
                                          .01e-3, 1e-3)
    f = time_functions.pulse(v1=v1n,
                             v2=v2n,
                             td=tdn,
                             tr=trn,
                             pw=pwn,
                             tf=tfn,
                             per=pern)
    FS = sympy.lambdify(
        time,
        sympy.N(
            F.subs(
                dict(v1=v1n, v2=v2n, td=tdn, tr=trn, pw=pwn, tf=tfn,
                     per=pern))))
    t = np.linspace(0, 1e-3, 1e3)
    for ti in t:
        assert np.allclose(f(ti), float(FS(ti)), rtol=1e-4)
def test_pulse():
    """Test time_functions.pulse"""
    # as we cannot define a periodic function in sympy, we will only check one
    # period
    v1, v2, td, tr, pw, tf, per, time = sympy.symbols('v1, v2, td, tr, pw, ' +
                                                      'tf, per, time')
    F = v1*Heaviside(td - time) + \
        ((v2 - v1)/tr*time + v1 +(v1 - v2)*td/tr)* \
        Heaviside(time - td)*Heaviside(td + tr - time) + \
        v2*Heaviside(time - td - tr)*Heaviside(td + tr + pw - time) + \
        ((v1 - v2)/tf*time + v2 + (-v1*(pw + td + tr) + v2*(pw + td + tr))/tf)* \
        Heaviside(time - td - tr - pw)*Heaviside(td + tr + pw + tf - time) + \
        v1*Heaviside(time - td - tr - pw - tf)

    v1n, v2n, tdn, trn, pwn, tfn, pern = (-2, 3, 0.1e-3, 0.05e-3, .5e-3,
                                          .01e-3, 1e-3)
    f = time_functions.pulse(
        v1=v1n, v2=v2n, td=tdn, tr=trn, pw=pwn, tf=tfn, per=pern)
    FS = sympy.lambdify(
        time,
        sympy.N(
            F.subs(
                dict(v1=v1n, v2=v2n, td=tdn, tr=trn, pw=pwn, tf=tfn,
                     per=pern))))
    t = np.linspace(0, 1e-3, 1e3)
    for ti in t:
        assert np.allclose(f(ti), float(FS(ti)), rtol=1e-4)
示例#3
0
def test():
    """Test pulse and sin API"""
    step = time_functions.pulse(v1=0, v2=1, td=500e-9, tr=1e-12, pw=1, tf=1e-12, per=2)
    damped_sin = time_functions.sin(vo=0, va=1, td=500e-9, freq=15e3, theta=5e3, phi=90.)
    exp = time_functions.exp(v1=.5, v2=-.05, td1=0, tau1=20e-6, td2=400e-6, tau2=20e-6)

    mycircuit = circuit.Circuit(title="Butterworth Example circuit", filename=None)

    gnd = mycircuit.get_ground_node()

    mycircuit.add_resistor(part_id="R1", n1="n1", n2="n2", value=600)
    mycircuit.add_inductor(part_id="L1", n1="n2", n2="n3", value=15.24e-3)
    mycircuit.add_capacitor(part_id="C1", n1="n3", n2=gnd, value=119.37e-9)
    mycircuit.add_inductor(part_id="L2", n1="n3", n2="n4", value=61.86e-3)
    mycircuit.add_capacitor(part_id="C2", n1="n4", n2=gnd, value=155.12e-9)
    mycircuit.add_resistor(part_id="R2", n1="n4", n2=gnd, value=1.2e3)

    mycircuit.add_vsource("V1", n1="n1", n2='n5', dc_value=3.3333, ac_value=.33333, function=step)
    mycircuit.add_vsource("V2", n1="n5", n2='n6', dc_value=3.3333, ac_value=.33333, function=damped_sin)
    mycircuit.add_vsource("V3", n1="n6", n2=gnd, dc_value=3.3333, ac_value=.33333, function=exp)

    op_analysis = ahkab.new_op(outfile='time_functions')
    ac_analysis = ahkab.new_ac(start=1e3, stop=1e5, points=100, outfile='time_functions')
    tran_analysis = ahkab.new_tran(tstart=0, tstop=1.2e-3, tstep=1e-6, x0=None, outfile='time_functions')

    testbench = testing.APITest('time_functions', mycircuit,
                                [op_analysis, ac_analysis, tran_analysis],
                                skip_on_travis=True, er=1e-3, ea=1e-5)
    testbench.setUp()
    testbench.test()

    if cli:
        r = ahkab.run(mycircuit, an_list=[op_analysis, ac_analysis, tran_analysis])
        fig = plt.figure()
        plt.title(mycircuit.title + " - TRAN Simulation")
        plt.plot(r['tran']['T'], r['tran']['VN1'], label="Input voltage")
        plt.hold(True)
        plt.plot(r['tran']['T'], r['tran']['VN4'], label="output voltage")
        plt.legend()
        plt.hold(False)
        plt.grid(True)
        #plt.ylim([0,1.2])
        plt.ylabel('Step response')
        plt.xlabel('Time [s]')
        fig.savefig('tran_plot.png')

        fig = plt.figure()
        plt.subplot(211)
        plt.semilogx(r['ac']['w'], np.abs(r['ac']['Vn4']), 'o-')
        plt.ylabel('abs(V(n4)) [V]')
        plt.title(mycircuit.title + " - AC Simulation")
        plt.subplot(212)
        plt.grid(True)
        plt.semilogx(r['ac']['w'], np.angle(r['ac']['Vn4']), 'o-')
        plt.xlabel('Angular frequency [rad/s]')
        plt.ylabel('arg(V(n4)) [rad]')
        fig.savefig('ac_plot.png')
    else:
        testbench.tearDown()
示例#4
0
文件: script.py 项目: leoli1119/ahkab
import numpy as np
import ahkab
from ahkab import ahkab, circuit, time_functions

mycircuit = circuit.Circuit(title="Butterworth Example circuit", filename=None)

gnd = mycircuit.get_ground_node()

mycircuit.add_resistor("R1", n1="n1", n2="n2", value=600)
mycircuit.add_inductor("L1", n1="n2", n2="n3", value=15.24e-3)
mycircuit.add_capacitor("C1", n1="n3", n2=gnd, value=119.37e-9)
mycircuit.add_inductor("L2", n1="n3", n2="n4", value=61.86e-3)
mycircuit.add_capacitor("C2", n1="n4", n2=gnd, value=155.12e-9)
mycircuit.add_resistor("R2", n1="n4", n2=gnd, value=1.2e3)

voltage_step = time_functions.pulse(v1=0, v2=1, td=500e-9, tr=1e-12, pw=1, tf=1e-12, per=2)
mycircuit.add_vsource("V1", n1="n1", n2=gnd, dc_value=5, ac_value=1, function=voltage_step)

print mycircuit

op_analysis = ahkab.new_op()
ac_analysis = ahkab.new_ac(start=1e3, stop=1e5, points=100)
tran_analysis = ahkab.new_tran(tstart=0, tstop=1.2e-3, tstep=1e-6, x0=None)

r = ahkab.run(mycircuit, an_list=[op_analysis, ac_analysis, tran_analysis])

import pylab

fig = pylab.figure()
pylab.title(mycircuit.title + " - TRAN Simulation")
pylab.plot(r['tran']['T'], r['tran']['VN1'], label="Input voltage")
示例#5
0
mycircuit = circuit.Circuit(title="Butterworth Example circuit", filename=None)

gnd = mycircuit.get_ground_node()

mycircuit.add_resistor("R1", n1="n1", n2="n2", value=600)
mycircuit.add_inductor("L1", n1="n2", n2="n3", value=15.24e-3)
mycircuit.add_capacitor("C1", n1="n3", n2=gnd, value=119.37e-9)
mycircuit.add_inductor("L2", n1="n3", n2="n4", value=61.86e-3)
mycircuit.add_capacitor("C2", n1="n4", n2=gnd, value=155.12e-9)
mycircuit.add_resistor("R2", n1="n4", n2=gnd, value=1.2e3)

voltage_step = time_functions.pulse(v1=0,
                                    v2=1,
                                    td=500e-9,
                                    tr=1e-12,
                                    pw=1,
                                    tf=1e-12,
                                    per=2)
mycircuit.add_vsource("V1",
                      n1="n1",
                      n2=gnd,
                      dc_value=5,
                      ac_value=1,
                      function=voltage_step)

print mycircuit

op_analysis = ahkab.new_op()
ac_analysis = ahkab.new_ac(start=1e3, stop=1e5, points=100)
tran_analysis = ahkab.new_tran(tstart=0, tstop=1.2e-3, tstep=1e-6, x0=None)
示例#6
0
 elif V_DC_pattern:
     V_Name = V_DC_pattern.group(1)
     V_Np = V_DC_pattern.group(2)
     V_Nn = V_DC_pattern.group(3)
     V_value = unit_transform(V_DC_pattern.group(7))
     mycircuit.add_vsource(V_Name, n1=V_Np, n2=V_Nn, dc_value=V_value)
 elif V_AC_pattern:
     V_Name = V_AC_pattern.group(1)
     V_Np = V_AC_pattern.group(2)
     V_Nn = V_AC_pattern.group(3)
     V_value = unit_transform(V_AC_pattern.group(7))
     V_AC_phase = V_AC_pattern.group(10)
     mycircuit.add_vsource(V_Name, n1=V_Np, n2=V_Nn, ac_value=V_value)
 elif V_PULSE_pattern:
     V_PULSE = PulseVoltageSource(line)
     v_pulse = time_functions.pulse(v1=V_PULSE.voltage_low, v2=V_PULSE.voltage_high, td=V_PULSE.delay,
                                    per=V_PULSE.period, tr=V_PULSE.rise, tf=V_PULSE.fall, pw=V_PULSE.width)
     mycircuit.add_vsource(V_PULSE.name, n1=V_PULSE.Np, n2=V_PULSE.Nn, function=v_pulse)
 elif I_DC_pattern:
     I_Name = I_DC_pattern.group(1)
     I_Np = I_DC_pattern.group(2)
     I_Nn = I_DC_pattern.group(3)
     I_value = unit_transform(I_DC_pattern.group(7))
     mycircuit.add_isource(I_Name, n1=I_Np, n2=I_Nn, dc_value=I_value)
 elif I_AC_pattern:
     I_Name = I_AC_pattern.group(1)
     I_Np = I_AC_pattern.group(2)
     I_Nn = I_AC_pattern.group(3)
     I_value = unit_transform(I_AC_pattern.group(7))
     mycircuit.add_isource(I_Name, n1=I_Np, n2=I_Nn, ac_value=I_value)
 elif E_pattern:
     E_name = E_pattern.group(1)
示例#7
0
                           n2="n%s" % (3 * i + 2),
                           value=Rs)
    mycircuit.add_capacitor("Cj%s" % (i),
                            n1="n%s" % (3 * i + 2),
                            n2=gnd,
                            value=Cj)
    # mycircuit.add_resistor("Rsub%s" % (i), n1="n%s" %(4*i-1), n2="n%s" %(4*i), value=20000.28)
    # mycircuit.add_capacitor("Cox%s" % (i), n1="n%s" %(4*i), n2=gnd, value=4.09e-20)

mycircuit.add_resistor("Rt", n1="n%s" % (3 * i + 1), n2=gnd, value=25.0)

# voltage_step = time_functions.pulse(v1=0, v2=1, td=500e-9, tr=1e-12, pw=1, tf=1e-12, per=2)
voltage_step = time_functions.pulse(v1=0,
                                    v2=1,
                                    td=0.1e-9,
                                    tr=1e-12,
                                    pw=0.5e-9,
                                    tf=1e-12,
                                    per=1)
# voltage_step = time_functions.pulse(v1=0, v2=1, td=1e-12, tr=1e-15, pw=1, tf=1e-15, per=1)
mycircuit.add_vsource("V1",
                      n1="n1",
                      n2=gnd,
                      dc_value=0,
                      ac_value=0.1,
                      function=voltage_step)

print mycircuit

op_analysis = ahkab.new_op()
ac_analysis = ahkab.new_ac(start=1e8, stop=5e10, points=1e2)