pkg load signal
pkg load financial
addpath('./octave/')
sig = load_samples ('keyfob.bin', 'bladeRF');

# Don't really need FIR
# b = fir1(200, 848e3/Fs, 'low');
# plot(abs(sig));

# Part of it
sig_pkt = sig(2.648e6:2.838e6);
plot(abs(sig_pkt));
plot(abs(fftshift(fft(sig_pkt))));

phase = unwrap(angle(sig_pkt));
plot(phase);
p = movavg(phase, 100, 100);
freq = p(2:end) - p(1:end-1);
subplot(2, 1, 1);
plot(freq);

# digitize
dig = f > .96;
subplot(2, 1, 2);
plot(dig);
axis([0000 140000 -.5 1.5]);