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()




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()

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

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()