|
|
|
@ -1,25 +1,13 @@
|
|
|
|
#!/usr/bin/env python3
|
|
|
|
|
|
|
|
"""
|
|
|
|
|
|
|
|
Examples:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
HackRF:
|
|
|
|
|
|
|
|
./.venv-sdr/bin/python read_energy_wide.py \
|
|
|
|
|
|
|
|
--backend hackrf \
|
|
|
|
|
|
|
|
--serial 0000000000000000a18c63dc2a83b813 \
|
|
|
|
|
|
|
|
--sample-rate 20000000 \
|
|
|
|
|
|
|
|
--base 6000 \
|
|
|
|
|
|
|
|
--roof 5700 \
|
|
|
|
|
|
|
|
--step 20
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Neptune / Pluto over IIO:
|
|
|
|
|
|
|
|
./.venv-sdr/bin/python read_energy_wide.py \
|
|
|
|
|
|
|
|
--backend iio \
|
|
|
|
|
|
|
|
--uri ip:192.168.2.1 \
|
|
|
|
|
|
|
|
--base 2402 \
|
|
|
|
|
|
|
|
--roof 2398 \
|
|
|
|
|
|
|
|
--step 1
|
|
|
|
|
|
|
|
"""
|
|
|
|
"""
|
|
|
|
|
|
|
|
./.venv-sdr/bin/python read_energy_wide.py \
|
|
|
|
|
|
|
|
--serial 0000000000000000a18c63dc2a83b813 \
|
|
|
|
|
|
|
|
--sample-rate 20000000 \
|
|
|
|
|
|
|
|
--base 6000 \
|
|
|
|
|
|
|
|
--roof 5700 \
|
|
|
|
|
|
|
|
--step 20
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
"""
|
|
|
|
|
|
|
|
#!/usr/bin/env python3
|
|
|
|
import argparse
|
|
|
|
import argparse
|
|
|
|
import math
|
|
|
|
import math
|
|
|
|
import re
|
|
|
|
import re
|
|
|
|
@ -28,7 +16,7 @@ import subprocess
|
|
|
|
import sys
|
|
|
|
import sys
|
|
|
|
import time
|
|
|
|
import time
|
|
|
|
from dataclasses import dataclass
|
|
|
|
from dataclasses import dataclass
|
|
|
|
from typing import Dict, List, Optional, Sequence, Tuple
|
|
|
|
from typing import Dict, List, Optional, Tuple
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
import numpy as np
|
|
|
|
import numpy as np
|
|
|
|
@ -45,7 +33,6 @@ except Exception as exc:
|
|
|
|
sys.exit(1)
|
|
|
|
sys.exit(1)
|
|
|
|
|
|
|
|
|
|
|
|
EPS = 1e-20
|
|
|
|
EPS = 1e-20
|
|
|
|
IIO_MIN_SAMPLE_RATE = 2083333
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@dataclass
|
|
|
|
@dataclass
|
|
|
|
@ -66,10 +53,10 @@ class ScanWindow:
|
|
|
|
pass_no: int = 0
|
|
|
|
pass_no: int = 0
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
class HackRfWideProbe(gr.top_block):
|
|
|
|
class WideProbeTop(gr.top_block):
|
|
|
|
def __init__(
|
|
|
|
def __init__(
|
|
|
|
self,
|
|
|
|
self,
|
|
|
|
serial: str,
|
|
|
|
index: int,
|
|
|
|
center_freq_hz: float,
|
|
|
|
center_freq_hz: float,
|
|
|
|
sample_rate: float,
|
|
|
|
sample_rate: float,
|
|
|
|
vec_len: int,
|
|
|
|
vec_len: int,
|
|
|
|
@ -80,7 +67,7 @@ class HackRfWideProbe(gr.top_block):
|
|
|
|
super().__init__("hackrf_energy_wide_probe")
|
|
|
|
super().__init__("hackrf_energy_wide_probe")
|
|
|
|
self.probe = blocks.probe_signal_vc(vec_len)
|
|
|
|
self.probe = blocks.probe_signal_vc(vec_len)
|
|
|
|
self.stream_to_vec = blocks.stream_to_vector(gr.sizeof_gr_complex * 1, vec_len)
|
|
|
|
self.stream_to_vec = blocks.stream_to_vector(gr.sizeof_gr_complex * 1, vec_len)
|
|
|
|
self.src = osmosdr.source(args=f"numchan=1 hackrf={serial}")
|
|
|
|
self.src = osmosdr.source(args=f"numchan=1 hackrf={index}")
|
|
|
|
self.src.set_time_unknown_pps(osmosdr.time_spec_t())
|
|
|
|
self.src.set_time_unknown_pps(osmosdr.time_spec_t())
|
|
|
|
self.src.set_sample_rate(sample_rate)
|
|
|
|
self.src.set_sample_rate(sample_rate)
|
|
|
|
self.src.set_center_freq(center_freq_hz, 0)
|
|
|
|
self.src.set_center_freq(center_freq_hz, 0)
|
|
|
|
@ -95,8 +82,8 @@ class HackRfWideProbe(gr.top_block):
|
|
|
|
for fn, val in (("set_gain", gain), ("set_if_gain", if_gain), ("set_bb_gain", bb_gain)):
|
|
|
|
for fn, val in (("set_gain", gain), ("set_if_gain", if_gain), ("set_bb_gain", bb_gain)):
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
getattr(self.src, fn)(val, 0)
|
|
|
|
getattr(self.src, fn)(val, 0)
|
|
|
|
except Exception as exc:
|
|
|
|
except Exception:
|
|
|
|
raise RuntimeError(f"failed to set {fn}={val}: {exc}") from exc
|
|
|
|
pass
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
self.src.set_bandwidth(0, 0)
|
|
|
|
self.src.set_bandwidth(0, 0)
|
|
|
|
except Exception:
|
|
|
|
except Exception:
|
|
|
|
@ -108,7 +95,7 @@ class HackRfWideProbe(gr.top_block):
|
|
|
|
self.connect((self.src, 0), (self.stream_to_vec, 0))
|
|
|
|
self.connect((self.src, 0), (self.stream_to_vec, 0))
|
|
|
|
self.connect((self.stream_to_vec, 0), (self.probe, 0))
|
|
|
|
self.connect((self.stream_to_vec, 0), (self.probe, 0))
|
|
|
|
|
|
|
|
|
|
|
|
def tune(self, freq_hz: float):
|
|
|
|
def tune(self, freq_hz: float) -> None:
|
|
|
|
self.src.set_center_freq(freq_hz, 0)
|
|
|
|
self.src.set_center_freq(freq_hz, 0)
|
|
|
|
|
|
|
|
|
|
|
|
def read_metrics(self) -> Tuple[float, float, float, int]:
|
|
|
|
def read_metrics(self) -> Tuple[float, float, float, int]:
|
|
|
|
@ -152,154 +139,374 @@ class HackRfWideProbe(gr.top_block):
|
|
|
|
return rms, power_lin, dbfs, samples
|
|
|
|
return rms, power_lin, dbfs, samples
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
class IioWideProbe:
|
|
|
|
def parse_hackrf_info() -> Dict[str, int]:
|
|
|
|
def __init__(self, args: argparse.Namespace):
|
|
|
|
try:
|
|
|
|
self.args = args
|
|
|
|
proc = subprocess.run(["hackrf_info"], capture_output=True, text=True, timeout=15)
|
|
|
|
self._static_signature: Optional[Tuple[Optional[float], Optional[float], str, Optional[float], str]] = None
|
|
|
|
except FileNotFoundError:
|
|
|
|
self._last_freq_hz: Optional[int] = None
|
|
|
|
raise RuntimeError("hackrf_info not found")
|
|
|
|
self._input_channels = [args.iio_i_channel]
|
|
|
|
except subprocess.TimeoutExpired:
|
|
|
|
if args.iio_q_channel not in self._input_channels:
|
|
|
|
raise RuntimeError("hackrf_info timeout")
|
|
|
|
self._input_channels.append(args.iio_q_channel)
|
|
|
|
text = (proc.stdout or "") + "\n" + (proc.stderr or "")
|
|
|
|
|
|
|
|
out: Dict[str, int] = {}
|
|
|
|
def _run(self, cmd: Sequence[str], binary: bool = False) -> str | bytes:
|
|
|
|
cur_idx: Optional[int] = None
|
|
|
|
proc = subprocess.run(list(cmd), capture_output=True, text=not binary)
|
|
|
|
for line in text.splitlines():
|
|
|
|
if proc.returncode != 0:
|
|
|
|
m = re.search(r"^Index:\s*(\d+)", line)
|
|
|
|
stderr = proc.stderr if binary else (proc.stderr or "")
|
|
|
|
if m:
|
|
|
|
stdout = "" if binary else (proc.stdout or "")
|
|
|
|
cur_idx = int(m.group(1))
|
|
|
|
details = stderr.strip() or stdout.strip() or f"exit code {proc.returncode}"
|
|
|
|
continue
|
|
|
|
raise RuntimeError(details)
|
|
|
|
m = re.search(r"^Serial number:\s*([0-9a-fA-F]+)", line)
|
|
|
|
return proc.stdout
|
|
|
|
if m and cur_idx is not None:
|
|
|
|
|
|
|
|
out[m.group(1).lower()] = cur_idx
|
|
|
|
def _run_binary(self, cmd: Sequence[str]) -> bytes:
|
|
|
|
if not out:
|
|
|
|
proc = subprocess.run(list(cmd), capture_output=True)
|
|
|
|
raise RuntimeError("no devices parsed from hackrf_info")
|
|
|
|
if proc.returncode != 0:
|
|
|
|
return out
|
|
|
|
details = proc.stderr.decode("utf-8", errors="ignore").strip() or proc.stdout.decode("utf-8", errors="ignore").strip()
|
|
|
|
|
|
|
|
raise RuntimeError(details or f"exit code {proc.returncode}")
|
|
|
|
|
|
|
|
return proc.stdout
|
|
|
|
def fmt(value: Optional[float], spec: str) -> str:
|
|
|
|
|
|
|
|
return "-" if value is None else format(value, spec)
|
|
|
|
def _read_input_attr(self, channel: str, attr: str) -> str:
|
|
|
|
|
|
|
|
out = self._run(
|
|
|
|
|
|
|
|
[
|
|
|
|
def build_windows(base_mhz: float, roof_mhz: float, step_mhz: float) -> List[ScanWindow]:
|
|
|
|
"iio_attr",
|
|
|
|
if step_mhz <= 0:
|
|
|
|
"-u",
|
|
|
|
raise ValueError("step must be > 0")
|
|
|
|
self.args.uri,
|
|
|
|
if base_mhz == roof_mhz:
|
|
|
|
"-i",
|
|
|
|
raise ValueError("base and roof must be different")
|
|
|
|
"-c",
|
|
|
|
|
|
|
|
self.args.iio_phy_device,
|
|
|
|
direction = -1.0 if roof_mhz < base_mhz else 1.0
|
|
|
|
channel,
|
|
|
|
edge = base_mhz
|
|
|
|
attr,
|
|
|
|
seq = 1
|
|
|
|
]
|
|
|
|
windows: List[ScanWindow] = []
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
while True:
|
|
|
|
|
|
|
|
next_edge = edge + direction * step_mhz
|
|
|
|
|
|
|
|
if direction < 0 and next_edge < roof_mhz:
|
|
|
|
|
|
|
|
next_edge = roof_mhz
|
|
|
|
|
|
|
|
if direction > 0 and next_edge > roof_mhz:
|
|
|
|
|
|
|
|
next_edge = roof_mhz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
low_mhz = min(edge, next_edge)
|
|
|
|
|
|
|
|
high_mhz = max(edge, next_edge)
|
|
|
|
|
|
|
|
center_mhz = (low_mhz + high_mhz) / 2.0
|
|
|
|
|
|
|
|
windows.append(
|
|
|
|
|
|
|
|
ScanWindow(
|
|
|
|
|
|
|
|
seq=seq,
|
|
|
|
|
|
|
|
start_mhz=edge,
|
|
|
|
|
|
|
|
end_mhz=next_edge,
|
|
|
|
|
|
|
|
low_mhz=low_mhz,
|
|
|
|
|
|
|
|
high_mhz=high_mhz,
|
|
|
|
|
|
|
|
center_mhz=center_mhz,
|
|
|
|
|
|
|
|
)
|
|
|
|
)
|
|
|
|
)
|
|
|
|
return str(out).strip()
|
|
|
|
|
|
|
|
|
|
|
|
if next_edge == roof_mhz:
|
|
|
|
def _set_input_attr(self, channel: str, attr: str, value: str):
|
|
|
|
break
|
|
|
|
self._run(
|
|
|
|
edge = next_edge
|
|
|
|
[
|
|
|
|
seq += 1
|
|
|
|
"iio_attr",
|
|
|
|
|
|
|
|
"-u",
|
|
|
|
return windows
|
|
|
|
self.args.uri,
|
|
|
|
|
|
|
|
"-q",
|
|
|
|
|
|
|
|
"-i",
|
|
|
|
def render(
|
|
|
|
"-c",
|
|
|
|
windows: List[ScanWindow],
|
|
|
|
self.args.iio_phy_device,
|
|
|
|
serial: str,
|
|
|
|
channel,
|
|
|
|
index: int,
|
|
|
|
attr,
|
|
|
|
sample_rate: float,
|
|
|
|
value,
|
|
|
|
base_mhz: float,
|
|
|
|
]
|
|
|
|
roof_mhz: float,
|
|
|
|
|
|
|
|
step_mhz: float,
|
|
|
|
|
|
|
|
started_at: float,
|
|
|
|
|
|
|
|
pass_no: int,
|
|
|
|
|
|
|
|
current_seq: int,
|
|
|
|
|
|
|
|
) -> None:
|
|
|
|
|
|
|
|
now = time.time()
|
|
|
|
|
|
|
|
capture_bw_mhz = sample_rate / 1e6
|
|
|
|
|
|
|
|
current_row = next((row for row in windows if row.seq == current_seq), None)
|
|
|
|
|
|
|
|
best_row = max(
|
|
|
|
|
|
|
|
(row for row in windows if row.status == "OK" and row.dbfs is not None),
|
|
|
|
|
|
|
|
key=lambda row: row.dbfs if row.dbfs is not None else float("-inf"),
|
|
|
|
|
|
|
|
default=None,
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
print("\x1b[2J\x1b[H", end="")
|
|
|
|
|
|
|
|
print("HackRF Wide Energy Monitor (relative power: RMS / linear / dBFS)")
|
|
|
|
|
|
|
|
print(
|
|
|
|
|
|
|
|
f"serial: {serial} | idx: {index} | sample-rate: {capture_bw_mhz:.3f} MHz | "
|
|
|
|
|
|
|
|
f"scan: {base_mhz:.3f}->{roof_mhz:.3f} MHz step {step_mhz:.3f} MHz | "
|
|
|
|
|
|
|
|
f"pass: {pass_no} | uptime: {int(now-started_at)}s | {time.strftime('%Y-%m-%d %H:%M:%S')}"
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
print()
|
|
|
|
|
|
|
|
header = (
|
|
|
|
|
|
|
|
f"{'cur':>3} {'seq':>3} {'window MHz':>23} {'center':>9} {'status':>8} "
|
|
|
|
|
|
|
|
f"{'dBFS':>9} {'rms':>10} {'power':>12} {'N':>5} {'age':>5} error"
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
print(header)
|
|
|
|
|
|
|
|
print("-" * len(header))
|
|
|
|
|
|
|
|
for row in windows:
|
|
|
|
|
|
|
|
age = "-" if row.updated_at <= 0 else f"{(now-row.updated_at):.1f}"
|
|
|
|
|
|
|
|
err = row.error
|
|
|
|
|
|
|
|
if len(err) > 50:
|
|
|
|
|
|
|
|
err = err[:47] + "..."
|
|
|
|
|
|
|
|
marker = ">>>" if row.seq == current_seq else ""
|
|
|
|
|
|
|
|
print(
|
|
|
|
|
|
|
|
f"{marker:>3} {row.seq:>3} "
|
|
|
|
|
|
|
|
f"{f'{row.high_mhz:.3f}-{row.low_mhz:.3f}':>23} {row.center_mhz:>9.3f} {row.status:>8} "
|
|
|
|
|
|
|
|
f"{fmt(row.dbfs, '.2f'):>9} {fmt(row.rms, '.6f'):>10} {fmt(row.power_lin, '.8f'):>12} "
|
|
|
|
|
|
|
|
f"{row.samples:>5} {age:>5} {err}"
|
|
|
|
)
|
|
|
|
)
|
|
|
|
|
|
|
|
print()
|
|
|
|
|
|
|
|
if best_row is not None:
|
|
|
|
|
|
|
|
best_age = "-" if best_row.updated_at <= 0 else f"{(now-best_row.updated_at):.1f}"
|
|
|
|
|
|
|
|
print(
|
|
|
|
|
|
|
|
f"{'':>3} {'MAX':>3} "
|
|
|
|
|
|
|
|
f"{f'{best_row.high_mhz:.3f}-{best_row.low_mhz:.3f}':>23} {best_row.center_mhz:>9.3f} {best_row.status:>8} "
|
|
|
|
|
|
|
|
f"{fmt(best_row.dbfs, '.2f'):>9} {fmt(best_row.rms, '.6f'):>10} {fmt(best_row.power_lin, '.8f'):>12} "
|
|
|
|
|
|
|
|
f"{best_row.samples:>5} {best_age:>5} pass={best_row.pass_no}"
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
elif current_row is not None:
|
|
|
|
|
|
|
|
current_age = "-" if current_row.updated_at <= 0 else f"{(now-current_row.updated_at):.1f}"
|
|
|
|
|
|
|
|
print(
|
|
|
|
|
|
|
|
f"{'':>3} {'MAX':>3} "
|
|
|
|
|
|
|
|
f"{f'{current_row.high_mhz:.3f}-{current_row.low_mhz:.3f}':>23} {current_row.center_mhz:>9.3f} {'INIT':>8} "
|
|
|
|
|
|
|
|
f"{fmt(None, '.2f'):>9} {fmt(None, '.6f'):>10} {fmt(None, '.8f'):>12} "
|
|
|
|
|
|
|
|
f"{0:>5} {current_age:>5} no successful windows yet"
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
print("Ctrl+C to stop. Window width equals step; sample-rate must be >= step to cover each window.")
|
|
|
|
|
|
|
|
sys.stdout.flush()
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def build_parser() -> argparse.ArgumentParser:
|
|
|
|
|
|
|
|
parser = argparse.ArgumentParser(description="Retune one HackRF across a wide frequency range and measure energy")
|
|
|
|
|
|
|
|
parser.add_argument("--serial", required=True, help="HackRF serial number from hackrf_info")
|
|
|
|
|
|
|
|
parser.add_argument("--sample-rate", type=float, required=True, help="Sample rate in Hz")
|
|
|
|
|
|
|
|
parser.add_argument("--base", type=float, required=True, help="Scan start edge in MHz")
|
|
|
|
|
|
|
|
parser.add_argument("--roof", type=float, required=True, help="Scan end edge in MHz")
|
|
|
|
|
|
|
|
parser.add_argument("--step", type=float, required=True, help="Window width / retune step in MHz")
|
|
|
|
|
|
|
|
parser.add_argument("--vec-len", type=int, default=4096, help="Probe vector length")
|
|
|
|
|
|
|
|
parser.add_argument("--settle", type=float, default=0.12, help="Wait time after retune before reading (s)")
|
|
|
|
|
|
|
|
parser.add_argument("--avg-reads", type=int, default=3, help="How many probe reads to average per window")
|
|
|
|
|
|
|
|
parser.add_argument("--pause-between-reads", type=float, default=0.02, help="Pause between averaged reads (s)")
|
|
|
|
|
|
|
|
parser.add_argument("--passes", type=int, default=0, help="Number of sweep passes, 0 means infinite")
|
|
|
|
|
|
|
|
parser.add_argument("--gain", type=float, default=16.0, help="General gain")
|
|
|
|
|
|
|
|
parser.add_argument("--if-gain", type=float, default=16.0, help="IF gain")
|
|
|
|
|
|
|
|
parser.add_argument("--bb-gain", type=float, default=16.0, help="BB gain")
|
|
|
|
|
|
|
|
return parser
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _set_output_attr(self, channel: str, attr: str, value: str):
|
|
|
|
def main() -> int:
|
|
|
|
self._run(
|
|
|
|
args = build_parser().parse_args()
|
|
|
|
[
|
|
|
|
serial = args.serial.lower()
|
|
|
|
"iio_attr",
|
|
|
|
|
|
|
|
"-u",
|
|
|
|
try:
|
|
|
|
self.args.uri,
|
|
|
|
windows = build_windows(args.base, args.roof, args.step)
|
|
|
|
"-q",
|
|
|
|
except ValueError as exc:
|
|
|
|
"-o",
|
|
|
|
print(f"invalid scan range: {exc}", file=sys.stderr)
|
|
|
|
"-c",
|
|
|
|
return 2
|
|
|
|
self.args.iio_phy_device,
|
|
|
|
|
|
|
|
channel,
|
|
|
|
step_hz = args.step * 1e6
|
|
|
|
attr,
|
|
|
|
if args.sample_rate < step_hz:
|
|
|
|
value,
|
|
|
|
print(
|
|
|
|
]
|
|
|
|
f"sample-rate {args.sample_rate:.0f} Hz is smaller than step window {step_hz:.0f} Hz; "
|
|
|
|
|
|
|
|
"this would leave gaps in the scan",
|
|
|
|
|
|
|
|
file=sys.stderr,
|
|
|
|
)
|
|
|
|
)
|
|
|
|
|
|
|
|
return 2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
serial_to_index = parse_hackrf_info()
|
|
|
|
|
|
|
|
except Exception as exc:
|
|
|
|
|
|
|
|
print(f"hackrf discovery failed: {exc}", file=sys.stderr)
|
|
|
|
|
|
|
|
return 3
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
index = serial_to_index.get(serial)
|
|
|
|
|
|
|
|
if index is None:
|
|
|
|
|
|
|
|
print(f"serial {serial} not found in hackrf_info", file=sys.stderr)
|
|
|
|
|
|
|
|
print("available serials:", file=sys.stderr)
|
|
|
|
|
|
|
|
for item_serial, item_index in sorted(serial_to_index.items(), key=lambda item: item[1]):
|
|
|
|
|
|
|
|
print(f" idx={item_index} serial={item_serial}", file=sys.stderr)
|
|
|
|
|
|
|
|
return 4
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
stop_requested = False
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def on_signal(signum, frame):
|
|
|
|
|
|
|
|
nonlocal stop_requested
|
|
|
|
|
|
|
|
stop_requested = True
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
signal.signal(signal.SIGINT, on_signal)
|
|
|
|
|
|
|
|
signal.signal(signal.SIGTERM, on_signal)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
probe: Optional[WideProbeTop] = None
|
|
|
|
|
|
|
|
started_at = time.time()
|
|
|
|
|
|
|
|
pass_no = 0
|
|
|
|
|
|
|
|
current_seq = windows[0].seq
|
|
|
|
|
|
|
|
|
|
|
|
def get_effective_sample_rate(self) -> float:
|
|
|
|
try:
|
|
|
|
if self.args.sample_rate is not None:
|
|
|
|
probe = WideProbeTop(
|
|
|
|
return float(self.args.sample_rate)
|
|
|
|
index=index,
|
|
|
|
raw = self._read_input_attr(self.args.iio_i_channel, "sampling_frequency")
|
|
|
|
center_freq_hz=windows[0].center_mhz * 1e6,
|
|
|
|
return float(raw)
|
|
|
|
sample_rate=args.sample_rate,
|
|
|
|
|
|
|
|
vec_len=args.vec_len,
|
|
|
|
def ensure_configured(self):
|
|
|
|
gain=args.gain,
|
|
|
|
sample_rate = None
|
|
|
|
if_gain=args.if_gain,
|
|
|
|
if self.args.sample_rate is not None:
|
|
|
|
bb_gain=args.bb_gain,
|
|
|
|
sample_rate = int(round(max(float(self.args.sample_rate), IIO_MIN_SAMPLE_RATE)))
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
bandwidth = None
|
|
|
|
|
|
|
|
if self.args.bandwidth is not None:
|
|
|
|
|
|
|
|
bandwidth = int(round(max(float(self.args.bandwidth), 200000.0)))
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
signature = (
|
|
|
|
|
|
|
|
None if sample_rate is None else float(sample_rate),
|
|
|
|
|
|
|
|
None if bandwidth is None else float(bandwidth),
|
|
|
|
|
|
|
|
self.args.iio_gain_mode,
|
|
|
|
|
|
|
|
self.args.iio_hardwaregain,
|
|
|
|
|
|
|
|
self.args.iio_port_select,
|
|
|
|
|
|
|
|
)
|
|
|
|
)
|
|
|
|
if signature == self._static_signature:
|
|
|
|
probe.start()
|
|
|
|
return
|
|
|
|
time.sleep(max(args.settle, 0.12))
|
|
|
|
|
|
|
|
|
|
|
|
for channel in self._input_channels:
|
|
|
|
while not stop_requested:
|
|
|
|
if sample_rate is not None:
|
|
|
|
pass_no += 1
|
|
|
|
self._set_input_attr(channel, "sampling_frequency", str(sample_rate))
|
|
|
|
for row in windows:
|
|
|
|
if bandwidth is not None:
|
|
|
|
if stop_requested:
|
|
|
|
self._set_input_attr(channel, "rf_bandwidth", str(bandwidth))
|
|
|
|
break
|
|
|
|
if self.args.iio_port_select:
|
|
|
|
current_seq = row.seq
|
|
|
|
self._set_input_attr(channel, "rf_port_select", self.args.iio_port_select)
|
|
|
|
try:
|
|
|
|
if self.args.iio_gain_mode:
|
|
|
|
probe.tune(row.center_mhz * 1e6)
|
|
|
|
self._set_input_attr(channel, "gain_control_mode", self.args.iio_gain_mode)
|
|
|
|
rms, power_lin, dbfs, samples = probe.read_window(
|
|
|
|
if self.args.iio_gain_mode == "manual" and self.args.iio_hardwaregain is not None:
|
|
|
|
settle=args.settle,
|
|
|
|
self._set_input_attr(channel, "hardwaregain", f"{self.args.iio_hardwaregain:.6f}")
|
|
|
|
avg_reads=args.avg_reads,
|
|
|
|
|
|
|
|
pause_between_reads=args.pause_between_reads,
|
|
|
|
self._static_signature = signature
|
|
|
|
)
|
|
|
|
|
|
|
|
row.status = "OK"
|
|
|
|
def tune(self, freq_hz: float):
|
|
|
|
row.rms = rms
|
|
|
|
target = int(round(freq_hz))
|
|
|
|
row.power_lin = power_lin
|
|
|
|
if self._last_freq_hz == target:
|
|
|
|
row.dbfs = dbfs
|
|
|
|
return
|
|
|
|
row.samples = samples
|
|
|
|
self._set_output_attr(self.args.iio_lo_channel, "frequency", str(target))
|
|
|
|
row.error = ""
|
|
|
|
self._last_freq_hz = target
|
|
|
|
row.updated_at = time.time()
|
|
|
|
|
|
|
|
row.pass_no = pass_no
|
|
|
|
def read_metrics(self, samples_per_read: int) -> Tuple[float, float, float, int]:
|
|
|
|
except Exception as exc:
|
|
|
|
cmd = [
|
|
|
|
row.status = "ERR"
|
|
|
|
"iio_readdev",
|
|
|
|
row.error = str(exc)
|
|
|
|
"-u",
|
|
|
|
row.updated_at = time.time()
|
|
|
|
self.args.uri,
|
|
|
|
render(
|
|
|
|
"-T",
|
|
|
|
windows=windows,
|
|
|
|
str(int(self.args.timeout_ms)),
|
|
|
|
serial=serial,
|
|
|
|
"-b",
|
|
|
|
index=index,
|
|
|
|
str(max(4, int(samples_per_read))),
|
|
|
|
sample_rate=args.sample_rate,
|
|
|
|
"-s",
|
|
|
|
base_mhz=args.base,
|
|
|
|
str(max(4, int(samples_per_read))),
|
|
|
|
roof_mhz=args.roof,
|
|
|
|
self.args.iio_device,
|
|
|
|
step_mhz=args.step,
|
|
|
|
self.args.iio_i_channel,
|
|
|
|
started_at=started_at,
|
|
|
|
self.args.iio_q_channel,
|
|
|
|
pass_no=pass_no,
|
|
|
|
]
|
|
|
|
current_seq=current_seq,
|
|
|
|
raw = self._run_binary(cmd)
|
|
|
|
)
|
|
|
|
values = np.frombuffer(raw, dtype="<i2")
|
|
|
|
if args.passes > 0 and pass_no >= args.passes:
|
|
|
|
if values.size == 0:
|
|
|
|
break
|
|
|
|
|
|
|
|
except Exception as exc:
|
|
|
|
|
|
|
|
print(f"scanner failed: {exc}", file=sys.stderr)
|
|
|
|
|
|
|
|
return 5
|
|
|
|
|
|
|
|
finally:
|
|
|
|
|
|
|
|
if probe is not None:
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
probe.stop()
|
|
|
|
|
|
|
|
probe.wait()
|
|
|
|
|
|
|
|
except Exception:
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
return 0
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
if __name__ == "__main__":
|
|
|
|
|
|
|
|
raise SystemExit(main())
|
|
|
|
|
|
|
|
#!/usr/bin/env python3
|
|
|
|
|
|
|
|
import argparse
|
|
|
|
|
|
|
|
import math
|
|
|
|
|
|
|
|
import re
|
|
|
|
|
|
|
|
import signal
|
|
|
|
|
|
|
|
import subprocess
|
|
|
|
|
|
|
|
import sys
|
|
|
|
|
|
|
|
import time
|
|
|
|
|
|
|
|
from dataclasses import dataclass
|
|
|
|
|
|
|
|
from typing import Dict, List, Optional, Tuple
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
import numpy as np
|
|
|
|
|
|
|
|
except Exception as exc:
|
|
|
|
|
|
|
|
print(f"numpy import failed: {exc}", file=sys.stderr)
|
|
|
|
|
|
|
|
sys.exit(1)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
from gnuradio import blocks, gr
|
|
|
|
|
|
|
|
import osmosdr
|
|
|
|
|
|
|
|
except Exception as exc:
|
|
|
|
|
|
|
|
print(f"gnuradio/osmosdr import failed: {exc}", file=sys.stderr)
|
|
|
|
|
|
|
|
print("Run with the SDR venv, e.g. .venv-sdr/bin/python read_energy_wide.py", file=sys.stderr)
|
|
|
|
|
|
|
|
sys.exit(1)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EPS = 1e-20
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@dataclass
|
|
|
|
|
|
|
|
class ScanWindow:
|
|
|
|
|
|
|
|
seq: int
|
|
|
|
|
|
|
|
start_mhz: float
|
|
|
|
|
|
|
|
end_mhz: float
|
|
|
|
|
|
|
|
low_mhz: float
|
|
|
|
|
|
|
|
high_mhz: float
|
|
|
|
|
|
|
|
center_mhz: float
|
|
|
|
|
|
|
|
status: str = "INIT"
|
|
|
|
|
|
|
|
rms: Optional[float] = None
|
|
|
|
|
|
|
|
power_lin: Optional[float] = None
|
|
|
|
|
|
|
|
dbfs: Optional[float] = None
|
|
|
|
|
|
|
|
samples: int = 0
|
|
|
|
|
|
|
|
updated_at: float = 0.0
|
|
|
|
|
|
|
|
error: str = ""
|
|
|
|
|
|
|
|
pass_no: int = 0
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
class WideProbeTop(gr.top_block):
|
|
|
|
|
|
|
|
def __init__(
|
|
|
|
|
|
|
|
self,
|
|
|
|
|
|
|
|
index: int,
|
|
|
|
|
|
|
|
center_freq_hz: float,
|
|
|
|
|
|
|
|
sample_rate: float,
|
|
|
|
|
|
|
|
vec_len: int,
|
|
|
|
|
|
|
|
gain: float,
|
|
|
|
|
|
|
|
if_gain: float,
|
|
|
|
|
|
|
|
bb_gain: float,
|
|
|
|
|
|
|
|
):
|
|
|
|
|
|
|
|
super().__init__("hackrf_energy_wide_probe")
|
|
|
|
|
|
|
|
self.probe = blocks.probe_signal_vc(vec_len)
|
|
|
|
|
|
|
|
self.stream_to_vec = blocks.stream_to_vector(gr.sizeof_gr_complex * 1, vec_len)
|
|
|
|
|
|
|
|
self.src = osmosdr.source(args=f"numchan=1 hackrf={index}")
|
|
|
|
|
|
|
|
self.src.set_time_unknown_pps(osmosdr.time_spec_t())
|
|
|
|
|
|
|
|
self.src.set_sample_rate(sample_rate)
|
|
|
|
|
|
|
|
self.src.set_center_freq(center_freq_hz, 0)
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
self.src.set_freq_corr(0, 0)
|
|
|
|
|
|
|
|
except Exception:
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
self.src.set_gain_mode(False, 0)
|
|
|
|
|
|
|
|
except Exception:
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
for fn, val in (("set_gain", gain), ("set_if_gain", if_gain), ("set_bb_gain", bb_gain)):
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
getattr(self.src, fn)(val, 0)
|
|
|
|
|
|
|
|
except Exception:
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
self.src.set_bandwidth(0, 0)
|
|
|
|
|
|
|
|
except Exception:
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
self.src.set_antenna("", 0)
|
|
|
|
|
|
|
|
except Exception:
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
self.connect((self.src, 0), (self.stream_to_vec, 0))
|
|
|
|
|
|
|
|
self.connect((self.stream_to_vec, 0), (self.probe, 0))
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def tune(self, freq_hz: float) -> None:
|
|
|
|
|
|
|
|
self.src.set_center_freq(freq_hz, 0)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def read_metrics(self) -> Tuple[float, float, float, int]:
|
|
|
|
|
|
|
|
arr = np.asarray(self.probe.level(), dtype=np.complex64)
|
|
|
|
|
|
|
|
if arr.size == 0:
|
|
|
|
raise RuntimeError("no samples")
|
|
|
|
raise RuntimeError("no samples")
|
|
|
|
if values.size % 2 != 0:
|
|
|
|
power_lin = float(np.mean(arr.real * arr.real + arr.imag * arr.imag))
|
|
|
|
raise RuntimeError(f"unexpected IQ sample payload: {values.size} int16 values")
|
|
|
|
|
|
|
|
iq = values.reshape(-1, 2).astype(np.float32, copy=False)
|
|
|
|
|
|
|
|
i = iq[:, 0]
|
|
|
|
|
|
|
|
q = iq[:, 1]
|
|
|
|
|
|
|
|
power_lin = float(np.mean(i * i + q * q))
|
|
|
|
|
|
|
|
rms = math.sqrt(max(power_lin, 0.0))
|
|
|
|
rms = math.sqrt(max(power_lin, 0.0))
|
|
|
|
dbfs = 10.0 * math.log10(max(power_lin, EPS))
|
|
|
|
dbfs = 10.0 * math.log10(max(power_lin, EPS))
|
|
|
|
return rms, power_lin, dbfs, int(iq.shape[0])
|
|
|
|
return rms, power_lin, dbfs, int(arr.size)
|
|
|
|
|
|
|
|
|
|
|
|
def read_window(self, settle: float, avg_reads: int, pause_between_reads: float, samples_per_read: int) -> Tuple[float, float, float, int]:
|
|
|
|
def read_window(self, settle: float, avg_reads: int, pause_between_reads: float) -> Tuple[float, float, float, int]:
|
|
|
|
if settle > 0:
|
|
|
|
if settle > 0:
|
|
|
|
time.sleep(settle)
|
|
|
|
time.sleep(settle)
|
|
|
|
|
|
|
|
|
|
|
|
@ -309,10 +516,10 @@ class IioWideProbe:
|
|
|
|
last_error: Optional[Exception] = None
|
|
|
|
last_error: Optional[Exception] = None
|
|
|
|
|
|
|
|
|
|
|
|
for idx in range(read_count):
|
|
|
|
for idx in range(read_count):
|
|
|
|
deadline = time.time() + max(1.0, self.args.timeout_ms / 1000.0)
|
|
|
|
deadline = time.time() + 1.0
|
|
|
|
while True:
|
|
|
|
while True:
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
_, power_lin, _, samples = self.read_metrics(samples_per_read)
|
|
|
|
_, power_lin, _, samples = self.read_metrics()
|
|
|
|
powers.append(power_lin)
|
|
|
|
powers.append(power_lin)
|
|
|
|
sample_sizes.append(samples)
|
|
|
|
sample_sizes.append(samples)
|
|
|
|
break
|
|
|
|
break
|
|
|
|
@ -334,21 +541,21 @@ class IioWideProbe:
|
|
|
|
def parse_hackrf_info() -> Dict[str, int]:
|
|
|
|
def parse_hackrf_info() -> Dict[str, int]:
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
proc = subprocess.run(["hackrf_info"], capture_output=True, text=True, timeout=15)
|
|
|
|
proc = subprocess.run(["hackrf_info"], capture_output=True, text=True, timeout=15)
|
|
|
|
except FileNotFoundError as exc:
|
|
|
|
except FileNotFoundError:
|
|
|
|
raise RuntimeError("hackrf_info not found") from exc
|
|
|
|
raise RuntimeError("hackrf_info not found")
|
|
|
|
except subprocess.TimeoutExpired as exc:
|
|
|
|
except subprocess.TimeoutExpired:
|
|
|
|
raise RuntimeError("hackrf_info timeout") from exc
|
|
|
|
raise RuntimeError("hackrf_info timeout")
|
|
|
|
text = (proc.stdout or "") + "\n" + (proc.stderr or "")
|
|
|
|
text = (proc.stdout or "") + "\n" + (proc.stderr or "")
|
|
|
|
out: Dict[str, int] = {}
|
|
|
|
out: Dict[str, int] = {}
|
|
|
|
cur_idx: Optional[int] = None
|
|
|
|
cur_idx: Optional[int] = None
|
|
|
|
for line in text.splitlines():
|
|
|
|
for line in text.splitlines():
|
|
|
|
match = re.search(r"^Index:\s*(\d+)", line)
|
|
|
|
m = re.search(r"^Index:\s*(\d+)", line)
|
|
|
|
if match:
|
|
|
|
if m:
|
|
|
|
cur_idx = int(match.group(1))
|
|
|
|
cur_idx = int(m.group(1))
|
|
|
|
continue
|
|
|
|
continue
|
|
|
|
match = re.search(r"^Serial number:\s*([0-9a-fA-F]+)", line)
|
|
|
|
m = re.search(r"^Serial number:\s*([0-9a-fA-F]+)", line)
|
|
|
|
if match and cur_idx is not None:
|
|
|
|
if m and cur_idx is not None:
|
|
|
|
out[match.group(1).lower()] = cur_idx
|
|
|
|
out[m.group(1).lower()] = cur_idx
|
|
|
|
if not out:
|
|
|
|
if not out:
|
|
|
|
raise RuntimeError("no devices parsed from hackrf_info")
|
|
|
|
raise RuntimeError("no devices parsed from hackrf_info")
|
|
|
|
return out
|
|
|
|
return out
|
|
|
|
@ -400,9 +607,8 @@ def build_windows(base_mhz: float, roof_mhz: float, step_mhz: float) -> List[Sca
|
|
|
|
|
|
|
|
|
|
|
|
def render(
|
|
|
|
def render(
|
|
|
|
windows: List[ScanWindow],
|
|
|
|
windows: List[ScanWindow],
|
|
|
|
backend: str,
|
|
|
|
serial: str,
|
|
|
|
device_ref: str,
|
|
|
|
index: int,
|
|
|
|
index_ref: str,
|
|
|
|
|
|
|
|
sample_rate: float,
|
|
|
|
sample_rate: float,
|
|
|
|
base_mhz: float,
|
|
|
|
base_mhz: float,
|
|
|
|
roof_mhz: float,
|
|
|
|
roof_mhz: float,
|
|
|
|
@ -420,93 +626,73 @@ def render(
|
|
|
|
default=None,
|
|
|
|
default=None,
|
|
|
|
)
|
|
|
|
)
|
|
|
|
print("\x1b[2J\x1b[H", end="")
|
|
|
|
print("\x1b[2J\x1b[H", end="")
|
|
|
|
print("Wide Energy Monitor (relative power: RMS / linear / dBFS)")
|
|
|
|
print("HackRF Wide Energy Monitor (relative power: RMS / linear / dBFS)")
|
|
|
|
print(
|
|
|
|
print(
|
|
|
|
f"backend: {backend} | device: {device_ref} | idx: {index_ref} | sample-rate: {capture_bw_mhz:.3f} MHz | "
|
|
|
|
f"serial: {serial} | idx: {index} | sample-rate: {capture_bw_mhz:.3f} MHz | "
|
|
|
|
f"scan: {base_mhz:.3f}->{roof_mhz:.3f} MHz step {step_mhz:.3f} MHz | "
|
|
|
|
f"scan: {base_mhz:.3f}->{roof_mhz:.3f} MHz step {step_mhz:.3f} MHz | "
|
|
|
|
f"pass: {pass_no} | uptime: {int(now-started_at)}s | {time.strftime('%Y-%m-%d %H:%M:%S')}"
|
|
|
|
f"pass: {pass_no} | uptime: {int(now-started_at)}s | {time.strftime('%Y-%m-%d %H:%M:%S')}"
|
|
|
|
)
|
|
|
|
)
|
|
|
|
print()
|
|
|
|
print()
|
|
|
|
header = (
|
|
|
|
header = (
|
|
|
|
f"{'cur':>3} {'seq':>3} {'window MHz':>23} {'center':>9} {'status':>8} "
|
|
|
|
f"{'cur':>3} {'seq':>3} {'window MHz':>23} {'center':>9} {'status':>8} "
|
|
|
|
f"{'dBFS':>9} {'rms':>10} {'power':>12} {'N':>6} {'age':>5} error"
|
|
|
|
f"{'dBFS':>9} {'rms':>10} {'power':>12} {'N':>5} {'age':>5} error"
|
|
|
|
)
|
|
|
|
)
|
|
|
|
print(header)
|
|
|
|
print(header)
|
|
|
|
print("-" * len(header))
|
|
|
|
print("-" * len(header))
|
|
|
|
for row in windows:
|
|
|
|
for row in windows:
|
|
|
|
age = "-" if row.updated_at <= 0 else f"{(now - row.updated_at):.1f}"
|
|
|
|
age = "-" if row.updated_at <= 0 else f"{(now-row.updated_at):.1f}"
|
|
|
|
err = row.error
|
|
|
|
err = row.error
|
|
|
|
if len(err) > 56:
|
|
|
|
if len(err) > 50:
|
|
|
|
err = err[:53] + "..."
|
|
|
|
err = err[:47] + "..."
|
|
|
|
marker = ">>>" if row.seq == current_seq else ""
|
|
|
|
marker = ">>>" if row.seq == current_seq else ""
|
|
|
|
print(
|
|
|
|
print(
|
|
|
|
f"{marker:>3} {row.seq:>3} "
|
|
|
|
f"{marker:>3} {row.seq:>3} "
|
|
|
|
f"{f'{row.high_mhz:.3f}-{row.low_mhz:.3f}':>23} {row.center_mhz:>9.3f} {row.status:>8} "
|
|
|
|
f"{f'{row.high_mhz:.3f}-{row.low_mhz:.3f}':>23} {row.center_mhz:>9.3f} {row.status:>8} "
|
|
|
|
f"{fmt(row.dbfs, '.2f'):>9} {fmt(row.rms, '.6f'):>10} {fmt(row.power_lin, '.8f'):>12} "
|
|
|
|
f"{fmt(row.dbfs, '.2f'):>9} {fmt(row.rms, '.6f'):>10} {fmt(row.power_lin, '.8f'):>12} "
|
|
|
|
f"{row.samples:>6} {age:>5} {err}"
|
|
|
|
f"{row.samples:>5} {age:>5} {err}"
|
|
|
|
)
|
|
|
|
)
|
|
|
|
print()
|
|
|
|
print()
|
|
|
|
if best_row is not None:
|
|
|
|
if best_row is not None:
|
|
|
|
best_age = "-" if best_row.updated_at <= 0 else f"{(now - best_row.updated_at):.1f}"
|
|
|
|
best_age = "-" if best_row.updated_at <= 0 else f"{(now-best_row.updated_at):.1f}"
|
|
|
|
print(
|
|
|
|
print(
|
|
|
|
f"{'':>3} {'MAX':>3} "
|
|
|
|
f"{'':>3} {'MAX':>3} "
|
|
|
|
f"{f'{best_row.high_mhz:.3f}-{best_row.low_mhz:.3f}':>23} {best_row.center_mhz:>9.3f} {best_row.status:>8} "
|
|
|
|
f"{f'{best_row.high_mhz:.3f}-{best_row.low_mhz:.3f}':>23} {best_row.center_mhz:>9.3f} {best_row.status:>8} "
|
|
|
|
f"{fmt(best_row.dbfs, '.2f'):>9} {fmt(best_row.rms, '.6f'):>10} {fmt(best_row.power_lin, '.8f'):>12} "
|
|
|
|
f"{fmt(best_row.dbfs, '.2f'):>9} {fmt(best_row.rms, '.6f'):>10} {fmt(best_row.power_lin, '.8f'):>12} "
|
|
|
|
f"{best_row.samples:>6} {best_age:>5} pass={best_row.pass_no}"
|
|
|
|
f"{best_row.samples:>5} {best_age:>5} pass={best_row.pass_no}"
|
|
|
|
)
|
|
|
|
)
|
|
|
|
elif current_row is not None:
|
|
|
|
elif current_row is not None:
|
|
|
|
current_age = "-" if current_row.updated_at <= 0 else f"{(now - current_row.updated_at):.1f}"
|
|
|
|
current_age = "-" if current_row.updated_at <= 0 else f"{(now-current_row.updated_at):.1f}"
|
|
|
|
print(
|
|
|
|
print(
|
|
|
|
f"{'':>3} {'MAX':>3} "
|
|
|
|
f"{'':>3} {'MAX':>3} "
|
|
|
|
f"{f'{current_row.high_mhz:.3f}-{current_row.low_mhz:.3f}':>23} {current_row.center_mhz:>9.3f} {'INIT':>8} "
|
|
|
|
f"{f'{current_row.high_mhz:.3f}-{current_row.low_mhz:.3f}':>23} {current_row.center_mhz:>9.3f} {'INIT':>8} "
|
|
|
|
f"{fmt(None, '.2f'):>9} {fmt(None, '.6f'):>10} {fmt(None, '.8f'):>12} "
|
|
|
|
f"{fmt(None, '.2f'):>9} {fmt(None, '.6f'):>10} {fmt(None, '.8f'):>12} "
|
|
|
|
f"{0:>6} {current_age:>5} no successful windows yet"
|
|
|
|
f"{0:>5} {current_age:>5} no successful windows yet"
|
|
|
|
)
|
|
|
|
)
|
|
|
|
print("Ctrl+C to stop. Window width equals step; sample-rate must be >= step to cover each window.")
|
|
|
|
print("Ctrl+C to stop. Window width equals step; sample-rate must be >= step to cover each window.")
|
|
|
|
sys.stdout.flush()
|
|
|
|
sys.stdout.flush()
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def build_parser() -> argparse.ArgumentParser:
|
|
|
|
def build_parser() -> argparse.ArgumentParser:
|
|
|
|
parser = argparse.ArgumentParser(description="Retune one SDR across a wide frequency range and measure energy")
|
|
|
|
parser = argparse.ArgumentParser(description="Retune one HackRF across a wide frequency range and measure energy")
|
|
|
|
parser.add_argument("--backend", choices=("hackrf", "iio"), default="hackrf", help="SDR backend")
|
|
|
|
parser.add_argument("--serial", required=True, help="HackRF serial number from hackrf_info")
|
|
|
|
parser.add_argument("--serial", help="HackRF serial number from hackrf_info")
|
|
|
|
parser.add_argument("--sample-rate", type=float, required=True, help="Sample rate in Hz")
|
|
|
|
parser.add_argument(
|
|
|
|
|
|
|
|
"--sample-rate",
|
|
|
|
|
|
|
|
type=float,
|
|
|
|
|
|
|
|
default=None,
|
|
|
|
|
|
|
|
help="Sample rate in Hz (required for hackrf, optional for iio to keep device setting)",
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
parser.add_argument("--bandwidth", type=float, default=None, help="RF bandwidth in Hz (iio only, optional)")
|
|
|
|
|
|
|
|
parser.add_argument("--base", type=float, required=True, help="Scan start edge in MHz")
|
|
|
|
parser.add_argument("--base", type=float, required=True, help="Scan start edge in MHz")
|
|
|
|
parser.add_argument("--roof", type=float, required=True, help="Scan end edge in MHz")
|
|
|
|
parser.add_argument("--roof", type=float, required=True, help="Scan end edge in MHz")
|
|
|
|
parser.add_argument("--step", type=float, required=True, help="Window width / retune step in MHz")
|
|
|
|
parser.add_argument("--step", type=float, required=True, help="Window width / retune step in MHz")
|
|
|
|
parser.add_argument("--vec-len", type=int, default=4096, help="Probe vector length / samples per read")
|
|
|
|
parser.add_argument("--vec-len", type=int, default=4096, help="Probe vector length")
|
|
|
|
parser.add_argument("--settle", type=float, default=0.12, help="Wait time after retune before reading (s)")
|
|
|
|
parser.add_argument("--settle", type=float, default=0.12, help="Wait time after retune before reading (s)")
|
|
|
|
parser.add_argument("--avg-reads", type=int, default=3, help="How many reads to average per window")
|
|
|
|
parser.add_argument("--avg-reads", type=int, default=3, help="How many probe reads to average per window")
|
|
|
|
parser.add_argument("--pause-between-reads", type=float, default=0.02, help="Pause between averaged reads (s)")
|
|
|
|
parser.add_argument("--pause-between-reads", type=float, default=0.02, help="Pause between averaged reads (s)")
|
|
|
|
parser.add_argument("--gain", type=float, default=16.0, help="General gain for HackRF")
|
|
|
|
parser.add_argument("--passes", type=int, default=0, help="Number of sweep passes, 0 means infinite")
|
|
|
|
parser.add_argument("--if-gain", type=float, default=16.0, help="IF gain for HackRF")
|
|
|
|
parser.add_argument("--gain", type=float, default=16.0, help="General gain")
|
|
|
|
parser.add_argument("--bb-gain", type=float, default=16.0, help="BB gain for HackRF")
|
|
|
|
parser.add_argument("--if-gain", type=float, default=16.0, help="IF gain")
|
|
|
|
parser.add_argument("--uri", default="ip:192.168.2.1", help="IIO URI, e.g. ip:192.168.2.1")
|
|
|
|
parser.add_argument("--bb-gain", type=float, default=16.0, help="BB gain")
|
|
|
|
parser.add_argument("--iio-device", default="cf-ad9361-lpc", help="IIO RX buffer device")
|
|
|
|
|
|
|
|
parser.add_argument("--iio-phy-device", default="ad9361-phy", help="IIO PHY device")
|
|
|
|
|
|
|
|
parser.add_argument("--iio-i-channel", default="voltage0", help="IIO I channel")
|
|
|
|
|
|
|
|
parser.add_argument("--iio-q-channel", default="voltage1", help="IIO Q channel")
|
|
|
|
|
|
|
|
parser.add_argument("--iio-lo-channel", default="altvoltage0", help="IIO LO channel for RX frequency")
|
|
|
|
|
|
|
|
parser.add_argument("--iio-port-select", default="A_BALANCED", help="IIO rf_port_select value")
|
|
|
|
|
|
|
|
parser.add_argument(
|
|
|
|
|
|
|
|
"--iio-gain-mode",
|
|
|
|
|
|
|
|
choices=("manual", "fast_attack", "slow_attack", "hybrid"),
|
|
|
|
|
|
|
|
default="slow_attack",
|
|
|
|
|
|
|
|
help="IIO gain control mode",
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
parser.add_argument("--iio-hardwaregain", type=float, default=None, help="IIO hardware gain in dB for manual mode")
|
|
|
|
|
|
|
|
parser.add_argument("--timeout-ms", type=int, default=4000, help="IIO read timeout in milliseconds")
|
|
|
|
|
|
|
|
return parser
|
|
|
|
return parser
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def main() -> int:
|
|
|
|
def main() -> int:
|
|
|
|
args = build_parser().parse_args()
|
|
|
|
args = build_parser().parse_args()
|
|
|
|
|
|
|
|
serial = args.serial.lower()
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
windows = build_windows(args.base, args.roof, args.step)
|
|
|
|
windows = build_windows(args.base, args.roof, args.step)
|
|
|
|
@ -514,72 +700,28 @@ def main() -> int:
|
|
|
|
print(f"invalid scan range: {exc}", file=sys.stderr)
|
|
|
|
print(f"invalid scan range: {exc}", file=sys.stderr)
|
|
|
|
return 2
|
|
|
|
return 2
|
|
|
|
|
|
|
|
|
|
|
|
probe: Optional[HackRfWideProbe | IioWideProbe] = None
|
|
|
|
step_hz = args.step * 1e6
|
|
|
|
device_ref = ""
|
|
|
|
if args.sample_rate < step_hz:
|
|
|
|
index_ref = "-"
|
|
|
|
print(
|
|
|
|
|
|
|
|
f"sample-rate {args.sample_rate:.0f} Hz is smaller than step window {step_hz:.0f} Hz; "
|
|
|
|
if args.backend == "hackrf":
|
|
|
|
"this would leave gaps in the scan",
|
|
|
|
if not args.serial:
|
|
|
|
file=sys.stderr,
|
|
|
|
print("--serial is required for --backend hackrf", file=sys.stderr)
|
|
|
|
|
|
|
|
return 2
|
|
|
|
|
|
|
|
if args.sample_rate is None:
|
|
|
|
|
|
|
|
print("--sample-rate is required for --backend hackrf", file=sys.stderr)
|
|
|
|
|
|
|
|
return 2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
serial = args.serial.lower()
|
|
|
|
|
|
|
|
step_hz = args.step * 1e6
|
|
|
|
|
|
|
|
if args.sample_rate < step_hz:
|
|
|
|
|
|
|
|
print(
|
|
|
|
|
|
|
|
f"sample-rate {args.sample_rate:.0f} Hz is smaller than step window {step_hz:.0f} Hz; this would leave gaps in the scan",
|
|
|
|
|
|
|
|
file=sys.stderr,
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
return 2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
|
|
serial_to_index = parse_hackrf_info()
|
|
|
|
|
|
|
|
except Exception as exc:
|
|
|
|
|
|
|
|
print(f"hackrf discovery failed: {exc}", file=sys.stderr)
|
|
|
|
|
|
|
|
return 3
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
index = serial_to_index.get(serial)
|
|
|
|
|
|
|
|
if index is None:
|
|
|
|
|
|
|
|
print(f"serial {serial} not found in hackrf_info", file=sys.stderr)
|
|
|
|
|
|
|
|
print("available serials:", file=sys.stderr)
|
|
|
|
|
|
|
|
for item_serial, item_index in sorted(serial_to_index.items(), key=lambda item: item[1]):
|
|
|
|
|
|
|
|
print(f" idx={item_index} serial={item_serial}", file=sys.stderr)
|
|
|
|
|
|
|
|
return 4
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
probe = HackRfWideProbe(
|
|
|
|
|
|
|
|
serial=serial,
|
|
|
|
|
|
|
|
center_freq_hz=windows[0].center_mhz * 1e6,
|
|
|
|
|
|
|
|
sample_rate=args.sample_rate,
|
|
|
|
|
|
|
|
vec_len=args.vec_len,
|
|
|
|
|
|
|
|
gain=args.gain,
|
|
|
|
|
|
|
|
if_gain=args.if_gain,
|
|
|
|
|
|
|
|
bb_gain=args.bb_gain,
|
|
|
|
|
|
|
|
)
|
|
|
|
)
|
|
|
|
effective_sample_rate = float(args.sample_rate)
|
|
|
|
return 2
|
|
|
|
device_ref = serial
|
|
|
|
|
|
|
|
index_ref = str(index)
|
|
|
|
try:
|
|
|
|
else:
|
|
|
|
serial_to_index = parse_hackrf_info()
|
|
|
|
probe = IioWideProbe(args)
|
|
|
|
except Exception as exc:
|
|
|
|
try:
|
|
|
|
print(f"hackrf discovery failed: {exc}", file=sys.stderr)
|
|
|
|
probe.ensure_configured()
|
|
|
|
return 3
|
|
|
|
effective_sample_rate = probe.get_effective_sample_rate()
|
|
|
|
|
|
|
|
except Exception as exc:
|
|
|
|
|
|
|
|
print(f"iio setup failed: {exc}", file=sys.stderr)
|
|
|
|
|
|
|
|
return 3
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
step_hz = args.step * 1e6
|
|
|
|
|
|
|
|
if effective_sample_rate < step_hz:
|
|
|
|
|
|
|
|
print(
|
|
|
|
|
|
|
|
f"effective sample-rate {effective_sample_rate:.0f} Hz is smaller than step window {step_hz:.0f} Hz; this would leave gaps in the scan",
|
|
|
|
|
|
|
|
file=sys.stderr,
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
return 2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
device_ref = args.uri
|
|
|
|
index = serial_to_index.get(serial)
|
|
|
|
index_ref = "iio"
|
|
|
|
if index is None:
|
|
|
|
|
|
|
|
print(f"serial {serial} not found in hackrf_info", file=sys.stderr)
|
|
|
|
|
|
|
|
print("available serials:", file=sys.stderr)
|
|
|
|
|
|
|
|
for item_serial, item_index in sorted(serial_to_index.items(), key=lambda item: item[1]):
|
|
|
|
|
|
|
|
print(f" idx={item_index} serial={item_serial}", file=sys.stderr)
|
|
|
|
|
|
|
|
return 4
|
|
|
|
|
|
|
|
|
|
|
|
stop_requested = False
|
|
|
|
stop_requested = False
|
|
|
|
|
|
|
|
|
|
|
|
@ -590,14 +732,23 @@ def main() -> int:
|
|
|
|
signal.signal(signal.SIGINT, on_signal)
|
|
|
|
signal.signal(signal.SIGINT, on_signal)
|
|
|
|
signal.signal(signal.SIGTERM, on_signal)
|
|
|
|
signal.signal(signal.SIGTERM, on_signal)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
probe: Optional[WideProbeTop] = None
|
|
|
|
started_at = time.time()
|
|
|
|
started_at = time.time()
|
|
|
|
pass_no = 0
|
|
|
|
pass_no = 0
|
|
|
|
current_seq = windows[0].seq
|
|
|
|
current_seq = windows[0].seq
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
if isinstance(probe, HackRfWideProbe):
|
|
|
|
probe = WideProbeTop(
|
|
|
|
probe.start()
|
|
|
|
index=index,
|
|
|
|
time.sleep(max(args.settle, 0.12))
|
|
|
|
center_freq_hz=windows[0].center_mhz * 1e6,
|
|
|
|
|
|
|
|
sample_rate=args.sample_rate,
|
|
|
|
|
|
|
|
vec_len=args.vec_len,
|
|
|
|
|
|
|
|
gain=args.gain,
|
|
|
|
|
|
|
|
if_gain=args.if_gain,
|
|
|
|
|
|
|
|
bb_gain=args.bb_gain,
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
probe.start()
|
|
|
|
|
|
|
|
time.sleep(max(args.settle, 0.12))
|
|
|
|
|
|
|
|
|
|
|
|
while not stop_requested:
|
|
|
|
while not stop_requested:
|
|
|
|
pass_no += 1
|
|
|
|
pass_no += 1
|
|
|
|
@ -607,19 +758,11 @@ def main() -> int:
|
|
|
|
current_seq = row.seq
|
|
|
|
current_seq = row.seq
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
probe.tune(row.center_mhz * 1e6)
|
|
|
|
probe.tune(row.center_mhz * 1e6)
|
|
|
|
if isinstance(probe, HackRfWideProbe):
|
|
|
|
rms, power_lin, dbfs, samples = probe.read_window(
|
|
|
|
rms, power_lin, dbfs, samples = probe.read_window(
|
|
|
|
settle=args.settle,
|
|
|
|
settle=args.settle,
|
|
|
|
avg_reads=args.avg_reads,
|
|
|
|
avg_reads=args.avg_reads,
|
|
|
|
pause_between_reads=args.pause_between_reads,
|
|
|
|
pause_between_reads=args.pause_between_reads,
|
|
|
|
)
|
|
|
|
)
|
|
|
|
|
|
|
|
else:
|
|
|
|
|
|
|
|
rms, power_lin, dbfs, samples = probe.read_window(
|
|
|
|
|
|
|
|
settle=args.settle,
|
|
|
|
|
|
|
|
avg_reads=args.avg_reads,
|
|
|
|
|
|
|
|
pause_between_reads=args.pause_between_reads,
|
|
|
|
|
|
|
|
samples_per_read=args.vec_len,
|
|
|
|
|
|
|
|
)
|
|
|
|
|
|
|
|
row.status = "OK"
|
|
|
|
row.status = "OK"
|
|
|
|
row.rms = rms
|
|
|
|
row.rms = rms
|
|
|
|
row.power_lin = power_lin
|
|
|
|
row.power_lin = power_lin
|
|
|
|
@ -632,13 +775,11 @@ def main() -> int:
|
|
|
|
row.status = "ERR"
|
|
|
|
row.status = "ERR"
|
|
|
|
row.error = str(exc)
|
|
|
|
row.error = str(exc)
|
|
|
|
row.updated_at = time.time()
|
|
|
|
row.updated_at = time.time()
|
|
|
|
|
|
|
|
|
|
|
|
render(
|
|
|
|
render(
|
|
|
|
windows=windows,
|
|
|
|
windows=windows,
|
|
|
|
backend=args.backend,
|
|
|
|
serial=serial,
|
|
|
|
device_ref=device_ref,
|
|
|
|
index=index,
|
|
|
|
index_ref=index_ref,
|
|
|
|
sample_rate=args.sample_rate,
|
|
|
|
sample_rate=effective_sample_rate,
|
|
|
|
|
|
|
|
base_mhz=args.base,
|
|
|
|
base_mhz=args.base,
|
|
|
|
roof_mhz=args.roof,
|
|
|
|
roof_mhz=args.roof,
|
|
|
|
step_mhz=args.step,
|
|
|
|
step_mhz=args.step,
|
|
|
|
@ -646,12 +787,13 @@ def main() -> int:
|
|
|
|
pass_no=pass_no,
|
|
|
|
pass_no=pass_no,
|
|
|
|
current_seq=current_seq,
|
|
|
|
current_seq=current_seq,
|
|
|
|
)
|
|
|
|
)
|
|
|
|
|
|
|
|
if args.passes > 0 and pass_no >= args.passes:
|
|
|
|
|
|
|
|
break
|
|
|
|
except Exception as exc:
|
|
|
|
except Exception as exc:
|
|
|
|
print(f"scanner failed: {exc}", file=sys.stderr)
|
|
|
|
print(f"scanner failed: {exc}", file=sys.stderr)
|
|
|
|
return 5
|
|
|
|
return 5
|
|
|
|
finally:
|
|
|
|
finally:
|
|
|
|
if isinstance(probe, HackRfWideProbe):
|
|
|
|
if probe is not None:
|
|
|
|
try:
|
|
|
|
try:
|
|
|
|
probe.stop()
|
|
|
|
probe.stop()
|
|
|
|
probe.wait()
|
|
|
|
probe.wait()
|
|
|
|
|
Sergey Revyakin