blob: b2fbedc465be2fb39b7f4608234202fe7a8e6b15 (
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
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
|
#!/usr/bin/env python3
#
# Plot fsk_demod statistic outputs.
#
# Copyright (C) 2018 Mark Jessop <[email protected]>
# Released under GNU GPL v3 or later
#
import json
import os
import sys
import time
import subprocess
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
_filename = sys.argv[1]
if len(sys.argv)>2:
_sps = int(sys.argv[2])
_x_label = "Seconds"
else:
_sps = 1
_x_label = "Sample Number"
_ebno = []
_fest1 = []
_fest2 = []
_fest3 = []
_fest4 = []
_ppm = []
_time = []
_f = open(_filename,'r')
for _line in _f:
if _line[0] != '{':
continue
try:
_data = json.loads(_line)
except Exception as e:
#print("Line parsing error: %s" % str(e))
continue
_ebno.append(_data['EbNodB'])
_fest1.append(_data['f1_est'])
_fest2.append(_data['f2_est'])
if 'f3_est' in _data:
_fest3.append(_data['f3_est'])
_fest4.append(_data['f4_est'])
_ppm.append(_data['ppm'])
if _time == []:
_time = [0]
else:
_time.append(_time[-1]+1.0/_sps)
_ebno_max = pd.Series(_ebno).rolling(10).max().dropna().tolist()
plt.figure()
plt.plot(_time[:len(_ebno_max)],_ebno_max)
plt.xlabel(_x_label)
plt.ylabel("Eb/N0 (dB)")
plt.title("Eb/N0")
plt.figure()
plt.plot(_time,_fest1, label="f1 est")
plt.plot(_time,_fest2, label="f2 est")
if len(_fest3) > 0:
plt.plot(_time,_fest3, label="f3 est")
plt.plot(_time,_fest4, label="f4 est")
plt.legend()
plt.xlabel(_x_label)
plt.ylabel("Frequency (Hz)")
plt.title("Frequency Estimator Outputs")
plt.figure()
plt.plot(_time,_ppm)
plt.xlabel(_x_label)
plt.ylabel("PPM")
plt.title("Demod PPM Estimate")
plt.show()
|
