Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import numpy as np
- import pylab as plt
- data = np.loadtxt("Measurements.csv",delimiter=",")
- phis = data[:,0]
- E0 = data[:,12]
- E1 = data[:,13]
- E2 = data[:,14]
- E3 = data[:,15]
- E4 = data[:,16]
- phipts = np.linspace(0,2*np.pi,300)
- ticks = [r'$0$', r'$\frac{\pi}{4}$', r'$\frac{\pi}{2}$', r'$\frac{3\pi}{4}$', r'$\pi$',
- r'$\frac{5\pi}{4}$', r'$\frac{3\pi}{2}$', r'$\frac{7\pi}{4}$', r'$2\pi$']
- plt.figure(figsize=(11,9))
- plt.plot(phis,E0,'o')
- plt.plot(phipts, 0.5*(1-np.abs(np.cos(phipts)) ) )
- plt.xlabel(r'$\varphi$',fontsize=20)
- plt.ylabel(r'$E_0$',fontsize=20)
- plt.xticks(np.arange(0,2*np.pi+0.1,step=np.pi/4),ticks)
- plt.tick_params(axis='both', which='major', labelsize=24, width=2.5, length=10)
- plt.savefig('E0.pdf')
- plt.figure(figsize=(11,9))
- # plt.tick_params(axis='both', which='major', labelsize=14, width=2.5, length=10)
- plt.plot(phis,E1,'o')
- plt.plot(phipts, 0.5*(1-np.abs(np.cos(phipts))**3 ) )
- plt.xlabel(r'$\varphi$',fontsize=20)
- plt.ylabel(r'$E_1$',fontsize=20)
- plt.xticks(np.arange(0,2*np.pi+0.1,step=np.pi/4),ticks)
- plt.tick_params(axis='both', which='major', labelsize=24, width=2.5, length=10)
- plt.savefig('E1.pdf')
- plt.figure(figsize=(11,9))
- plt.plot(phis,E2,'o')
- plt.plot(phipts, 0.5*(1-np.abs(np.cos(phipts))**3 ) )
- plt.xlabel(r'$\varphi$',fontsize=20)
- plt.ylabel(r'$E_2$',fontsize=20)
- plt.xticks(np.arange(0,2*np.pi+0.1,step=np.pi/4),ticks)
- plt.tick_params(axis='both', which='major', labelsize=24, width=2.5, length=10)
- plt.savefig('E2.pdf')
- plt.figure(figsize=(11,9))
- plt.plot(phis,E3,'o')
- plt.plot(phipts, 0.5*(1-np.abs(np.cos(phipts))**2 ) )
- plt.xlabel(r'$\varphi$',fontsize=20)
- plt.ylabel(r'$E_3$',fontsize=20)
- plt.xticks(np.arange(0,2*np.pi+0.1,step=np.pi/4),ticks)
- plt.tick_params(axis='both', which='major', labelsize=24, width=2.5, length=10)
- plt.savefig('E3.pdf')
- plt.figure(figsize=(11,9))
- plt.plot(phis,E4,'o')
- plt.plot(phipts, 0.5*(1-np.abs(np.cos(phipts)) ) )
- plt.xlabel(r'$\varphi$',fontsize=20)
- plt.ylabel(r'$E_4$',fontsize=20)
- plt.xticks(np.arange(0,2*np.pi+0.1,step=np.pi/4),ticks)
- plt.tick_params(axis='both', which='major', labelsize=24, width=2.5, length=10)
- plt.savefig('E4.pdf')
- # plt.show()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement