from IPython.display import Image
Image(filename =r'switching_1_fig_1.png', width=320)

import numpy as np
import matplotlib.pyplot as plt 
import gseim_calc as calc
from setsize import set_size

I1 = 5.0
V1 = 15.0 # to be changed by user

t_start = 0.0
t1 = 2e-6
t_on = 1e-6
t2 = t1 + t_on
t_end = 5e-6

n_div = 2000
t = np.linspace(t_start, t_end, (n_div+1))

slope_V = V1/t_on
slope_I = I1/t_on

l_I = []
l_V = []
for i, t0 in enumerate(t):
    if t0 < t1:
        I = 0.0
        V = V1
    elif t0 < t2:
        delt = t0-t1
        I = slope_I*delt
        V = V1 - slope_V*delt
    else:
        I = I1
        V = 0.0
    l_I.append(I)
    l_V.append(V)

np_I = np.array(l_I)
np_V = np.array(l_V)
P = np_I*np_V

Pmax = np.max(P)

plot_min_P, plot_max_P = calc.delta_plot_1(P   , 0.2)
plot_min_V, plot_max_V = calc.delta_plot_1(np_V, 0.2)
plot_min_I, plot_max_I = calc.delta_plot_1(np_I, 0.2)

print('peak power dissipation:', "%11.4E"%Pmax)

color1='blue'
color2='red'
color3='green'

fig, ax = plt.subplots(3, sharex=False)
plt.subplots_adjust(wspace=0, hspace=0.0)

set_size(6.5, 5, ax[0])

for i in range(3):
    ax[i].set_xlim(left=0.0, right=t_end*1e6)
    ax[i].grid(color='#CCCCCC', linestyle='solid', linewidth=0.5)

ax[0].set_ylim(bottom=plot_min_V, top=plot_max_V)
ax[1].set_ylim(bottom=plot_min_I, top=plot_max_I)
ax[2].set_ylim(bottom=plot_min_P, top=plot_max_P)

ax[0].set_ylabel(r'$v$', fontsize=12)
ax[1].set_ylabel(r'$i$', fontsize=12)
ax[2].set_ylabel(r'$p$', fontsize=12)

ax[0].tick_params(labelbottom=False)
ax[1].tick_params(labelbottom=False)

ax[0].plot(t*1e6, np_V, color=color1, linewidth=1.0, label="$v$")
ax[1].plot(t*1e6, np_I, color=color2, linewidth=1.0, label="$i$")
ax[2].plot(t*1e6, P   , color=color3, linewidth=1.0, label="$p$")

ax[2].set_xlabel(r'time ($\mu$sec)', fontsize=12)

#plt.tight_layout()
plt.show()

peak power dissipation:  1.8750E+01

Switching