Commit 27608a2ff3e4db84f1b59d7b91e5098297b16d9a
1 parent
9afb8a614f
Exists in
master
replace 4 spaces by tabs
Showing 1 changed file with 7 additions and 7 deletions Inline Diff
ringdown.py
| import matplotlib.pyplot as plt | 1 | 1 | import matplotlib.pyplot as plt | |
| import numpy as np | 2 | 2 | import numpy as np | |
| import csv, glob, sys | 3 | 3 | import csv, glob, sys | |
| from scipy.optimize import curve_fit | 4 | 4 | from scipy.optimize import curve_fit | |
| 5 | 5 | |||
| list_files = (glob.glob(sys.argv[1])) | 6 | 6 | list_files = (glob.glob(sys.argv[1])) | |
| 7 | 7 | |||
| data = [] | 8 | 8 | data = [] | |
| 9 | 9 | |||
| for f in list_files: | 10 | 10 | for f in list_files: | |
| data_iter = csv.reader(open(f, 'r'), delimiter = ',', quotechar = '"') | 11 | 11 | data_iter = csv.reader(open(f, 'r'), delimiter = ',', quotechar = '"') | |
| for i in range(2): | 12 | 12 | for i in range(2): | |
| data_iter.next() | 13 | 13 | data_iter.next() | |
| #temp_data = [value for value in data_iter] #if ends with a number | 14 | 14 | #temp_data = [value for value in data_iter] #if ends with a number | |
| temp_data = [value[:-1] for value in data_iter]#if ends with a comma | 15 | 15 | temp_data = [value[:-1] for value in data_iter]#if ends with a comma | |
| data.extend(temp_data) | 16 | 16 | data.extend(temp_data) | |
| 17 | 17 | |||
| data = np.asarray(data, dtype = float) | 18 | 18 | data = np.asarray(data, dtype = float) | |
| 19 | 19 | |||
| del(temp_data, list_files, value, f) | 20 | 20 | del(temp_data, list_files, value, f) | |
| 21 | 21 | |||
| plt.subplot(111) | 22 | 22 | plt.subplot(111) | |
| plt.clf() | 23 | 23 | plt.clf() | |
| 24 | 24 | |||
| plt.plot(data[data[:,0]>=0,0], data[data[:,0]>=0,1], label ='mes') | 25 | 25 | plt.plot(data[data[:,0]>=0,0], data[data[:,0]>=0,1], label ='mes') | |
| 26 | 26 | |||
| def func(t, tau, A, w , a, b): | 27 | 27 | def func(t, tau, A, w , a, b): | |
| return A * np.exp(-t/tau) * np.sin((w+a*t)*t) + b | 28 | 28 | return A * np.exp(-t/tau) * np.sin((w+a*t)*t) + b | |
| 29 | 29 | |||
| popt, pcov = curve_fit(func, data[data[:,0]>=0,0], data[data[:,0]>=0,1], p0 = [1e-5, 1e-2, 1e5, 1e10, 1e-4], maxfev=10000) | 30 | 30 | popt, pcov = curve_fit(func, data[data[:,0]>=0,0], data[data[:,0]>=0,1], p0 = [1e-5, 1e-2, 1e5, 1e10, 1e-4], maxfev=10000) | |
| yfit = func(data[data[:,0]>=0,0], *popt) | 31 | 31 | yfit = func(data[data[:,0]>=0,0], *popt) | |
| 32 | 32 | |||
| tau = float(popt[0]) | 33 | 33 | tau = float(popt[0]) | |
| Q = np.pi*299792458*tau/(2.*140e-3) | 34 | 34 | Q = np.pi*299792458*tau/(2.*140e-3) | |
| 35 | 35 | |||
| #print(np.sqrt(np.diag(pcov))) | 36 | 36 | #print(np.sqrt(np.diag(pcov))) | |
| print(tau, Q) | 37 | 37 | print(tau, Q) | |
| 38 | 38 | |||
| plt.plot(data[data[:,0]>=0,0], yfit, label ='fit') | 39 | 39 | plt.plot(data[data[:,0]>=0,0], yfit, label ='fit') |