2015-08-27 12:59:03 -04:00
|
|
|
/*
|
|
|
|
* Copyright (C) 2015 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 "dsp_squelch.hpp"
|
|
|
|
|
|
|
|
#include <cstdint>
|
|
|
|
#include <array>
|
|
|
|
|
2016-01-31 03:34:24 -05:00
|
|
|
bool FMSquelch::execute(const buffer_f32_t& audio) {
|
2023-05-18 16:16:05 -04:00
|
|
|
if (threshold_squared == 0.0f) {
|
|
|
|
return true;
|
|
|
|
}
|
2016-01-31 03:34:24 -05:00
|
|
|
|
2023-09-08 13:41:09 -04:00
|
|
|
// TODO: alloca temp buffer, assert audio.count
|
|
|
|
std::array<float, 32> squelch_energy_buffer;
|
|
|
|
const buffer_f32_t squelch_energy{squelch_energy_buffer.data(), audio.count};
|
2023-05-18 16:16:05 -04:00
|
|
|
non_audio_hpf.execute(audio, squelch_energy);
|
2015-08-27 12:59:03 -04:00
|
|
|
|
2023-08-31 02:05:49 -04:00
|
|
|
// "Non-audio" implies "noise" here. Find the loudest noise sample.
|
2023-05-18 16:16:05 -04:00
|
|
|
float non_audio_max_squared = 0;
|
2023-09-08 13:41:09 -04:00
|
|
|
for (size_t i = 0; i < squelch_energy.count; ++i) {
|
|
|
|
auto sample = squelch_energy.p[i];
|
|
|
|
float sample_squared = sample * sample;
|
|
|
|
|
|
|
|
if (sample_squared > non_audio_max_squared)
|
2023-05-18 16:16:05 -04:00
|
|
|
non_audio_max_squared = sample_squared;
|
|
|
|
}
|
2015-08-27 12:59:03 -04:00
|
|
|
|
2023-08-31 02:05:49 -04:00
|
|
|
// Is the noise less than the threshold?
|
2023-05-18 16:16:05 -04:00
|
|
|
return (non_audio_max_squared < threshold_squared);
|
2016-01-31 03:34:24 -05:00
|
|
|
}
|
|
|
|
|
|
|
|
void FMSquelch::set_threshold(const float new_value) {
|
2023-05-18 16:16:05 -04:00
|
|
|
threshold_squared = new_value * new_value;
|
2015-08-27 12:59:03 -04:00
|
|
|
}
|
2023-08-31 02:05:49 -04:00
|
|
|
|
|
|
|
bool FMSquelch::enabled() const {
|
|
|
|
return threshold_squared > 0.0;
|
|
|
|
}
|