hcui91 / e7-lrc

This repository contained codes for the 2nd year experiment "Study of an Oscillating Electrical System Driven into Resonance". It would import data from csv file and output plotting with phase difference and angle. Please be aware of possible plagiarism check if using the code without any modification.

E7-LRC

This repository contained codes for the 2nd year experiment Study of an Oscillating Electrical System Driven into Resonance. It would import data from csv file and output plotting with phase difference and angle.

Please be aware of possible plagiarism check if using the code without any modification.

From experiment to this code

From your RLC circuit

You will need to change the variables from line 13

##### Theoritical value #####

L=18e-3         #Inductance of the inductor in Henry (H)

C=1800e-12      #Capacitance of the capacitor in Faraday (F)

R=820           #Resistance of the resistor in Ohm

Note: 18e-3 means 18*10^{-3}

From your experiment data

I would recommand the arrangement of dataset below in a comma delimiated file.

#frequency(Hz)	peak-to-peak voltage(volt)	unc voltage (volt)	y0	y1
10000	2.30	0.05	2.30	0.50

But feel free to change the order of columns, and update it in the code line 28&29. Change the order of variables.

freq, volt, delta_volt, y0, y1, delta_y = np.loadtxt("low_damp.csv",delimiter=",",unpack=True)   # Change file name and delimiter if required
freq_1, volt_1, delta_volt_1, y0_1, y1_1 = np.loadtxt("high_damp.csv",delimiter=",",unpack=True)

Tricks in coding

Better way to display axis

In the plot of phase angle, it is more convenient to display y-axis in unit of pi. On line 93, I used plt.yticks to modify the display of axis.

plt.yticks(np.arange(0,np.pi+0.1,step=0.25*np.pi), ('0',r'$\frac{1}{4} \pi$',r'$\frac{1}{2} \pi$',r'$\frac{3}{4} \pi$',r'$\pi$'))

It will devide pi into 4 ticks with step 0.25*pi, and display them on the axis. Text type of display is used to render \pi in Tex (matplotlib supports inline math with r'' in text contents).