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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
|
import sys
import traceback
from datetime import datetime
import u3
import u6
import ue9
# MAX_REQUESTS is the number of packets to be read.
MAX_REQUESTS = 75
# SCAN_FREQUENCY is the scan frequency of stream mode in Hz
SCAN_FREQUENCY = 5000
d = None
###############################################################################
# U3
# Uncomment these lines to stream from a U3
###############################################################################
'''
# At high frequencies ( >5 kHz), the number of samples will be MAX_REQUESTS
# times 48 (packets per request) times 25 (samples per packet).
d = u3.U3()
# To learn the if the U3 is an HV
d.configU3()
# For applying the proper calibration to readings.
d.getCalibrationData()
# Set the FIO0 and FIO1 to Analog (d3 = b00000011)
d.configIO(FIOAnalog=3)
print("Configuring U3 stream")
d.streamConfig(NumChannels=2, PChannels=[0, 1], NChannels=[31, 31], Resolution=3, ScanFrequency=SCAN_FREQUENCY)
'''
###############################################################################
# U6
# Uncomment these lines to stream from a U6
###############################################################################
'''
# At high frequencies ( >5 kHz), the number of samples will be MAX_REQUESTS
# times 48 (packets per request) times 25 (samples per packet).
d = u6.U6()
# For applying the proper calibration to readings.
d.getCalibrationData()
print("Configuring U6 stream")
d.streamConfig(NumChannels=2, ChannelNumbers=[0, 1], ChannelOptions=[0, 0], SettlingFactor=1, ResolutionIndex=1, ScanFrequency=SCAN_FREQUENCY)
'''
###############################################################################
# UE9
# Uncomment these lines to stream from a UE9
###############################################################################
# At 96 Hz or higher frequencies, the number of samples will be MAX_REQUESTS
# times 8 (packets per request) times 16 (samples per packet).
# Currently over ethernet packets per request is 1.
d = ue9.UE9()
#d = ue9.UE9(ethernet=True, ipAddress="192.168.1.209") # Over TCP/ethernet connect to UE9 with IP address 192.168.1.209
# For applying the proper calibration to readings.
c=d.getCalibrationData()
# by evmik
# fixing missing slope for gain '0'
c['AINSlopes']['0']= 0.0000775030
print("Configuring UE9 stream")
d.streamConfig(NumChannels=2, ChannelNumbers=[0, 1], ChannelOptions=[0, 0], SettlingTime=0, Resolution=12, ScanFrequency=SCAN_FREQUENCY)
'''
'''
if d is None:
print("""Configure a device first.
Please open streamTest.py in a text editor and uncomment the lines for your device.
Exiting...""")
sys.exit(0)
try:
print("Start stream")
d.streamStart()
start = datetime.now()
print("Start time is %s" % start)
missed = 0
dataCount = 0
packetCount = 0
for r in d.streamData():
if r is not None:
# Our stop condition
if dataCount >= MAX_REQUESTS:
break
if r["errors"] != 0:
print("Errors counted: %s ; %s" % (r["errors"], datetime.now()))
if r["numPackets"] != d.packetsPerRequest:
print("----- UNDERFLOW : %s ; %s" %
(r["numPackets"], datetime.now()))
if r["missed"] != 0:
missed += r['missed']
print("+++ Missed %s" % r["missed"])
# Comment out these prints and do something with r
print("Average of %s AIN0, %s AIN1 readings: %s, %s" %
(len(r["AIN0"]), len(r["AIN1"]), sum(r["AIN0"])/len(r["AIN0"]), sum(r["AIN1"])/len(r["AIN1"])))
dataCount += 1
packetCount += r['numPackets']
else:
# Got no data back from our read.
# This only happens if your stream isn't faster than the USB read
# timeout, ~1 sec.
print("No data ; %s" % datetime.now())
except:
print("".join(i for i in traceback.format_exc()))
finally:
stop = datetime.now()
d.streamStop()
print("Stream stopped.\n")
d.close()
sampleTotal = packetCount * d.streamSamplesPerPacket
scanTotal = sampleTotal / 2 # sampleTotal / NumChannels
print("%s requests with %s packets per request with %s samples per packet = %s samples total." %
(dataCount, (float(packetCount)/dataCount), d.streamSamplesPerPacket, sampleTotal))
print("%s samples were lost due to errors." % missed)
sampleTotal -= missed
print("Adjusted number of samples = %s" % sampleTotal)
runTime = (stop-start).seconds + float((stop-start).microseconds)/1000000
print("The experiment took %s seconds." % runTime)
print("Actual Scan Rate = %s Hz" % SCAN_FREQUENCY)
print("Timed Scan Rate = %s scans / %s seconds = %s Hz" %
(scanTotal, runTime, float(scanTotal)/runTime))
print("Timed Sample Rate = %s samples / %s seconds = %s Hz" %
(sampleTotal, runTime, float(sampleTotal)/runTime))
|
