Phasors
In the circuit given below, $R_1=R_2=R$, $C_1=C_2=C$, $V_1=10\angle {0}$, $V_2=V\angle {\phi}$.- Derive an expression for $\phi$ in terms of $R$, $C$, and $f$ (frequency).
- Find the condition for $\phi=-45^{\circ}$.
- Plot $\phi$ versus $T\,(=R\,C)$, with $\omega = 2\,\pi \times 50$ rad/s. (Keep $C=1\,$mF, and vary $R$.)
- Plot $\phi$ versus $\omega$ for $R=1\,\Omega$, $C=1\,$mF.
In [1]:
from IPython.display import Image
Image(filename =r'phasor_8_fig_1.png', width=250)
Out[1]:
In [2]:
# run this cell to view the circuit file.
%pycat phasor_8_orig.in
In [3]:
import numpy as np
import gseim_calc as calc
import os
import dos_unix
import matplotlib.pyplot as plt
from matplotlib.ticker import (MultipleLocator, AutoMinorLocator)
from setsize import set_size
c = 1.0e-3
s_C = "%11.4E"%c
s_f_hz = '50'
# first two args in logspace are powers of 10
# we will vary R from 0.1 ohm to 10k.
r_arr = np.logspace(-2.0, 3.0, 100, endpoint=True)
l_phi = []
l_T = []
slv = calc.slv("phasor_8_orig.in")
i_slv = 0
i_out = 0
filename = slv.l_filename_all[i_slv][i_out]
col_phi = slv.get_index(i_slv,i_out,"phase_of_V2_ac")
for r in r_arr:
s_R = "%11.4E"%r
l = [
('$R', s_R),
('$C', s_C),
('$frequency', s_f_hz),
]
print('R:', s_R)
calc.replace_strings_1("phasor_8_orig.in", "phasor_8.in", l)
# uncomment for windows:
#dos_unix.d2u("phasor_8.in")
os.system('run_gseim phasor_8.in')
print('r =', "%9.4f"%r, 'ran successfully')
u = np.loadtxt(filename)
phi = u[col_phi]
l_phi.append(phi)
l_T.append(r*c)
R: 1.0000E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0100 ran successfully R: 1.1233E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0112 ran successfully R: 1.2619E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0126 ran successfully R: 1.4175E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0142 ran successfully R: 1.5923E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0159 ran successfully R: 1.7886E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0179 ran successfully R: 2.0092E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0201 ran successfully R: 2.2570E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0226 ran successfully R: 2.5354E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0254 ran successfully R: 2.8480E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0285 ran successfully R: 3.1993E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0320 ran successfully R: 3.5938E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0359 ran successfully R: 4.0370E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0404 ran successfully R: 4.5349E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0453 ran successfully R: 5.0941E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0509 ran successfully R: 5.7224E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0572 ran successfully R: 6.4281E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0643 ran successfully R: 7.2208E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0722 ran successfully R: 8.1113E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0811 ran successfully R: 9.1116E-02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.0911 ran successfully R: 1.0235E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.1024 ran successfully R: 1.1498E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.1150 ran successfully R: 1.2915E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.1292 ran successfully R: 1.4508E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.1451 ran successfully R: 1.6298E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.1630 ran successfully R: 1.8307E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.1831 ran successfully R: 2.0565E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.2057 ran successfully R: 2.3101E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.2310 ran successfully R: 2.5950E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.2595 ran successfully R: 2.9151E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.2915 ran successfully R: 3.2745E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.3275 ran successfully R: 3.6784E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.3678 ran successfully R: 4.1320E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.4132 ran successfully R: 4.6416E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.4642 ran successfully R: 5.2140E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.5214 ran successfully R: 5.8570E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.5857 ran successfully R: 6.5793E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.6579 ran successfully R: 7.3907E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.7391 ran successfully R: 8.3022E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.8302 ran successfully R: 9.3260E-01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 0.9326 ran successfully R: 1.0476E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 1.0476 ran successfully R: 1.1768E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 1.1768 ran successfully R: 1.3219E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 1.3219 ran successfully R: 1.4850E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 1.4850 ran successfully R: 1.6681E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 1.6681 ran successfully R: 1.8738E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 1.8738 ran successfully R: 2.1049E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 2.1049 ran successfully R: 2.3645E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 2.3645 ran successfully R: 2.6561E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 2.6561 ran successfully R: 2.9836E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 2.9836 ran successfully R: 3.3516E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 3.3516 ran successfully R: 3.7649E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 3.7649 ran successfully R: 4.2292E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 4.2292 ran successfully R: 4.7508E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 4.7508 ran successfully R: 5.3367E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 5.3367 ran successfully R: 5.9948E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 5.9948 ran successfully R: 6.7342E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 6.7342 ran successfully R: 7.5646E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 7.5646 ran successfully R: 8.4975E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 8.4975 ran successfully R: 9.5455E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 9.5455 ran successfully R: 1.0723E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 10.7227 ran successfully R: 1.2045E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 12.0450 ran successfully R: 1.3530E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 13.5305 ran successfully R: 1.5199E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 15.1991 ran successfully R: 1.7074E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 17.0735 ran successfully R: 1.9179E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 19.1791 ran successfully R: 2.1544E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 21.5443 ran successfully R: 2.4201E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 24.2013 ran successfully R: 2.7186E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 27.1859 ran successfully R: 3.0539E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 30.5386 ran successfully R: 3.4305E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 34.3047 ran successfully R: 3.8535E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 38.5353 ran successfully R: 4.3288E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 43.2876 ran successfully R: 4.8626E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 48.6260 ran successfully R: 5.4623E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 54.6228 ran successfully R: 6.1359E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 61.3591 ran successfully R: 6.8926E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 68.9261 ran successfully R: 7.7426E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 77.4264 ran successfully R: 8.6975E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 86.9749 ran successfully R: 9.7701E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 97.7010 ran successfully R: 1.0975E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 109.7499 ran successfully R: 1.2328E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 123.2847 ran successfully R: 1.3849E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 138.4886 ran successfully R: 1.5557E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 155.5676 ran successfully R: 1.7475E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 174.7528 ran successfully R: 1.9630E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 196.3041 ran successfully R: 2.2051E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 220.5131 ran successfully R: 2.4771E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 247.7076 ran successfully R: 2.7826E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 278.2559 ran successfully R: 3.1257E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 312.5716 ran successfully R: 3.5112E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 351.1192 ran successfully R: 3.9442E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 394.4206 ran successfully R: 4.4306E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 443.0621 ran successfully R: 4.9770E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 497.7024 ran successfully R: 5.5908E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 559.0810 ran successfully R: 6.2803E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 628.0291 ran successfully R: 7.0548E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 705.4802 ran successfully R: 7.9248E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 792.4829 ran successfully R: 8.9022E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 890.2151 ran successfully R: 1.0000E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. r = 1000.0000 ran successfully
In [4]:
fig, ax = plt.subplots()
set_size(5.5, 3, ax)
plt.xscale('log')
plt.grid(color='#CCCCCC', linestyle='solid', linewidth=0.5)
ax.set_xlabel('$T$',fontsize=12)
ax.set_ylabel('$\phi$',fontsize=12)
ax.plot(l_T, l_phi, color="blue", linewidth=1.0)
plt.tight_layout()
plt.show()
In [5]:
import numpy as np
import gseim_calc as calc
import os
import dos_unix
import matplotlib.pyplot as plt
from matplotlib.ticker import (MultipleLocator, AutoMinorLocator)
from setsize import set_size
s_C = "1m"
s_R = "1"
# first two args in logspace are powers of 10
f_hz_arr = np.logspace(0.0, 4.0, 100, endpoint=True)
l_phi = []
l_f_hz = []
slv = calc.slv("phasor_8_orig.in")
i_slv = 0
i_out = 0
filename = slv.l_filename_all[i_slv][i_out]
col_phi = slv.get_index(i_slv,i_out,"phase_of_V2_ac")
for f_hz in f_hz_arr:
s_f_hz = "%11.4E"%f_hz
l = [
('$R', s_R),
('$C', s_C),
('$frequency', s_f_hz),
]
print('f_hz:', s_f_hz)
calc.replace_strings_1("phasor_8_orig.in", "phasor_8.in", l)
# uncomment for windows:
#dos_unix.d2u("phasor_8.in")
os.system('run_gseim phasor_8.in')
u = np.loadtxt(filename)
phi = u[col_phi]
l_phi.append(phi)
l_f_hz.append(f_hz)
f_hz: 1.0000E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.0975E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.2045E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.3219E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.4508E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.5923E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.7475E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.9179E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.1049E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.3101E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.5354E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.7826E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.0539E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.3516E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.6784E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.0370E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.4306E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.8626E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.3367E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.8570E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 6.4281E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 7.0548E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 7.7426E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 8.4975E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 9.3260E+00 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.0235E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.1233E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.2328E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.3530E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.4850E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.6298E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.7886E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.9630E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.1544E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.3645E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.5950E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.8480E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.1257E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.4305E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.7649E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.1320E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.5349E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.9770E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.4623E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.9948E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 6.5793E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 7.2208E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 7.9248E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 8.6975E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 9.5455E+01 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.0476E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.1498E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.2619E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.3849E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.5199E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.6681E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.8307E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.0092E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.2051E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.4201E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.6561E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.9151E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.1993E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.5112E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.8535E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.2292E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.6416E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.0941E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.5908E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 6.1359E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 6.7342E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 7.3907E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 8.1113E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 8.9022E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 9.7701E+02 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.0723E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.1768E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.2915E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.4175E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.5557E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.7074E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.8738E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.0565E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.2570E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.4771E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.7186E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 2.9836E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.2745E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.5938E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 3.9442E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.3288E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 4.7508E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.2140E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 5.7224E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 6.2803E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 6.8926E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 7.5646E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 8.3022E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 9.1116E+03 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed. f_hz: 1.0000E+04 Circuit: filename = phasor_8.in main: i_solve = 0 GSEIM: Program completed.
In [6]:
fig, ax = plt.subplots()
set_size(5.5, 3, ax)
plt.xscale('log')
plt.grid(color='#CCCCCC', linestyle='solid', linewidth=0.5)
ax.set_xlabel('$f$',fontsize=12)
ax.set_ylabel('$\phi$',fontsize=12)
ax.plot(l_f_hz, l_phi, color="blue", linewidth=1.0)
plt.tight_layout()
plt.show()
This notebook was contributed by Prof. Nakul Narayanan K, Govt. Engineering College, Thrissur. He may be contacted at nakul@gectcr.ac.in.
In [ ]: