import numpy as np
import math
import matplotlib.pyplot as plt
L=100e-3
C=100e-6
s=np.arange(0.1,7,0.051)
freq=pow(10,s)
XL=(2*np.pi*freq)*L
XC=-1/(2*np.pi*freq*C)
plt.figure(figsize=(15,10))
plt.title('Xl+XC',fontsize=18)
plt.ylabel('XL+XC / $\Omega $',fontsize=18)
plt.xlabel('Frekvenca / Hz',fontsize=18)
plt.tick_params(labelsize=18)
plt.loglog(freq,abs(XL+XC))
plt.grid()
plt.show()
Kako bi se mogla spreminjati induktivnost s frekvenco, da bi imeli konstanten XL+XC?
K=5
L=(1/(2*np.pi*freq*C)-K)/(2*np.pi*freq)
XL=(2*np.pi*freq)*L
XC=-1/(2*np.pi*freq*C)
plt.figure(figsize=(15,10))
plt.title('Želena induktivnost',fontsize=18)
plt.ylabel('L / H',fontsize=18)
plt.xlabel('Frekvenca / Hz',fontsize=18)
plt.tick_params(labelsize=18)
# plt.ylim(0,2)
plt.xlim(1,1000)
plt.loglog(freq,L)
plt.grid()
plt.show()
Rezultat kaže, da bi bio težko doseči tako spreminjanje induktivnosti s frekvenco
K=0
L=60-3
#s=np.arange(0.01,3,0.1)
s=np.linspace(0.01,2,8)
freq=pow(10,s)
C=1 /(K+2*np.pi*freq*L)/(2*np.pi*freq)
XL=(2*np.pi*freq)*L
XC=-1/(2*np.pi*freq*C)
plt.figure(figsize=(15,10))
plt.title('Želena kapacitivnost',fontsize=18)
plt.ylabel('C / H',fontsize=18)
plt.xlabel('Frekvenca / Hz',fontsize=18)
plt.tick_params(labelsize=18)
plt.ylim(1e-2,1e-8)
plt.xlim(1,1000)
plt.loglog(freq,C,'o-')
plt.grid()
C
array([4.24390260e-04, 1.14602470e-04, 3.09472847e-05, 8.35701385e-06, 2.25673047e-06, 6.09408159e-07, 1.64564758e-07, 4.44391156e-08])
Koliko bo XL+XC, če znotraj frekvenčnega območja lahko vključujem 8 kondenzatorjev? Kolikšne vrednosti morajo imeti?
L=60e-3
s=np.linspace(0.01,2,8)
freq1=[]
X=[]
for i in range(1,8):
s1=np.linspace(s[i-1],s[i],10)
freq=pow(10,s1)
freq1=np.concatenate([freq1,freq])
XL=(2*np.pi*freq)*L
XC=-0.0005*1/(2*np.pi*freq*C[i])
X=np.concatenate([X,XL+XC])
print('Copti = ', C[i]/0.005, 'za frekvence med ',freq[0],'in ',freq[9])
print(freq1,X)
plt.figure(figsize=(15,10))
plt.title('Xl+XC',fontsize=18)
plt.ylabel('XL+XC / $\Omega $',fontsize=18)
plt.xlabel('Frekvenca / Hz',fontsize=18)
plt.tick_params(labelsize=18)
#plt.ylim(1,1000)
plt.semilogx(freq1,(X),'.-')
plt.grid()
plt.show()
Copti = 0.0229204940560007 za frekvence med 1.023292992280754 in 1.9691813515545467 Copti = 0.006189456943298567 za frekvence med 1.9691813515545467 in 3.789408531634211 Copti = 0.0016714027699118136 za frekvence med 3.789408531634211 in 7.292176014304688 Copti = 0.00045134609463493377 za frekvence med 7.292176014304688 in 14.032752230245311 Copti = 0.00012188163188993342 za frekvence med 14.032752230245311 in 27.00402935545858 Copti = 3.291295165446955e-05 za frekvence med 27.00402935545858 in 51.96540133151923 Copti = 8.887823126520855e-06 za frekvence med 51.96540133151923 in 100.0 [ 1.02329299 1.10049304 1.18351728 1.27280509 1.36882903 1.47209727 1.58315636 1.70259405 1.83104245 1.96918135 1.96918135 2.11774183 2.27751013 2.44933179 2.63411615 2.83284115 3.04655851 3.27639931 3.52357992 3.78940853 3.78940853 4.07529198 4.38274327 4.71338954 5.06898067 5.45139858 5.86266715 6.30496296 6.78062679 7.29217601 7.29217601 7.84231792 8.43396406 9.07024562 9.75453001 10.49043871 11.2818664 12.13300157 13.04834874 14.03275223 14.03275223 15.09142184 16.22996041 17.45439347 18.77120115 20.18735243 21.71034209 23.34823029 25.10968531 27.00402936 27.00402936 29.04128795 31.23224297 33.58848969 36.12249817 38.84767926 41.77845553 44.93033767 48.32000651 51.96540133 51.96540133 55.88581481 60.10199512 64.63625572 69.51259347 74.75681562 80.39667637 86.46202381 92.98495782 100. ] [-2.92799977e-01 -2.16094107e-01 -1.40531881e-01 -6.57133973e-02 8.75730834e-03 8.32743606e-02 1.58232129e-01 2.34027316e-01 3.11061056e-01 3.89741037e-01 -5.63451776e-01 -4.15842274e-01 -2.70433552e-01 -1.26456056e-01 1.68521903e-02 1.60249625e-01 3.04495155e-01 4.50352177e-01 5.98592617e-01 7.50001014e-01 -1.08428254e+00 -8.00229121e-01 -5.20410783e-01 -2.43346637e-01 3.24296357e-02 3.08377537e-01 5.85957478e-01 8.66638505e-01 1.15190608e+00 1.44326994e+00 -2.08654704e+00 -1.53992676e+00 -1.00145630e+00 -4.68285880e-01 6.24062066e-02 5.93428568e-01 1.12759156e+00 1.66772214e+00 2.21667889e+00 2.77736705e+00 -4.01526206e+00 -2.96336932e+00 -1.92715974e+00 -9.01149357e-01 1.20091840e-01 1.14196860e+00 2.16988906e+00 3.20929330e+00 4.26568222e+00 5.34464660e+00 -7.72679890e+00 -5.70258142e+00 -3.70854394e+00 -1.73413335e+00 2.31099610e-01 2.19755563e+00 4.17564187e+00 6.17582703e+00 8.20869678e+00 1.02850098e+01 -1.48691220e+01 -1.09738042e+01 -7.13656365e+00 -3.33709220e+00 4.44718224e-01 4.22888225e+00 8.03542696e+00 1.18844979e+01 1.57964657e+01 1.97920337e+01]