Commit ee82fa8490e3105729eaa4ca5549ba0c896abddc
1 parent
84d931fb3d
Exists in
master
.
Showing 1 changed file with 0 additions and 45 deletions Inline Diff
RingDown.py
| import matplotlib.pyplot as plt | 1 | File was deleted | ||
| import numpy as np | 2 | |||
| import csv, glob, sys | 3 | |||
| from scipy.optimize import curve_fit | 4 | |||
| 5 | ||||
| list_files = (glob.glob(sys.argv[1])) | 6 | |||
| 7 | ||||
| data = [] | 8 | |||
| 9 | ||||
| for f in list_files: | 10 | |||
| data_iter = csv.reader(open(f, 'r'), delimiter = ',', quotechar = '"') | 11 | |||
| for i in range(2): | 12 | |||
| data_iter.next() | 13 | |||
| #temp_data = [value for value in data_iter] #if ends with a number | 14 | |||
| temp_data = [value[:-1] for value in data_iter]#if ends with a comma | 15 | |||
| data.extend(temp_data) | 16 | |||
| 17 | ||||
| data = np.asarray(data, dtype = float) | 18 | |||
| 19 | ||||
| del(temp_data, list_files, value, f) | 20 | |||
| 21 | ||||
| plt.subplot(111) | 22 | |||
| plt.clf() | 23 | |||
| 24 | ||||
| plt.plot(data[data[:,0]>=0,0], data[data[:,0]>=0,1], label ='mes') | 25 | |||
| 26 | ||||
| def func(t, tau, A, w , a, b): | 27 | |||
| return A * np.exp(-t/tau) * np.sin((w+a*t)*t) + b | 28 | |||
| 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 | |||
| yfit = func(data[data[:,0]>=0,0], *popt) | 31 | |||
| 32 |