#!/usr/bin/python3
print('*** EVM(podkrmiljenje oddajnika) - S53MV 19.04.2023 ***')

import math
pmin=0		#začetni P
pmax=16		#končni P
emin=1e-5	#zacetni EVM
emax=100	#koncni EVM

p=pmin
kp=0.1
bb=[]
ee=[]
while p<=pmax:
	bb+=[p]
	pp=10**(p/10)
	e=100*math.sqrt((pp+1)*math.exp(-pp))
	ee+=[e]
	p+=kp

import matplotlib.pyplot as plt
fig=plt.figure(figsize=(10,10))
plt.rcParams.update({'font.size': 16})
plt.semilogy([pmin,pmax],[1,1],'b-')
plt.semilogy(bb,ee,'r-',)
plt.axis([pmin,pmax,emin,emax])
plt.xlabel('Podkrmiljenje OFDM oddajnika P1dB/<P> [dB]')
plt.ylabel('Napaka EVM sqrt(<Pe>/<P>) [%]')
plt.grid(c='k',linestyle='--')
fig.savefig('evm.svg', transparent=True)
plt.show()

print('*** Konec ***')