HackRF采集调频广播基带数据并使用Ocatve解调播放

clc;
clear all;
pkg load signal;

filename = ‘FM93_8Msps.wav‘
fid = fopen(filename, ‘r‘);
[data] = fread(fid, inf, ‘int8‘);
fclose(fid);

I = data(1: 2: end);
Q = data(2: 2: end);

B = fir1(131, 0.02);
A = [1];

I = filter(B, A, I);
Q = filter(B, A, Q);

I = I(1: 50: end);
Q = Q(1: 50: end);

% 160Ksps
S = I .+ Q*1i;

% Demod
d = S(1: end - 1) .* conj(S(2: end));

s = atan2(imag(d), real(d))/pi;

B = fir1(51, 0.1);
A = [1];

s = filter(B, A, s);
s = s(1: 10: end);

% 16Ksps
sound(s, 16000);