2015-08-27 16:19:34 -04:00
|
|
|
/*
|
|
|
|
* Copyright (C) 2014 Jared Boone, ShareBrained Technology, Inc.
|
|
|
|
*
|
|
|
|
* This file is part of PortaPack.
|
|
|
|
*
|
|
|
|
* This program is free software; you can redistribute it and/or modify
|
|
|
|
* it under the terms of the GNU General Public License as published by
|
|
|
|
* the Free Software Foundation; either version 2, or (at your option)
|
|
|
|
* any later version.
|
|
|
|
*
|
|
|
|
* This program is distributed in the hope that it will be useful,
|
|
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
|
|
* GNU General Public License for more details.
|
|
|
|
*
|
|
|
|
* You should have received a copy of the GNU General Public License
|
|
|
|
* along with this program; see the file COPYING. If not, write to
|
|
|
|
* the Free Software Foundation, Inc., 51 Franklin Street,
|
|
|
|
* Boston, MA 02110-1301, USA.
|
|
|
|
*/
|
|
|
|
|
|
|
|
#include "proc_wfm_audio.hpp"
|
|
|
|
|
|
|
|
#include <cstdint>
|
|
|
|
|
2016-01-02 13:35:23 -05:00
|
|
|
WidebandFMAudio::WidebandFMAudio() {
|
|
|
|
constexpr size_t baseband_fs = 3072000;
|
|
|
|
|
|
|
|
constexpr size_t decim_0_input_fs = baseband_fs;
|
|
|
|
constexpr size_t decim_0_decimation_factor = 4;
|
|
|
|
constexpr size_t decim_0_output_fs = decim_0_input_fs / decim_0_decimation_factor;
|
|
|
|
|
|
|
|
constexpr size_t decim_1_input_fs = decim_0_output_fs;
|
|
|
|
constexpr size_t decim_1_decimation_factor = 2;
|
|
|
|
constexpr size_t decim_1_output_fs = decim_1_input_fs / decim_1_decimation_factor;
|
|
|
|
|
|
|
|
constexpr size_t demod_input_fs = decim_1_output_fs;
|
|
|
|
|
|
|
|
decim_0.configure(taps_200k_wfm_decim_0.taps, 33554432);
|
|
|
|
decim_1.configure(taps_200k_wfm_decim_1.taps, 131072);
|
|
|
|
demod.configure(demod_input_fs, 75000);
|
|
|
|
audio_filter.configure(taps_64_lp_156_198.taps);
|
|
|
|
}
|
|
|
|
|
2015-12-10 15:36:12 -05:00
|
|
|
void WidebandFMAudio::execute(const buffer_c8_t& buffer) {
|
2016-01-02 13:57:47 -05:00
|
|
|
std::array<complex16_t, 512> dst;
|
2016-01-02 13:35:23 -05:00
|
|
|
const buffer_c16_t dst_buffer {
|
|
|
|
dst.data(),
|
|
|
|
dst.size()
|
|
|
|
};
|
|
|
|
|
|
|
|
const auto decim_0_out = decim_0.execute(buffer, dst_buffer);
|
|
|
|
const auto decim_1_out = decim_1.execute(decim_0_out, dst_buffer);
|
|
|
|
const auto decimator_out = decim_1_out;
|
2015-08-27 16:19:34 -04:00
|
|
|
|
|
|
|
const buffer_s16_t work_audio_buffer {
|
|
|
|
(int16_t*)decimator_out.p,
|
|
|
|
sizeof(*decimator_out.p) * decimator_out.count
|
|
|
|
};
|
|
|
|
|
|
|
|
auto channel = decimator_out;
|
|
|
|
|
|
|
|
// TODO: Feed channel_stats post-decimation data?
|
|
|
|
feed_channel_stats(channel);
|
|
|
|
//feed_channel_spectrum(channel);
|
|
|
|
|
2016-01-02 13:35:23 -05:00
|
|
|
/* 384kHz complex<int16_t>[256]
|
2015-08-27 16:19:34 -04:00
|
|
|
* -> FM demodulation
|
2016-01-02 13:35:23 -05:00
|
|
|
* -> 384kHz int16_t[256] */
|
2015-08-27 16:19:34 -04:00
|
|
|
/* TODO: To improve adjacent channel rejection, implement complex channel filter:
|
|
|
|
* pass < +/- 100kHz, stop > +/- 200kHz
|
|
|
|
*/
|
|
|
|
|
|
|
|
auto audio_oversampled = demod.execute(decimator_out, work_audio_buffer);
|
|
|
|
|
|
|
|
/* 384kHz int16_t[256]
|
|
|
|
* -> 4th order CIC decimation by 2, gain of 1
|
|
|
|
* -> 192kHz int16_t[128] */
|
2016-01-02 13:35:23 -05:00
|
|
|
auto audio_4fs = audio_dec_1.execute(audio_oversampled, work_audio_buffer);
|
2015-08-27 16:19:34 -04:00
|
|
|
|
|
|
|
/* 192kHz int16_t[128]
|
|
|
|
* -> 4th order CIC decimation by 2, gain of 1
|
|
|
|
* -> 96kHz int16_t[64] */
|
2016-01-02 13:35:23 -05:00
|
|
|
auto audio_2fs = audio_dec_2.execute(audio_4fs, work_audio_buffer);
|
2015-08-27 16:19:34 -04:00
|
|
|
|
|
|
|
/* 96kHz int16_t[64]
|
|
|
|
* -> FIR filter, <15kHz (0.156fs) pass, >19kHz (0.198fs) stop, gain of 1
|
|
|
|
* -> 48kHz int16_t[32] */
|
|
|
|
auto audio = audio_filter.execute(audio_2fs, work_audio_buffer);
|
|
|
|
|
|
|
|
/* -> 48kHz int16_t[32] */
|
|
|
|
audio_hpf.execute_in_place(audio);
|
2016-01-01 23:56:30 -05:00
|
|
|
audio_deemph.execute_in_place(audio);
|
|
|
|
|
2015-08-27 16:19:34 -04:00
|
|
|
fill_audio_buffer(audio);
|
|
|
|
}
|