blob: 3d46a9e753e92524cfa1b1b08df2ccf6e63e25ac (
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
|
import asyncio
import pyaudio
import numpy as np
import scipy.signal as signal
from rtlsdr import RtlSdr
pa = pyaudio.PyAudio()
stream = pa.open(format=pyaudio.paInt16, channels=1, rate=50000, output=True)
async def streaming():
sdr = RtlSdr()
sdr.sample_rate = 1.2e6
sdr.center_freq = 91.8e6
prev = None
stream.write(np.array([0] * 100000).astype('int16'))
stream.start_stream()
async for samples in sdr.stream():
samples = np.array(samples).astype('complex64')
samples = signal.decimate(samples, int(1.2e6 / 200e3))
if prev is not None:
samples = np.insert(samples, 0, prev)
prev = samples[-1]
samples = np.angle(samples[1:] * np.conj(samples[:-1]))
x = np.exp(-1 / (200e3 * 75e-6))
samples = signal.lfilter([1-x], [1,-x], samples)
samples = signal.decimate(samples, int(200e3 / 50e3))
samples *= 10000
stream.write(samples.astype('int16'))
await sdr.stop()
sdr.close()
loop = asyncio.get_event_loop()
loop.run_until_complete(streaming())
|
