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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
|
#
# Copyright (c) 2024 Analog Devices Inc.
#
# This file is part of libm2k
# (see http://www.github.com/analogdevicesinc/libm2k).
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation, either version 2.1 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
from scipy.signal import find_peaks
import matplotlib.pyplot as plt
import libm2k
from helpers import get_result_files, plot_to_file
from open_context import ctx_timeout, ctx, ain, aout, trig
gen_reports = True
def test_dac_artifact(ain, aout, trig, chn):
# check for DAC artifacts at slow sample rates
file_name, dir_name, csv_path = get_result_files(gen_reports)
buffer_size = 4000
sampling_frequency_in = 1000000
sampling_frequency_out = 750000
# setup adc
ain.setSampleRate(sampling_frequency_in)
ain.setRange(chn, libm2k.PLUS_MINUS_2_5V)
ain.enableChannel(chn, True)
# setup trigger
trig.setAnalogSource(chn)
trig.setAnalogCondition(chn, libm2k.RISING_EDGE_ANALOG)
trig.setAnalogLevel(chn, 0.5) # Set trigger level at 0.5
trig.setAnalogDelay(0) # Trigger is centered
trig.setAnalogMode(chn, libm2k.ANALOG)
# setup dac
aout.stop()
aout.setCyclic(False)
aout.enableChannel(chn, True)
aout.setSampleRate(chn, sampling_frequency_out)
buffer = generate_clock_signal()
ain.startAcquisition(buffer_size)
aout.push(chn, buffer)
data = ain.getSamples(buffer_size)[chn]
peaks = find_peaks(data, height=2.5)
if peaks == 0:
dac_artifact = 0
else:
dac_artifact = 1
if gen_reports:
plot_to_file('Spikes on chn '+ str(chn), data, dir_name, 'spikes_ch'+str(chn)+'.png')
return dac_artifact
def generate_clock_signal():
buffer = []
for i in range(256):
buffer.append(1)
for i in range(256):
buffer.append(0)
return buffer
|
