1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
#!/usr/bin/env python
import numpy as np
import matplotlib.pyplot as plt
dat = np.loadtxt("tone.csv", dtype=np.int32, delimiter=",")
# columns :"freq, threshold, num_samples, detector_output, avg_rms"
fig, ax = plt.subplots()
plt.scatter(dat[...,0], dat[...,1] / dat[...,4], c=dat[...,3].astype(np.float32))
plt.xlabel("freq")
plt.ylabel("threshold over rms")
plt.show()
if 0:
thresholds = [200, 2000, 4000, 8000, 12000]
fig, ax = plt.subplots()
for th in thresholds:
dat_th = dat[dat[...,1] == th]
dat_th[...,1]
ax.plot(dat_th[...,0], dat_th[...,3], label="{}".format(th))
legend = ax.legend(loc='upper left', shadow=True)
plt.show()
|
