filter = 4;
f_mod = 175e3;
snr = 10;


ais_settings();

training_sequence = [0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1];
U = [training_sequence, 0 1 1 1 1 1 1 0 1 1 0 0 1 1 0 1 0 1 1 1 1 1 1 0];


U = 2*(U-0.5);


S_x = ais_modulate(U, f_mod(k));

S_y = ais_awgn(S_x, snr);
I_y = real(S_y);
Q_y = imag(S_y);

[I_bb, Q_bb] = ais_chansel(I_y, Q_y, 'AIS1', filter);

[f_est1, out1] = ais_freqdet_arctan(I_bb, Q_bb, 300, 12);

figure(1); plot(out1)