mirror of
https://github.com/eried/portapack-mayhem.git
synced 2024-10-01 01:26:06 -04:00
106 lines
3.3 KiB
C++
106 lines
3.3 KiB
C++
/*
|
|
* 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 "spectrum_collector.hpp"
|
|
|
|
#include "dsp_fft.hpp"
|
|
|
|
#include "utility.hpp"
|
|
#include "event_m4.hpp"
|
|
#include "portapack_shared_memory.hpp"
|
|
|
|
#include <algorithm>
|
|
|
|
void SpectrumCollector::start() {
|
|
streaming = true;
|
|
ChannelSpectrumConfigMessage message { &fifo };
|
|
shared_memory.application_queue.push(message);
|
|
}
|
|
|
|
void SpectrumCollector::stop() {
|
|
streaming = false;
|
|
fifo.reset_in();
|
|
}
|
|
|
|
void SpectrumCollector::set_decimation_factor(
|
|
const size_t decimation_factor
|
|
) {
|
|
channel_spectrum_decimator.set_factor(decimation_factor);
|
|
}
|
|
|
|
/* TODO: Refactor to register task with idle thread?
|
|
* It's sad that the idle thread has to call all the way back here just to
|
|
* perform the deferred task on the buffer of data we prepared.
|
|
*/
|
|
|
|
void SpectrumCollector::feed(
|
|
const buffer_c16_t& channel,
|
|
const uint32_t filter_pass_frequency,
|
|
const uint32_t filter_stop_frequency
|
|
) {
|
|
// Called from baseband processing thread.
|
|
channel_filter_pass_frequency = filter_pass_frequency;
|
|
channel_filter_stop_frequency = filter_stop_frequency;
|
|
|
|
channel_spectrum_decimator.feed(
|
|
channel,
|
|
[this](const buffer_c16_t& data) {
|
|
this->post_message(data);
|
|
}
|
|
);
|
|
}
|
|
|
|
void SpectrumCollector::post_message(const buffer_c16_t& data) {
|
|
// Called from baseband processing thread.
|
|
if( streaming && !channel_spectrum_request_update ) {
|
|
fft_swap(data, channel_spectrum);
|
|
channel_spectrum_sampling_rate = data.sampling_rate;
|
|
channel_spectrum_request_update = true;
|
|
events_flag(EVT_MASK_SPECTRUM);
|
|
}
|
|
}
|
|
|
|
void SpectrumCollector::update() {
|
|
// Called from idle thread (after EVT_MASK_SPECTRUM is flagged)
|
|
if( streaming && channel_spectrum_request_update ) {
|
|
/* Decimated buffer is full. Compute spectrum. */
|
|
fft_c_preswapped(channel_spectrum);
|
|
|
|
ChannelSpectrum spectrum;
|
|
spectrum.sampling_rate = channel_spectrum_sampling_rate;
|
|
spectrum.channel_filter_pass_frequency = channel_filter_pass_frequency;
|
|
spectrum.channel_filter_stop_frequency = channel_filter_stop_frequency;
|
|
for(size_t i=0; i<spectrum.db.size(); i++) {
|
|
// Three point Hamming window.
|
|
const auto corrected_sample = channel_spectrum[i] * 0.54f
|
|
+ (channel_spectrum[(i-1) & 0xff] + channel_spectrum[(i+1) & 0xff]) * -0.23f;
|
|
const auto mag2 = magnitude_squared(corrected_sample);
|
|
const float db = complex16_mag_squared_to_dbv_norm(mag2);
|
|
constexpr float mag_scale = 5.0f;
|
|
const unsigned int v = (db * mag_scale) + 255.0f;
|
|
spectrum.db[i] = std::max(0U, std::min(255U, v));
|
|
}
|
|
fifo.in(spectrum);
|
|
}
|
|
|
|
channel_spectrum_request_update = false;
|
|
}
|