mirror of
https://github.com/eried/portapack-mayhem.git
synced 2025-02-25 09:01:12 -05:00
032168234a
* Adding_TX_BW_user_selection_to_the_Mic_SSB_mode * TX_SSB_BW_adjustm_2k_3k_4k * init_value_correction * Comment_correction * Leaving_only_OK_txbw_SSB_2K_3K * Limti_SSB_TX_BW_2K_3K
87 lines
2.3 KiB
C++
87 lines
2.3 KiB
C++
/*
|
|
* Copyright (C) 2020 Belousov Oleg
|
|
*
|
|
* 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_hilbert.hpp"
|
|
#include "dsp_sos_config.hpp"
|
|
|
|
namespace dsp {
|
|
|
|
HilbertTransform::HilbertTransform() {
|
|
n = 0;
|
|
|
|
sos_input.configure(half_band_lpf_config);
|
|
sos_i.configure(half_band_lpf_config);
|
|
sos_q.configure(half_band_lpf_config);
|
|
}
|
|
|
|
void HilbertTransform::execute(float in, float& out_i, float& out_q) {
|
|
float a = 0, b = 0;
|
|
|
|
float in_filtered = sos_input.execute(in) * 1.0f; // Anti-aliasing LPF at fs/4 mic audio filter front-end.
|
|
// it is exactly matching with the usefull mic BW_cut_off =fs/4 of the following Hilbert Transf. implementation
|
|
// Ex. to TX SSB 3khz BW (3KHZ =fs/4) ,we need to run Hiblert_Trf fs = 12khz . And our anti-aliasing filter is also fcut =fs/4 when we use a half_band_LPF .
|
|
|
|
switch (n) {
|
|
case 0:
|
|
a = in_filtered;
|
|
b = 0;
|
|
break;
|
|
case 1:
|
|
a = 0;
|
|
b = -in_filtered;
|
|
break;
|
|
case 2:
|
|
a = -in_filtered;
|
|
b = 0;
|
|
break;
|
|
case 3:
|
|
a = 0;
|
|
b = in_filtered;
|
|
break;
|
|
}
|
|
|
|
float i = sos_i.execute(a) * 2.0f;
|
|
float q = sos_q.execute(b) * 2.0f;
|
|
|
|
switch (n) {
|
|
case 0:
|
|
out_i = i;
|
|
out_q = q;
|
|
break;
|
|
case 1:
|
|
out_i = -q;
|
|
out_q = i;
|
|
break;
|
|
case 2:
|
|
out_i = -i;
|
|
out_q = -q;
|
|
break;
|
|
case 3:
|
|
out_i = q;
|
|
out_q = -i;
|
|
break;
|
|
}
|
|
|
|
n = (n + 1) % 4;
|
|
}
|
|
|
|
} /* namespace dsp */
|