1
0
mirror of https://github.com/f4exb/sdrangel.git synced 2024-11-05 16:31:15 -05:00
sdrangel/plugins/channelrx/demoddatv/datvdemod.cpp

986 lines
29 KiB
C++
Raw Normal View History

2018-02-22 16:52:49 -05:00
///////////////////////////////////////////////////////////////////////////////////
// Copyright (C) 2018 F4HKW //
// for F4EXB / SDRAngel //
// using LeanSDR Framework (C) 2016 F4DAV //
// //
// 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/>. //
///////////////////////////////////////////////////////////////////////////////////
#include "datvdemod.h"
#include <QTime>
#include <QDebug>
#include <stdio.h>
#include <complex.h>
#include "audio/audiooutput.h"
#include "dsp/dspengine.h"
#include "dsp/downchannelizer.h"
#include "dsp/threadedbasebandsamplesink.h"
#include "device/devicesourceapi.h"
const QString DATVDemod::m_channelIdURI = "sdrangel.channel.demoddatv";
const QString DATVDemod::m_channelId = "DATVDemod";
MESSAGE_CLASS_DEFINITION(DATVDemod::MsgConfigureDATVDemod, Message)
MESSAGE_CLASS_DEFINITION(DATVDemod::MsgConfigureChannelizer, Message)
DATVDemod::DATVDemod(DeviceSourceAPI *deviceAPI) :
2018-02-25 19:04:45 -05:00
ChannelSinkAPI(m_channelIdURI),
m_blnNeedConfigUpdate(false),
m_deviceAPI(deviceAPI),
m_objRegisteredDATVScreen(NULL),
m_objRegisteredVideoRender(NULL),
m_objVideoStream(NULL),
m_objRenderThread(NULL),
m_blnRenderingVideo(false),
m_blnStartStopVideo(false),
m_enmModulation(BPSK /*DATV_FM1*/),
m_objSettingsMutex(QMutex::NonRecursive)
2018-02-22 16:52:49 -05:00
{
setObjectName("DATVDemod");
//*************** DATV PARAMETERS ***************
2018-02-25 19:04:45 -05:00
m_blnInitialized=false;
CleanUpDATVFramework(false);
2018-02-22 16:52:49 -05:00
m_objVideoStream = new DATVideostream();
m_objRFFilter = new fftfilt(-256000.0 / 1024000.0, 256000.0 / 1024000.0, rfFilterFftLength);
m_channelizer = new DownChannelizer(this);
m_threadedChannelizer = new ThreadedBasebandSampleSink(m_channelizer, this);
m_deviceAPI->addThreadedSink(m_threadedChannelizer);
m_deviceAPI->addChannelAPI(this);
connect(m_channelizer, SIGNAL(inputSampleRateChanged()), this, SLOT(channelSampleRateChanged()));
2018-02-25 19:04:45 -05:00
//To setup correct Sample Rate
channelSampleRateChanged();
2018-02-22 16:52:49 -05:00
}
DATVDemod::~DATVDemod()
{
2018-02-25 19:04:45 -05:00
m_blnInitialized=false;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
if(m_objVideoStream!=NULL)
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
//Immediately exit from DATVideoStream if waiting for data before killing thread
m_objVideoStream->ThreadTimeOut=0;
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
if(m_objRenderThread!=NULL)
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
if(m_objRenderThread->isRunning())
{
m_objRenderThread->stopRendering();
m_objRenderThread->quit();
}
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
m_objRenderThread->wait(2000);
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
CleanUpDATVFramework(true);
2018-02-22 16:52:49 -05:00
m_deviceAPI->removeChannelAPI(this);
m_deviceAPI->removeThreadedSink(m_threadedChannelizer);
delete m_threadedChannelizer;
delete m_channelizer;
}
2018-02-25 19:04:45 -05:00
void DATVDemod::channelSampleRateChanged()
{
qDebug() << "DATVDemod::channelSampleRateChanged:"
<< " intMsps: " << m_channelizer->getInputSampleRate();
if(m_objRunning.intMsps!=m_channelizer->getInputSampleRate())
{
m_objRunning.intMsps = m_channelizer->getInputSampleRate();
m_objRunning.intSampleRate = m_objRunning.intMsps;
ApplySettings();
}
}
2018-02-22 16:52:49 -05:00
bool DATVDemod::SetDATVScreen(DATVScreen *objScreen)
{
m_objRegisteredDATVScreen = objScreen;
2018-02-24 18:07:08 -05:00
return true;
2018-02-22 16:52:49 -05:00
}
DATVideostream * DATVDemod::SetVideoRender(DATVideoRender *objScreen)
{
m_objRegisteredVideoRender = objScreen;
2018-02-25 19:04:45 -05:00
m_objRenderThread = new DATVideoRenderThread(m_objRegisteredVideoRender,m_objVideoStream);
2018-02-22 16:52:49 -05:00
return m_objVideoStream;
}
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
bool DATVDemod::PlayVideo(bool blnStartStop)
{
2018-02-25 19:04:45 -05:00
if(m_objVideoStream==NULL)
2018-02-22 16:52:49 -05:00
{
return false;
}
2018-02-25 19:04:45 -05:00
if(m_objRegisteredVideoRender==NULL)
2018-02-22 16:52:49 -05:00
{
return false;
}
2018-02-25 19:04:45 -05:00
if(m_objRenderThread==NULL)
2018-02-22 16:52:49 -05:00
{
return false;
}
2018-02-25 19:04:45 -05:00
if (m_blnStartStopVideo && !blnStartStop)
2018-02-22 16:52:49 -05:00
{
return true;
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
if(blnStartStop==true)
{
m_blnStartStopVideo=true;
}
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
if(m_objRenderThread->isRunning())
{
if(blnStartStop==true)
{
m_objRenderThread->stopRendering();
}
return true;
}
if(m_objVideoStream->bytesAvailable()>0)
{
m_objRenderThread->setStreamAndRenderer(m_objRegisteredVideoRender,m_objVideoStream);
m_objVideoStream->MultiThreaded=true;
m_objVideoStream->ThreadTimeOut=5000; //5000 ms
m_objRenderThread->start();
}
2018-02-22 16:52:49 -05:00
return true;
}
2018-02-25 19:04:45 -05:00
void DATVDemod::configure(MessageQueue* objMessageQueue,
int intRFBandwidth,
int intCenterFrequency,
dvb_version enmStandard,
DATVModulation enmModulation,
code_rate enmFEC,
int intSymbolRate,
int intNotchFilters,
bool blnAllowDrift,
bool blnFastLock,
dvb_sampler enmFilter,
bool blnHardMetric,
float fltRollOff,
bool blnViterbi,
int intExcursion)
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
Message* msgCmd = MsgConfigureDATVDemod::create(intRFBandwidth,intCenterFrequency,enmStandard, enmModulation, enmFEC, intSymbolRate, intNotchFilters, blnAllowDrift,blnFastLock,enmFilter,blnHardMetric,fltRollOff, blnViterbi,intExcursion);
2018-02-22 16:52:49 -05:00
objMessageQueue->push(msgCmd);
}
2018-02-25 19:04:45 -05:00
void DATVDemod::InitDATVParameters(int intMsps,
int intRFBandwidth,
int intCenterFrequency,
dvb_version enmStandard,
DATVModulation enmModulation,
code_rate enmFEC,
int intSampleRate,
int intSymbolRate,
int intNotchFilters,
bool blnAllowDrift,
bool blnFastLock,
dvb_sampler enmFilter,
bool blnHardMetric,
float fltRollOff,
bool blnViterbi,
int intExcursion)
2018-02-22 16:52:49 -05:00
{
Real fltLowCut;
Real fltHiCut;
m_objSettingsMutex.lock();
2018-02-25 19:04:45 -05:00
m_blnInitialized=false;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
//Bandpass filter shaping
fltLowCut = -((float)intRFBandwidth / 2.0) / (float)intMsps;
fltHiCut = ((float)intRFBandwidth / 2.0) / (float)intMsps;
2018-02-22 16:52:49 -05:00
m_objRFFilter->create_filter(fltLowCut, fltHiCut);
2018-02-25 19:04:45 -05:00
m_objNCO.setFreq(-(float)intCenterFrequency,(float)intMsps);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
//Config update
2018-02-22 16:52:49 -05:00
m_objRunning.intMsps = intMsps;
m_objRunning.intCenterFrequency = intCenterFrequency;
m_objRunning.intRFBandwidth = intRFBandwidth;
m_objRunning.enmStandard = enmStandard;
m_objRunning.enmModulation = enmModulation;
m_objRunning.enmFEC = enmFEC;
m_objRunning.intSampleRate = intSampleRate;
m_objRunning.intSymbolRate = intSymbolRate;
m_objRunning.intNotchFilters = intNotchFilters;
m_objRunning.blnAllowDrift = blnAllowDrift;
m_objRunning.blnFastLock = blnFastLock;
2018-02-25 19:04:45 -05:00
m_objRunning.enmFilter = enmFilter;
2018-02-22 16:52:49 -05:00
m_objRunning.blnHardMetric = blnHardMetric;
2018-02-25 19:04:45 -05:00
m_objRunning.fltRollOff = fltRollOff;
2018-02-22 16:52:49 -05:00
m_objRunning.blnViterbi = blnViterbi;
2018-02-25 19:04:45 -05:00
m_objRunning.intExcursion = intExcursion;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
m_blnInitialized=true;
2018-02-22 16:52:49 -05:00
m_objSettingsMutex.unlock();
2018-02-25 19:04:45 -05:00
m_blnNeedConfigUpdate=true;
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
void DATVDemod::CleanUpDATVFramework(bool blnRelease)
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
if(blnRelease==true)
{
if(m_objScheduler!=NULL)
{
m_objScheduler->shutdown();
delete m_objScheduler;
}
// NOTCH FILTER
if(r_auto_notch!=NULL) delete r_auto_notch;
if(p_autonotched!=NULL) delete p_autonotched;
// FREQUENCY CORRECTION : DEROTATOR
if(p_derot!=NULL) delete p_derot;
if(r_derot!=NULL) delete r_derot;
// CNR ESTIMATION
if(p_cnr!=NULL) delete p_cnr;
if(r_cnr!=NULL) delete r_cnr;
//FILTERING
if(r_resample!=NULL) delete r_resample;
if(p_resampled!=NULL) delete p_resampled;
if(coeffs!=NULL) delete coeffs;
// OUTPUT PREPROCESSED DATA
if(sampler!=NULL) delete sampler;
if(coeffs_sampler!=NULL) delete coeffs_sampler;
if(p_symbols!=NULL) delete p_symbols;
if(p_freq!=NULL) delete p_freq;
if(p_ss!=NULL) delete p_ss;
if(p_mer!=NULL) delete p_mer;
if(p_sampled!=NULL) delete p_sampled;
//DECIMATION
if(p_decimated!=NULL) delete p_decimated;
if(p_decim!=NULL) delete p_decim;
if(r_ppout!=NULL) delete r_ppout;
//GENERIC CONSTELLATION RECEIVER
if(m_objDemodulator!=NULL) delete m_objDemodulator;
//DECONVOLUTION AND SYNCHRONIZATION
if(p_bytes!=NULL) delete p_bytes;
if(r_deconv!=NULL) delete r_deconv;
if(r!=NULL) delete r;
if(p_descrambled!=NULL) delete p_descrambled;
if(p_frames!=NULL) delete p_frames;
if(r_etr192_descrambler!=NULL) delete r_etr192_descrambler;
if(r_sync!=NULL) delete r_sync;
if(p_mpegbytes!=NULL) delete p_mpegbytes;
if(p_lock!=NULL) delete p_lock;
if(p_locktime!=NULL) delete p_locktime;
if(r_sync_mpeg!=NULL) delete r_sync_mpeg;
// DEINTERLEAVING
if(p_rspackets!=NULL) delete p_rspackets;
if(r_deinter!=NULL) delete r_deinter;
if(p_vbitcount!=NULL) delete p_vbitcount;
if(p_verrcount!=NULL) delete p_verrcount;
if(p_rtspackets!=NULL) delete p_rtspackets;
if(r_rsdec!=NULL) delete r_rsdec;
//BER ESTIMATION
if(p_vber!=NULL) delete p_vber;
if(r_vber!=NULL) delete r_vber;
// DERANDOMIZATION
if(p_tspackets!=NULL) delete p_tspackets;
if(r_derand!=NULL) delete r_derand;
//OUTPUT : To remove
if(r_stdout!=NULL) delete r_stdout;
if(r_videoplayer!=NULL) delete r_videoplayer;
//CONSTELLATION
if(r_scope_symbols!=NULL) delete r_scope_symbols;
// INPUT
//if(p_rawiq!=NULL) delete p_rawiq;
//if(p_rawiq_writer!=NULL) delete p_rawiq_writer;
//if(p_preprocessed!=NULL) delete p_preprocessed;
}
m_objScheduler=NULL;
2018-02-22 16:52:49 -05:00
// INPUT
2018-02-25 19:04:45 -05:00
p_rawiq = NULL;
p_rawiq_writer = NULL;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
p_preprocessed = NULL;
2018-02-22 16:52:49 -05:00
// NOTCH FILTER
2018-02-25 19:04:45 -05:00
r_auto_notch = NULL;
p_autonotched = NULL;
2018-02-22 16:52:49 -05:00
// FREQUENCY CORRECTION : DEROTATOR
2018-02-25 19:04:45 -05:00
p_derot = NULL;
r_derot=NULL;
2018-02-22 16:52:49 -05:00
// CNR ESTIMATION
2018-02-25 19:04:45 -05:00
p_cnr = NULL;
r_cnr = NULL;
2018-02-22 16:52:49 -05:00
//FILTERING
2018-02-25 19:04:45 -05:00
r_resample = NULL;
p_resampled = NULL;
coeffs = NULL;
ncoeffs=0;
2018-02-22 16:52:49 -05:00
// OUTPUT PREPROCESSED DATA
2018-02-25 19:04:45 -05:00
sampler = NULL;
coeffs_sampler=NULL;
ncoeffs_sampler=0;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
p_symbols = NULL;
p_freq = NULL;
p_ss = NULL;
p_mer = NULL;
p_sampled = NULL;
2018-02-22 16:52:49 -05:00
//DECIMATION
2018-02-25 19:04:45 -05:00
p_decimated = NULL;
p_decim = NULL;
r_ppout = NULL;
2018-02-22 16:52:49 -05:00
//GENERIC CONSTELLATION RECEIVER
2018-02-25 19:04:45 -05:00
m_objDemodulator = NULL;
2018-02-22 16:52:49 -05:00
//DECONVOLUTION AND SYNCHRONIZATION
2018-02-25 19:04:45 -05:00
p_bytes=NULL;
r_deconv=NULL;
r = NULL;
p_descrambled = NULL;
p_frames = NULL;
r_etr192_descrambler = NULL;
r_sync = NULL;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
p_mpegbytes = NULL;
p_lock = NULL;
p_locktime = NULL;
r_sync_mpeg = NULL;
2018-02-22 16:52:49 -05:00
// DEINTERLEAVING
2018-02-25 19:04:45 -05:00
p_rspackets = NULL;
r_deinter = NULL;
p_vbitcount = NULL;
p_verrcount = NULL;
p_rtspackets = NULL;
r_rsdec = NULL;
2018-02-22 16:52:49 -05:00
//BER ESTIMATION
2018-02-25 19:04:45 -05:00
p_vber = NULL;
r_vber = NULL;
2018-02-22 16:52:49 -05:00
// DERANDOMIZATION
2018-02-25 19:04:45 -05:00
p_tspackets = NULL;
r_derand = NULL;
2018-02-22 16:52:49 -05:00
//OUTPUT : To remove
2018-02-25 19:04:45 -05:00
r_stdout = NULL;
r_videoplayer = NULL;
2018-02-22 16:52:49 -05:00
//CONSTELLATION
2018-02-25 19:04:45 -05:00
r_scope_symbols = NULL;
2018-02-22 16:52:49 -05:00
}
void DATVDemod::InitDATVFramework()
{
2018-02-25 19:04:45 -05:00
m_blnDVBInitialized=false;
m_lngReadIQ=0;
CleanUpDATVFramework(false);
qDebug() << "DATVDemod::InitDATVParameters:"
<< " - Msps: " << m_objRunning.intMsps
<< " - Sample Rate: " << m_objRunning.intSampleRate
<< " - Symbol Rate: " << m_objRunning.intSymbolRate
<< " - Modulation: " << m_objRunning.enmModulation
<< " - Notch Filters: " << m_objRunning.intNotchFilters
<< " - Allow Drift: " << m_objRunning.blnAllowDrift
<< " - Fast Lock: " << m_objRunning.blnFastLock
<< " - Filter: " << m_objRunning.enmFilter
<< " - HARD METRIC: " << m_objRunning.blnHardMetric
<< " - RollOff: " << m_objRunning.fltRollOff
<< " - Viterbi: " << m_objRunning.blnViterbi
<< " - Excursion: " << m_objRunning.intExcursion;
2018-02-22 16:52:49 -05:00
m_objCfg.standard = m_objRunning.enmStandard;
m_objCfg.fec = m_objRunning.enmFEC;
m_objCfg.Fs = (float) m_objRunning.intSampleRate;
m_objCfg.Fm = (float) m_objRunning.intSymbolRate;
m_objCfg.fastlock = m_objRunning.blnFastLock;
2018-02-25 19:04:45 -05:00
m_objCfg.sampler = m_objRunning.enmFilter;
m_objCfg.rolloff=m_objRunning.fltRollOff; //0...1
m_objCfg.rrc_rej=(float) m_objRunning.intExcursion; //dB
m_objCfg.rrc_steps=0; //auto
switch(m_objRunning.enmModulation)
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
case BPSK:
m_objCfg.constellation = cstln_lut<256>::BPSK;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case QPSK:
m_objCfg.constellation = cstln_lut<256>::QPSK;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case PSK8:
m_objCfg.constellation = cstln_lut<256>::PSK8;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case APSK16:
m_objCfg.constellation = cstln_lut<256>::APSK16;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case APSK32:
m_objCfg.constellation = cstln_lut<256>::APSK32;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case APSK64E:
m_objCfg.constellation = cstln_lut<256>::APSK64E;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case QAM16:
m_objCfg.constellation = cstln_lut<256>::QAM16;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case QAM64:
m_objCfg.constellation = cstln_lut<256>::QAM64;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
case QAM256:
m_objCfg.constellation = cstln_lut<256>::QAM256;
2018-02-22 16:52:49 -05:00
break;
2018-02-25 19:04:45 -05:00
default:
m_objCfg.constellation = cstln_lut<256>::BPSK;
2018-02-22 16:52:49 -05:00
break;
}
m_objCfg.allow_drift = m_objRunning.blnAllowDrift;
m_objCfg.anf = m_objRunning.intNotchFilters;
m_objCfg.hard_metric = m_objRunning.blnHardMetric;
2018-02-25 19:04:45 -05:00
m_objCfg.sampler = m_objRunning.enmFilter;
2018-02-22 16:52:49 -05:00
m_objCfg.viterbi = m_objRunning.blnViterbi;
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
// Min buffer size for baseband data
// scopes: 1024
// ss_estimator: 1024
// anf: 4096
// cstln_receiver: reads in chunks of 128+1
BUF_BASEBAND = 4096 * m_objCfg.buf_factor;
// Min buffer size for IQ symbols
// cstln_receiver: writes in chunks of 128/omega symbols (margin 128)
// deconv_sync: reads at least 64+32
// A larger buffer improves performance significantly.
BUF_SYMBOLS = 1024 * m_objCfg.buf_factor;
// Min buffer size for unsynchronized bytes
// deconv_sync: writes 32 bytes
// mpeg_sync: reads up to 204*scan_syncs = 1632 bytes
BUF_BYTES = 2048 * m_objCfg.buf_factor;
// Min buffer size for synchronized (but interleaved) bytes
// mpeg_sync: writes 1 rspacket
// deinterleaver: reads 17*11*12+204 = 2448 bytes
BUF_MPEGBYTES = 2448 * m_objCfg.buf_factor;
// Min buffer size for packets: 1
BUF_PACKETS = m_objCfg.buf_factor;
// Min buffer size for misc measurements: 1
BUF_SLOW = m_objCfg.buf_factor;
2018-02-25 19:04:45 -05:00
m_lngExpectedReadIQ = BUF_BASEBAND;
2018-02-22 16:52:49 -05:00
m_objScheduler = new scheduler();
//***************
p_rawiq = new pipebuf<cf32>(m_objScheduler, "rawiq", BUF_BASEBAND);
p_rawiq_writer = new pipewriter<cf32>(*p_rawiq);
p_preprocessed = p_rawiq;
// NOTCH FILTER
2018-02-25 19:04:45 -05:00
if ( m_objCfg.anf>0 )
2018-02-22 16:52:49 -05:00
{
p_autonotched = new pipebuf<cf32>(m_objScheduler, "autonotched", BUF_BASEBAND);
r_auto_notch = new auto_notch<f32>(m_objScheduler, *p_preprocessed, *p_autonotched, m_objCfg.anf, 0);
p_preprocessed = p_autonotched;
}
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
// FREQUENCY CORRECTION
2018-02-25 19:04:45 -05:00
//******** -> if ( m_objCfg.Fderot>0 )
2018-02-22 16:52:49 -05:00
// CNR ESTIMATION
p_cnr = new pipebuf<f32>(m_objScheduler, "cnr", BUF_SLOW);
2018-02-25 19:04:45 -05:00
if ( m_objCfg.cnr==true )
{
2018-02-25 19:04:45 -05:00
r_cnr = new cnr_fft<f32>(m_objScheduler, *p_preprocessed, *p_cnr, m_objCfg.Fm/m_objCfg.Fs);
2018-02-22 16:52:49 -05:00
r_cnr->decimation = decimation(m_objCfg.Fs, 1); // 1 Hz
}
// FILTERING
int decim = 1;
2018-02-25 19:04:45 -05:00
//******** -> if ( m_objCfg.resample )
2018-02-22 16:52:49 -05:00
// DECIMATION
// (Unless already done in resampler)
2018-02-25 19:04:45 -05:00
//******** -> if ( !m_objCfg.resample && m_objCfg.decim>1 )
2018-02-22 16:52:49 -05:00
//Resampling FS
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
// Generic constellation receiver
p_symbols = new pipebuf<softsymbol>(m_objScheduler, "PSK soft-symbols", BUF_SYMBOLS);
2018-02-25 19:04:45 -05:00
p_freq = new pipebuf<f32> (m_objScheduler, "freq", BUF_SLOW);
p_ss = new pipebuf<f32> (m_objScheduler, "SS", BUF_SLOW);
p_mer = new pipebuf<f32> (m_objScheduler, "MER", BUF_SLOW);
p_sampled = new pipebuf<cf32> (m_objScheduler, "PSK symbols", BUF_BASEBAND);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
switch ( m_objCfg.sampler )
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
case SAMP_NEAREST:
sampler = new nearest_sampler<float>();
break;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
case SAMP_LINEAR:
sampler = new linear_sampler<float>();
break;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
case SAMP_RRC:
{
2018-02-25 19:04:45 -05:00
if ( m_objCfg.rrc_steps == 0 )
{
2018-02-22 16:52:49 -05:00
// At least 64 discrete sampling points between symbols
2018-02-25 19:04:45 -05:00
m_objCfg.rrc_steps = max(1, (int)(64*m_objCfg.Fm / m_objCfg.Fs));
}
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
float Frrc = m_objCfg.Fs * m_objCfg.rrc_steps; // Sample freq of the RRC filter
float transition = (m_objCfg.Fm/2) * m_objCfg.rolloff;
int order = m_objCfg.rrc_rej * Frrc / (22*transition);
ncoeffs_sampler = filtergen::root_raised_cosine(order, m_objCfg.Fm/Frrc, m_objCfg.rolloff, &coeffs_sampler);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
sampler = new fir_sampler<float,float>(ncoeffs_sampler, coeffs_sampler, m_objCfg.rrc_steps);
break;
}
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
default:
qCritical("DATVDemod::InitDATVFramework: Interpolator not implemented");
return;
2018-02-22 16:52:49 -05:00
}
m_objDemodulator = new cstln_receiver<f32>(m_objScheduler, sampler, *p_preprocessed, *p_symbols, p_freq, p_ss, p_mer, p_sampled);
2018-02-25 19:04:45 -05:00
if ( m_objCfg.standard == DVB_S )
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
if ( m_objCfg.constellation != cstln_lut<256>::QPSK && m_objCfg.constellation != cstln_lut<256>::BPSK )
{
fprintf(stderr, "Warning: non-standard constellation for DVB-S\n");
}
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
if ( m_objCfg.standard == DVB_S2 )
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
// For DVB-S2 testing only.
// Constellation should be determined from PL signalling.
fprintf(stderr, "DVB-S2: Testing symbol sampler only.\n");
2018-02-22 16:52:49 -05:00
}
m_objDemodulator->cstln = make_dvbs2_constellation(m_objCfg.constellation, m_objCfg.fec);
2018-02-25 19:04:45 -05:00
if ( m_objCfg.hard_metric )
{
2018-02-25 19:04:45 -05:00
m_objDemodulator->cstln->harden();
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
m_objDemodulator->set_omega(m_objCfg.Fs/m_objCfg.Fm);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
//******** if ( m_objCfg.Ftune )
//{
// m_objDemodulator->set_freq(m_objCfg.Ftune/m_objCfg.Fs);
//}
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
if ( m_objCfg.allow_drift )
{
2018-02-25 19:04:45 -05:00
m_objDemodulator->set_allow_drift(true);
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
//******** -> if ( m_objCfg.viterbi )
if ( m_objCfg.viterbi )
{
2018-02-25 19:04:45 -05:00
m_objDemodulator->pll_adjustment /= 6;
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
m_objDemodulator->meas_decimation = decimation(m_objCfg.Fs, m_objCfg.Finfo);
// TRACKING FILTERS
2018-02-25 19:04:45 -05:00
if ( r_cnr )
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
r_cnr->freq_tap = &m_objDemodulator->freq_tap;
r_cnr->tap_multiplier = 1.0 / decim;
2018-02-22 16:52:49 -05:00
}
//constellation
2018-02-25 19:04:45 -05:00
m_objRegisteredDATVScreen->resizeDATVScreen(256,256);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
r_scope_symbols = new datvconstellation<f32>(m_objScheduler, *p_sampled, -128,128, NULL, m_objRegisteredDATVScreen);
2018-02-22 16:52:49 -05:00
r_scope_symbols->decimation = 1;
r_scope_symbols->cstln = &m_objDemodulator->cstln;
// DECONVOLUTION AND SYNCHRONIZATION
p_bytes = new pipebuf<u8>(m_objScheduler, "bytes", BUF_BYTES);
r_deconv = NULL;
2018-02-25 19:04:45 -05:00
//******** -> if ( m_objCfg.viterbi )
if ( m_objCfg.viterbi )
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
if ( m_objCfg.fec==FEC23 && (m_objDemodulator->cstln->nsymbols==4 || m_objDemodulator->cstln->nsymbols==64) )
{
m_objCfg.fec = FEC46;
}
2018-02-25 19:04:45 -05:00
//To uncomment -> Linking Problem : undefined symbol: _ZN7leansdr21viterbi_dec_interfaceIhhiiE6updateEPiS2_
r = new viterbi_sync(m_objScheduler, (*p_symbols), (*p_bytes), m_objDemodulator->cstln, m_objCfg.fec);
2018-02-25 19:04:45 -05:00
if ( m_objCfg.fastlock )
{
r->resync_period = 1;
}
2018-02-22 16:52:49 -05:00
}
else
{
r_deconv = make_deconvol_sync_simple(m_objScheduler, (*p_symbols), (*p_bytes), m_objCfg.fec);
r_deconv->fastlock = m_objCfg.fastlock;
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
//******* -> if ( m_objCfg.hdlc )
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
p_mpegbytes = new pipebuf<u8> (m_objScheduler, "mpegbytes", BUF_MPEGBYTES);
p_lock = new pipebuf<int> (m_objScheduler, "lock", BUF_SLOW);
p_locktime = new pipebuf<u32> (m_objScheduler, "locktime", BUF_PACKETS);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
r_sync_mpeg = new mpeg_sync<u8,0>(m_objScheduler, *p_bytes, *p_mpegbytes, r_deconv, p_lock, p_locktime);
r_sync_mpeg->fastlock = m_objCfg.fastlock;
2018-02-22 16:52:49 -05:00
// DEINTERLEAVING
2018-02-25 19:04:45 -05:00
p_rspackets = new pipebuf< rspacket<u8> >(m_objScheduler, "RS-enc packets", BUF_PACKETS);
2018-02-22 16:52:49 -05:00
r_deinter = new deinterleaver<u8>(m_objScheduler, *p_mpegbytes, *p_rspackets);
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
// REED-SOLOMON
p_vbitcount = new pipebuf<int>(m_objScheduler, "Bits processed", BUF_PACKETS);
p_verrcount = new pipebuf<int>(m_objScheduler, "Bits corrected", BUF_PACKETS);
p_rtspackets = new pipebuf<tspacket>(m_objScheduler, "rand TS packets", BUF_PACKETS);
2018-02-25 19:04:45 -05:00
r_rsdec = new rs_decoder<u8,0> (m_objScheduler, *p_rspackets, *p_rtspackets, p_vbitcount, p_verrcount);
2018-02-22 16:52:49 -05:00
// BER ESTIMATION
2018-02-25 19:04:45 -05:00
/*
p_vber = new pipebuf<float> (m_objScheduler, "VBER", BUF_SLOW);
r_vber = new rate_estimator<float> (m_objScheduler, *p_verrcount, *p_vbitcount, *p_vber);
r_vber->sample_size = m_objCfg.Fm/2; // About twice per second, depending on CR
2018-02-22 16:52:49 -05:00
// Require resolution better than 2E-5
2018-02-25 19:04:45 -05:00
if ( r_vber->sample_size < 50000 )
2018-02-22 16:52:49 -05:00
{
r_vber->sample_size = 50000;
}
2018-02-25 19:04:45 -05:00
*/
2018-02-22 16:52:49 -05:00
// DERANDOMIZATION
p_tspackets = new pipebuf<tspacket>(m_objScheduler, "TS packets", BUF_PACKETS);
r_derand = new derandomizer(m_objScheduler, *p_rtspackets, *p_tspackets);
2018-02-25 19:04:45 -05:00
// OUTPUT
2018-02-25 19:04:45 -05:00
r_videoplayer = new datvvideoplayer<tspacket>(m_objScheduler, *p_tspackets,m_objVideoStream);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
m_blnDVBInitialized=true;
2018-02-22 16:52:49 -05:00
}
2018-02-24 18:07:08 -05:00
void DATVDemod::feed(const SampleVector::const_iterator& begin, const SampleVector::const_iterator& end, bool firstOfBurst __attribute__((unused)))
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
qint16 * ptrBufferToRelease=NULL;
2018-02-22 16:52:49 -05:00
float fltI;
float fltQ;
cf32 objIQ;
//Complex objC;
fftfilt::cmplx *objRF;
int intRFOut;
2018-02-25 19:04:45 -05:00
//********** Bis repetita : Let's rock and roll buddy ! **********
2018-02-22 16:52:49 -05:00
#ifdef EXTENDED_DIRECT_SAMPLE
qint16 * ptrBuffer;
qint32 intLen;
//********** Reading direct samples **********
SampleVector::const_iterator it = begin;
intLen = it->intLen;
ptrBuffer = it->ptrBuffer;
ptrBufferToRelease = ptrBuffer;
++it;
for(qint32 intInd=0; intInd<intLen-1; intInd +=2)
{
fltI= ((qint32) (*ptrBuffer)) << 4;
ptrBuffer ++;
fltQ= ((qint32) (*ptrBuffer)) << 4;
ptrBuffer ++;
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
#else
for (SampleVector::const_iterator it = begin; it != end; ++it /* ++it **/)
{
fltI = it->real();
fltQ = it->imag();
#endif
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
//********** demodulation **********
2018-02-25 19:04:45 -05:00
if (m_blnNeedConfigUpdate)
2018-02-22 16:52:49 -05:00
{
2018-02-25 19:04:45 -05:00
m_objSettingsMutex.lock();
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
m_blnNeedConfigUpdate=false;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
InitDATVFramework();
m_objSettingsMutex.unlock();
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
//********** iq stream ****************
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
Complex objC(fltI,fltQ);
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
objC *= m_objNCO.nextIQ();
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
intRFOut = m_objRFFilter->runFilt(objC, &objRF); // filter RF before demod
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
for (int intI = 0 ; intI < intRFOut; intI++)
{
objIQ.re = objRF->real();
objIQ.im = objRF->imag();
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
objRF ++;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
if (m_blnDVBInitialized
&& (p_rawiq_writer!=NULL)
&& (m_objScheduler!=NULL))
{
2018-02-22 16:52:49 -05:00
p_rawiq_writer->write(objIQ);
m_lngReadIQ++;
2018-02-25 19:04:45 -05:00
//Leave +1 by safety
if((m_lngReadIQ+1)>=p_rawiq_writer->writable())
{
m_objScheduler->step();
m_lngReadIQ=0;
delete p_rawiq_writer;
p_rawiq_writer = new pipewriter<cf32>(*p_rawiq);
}
2018-02-22 16:52:49 -05:00
}
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
}
//********** demodulation **********
}
if(ptrBufferToRelease!=NULL)
{
delete ptrBufferToRelease;
}
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
}
void DATVDemod::start()
{
2018-02-25 19:04:45 -05:00
2018-02-22 16:52:49 -05:00
}
void DATVDemod::stop()
{
}
bool DATVDemod::handleMessage(const Message& cmd)
{
2018-02-25 19:04:45 -05:00
qDebug() << "DATVDemod::handleMessage";
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
if (MsgConfigureDATVDemod::match(cmd))
{
MsgConfigureDATVDemod& objCfg = (MsgConfigureDATVDemod&) cmd;
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
if((objCfg.m_objMsgConfig.blnAllowDrift != m_objRunning.blnAllowDrift)
|| (objCfg.m_objMsgConfig.intRFBandwidth != m_objRunning.intRFBandwidth)
|| (objCfg.m_objMsgConfig.intCenterFrequency != m_objRunning.intCenterFrequency)
|| (objCfg.m_objMsgConfig.blnFastLock != m_objRunning.blnFastLock)
|| (objCfg.m_objMsgConfig.blnHardMetric != m_objRunning.blnHardMetric)
|| (objCfg.m_objMsgConfig.enmFilter != m_objRunning.enmFilter)
|| (objCfg.m_objMsgConfig.fltRollOff != m_objRunning.fltRollOff)
|| (objCfg.m_objMsgConfig.blnViterbi != m_objRunning.blnViterbi)
|| (objCfg.m_objMsgConfig.enmFEC != m_objRunning.enmFEC)
|| (objCfg.m_objMsgConfig.enmModulation != m_objRunning.enmModulation)
|| (objCfg.m_objMsgConfig.enmStandard != m_objRunning.enmStandard)
|| (objCfg.m_objMsgConfig.intNotchFilters != m_objRunning.intNotchFilters)
|| (objCfg.m_objMsgConfig.intSymbolRate != m_objRunning.intSymbolRate)
|| (objCfg.m_objMsgConfig.intExcursion != m_objRunning.intExcursion))
{
2018-02-22 16:52:49 -05:00
m_objRunning.blnAllowDrift = objCfg.m_objMsgConfig.blnAllowDrift;
m_objRunning.blnFastLock = objCfg.m_objMsgConfig.blnFastLock;
m_objRunning.blnHardMetric = objCfg.m_objMsgConfig.blnHardMetric;
2018-02-25 19:04:45 -05:00
m_objRunning.enmFilter = objCfg.m_objMsgConfig.enmFilter;
m_objRunning.fltRollOff = objCfg.m_objMsgConfig.fltRollOff;
2018-02-22 16:52:49 -05:00
m_objRunning.blnViterbi = objCfg.m_objMsgConfig.blnViterbi;
m_objRunning.enmFEC = objCfg.m_objMsgConfig.enmFEC;
m_objRunning.enmModulation = objCfg.m_objMsgConfig.enmModulation;
m_objRunning.enmStandard = objCfg.m_objMsgConfig.enmStandard;
m_objRunning.intNotchFilters = objCfg.m_objMsgConfig.intNotchFilters;
m_objRunning.intSymbolRate = objCfg.m_objMsgConfig.intSymbolRate;
m_objRunning.intRFBandwidth = objCfg.m_objMsgConfig.intRFBandwidth;
m_objRunning.intCenterFrequency = objCfg.m_objMsgConfig.intCenterFrequency;
2018-02-25 19:04:45 -05:00
m_objRunning.intExcursion = objCfg.m_objMsgConfig.intExcursion;
qDebug() << "ATVDemod::handleMessage: MsgConfigureDATVDemod: sampleRate: " << m_objRunning.intMsps
<< " sampleRate: " << m_objRunning.intSampleRate;
2018-02-22 16:52:49 -05:00
ApplySettings();
2018-02-25 19:04:45 -05:00
}
2018-02-22 16:52:49 -05:00
2018-02-25 19:04:45 -05:00
return true;
}
else
{
return false;
}
2018-02-22 16:52:49 -05:00
}
void DATVDemod::ApplySettings()
{
2018-02-25 19:04:45 -05:00
if(m_objRunning.intMsps==0)
2018-02-22 16:52:49 -05:00
{
return;
}
2018-02-25 19:04:45 -05:00
InitDATVParameters(m_objRunning.intMsps,
m_objRunning.intRFBandwidth,
m_objRunning.intCenterFrequency,
m_objRunning.enmStandard,
m_objRunning.enmModulation,
m_objRunning.enmFEC,
m_objRunning.intSampleRate,
m_objRunning.intSymbolRate,
m_objRunning.intNotchFilters,
m_objRunning.blnAllowDrift,
m_objRunning.blnFastLock,
m_objRunning.enmFilter,
m_objRunning.blnHardMetric,
m_objRunning.fltRollOff,
m_objRunning.blnViterbi,
m_objRunning.intExcursion);
2018-02-22 16:52:49 -05:00
}
int DATVDemod::GetSampleRate()
{
return m_objRunning.intMsps;
}