1 |
#!/usr/bin/python |
1 |
#!/usr/bin/python |
2 |
# |
2 |
# |
3 |
# Sample of measiring and correct frequency with ACOUNTER02A |
3 |
# Sample of measuring and frequency correction with ACOUNTER02A |
4 |
|
4 |
|
5 |
import time |
5 |
import time |
6 |
import datetime |
6 |
import datetime |
7 |
import sys |
7 |
import sys |
8 |
from pymlab import config |
8 |
from pymlab import config |
9 |
|
9 |
|
10 |
#### Script Arguments ############################################### |
10 |
#### Script Arguments ############################################### |
11 |
|
11 |
|
12 |
if (len(sys.argv) > 3) or (len(sys.argv) < 2): |
12 |
if (len(sys.argv) > 3) or (len(sys.argv) < 2): |
13 |
sys.stderr.write("Invalid number of arguments.\n") |
13 |
sys.stderr.write("Invalid number of arguments.\n") |
14 |
sys.stderr.write("Usage: %s PORT_ADDRESS [REQUIERED_MHz]\n" % (sys.argv[0], )) |
14 |
sys.stderr.write("Usage: %s PORT_ADDRESS [REQUIERED_MHz]\n" % (sys.argv[0], )) |
15 |
sys.exit(1) |
15 |
sys.exit(1) |
16 |
|
16 |
|
17 |
port = eval(sys.argv[1]) |
17 |
port = eval(sys.argv[1]) |
18 |
#### Sensor Configuration ########################################### |
18 |
#### Sensor Configuration ########################################### |
19 |
|
19 |
|
20 |
cfg = config.Config( |
20 |
cfg = config.Config( |
- |
|
21 |
i2c = { |
21 |
port = port, |
22 |
"port": port, |
- |
|
23 |
}, |
22 |
bus = [ |
24 |
bus = [ |
23 |
{ |
25 |
{ |
24 |
"type": "i2chub", |
26 |
"type": "i2chub", |
25 |
"address": 0x70, |
27 |
"address": 0x70, |
26 |
"children": [ |
28 |
"children": [ |
27 |
{ "name":"counter", "type":"acount02", "channel": 2, }, |
29 |
{ "name":"counter", "type":"acount02", "channel": 2, }, |
28 |
{ "name":"clkgen", "type":"clkgen01", "channel": 5, }, |
30 |
{ "name":"clkgen", "type":"clkgen01", "channel": 5, }, |
29 |
], |
31 |
], |
30 |
}, |
32 |
}, |
31 |
], |
33 |
], |
32 |
) |
34 |
) |
33 |
cfg.initialize() |
35 |
cfg.initialize() |
34 |
|
36 |
|
35 |
print "Frequency counter readout example \r\n" |
37 |
print "RMDS Station frequency management test software \r\n" |
36 |
fcount = cfg.get_device("counter") |
38 |
fcount = cfg.get_device("counter") |
37 |
fgen = cfg.get_device("clkgen") |
39 |
fgen = cfg.get_device("clkgen") |
38 |
time.sleep(0.5) |
40 |
time.sleep(0.5) |
39 |
frequency = fcount.get_freq() |
41 |
frequency = fcount.get_freq() |
40 |
rfreq = fgen.get_rfreq() |
42 |
rfreq = fgen.get_rfreq() |
41 |
hsdiv = fgen.get_hs_div() |
43 |
hsdiv = fgen.get_hs_div() |
42 |
n1 = fgen.get_n1_div() |
44 |
n1 = fgen.get_n1_div() |
43 |
fdco = 0 |
45 |
#fdco = 0 |
44 |
fxtal = 0 |
46 |
#fxtal = 0 |
45 |
regs = [0, 0, 0] |
47 |
#regs = [0, 0, 0] |
46 |
|
48 |
|
47 |
|
49 |
|
48 |
#### Data Logging ################################################### |
50 |
#### Data Logging ################################################### |
49 |
|
51 |
|
50 |
try: |
52 |
try: |
51 |
with open("frequency.log", "a") as f: |
53 |
with open("frequency.log", "a") as f: |
52 |
while True: |
54 |
while True: |
53 |
now = datetime.datetime.now() |
55 |
now = datetime.datetime.now() |
54 |
if (now.second == 15) or (now.second == 35) or (now.second == 55): |
56 |
if (now.second == 15) or (now.second == 35) or (now.second == 55): |
55 |
frequency = fcount.get_freq() |
57 |
frequency = fcount.get_freq() |
56 |
if (len(sys.argv) == 3): |
58 |
if (len(sys.argv) == 3): |
57 |
regs = fgen.set_freq(frequency/1e6, float(eval(sys.argv[2]))) |
59 |
regs = fgen.set_freq(frequency/1e6, float(eval(sys.argv[2]))) |
58 |
now = datetime.datetime.now() |
60 |
now = datetime.datetime.now() |
59 |
|
61 |
|
60 |
rfreq = fgen.get_rfreq() |
62 |
rfreq = fgen.get_rfreq() |
61 |
hsdiv = fgen.get_hs_div() |
63 |
hsdiv = fgen.get_hs_div() |
62 |
n1 = fgen.get_n1_div() |
64 |
n1 = fgen.get_n1_div() |
63 |
fdco = (frequency/1e6) * hsdiv * n1 |
65 |
fdco = (frequency/1e6) * hsdiv * n1 |
64 |
fxtal = fdco / rfreq |
66 |
fxtal = fdco / rfreq |
65 |
|
67 |
|
66 |
sys.stdout.write("frequency: " + str(frequency) + " Hz Time: " + str(now.second)) |
68 |
sys.stdout.write("frequency: " + str(frequency) + " Hz Time: " + str(now.second)) |
67 |
sys.stdout.write(" RFREQ: " + str(rfreq) + " HSDIV: " + str(hsdiv) + " N1: " + str(n1)) |
69 |
sys.stdout.write(" RFREQ: " + str(rfreq) + " HSDIV: " + str(hsdiv) + " N1: " + str(n1)) |
68 |
sys.stdout.write(" fdco: " + str(fdco) + " fxtal: " + str(fxtal) + "\r") |
70 |
sys.stdout.write(" fdco: " + str(fdco) + " fxtal: " + str(fxtal) + "\r") |
69 |
f.write("%d\t%s\t%.3f\n" % (time.time(), datetime.datetime.now().isoformat(), frequency)) |
71 |
f.write("%d\t%s\t%.3f\n" % (time.time(), datetime.datetime.now().isoformat(), frequency)) |
70 |
|
72 |
|
71 |
sys.stdout.flush() |
73 |
sys.stdout.flush() |
72 |
time.sleep(0.9) |
74 |
time.sleep(0.9) |
73 |
except KeyboardInterrupt: |
75 |
except KeyboardInterrupt: |
74 |
sys.stdout.write("\r\n") |
76 |
sys.stdout.write("\r\n") |
75 |
sys.exit(0) |
77 |
sys.exit(0) |
76 |
f.close() |
78 |
f.close() |