aboutsummaryrefslogtreecommitdiffstats
path: root/tests/correlator/acproc.py
blob: 5f6ace39d9a58450e507d197391972463d3ff32a (plain)
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
#!/usr/bin/env python3

import numpy as np
import matplotlib.pyplot as plt
import acgen

# Parameters
# samples per symbol
sps = 4
# number of initial bytes (to ignore)
nzeros = 10
# length of the access code in bytes
aclen = 2

# Create samples
print("Modulating symbols")
acgen.main()

# Extract one sequence
print("Extracting symbol sequence")

# raw data
data = np.fromfile("acgen.dat", dtype=np.complex64)
plt.plot(data.real)
plt.plot(data.imag)
plt.title("Raw Data (time domain)")
plt.show()

# take only symbols
symbols = data[1::sps]
plt.plot(symbols.real, symbols.imag)
plt.title("Symbols only (constellation)")
plt.show()

# where ac symbols start, in symbols
ac_start = nzeros * 8
ac_end = ac_start + aclen * 8

ac = symbols[ac_start:ac_end]

fig, (ax1, ax2) = plt.subplots(2, 1)
fig.tight_layout()

ax1.plot(ac.real, ac.imag)
ax1.set_title("Symbols of Access Code (constellation)")

ax2.plot(ac.real, ".-")
ax2.plot(ac.imag, ".-")
ax2.set_title("Symbols of Access Code (time)")
plt.show()

# print the symbols
print(f"Generated {len(ac)} symbols from a {aclen} byte sequence")
print(list(ac))
print("Reversed symbols (for FIR filter)")
print(list(ac[::-1]))