sdrangel/plugins/channelrx/demodam/amdemod.h

222 wiersze
5.9 KiB
C
Czysty Zwykły widok Historia

2015-05-12 00:53:35 +00:00
///////////////////////////////////////////////////////////////////////////////////
// Copyright (C) 2015 Edouard Griffiths, F4EXB. //
// //
// 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 as version 3 of the License, or //
// //
// 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 V3 for more details. //
// //
// You should have received a copy of the GNU General Public License //
// along with this program. If not, see <http://www.gnu.org/licenses/>. //
///////////////////////////////////////////////////////////////////////////////////
#ifndef INCLUDE_AMDEMOD_H
#define INCLUDE_AMDEMOD_H
#include <dsp/basebandsamplesink.h>
#include <QMutex>
2015-05-12 00:53:35 +00:00
#include <vector>
#include "dsp/nco.h"
#include "dsp/interpolator.h"
#include "dsp/movingaverage.h"
2015-05-12 10:12:13 +00:00
#include "dsp/agc.h"
2015-05-12 00:53:35 +00:00
#include "audio/audiofifo.h"
#include "util/message.h"
class AMDemod : public BasebandSampleSink {
Q_OBJECT
2015-05-12 00:53:35 +00:00
public:
AMDemod();
2015-05-12 00:53:35 +00:00
~AMDemod();
void configure(MessageQueue* messageQueue, Real rfBandwidth, Real volume, Real squelch, bool audioMute);
2015-05-12 00:53:35 +00:00
virtual void feed(const SampleVector::const_iterator& begin, const SampleVector::const_iterator& end, bool po);
2015-08-17 06:29:34 +00:00
virtual void start();
virtual void stop();
virtual bool handleMessage(const Message& cmd);
2015-05-12 00:53:35 +00:00
Real getMagSq() const { return m_magsq; }
bool getSquelchOpen() const { return m_squelchOpen; }
void getMagSqLevels(Real& avg, Real& peak, int& nbSamples)
{
avg = m_magsqSum / m_magsqCount;
peak = m_magsqPeak;
nbSamples = m_magsqCount;
m_magsqSum = 0.0f;
m_magsqPeak = 0.0f;
m_magsqCount = 0;
}
2015-05-12 00:53:35 +00:00
private:
class MsgConfigureAMDemod : public Message {
MESSAGE_CLASS_DECLARATION
public:
Real getRFBandwidth() const { return m_rfBandwidth; }
Real getVolume() const { return m_volume; }
Real getSquelch() const { return m_squelch; }
2015-12-26 17:23:55 +00:00
bool getAudioMute() const { return m_audioMute; }
2015-05-12 00:53:35 +00:00
static MsgConfigureAMDemod* create(Real rfBandwidth, Real volume, Real squelch, bool audioMute)
2015-05-12 00:53:35 +00:00
{
return new MsgConfigureAMDemod(rfBandwidth, volume, squelch, audioMute);
2015-05-12 00:53:35 +00:00
}
private:
Real m_rfBandwidth;
Real m_volume;
Real m_squelch;
2015-12-26 17:23:55 +00:00
bool m_audioMute;
2015-05-12 00:53:35 +00:00
MsgConfigureAMDemod(Real rfBandwidth, Real volume, Real squelch, bool audioMute) :
2015-05-12 00:53:35 +00:00
Message(),
m_rfBandwidth(rfBandwidth),
m_volume(volume),
2015-12-26 17:23:55 +00:00
m_squelch(squelch),
m_audioMute(audioMute)
2015-05-12 00:53:35 +00:00
{ }
};
struct AudioSample {
qint16 l;
qint16 r;
};
typedef std::vector<AudioSample> AudioVector;
enum RateState {
RSInitialFill,
RSRunning
};
struct Config {
int m_inputSampleRate;
qint64 m_inputFrequencyOffset;
Real m_rfBandwidth;
Real m_squelch;
Real m_volume;
quint32 m_audioSampleRate;
2015-12-26 17:23:55 +00:00
bool m_audioMute;
2015-05-12 00:53:35 +00:00
Config() :
m_inputSampleRate(-1),
m_inputFrequencyOffset(0),
m_rfBandwidth(-1),
m_squelch(0),
m_volume(0),
2015-12-26 17:23:55 +00:00
m_audioSampleRate(0),
m_audioMute(false)
2015-05-12 00:53:35 +00:00
{ }
};
Config m_config;
Config m_running;
NCO m_nco;
Interpolator m_interpolator;
Real m_interpolatorDistance;
Real m_interpolatorDistanceRemain;
Real m_squelchLevel;
int m_squelchCount;
bool m_squelchOpen;
Real m_magsq;
Real m_magsqSum;
Real m_magsqPeak;
int m_magsqCount;
2015-05-12 00:53:35 +00:00
2015-06-20 07:28:57 +00:00
MovingAverage<Real> m_movingAverage;
2015-05-12 10:12:13 +00:00
SimpleAGC m_volumeAGC;
2015-05-12 00:53:35 +00:00
AudioVector m_audioBuffer;
uint m_audioBufferFill;
AudioFifo m_audioFifo;
QMutex m_settingsMutex;
2015-05-12 00:53:35 +00:00
void apply();
void processOneSample(Complex &ci)
{
Real magsq = ci.real() * ci.real() + ci.imag() * ci.imag();
magsq /= (1<<30);
m_movingAverage.feed(magsq);
m_magsq = m_movingAverage.average();
m_magsqSum += magsq;
if (magsq > m_magsqPeak)
{
m_magsqPeak = magsq;
}
m_magsqCount++;
if (m_magsq >= m_squelchLevel)
{
if (m_squelchCount <= m_running.m_audioSampleRate / 10)
{
m_squelchCount++;
}
}
else
{
if (m_squelchCount > 1)
{
m_squelchCount -= 2;
}
}
qint16 sample;
if ((m_squelchCount >= m_running.m_audioSampleRate / 20) && !m_running.m_audioMute)
{
Real demod = sqrt(magsq);
if (demod < -1)
{
demod = -1;
}
else if (demod > 1)
{
demod = 1;
}
m_volumeAGC.feed(demod);
Real attack = (m_squelchCount - (m_running.m_audioSampleRate / 20)) / (Real) (m_running.m_audioSampleRate / 20);
demod *= ((0.003 * attack) / m_volumeAGC.getValue());
demod *= m_running.m_volume;
sample = demod * 32700 * 16;
m_squelchOpen = true;
}
else
{
sample = 0;
m_squelchOpen = false;
}
m_audioBuffer[m_audioBufferFill].l = sample;
m_audioBuffer[m_audioBufferFill].r = sample;
++m_audioBufferFill;
if (m_audioBufferFill >= m_audioBuffer.size())
{
uint res = m_audioFifo.write((const quint8*)&m_audioBuffer[0], m_audioBufferFill, 10);
if (res != m_audioBufferFill)
{
qDebug("AMDemod::feed: %u/%u audio samples written", res, m_audioBufferFill);
}
m_audioBufferFill = 0;
}
}
2015-05-12 00:53:35 +00:00
};
#endif // INCLUDE_AMDEMOD_H