VE7UWU/sdrangel
1
0
Fork 0
mirror of https://github.com/f4exb/sdrangel.git synced 2024-11-15 12:51:49 -05:00
Code Issues Projects Releases Wiki Activity

Tx ph.2: Interpolator polyphase filter: add possibility to specify the number of taps per phase

Browse Source
This commit is contained in:
f4exb 2016-10-31 23:40:46 +01:00
parent 283742cada
commit 3173bc0b07
2 changed files with 49 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

46
sdrbase/dsp/interpolator.cpp
View File

@ -3,20 +3,36 @@
#include <vector> #include <vector>
#include "dsp/interpolator.h" #include "dsp/interpolator.h"
static std::vector<Real> createPolyphaseLowPass(
int phaseSteps, void Interpolator::createPolyphaseLowPass(
double gain, std::vector<Real>& taps,
double sampleRateHz, int phaseSteps,
double cutoffFreqHz, double gain,
double transitionWidthHz, double sampleRateHz,
double oobAttenuationdB) double cutoffFreqHz,
double transitionWidthHz,
double oobAttenuationdB)
{ {
int ntaps = (int)(oobAttenuationdB * sampleRateHz / (22.0 * transitionWidthHz)); double nbTapsPerPhase = (oobAttenuationdB * sampleRateHz) / (22.0 * transitionWidthHz * phaseSteps);
return createPolyphaseLowPass(taps, phaseSteps, gain, sampleRateHz, cutoffFreqHz, nbTapsPerPhase);
}
void Interpolator::createPolyphaseLowPass(
std::vector<Real>& taps,
int phaseSteps,
double gain,
double sampleRateHz,
double cutoffFreqHz,
double nbTapsPerPhase)
{
int ntaps = (int)(nbTapsPerPhase * phaseSteps);
qDebug("Interpolator::createPolyphaseLowPass: ntaps: %d", ntaps);
if((ntaps % 2) != 0) if((ntaps % 2) != 0)
ntaps++; ntaps++;
ntaps *= phaseSteps; ntaps *= phaseSteps;
std::vector<float> taps(ntaps); taps.resize(ntaps);
std::vector<float> window(ntaps); std::vector<float> window(ntaps);
for(int n = 0; n < ntaps; n++) for(int n = 0; n < ntaps; n++)
@ -37,8 +53,6 @@ static std::vector<Real> createPolyphaseLowPass(
for(int i = 0; i < ntaps; i++) for(int i = 0; i < ntaps; i++)
taps[i] *= gain; taps[i] *= gain;
return taps;
} }
Interpolator::Interpolator() : Interpolator::Interpolator() :
@ -52,17 +66,19 @@ Interpolator::~Interpolator()
free(); free();
} }
void Interpolator::create(int phaseSteps, double sampleRate, double cutoff) void Interpolator::create(int phaseSteps, double sampleRate, double cutoff, double nbTapsPerPhase)
{ {
free(); free();
std::vector<Real> taps = createPolyphaseLowPass( std::vector<Real> taps;
createPolyphaseLowPass(
taps,
phaseSteps, // number of polyphases phaseSteps, // number of polyphases
1.0, // gain 1.0, // gain
phaseSteps * sampleRate, // sampling frequency phaseSteps * sampleRate, // sampling frequency
cutoff, // hz beginning of transition band cutoff, // hz beginning of transition band
sampleRate / 5.0, // hz width of transition band nbTapsPerPhase);
20.0); // out of band attenuation
// init state // init state
m_ptr = 0; m_ptr = 0;

19
sdrbase/dsp/interpolator.h
View File

@ -16,7 +16,7 @@ public:
Interpolator(); Interpolator();
~Interpolator(); ~Interpolator();
void create(int phaseSteps, double sampleRate, double cutoff); void create(int phaseSteps, double sampleRate, double cutoff, double nbTapsPerPhase = 4.5);
void free(); void free();
// Original code allowed for upsampling, but was never used that way // Original code allowed for upsampling, but was never used that way
@ -85,6 +85,23 @@ private:
int m_phaseSteps; int m_phaseSteps;
int m_nTaps; int m_nTaps;
static void createPolyphaseLowPass(
std::vector<Real>& taps,
int phaseSteps,
double gain,
double sampleRateHz,
double cutoffFreqHz,
double transitionWidthHz,
double oobAttenuationdB);
static void createPolyphaseLowPass(
std::vector<Real>& taps,
int phaseSteps,
double gain,
double sampleRateHz,
double cutoffFreqHz,
double nbTapsPerPhase);
void createTaps(int nTaps, double sampleRate, double cutoff, std::vector<Real>* taps); void createTaps(int nTaps, double sampleRate, double cutoff, std::vector<Real>* taps);
void advanceFilter(const Complex& next) void advanceFilter(const Complex& next)