2014-12-16 18:27:32 -05:00
|
|
|
#include "CubicSDRDefs.h"
|
2015-02-22 01:01:28 -05:00
|
|
|
#include "DemodulatorThread.h"
|
2014-12-16 18:27:32 -05:00
|
|
|
#include <vector>
|
|
|
|
|
2015-02-22 01:01:28 -05:00
|
|
|
#include <cmath>
|
|
|
|
#ifndef M_PI
|
|
|
|
#define M_PI 3.14159265358979323846
|
|
|
|
#endif
|
|
|
|
|
2014-12-16 18:27:32 -05:00
|
|
|
#ifdef __APPLE__
|
|
|
|
#include <pthread.h>
|
|
|
|
#endif
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
DemodulatorThread::DemodulatorThread(DemodulatorThreadPostInputQueue* iqInputQueue, DemodulatorThreadControlCommandQueue *threadQueueControl,
|
2014-12-24 01:28:33 -05:00
|
|
|
DemodulatorThreadCommandQueue* threadQueueNotify) :
|
2015-01-03 17:07:39 -05:00
|
|
|
iqInputQueue(iqInputQueue), audioVisOutputQueue(NULL), audioOutputQueue(NULL), iqAutoGain(NULL), amOutputCeil(1), amOutputCeilMA(1), amOutputCeilMAA(
|
|
|
|
1), stereo(false), terminated(
|
|
|
|
false), demodulatorType(DEMOD_TYPE_FM), threadQueueNotify(threadQueueNotify), threadQueueControl(threadQueueControl), squelchLevel(0), signalLevel(
|
2015-03-15 21:02:26 -04:00
|
|
|
0), squelchEnabled(false), audioSampleRate(0) {
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
demodFM = freqdem_create(0.5);
|
2015-02-12 01:38:53 -05:00
|
|
|
demodAM_USB = ampmodem_create(0.5, 0.0, LIQUID_AMPMODEM_USB, 1);
|
|
|
|
demodAM_LSB = ampmodem_create(0.5, 0.0, LIQUID_AMPMODEM_LSB, 1);
|
2015-01-07 20:23:15 -05:00
|
|
|
demodAM_DSB = ampmodem_create(0.5, 0.0, LIQUID_AMPMODEM_DSB, 1);
|
2015-01-03 17:07:39 -05:00
|
|
|
demodAM_DSB_CSP = ampmodem_create(0.5, 0.0, LIQUID_AMPMODEM_DSB, 0);
|
|
|
|
demodAM = demodAM_DSB_CSP;
|
2015-06-03 05:23:16 -04:00
|
|
|
|
|
|
|
// advanced demodulators
|
2015-06-10 04:50:14 -04:00
|
|
|
// This could properly be done easier.
|
|
|
|
|
|
|
|
demodulatorCons = 2;
|
|
|
|
|
|
|
|
demodASK = modem_create(LIQUID_MODEM_ASK2);
|
|
|
|
demodASK2 = modem_create(LIQUID_MODEM_ASK2);
|
|
|
|
demodASK4 = modem_create(LIQUID_MODEM_ASK4);
|
|
|
|
demodASK8 = modem_create(LIQUID_MODEM_ASK8);
|
|
|
|
demodASK16 = modem_create(LIQUID_MODEM_ASK16);
|
|
|
|
demodASK32 = modem_create(LIQUID_MODEM_ASK32);
|
|
|
|
demodASK64 = modem_create(LIQUID_MODEM_ASK64);
|
|
|
|
demodASK128 = modem_create(LIQUID_MODEM_ASK128);
|
|
|
|
demodASK256 = modem_create(LIQUID_MODEM_ASK256);
|
|
|
|
|
2015-06-09 11:00:05 -04:00
|
|
|
demodAPSK = modem_create(LIQUID_MODEM_APSK256);
|
2015-06-03 05:23:16 -04:00
|
|
|
demodBPSK = modem_create(LIQUID_MODEM_BPSK);
|
2015-06-09 11:00:05 -04:00
|
|
|
demodDPSK = modem_create(LIQUID_MODEM_DPSK256);
|
|
|
|
demodPSK = modem_create(LIQUID_MODEM_PSK256);
|
2015-06-03 05:23:16 -04:00
|
|
|
demodOOK = modem_create(LIQUID_MODEM_OOK);
|
2015-06-09 11:00:05 -04:00
|
|
|
demodSQAM = modem_create(LIQUID_MODEM_SQAM128);
|
2015-06-03 05:23:16 -04:00
|
|
|
demodST = modem_create(LIQUID_MODEM_V29);
|
2015-06-09 11:00:05 -04:00
|
|
|
demodQAM = modem_create(LIQUID_MODEM_QAM256);
|
2015-06-03 05:23:16 -04:00
|
|
|
demodQPSK = modem_create(LIQUID_MODEM_QPSK);
|
2015-06-09 11:00:05 -04:00
|
|
|
|
|
|
|
currentDemodLock = false;
|
|
|
|
|
2014-12-16 18:27:32 -05:00
|
|
|
}
|
|
|
|
DemodulatorThread::~DemodulatorThread() {
|
|
|
|
}
|
|
|
|
|
|
|
|
#ifdef __APPLE__
|
|
|
|
void *DemodulatorThread::threadMain() {
|
|
|
|
#else
|
2014-12-16 21:30:03 -05:00
|
|
|
void DemodulatorThread::threadMain() {
|
2014-12-16 18:27:32 -05:00
|
|
|
#endif
|
2014-12-16 20:33:44 -05:00
|
|
|
#ifdef __APPLE__
|
2014-12-18 20:11:25 -05:00
|
|
|
pthread_t tID = pthread_self(); // ID of this thread
|
|
|
|
int priority = sched_get_priority_max( SCHED_FIFO )-1;
|
2014-12-18 21:39:32 -05:00
|
|
|
sched_param prio = {priority}; // scheduling priority of thread
|
2014-12-18 20:11:25 -05:00
|
|
|
pthread_setschedparam(tID, SCHED_FIFO, &prio);
|
2014-12-16 20:33:44 -05:00
|
|
|
#endif
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
msresamp_rrrf audioResampler = NULL;
|
|
|
|
msresamp_rrrf stereoResampler = NULL;
|
|
|
|
firfilt_rrrf firStereoLeft = NULL;
|
|
|
|
firfilt_rrrf firStereoRight = NULL;
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-02-12 01:38:53 -05:00
|
|
|
liquid_float_complex x, y, z[2];
|
|
|
|
float rz[2];
|
2014-12-27 12:23:09 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
firhilbf firStereoR2C = firhilbf_create(5, 60.0f);
|
|
|
|
firhilbf firStereoC2R = firhilbf_create(5, 60.0f);
|
2015-01-01 18:08:54 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
nco_crcf stereoShifter = nco_crcf_create(LIQUID_NCO);
|
|
|
|
double stereoShiftFrequency = 0;
|
2015-01-01 18:08:54 -05:00
|
|
|
|
2015-02-19 23:51:24 -05:00
|
|
|
// half band filter used for side-band elimination
|
|
|
|
resamp2_cccf ssbFilt = resamp2_cccf_create(12,-0.25f,60.0f);
|
2015-01-25 01:27:37 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
// Automatic IQ gain
|
|
|
|
iqAutoGain = agc_crcf_create();
|
|
|
|
agc_crcf_set_bandwidth(iqAutoGain, 0.9);
|
2014-12-21 16:08:32 -05:00
|
|
|
|
2015-01-22 21:32:32 -05:00
|
|
|
AudioThreadInput *ati_vis = new AudioThreadInput;
|
|
|
|
ati_vis->data.reserve(DEMOD_VIS_SIZE);
|
|
|
|
|
2014-12-16 21:30:03 -05:00
|
|
|
std::cout << "Demodulator thread started.." << std::endl;
|
2014-12-24 01:28:33 -05:00
|
|
|
|
2015-02-03 12:25:02 -05:00
|
|
|
switch (demodulatorType.load()) {
|
2015-02-02 20:10:55 -05:00
|
|
|
case DEMOD_TYPE_FM:
|
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_LSB:
|
2015-02-12 22:46:45 -05:00
|
|
|
demodAM = demodAM_LSB;
|
2015-02-02 20:10:55 -05:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_USB:
|
2015-02-12 22:46:45 -05:00
|
|
|
demodAM = demodAM_USB;
|
2015-02-02 20:10:55 -05:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_DSB:
|
|
|
|
demodAM = demodAM_DSB;
|
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_AM:
|
|
|
|
demodAM = demodAM_DSB_CSP;
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
|
2015-01-08 21:12:49 -05:00
|
|
|
terminated = false;
|
|
|
|
|
2014-12-16 21:30:03 -05:00
|
|
|
while (!terminated) {
|
2014-12-23 01:12:14 -05:00
|
|
|
DemodulatorThreadPostIQData *inp;
|
2015-01-03 17:07:39 -05:00
|
|
|
iqInputQueue->pop(inp);
|
2014-12-24 01:28:33 -05:00
|
|
|
std::lock_guard < std::mutex > lock(inp->m_mutex);
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2014-12-23 01:12:14 -05:00
|
|
|
int bufSize = inp->data.size();
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2014-12-16 21:30:03 -05:00
|
|
|
if (!bufSize) {
|
2014-12-24 01:28:33 -05:00
|
|
|
inp->decRefCount();
|
2014-12-16 21:30:03 -05:00
|
|
|
continue;
|
|
|
|
}
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
if (audioResampler == NULL) {
|
2015-01-03 17:07:39 -05:00
|
|
|
audioResampler = inp->audioResampler;
|
|
|
|
stereoResampler = inp->stereoResampler;
|
2015-03-15 21:02:26 -04:00
|
|
|
firStereoLeft = inp->firStereoLeft;
|
|
|
|
firStereoRight = inp->firStereoRight;
|
|
|
|
audioSampleRate = inp->audioSampleRate;
|
2015-01-04 13:20:31 -05:00
|
|
|
} else if (audioResampler != inp->audioResampler) {
|
2015-01-03 17:07:39 -05:00
|
|
|
msresamp_rrrf_destroy(audioResampler);
|
|
|
|
msresamp_rrrf_destroy(stereoResampler);
|
|
|
|
audioResampler = inp->audioResampler;
|
|
|
|
stereoResampler = inp->stereoResampler;
|
2015-03-15 21:02:26 -04:00
|
|
|
audioSampleRate = inp->audioSampleRate;
|
2015-01-01 03:48:32 -05:00
|
|
|
|
2015-01-07 21:25:35 -05:00
|
|
|
if (demodAM) {
|
|
|
|
ampmodem_reset(demodAM);
|
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
freqdem_reset(demodFM);
|
2014-12-16 21:30:03 -05:00
|
|
|
}
|
2014-12-16 20:33:44 -05:00
|
|
|
|
2015-03-15 21:02:26 -04:00
|
|
|
if (firStereoLeft != inp->firStereoLeft) {
|
|
|
|
if (firStereoLeft != NULL) {
|
|
|
|
firfilt_rrrf_destroy(firStereoLeft);
|
|
|
|
}
|
|
|
|
firStereoLeft = inp->firStereoLeft;
|
|
|
|
}
|
|
|
|
|
|
|
|
if (firStereoRight != inp->firStereoRight) {
|
|
|
|
if (firStereoRight != NULL) {
|
|
|
|
firfilt_rrrf_destroy(firStereoRight);
|
|
|
|
}
|
|
|
|
firStereoRight = inp->firStereoRight;
|
|
|
|
}
|
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
if (agcData.size() != bufSize) {
|
|
|
|
if (agcData.capacity() < bufSize) {
|
|
|
|
agcData.reserve(bufSize);
|
|
|
|
agcAMData.reserve(bufSize);
|
2014-12-24 03:03:34 -05:00
|
|
|
}
|
2015-01-04 13:20:31 -05:00
|
|
|
agcData.resize(bufSize);
|
|
|
|
agcAMData.resize(bufSize);
|
2014-12-24 03:03:34 -05:00
|
|
|
}
|
2014-12-16 21:30:03 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
double audio_resample_ratio = inp->audioResampleRatio;
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
if (demodOutputData.size() != bufSize) {
|
|
|
|
if (demodOutputData.capacity() < bufSize) {
|
|
|
|
demodOutputData.reserve(bufSize);
|
2014-12-24 03:03:34 -05:00
|
|
|
}
|
2015-01-04 13:20:31 -05:00
|
|
|
demodOutputData.resize(bufSize);
|
2014-12-24 03:03:34 -05:00
|
|
|
}
|
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
int audio_out_size = ceil((double) (bufSize) * audio_resample_ratio) + 512;
|
2014-12-27 12:45:21 -05:00
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
agc_crcf_execute_block(iqAutoGain, &(inp->data[0]), bufSize, &agcData[0]);
|
2015-01-01 03:48:32 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
float currentSignalLevel = 0;
|
2015-01-01 03:48:32 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
currentSignalLevel = ((60.0 / fabs(agc_crcf_get_rssi(iqAutoGain))) / 15.0 - signalLevel);
|
2015-01-01 03:48:32 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (agc_crcf_get_signal_level(iqAutoGain) > currentSignalLevel) {
|
|
|
|
currentSignalLevel = agc_crcf_get_signal_level(iqAutoGain);
|
2015-01-01 03:48:32 -05:00
|
|
|
}
|
|
|
|
|
2015-01-01 18:08:54 -05:00
|
|
|
if (demodulatorType == DEMOD_TYPE_FM) {
|
2015-01-04 13:20:31 -05:00
|
|
|
freqdem_demodulate_block(demodFM, &agcData[0], bufSize, &demodOutputData[0]);
|
2015-01-01 18:08:54 -05:00
|
|
|
} else {
|
|
|
|
float p;
|
2015-06-04 16:38:43 -04:00
|
|
|
unsigned int bitstream;
|
2015-01-03 19:03:16 -05:00
|
|
|
switch (demodulatorType.load()) {
|
2015-01-01 18:08:54 -05:00
|
|
|
case DEMOD_TYPE_LSB:
|
2015-02-19 23:51:24 -05:00
|
|
|
for (int i = 0; i < bufSize; i++) { // Reject upper band
|
|
|
|
resamp2_cccf_filter_execute(ssbFilt,inp->data[i],&x,&y);
|
|
|
|
ampmodem_demodulate(demodAM, x, &demodOutputData[i]);
|
2015-02-12 22:46:45 -05:00
|
|
|
}
|
2015-01-01 18:08:54 -05:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_USB:
|
2015-02-19 23:51:24 -05:00
|
|
|
for (int i = 0; i < bufSize; i++) { // Reject lower band
|
|
|
|
resamp2_cccf_filter_execute(ssbFilt,inp->data[i],&x,&y);
|
|
|
|
ampmodem_demodulate(demodAM, y, &demodOutputData[i]);
|
2015-02-12 22:46:45 -05:00
|
|
|
}
|
2015-01-01 18:08:54 -05:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_AM:
|
2015-01-07 20:23:15 -05:00
|
|
|
case DEMOD_TYPE_DSB:
|
2015-02-12 02:14:22 -05:00
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
ampmodem_demodulate(demodAM, inp->data[i], &demodOutputData[i]);
|
|
|
|
}
|
2015-01-01 18:08:54 -05:00
|
|
|
break;
|
2015-06-04 16:38:43 -04:00
|
|
|
case DEMOD_TYPE_ASK:
|
2015-06-10 04:50:14 -04:00
|
|
|
if(demodulatorCons == 2) {
|
|
|
|
demodASK = demodASK2;
|
|
|
|
}
|
2015-06-04 16:38:43 -04:00
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodASK, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-10 04:50:14 -04:00
|
|
|
updateDemodulatorLock(demodASK, 0.5f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_BPSK:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodBPSK, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodBPSK, 0.8f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_DPSK:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodDPSK, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodDPSK, 0.8f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_PSK:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodPSK, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodPSK, 0.8f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_OOK:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodOOK, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodOOK, 0.8f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_SQAM:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodSQAM, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodSQAM, 0.8f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_ST:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodST, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodST, 0.8f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_QAM:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodQAM, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodQAM, 0.5f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_QPSK:
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
modem_demodulate(demodQPSK, inp->data[i], &bitstream);
|
2015-06-09 11:00:05 -04:00
|
|
|
// std::cout << bitstream << std::endl;
|
2015-06-04 16:38:43 -04:00
|
|
|
}
|
2015-06-09 11:00:05 -04:00
|
|
|
updateDemodulatorLock(demodQPSK, 0.8f);
|
2015-06-04 16:38:43 -04:00
|
|
|
break;
|
2015-01-01 03:48:32 -05:00
|
|
|
}
|
|
|
|
|
2015-03-01 10:30:02 -05:00
|
|
|
amOutputCeilMA = amOutputCeilMA + (amOutputCeil - amOutputCeilMA) * 0.025;
|
|
|
|
amOutputCeilMAA = amOutputCeilMAA + (amOutputCeilMA - amOutputCeilMAA) * 0.025;
|
2015-01-01 03:48:32 -05:00
|
|
|
|
2015-02-12 02:14:22 -05:00
|
|
|
amOutputCeil = 0;
|
|
|
|
|
|
|
|
for (int i = 0; i < bufSize; i++) {
|
|
|
|
if (demodOutputData[i] > amOutputCeil) {
|
|
|
|
amOutputCeil = demodOutputData[i];
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
float gain = 0.5 / amOutputCeilMAA;
|
2015-01-01 03:48:32 -05:00
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
for (int i = 0; i < bufSize; i++) {
|
2015-01-03 17:07:39 -05:00
|
|
|
demodOutputData[i] *= gain;
|
2015-01-01 03:48:32 -05:00
|
|
|
}
|
|
|
|
}
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (audio_out_size != resampledOutputData.size()) {
|
|
|
|
if (resampledOutputData.capacity() < audio_out_size) {
|
|
|
|
resampledOutputData.reserve(audio_out_size);
|
2014-12-27 12:45:21 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
resampledOutputData.resize(audio_out_size);
|
2014-12-27 12:45:21 -05:00
|
|
|
}
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
unsigned int numAudioWritten;
|
2015-01-04 13:20:31 -05:00
|
|
|
msresamp_rrrf_execute(audioResampler, &demodOutputData[0], bufSize, &resampledOutputData[0], &numAudioWritten);
|
2014-12-27 12:45:21 -05:00
|
|
|
|
2014-12-26 20:58:42 -05:00
|
|
|
if (stereo) {
|
2015-01-04 13:20:31 -05:00
|
|
|
if (demodStereoData.size() != bufSize) {
|
|
|
|
if (demodStereoData.capacity() < bufSize) {
|
|
|
|
demodStereoData.reserve(bufSize);
|
2014-12-27 12:45:21 -05:00
|
|
|
}
|
2015-01-04 13:20:31 -05:00
|
|
|
demodStereoData.resize(bufSize);
|
2014-12-27 12:45:21 -05:00
|
|
|
}
|
2014-12-26 20:58:42 -05:00
|
|
|
|
2015-03-22 20:47:07 -04:00
|
|
|
double freq = (2.0 * M_PI) * ((double) 38000) / ((double) inp->sampleRate);
|
2014-12-26 20:58:42 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (stereoShiftFrequency != freq) {
|
|
|
|
nco_crcf_set_frequency(stereoShifter, freq);
|
|
|
|
stereoShiftFrequency = freq;
|
2014-12-27 12:23:09 -05:00
|
|
|
}
|
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
for (int i = 0; i < bufSize; i++) {
|
2015-01-03 17:07:39 -05:00
|
|
|
firhilbf_r2c_execute(firStereoR2C, demodOutputData[i], &x);
|
|
|
|
nco_crcf_mix_down(stereoShifter, x, &y);
|
|
|
|
nco_crcf_step(stereoShifter);
|
|
|
|
firhilbf_c2r_execute(firStereoC2R, y, &demodStereoData[i]);
|
2014-12-26 20:58:42 -05:00
|
|
|
}
|
2014-12-24 03:03:34 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (audio_out_size != resampledStereoData.size()) {
|
|
|
|
if (resampledStereoData.capacity() < audio_out_size) {
|
|
|
|
resampledStereoData.reserve(audio_out_size);
|
2014-12-27 12:45:21 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
resampledStereoData.resize(audio_out_size);
|
2014-12-24 03:03:34 -05:00
|
|
|
}
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-01-04 13:20:31 -05:00
|
|
|
msresamp_rrrf_execute(stereoResampler, &demodStereoData[0], bufSize, &resampledStereoData[0], &numAudioWritten);
|
2014-12-26 20:58:42 -05:00
|
|
|
}
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (currentSignalLevel > signalLevel) {
|
|
|
|
signalLevel = signalLevel + (currentSignalLevel - signalLevel) * 0.5;
|
2014-12-31 19:45:01 -05:00
|
|
|
} else {
|
2015-01-03 17:07:39 -05:00
|
|
|
signalLevel = signalLevel + (currentSignalLevel - signalLevel) * 0.05;
|
2014-12-31 19:45:01 -05:00
|
|
|
}
|
|
|
|
|
2014-12-27 15:04:43 -05:00
|
|
|
AudioThreadInput *ati = NULL;
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (audioOutputQueue != NULL) {
|
|
|
|
if (!squelchEnabled || (signalLevel >= squelchLevel)) {
|
2014-12-24 01:28:33 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
for (outputBuffersI = outputBuffers.begin(); outputBuffersI != outputBuffers.end(); outputBuffersI++) {
|
|
|
|
if ((*outputBuffersI)->getRefCount() <= 0) {
|
|
|
|
ati = (*outputBuffersI);
|
2014-12-24 01:28:33 -05:00
|
|
|
break;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
if (ati == NULL) {
|
|
|
|
ati = new AudioThreadInput;
|
2015-01-03 17:07:39 -05:00
|
|
|
outputBuffers.push_back(ati);
|
2014-12-24 01:28:33 -05:00
|
|
|
}
|
|
|
|
|
2015-03-15 21:02:26 -04:00
|
|
|
ati->sampleRate = audioSampleRate;
|
2014-12-24 01:28:33 -05:00
|
|
|
ati->setRefCount(1);
|
2014-12-26 20:58:42 -05:00
|
|
|
|
|
|
|
if (stereo) {
|
|
|
|
ati->channels = 2;
|
2015-01-03 17:07:39 -05:00
|
|
|
if (ati->data.capacity() < (numAudioWritten * 2)) {
|
|
|
|
ati->data.reserve(numAudioWritten * 2);
|
2014-12-27 12:45:21 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
ati->data.resize(numAudioWritten * 2);
|
|
|
|
for (int i = 0; i < numAudioWritten; i++) {
|
2014-12-27 15:04:43 -05:00
|
|
|
float l, r;
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
firfilt_rrrf_push(firStereoLeft, (resampledOutputData[i] - (resampledStereoData[i])));
|
|
|
|
firfilt_rrrf_execute(firStereoLeft, &l);
|
2014-12-27 15:04:43 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
firfilt_rrrf_push(firStereoRight, (resampledOutputData[i] + (resampledStereoData[i])));
|
|
|
|
firfilt_rrrf_execute(firStereoRight, &r);
|
2014-12-27 15:04:43 -05:00
|
|
|
|
|
|
|
ati->data[i * 2] = l;
|
|
|
|
ati->data[i * 2 + 1] = r;
|
2014-12-26 20:58:42 -05:00
|
|
|
}
|
|
|
|
} else {
|
|
|
|
ati->channels = 1;
|
2015-01-03 17:07:39 -05:00
|
|
|
ati->data.assign(resampledOutputData.begin(), resampledOutputData.begin() + numAudioWritten);
|
2014-12-26 20:58:42 -05:00
|
|
|
}
|
2014-12-24 01:28:33 -05:00
|
|
|
|
2015-02-08 18:59:08 -05:00
|
|
|
std::vector<float>::iterator data_i;
|
|
|
|
ati->peak = 0;
|
|
|
|
for (data_i = ati->data.begin(); data_i != ati->data.end(); data_i++) {
|
2015-02-12 02:14:22 -05:00
|
|
|
float p = fabs(*data_i);
|
|
|
|
if (p > ati->peak) {
|
2015-02-08 18:59:08 -05:00
|
|
|
ati->peak = p;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
audioOutputQueue->push(ati);
|
2014-12-21 16:08:32 -05:00
|
|
|
}
|
2014-12-16 21:30:03 -05:00
|
|
|
}
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (ati && audioVisOutputQueue != NULL && audioVisOutputQueue->empty()) {
|
2014-12-21 17:37:41 -05:00
|
|
|
|
2014-12-21 16:08:32 -05:00
|
|
|
int num_vis = DEMOD_VIS_SIZE;
|
2014-12-26 22:20:50 -05:00
|
|
|
if (stereo) {
|
2014-12-27 15:04:43 -05:00
|
|
|
ati_vis->channels = 2;
|
|
|
|
int stereoSize = ati->data.size();
|
2014-12-26 22:20:50 -05:00
|
|
|
if (stereoSize > DEMOD_VIS_SIZE) {
|
|
|
|
stereoSize = DEMOD_VIS_SIZE;
|
|
|
|
}
|
2015-01-22 21:32:32 -05:00
|
|
|
|
2014-12-26 22:20:50 -05:00
|
|
|
ati_vis->data.resize(stereoSize);
|
|
|
|
|
|
|
|
for (int i = 0; i < stereoSize / 2; i++) {
|
2015-01-01 03:48:32 -05:00
|
|
|
ati_vis->data[i] = ati->data[i * 2];
|
|
|
|
ati_vis->data[i + stereoSize / 2] = ati->data[i * 2 + 1];
|
2014-12-21 16:08:32 -05:00
|
|
|
}
|
2014-12-18 21:39:32 -05:00
|
|
|
} else {
|
2014-12-27 15:04:43 -05:00
|
|
|
ati_vis->channels = 1;
|
2015-01-04 13:20:31 -05:00
|
|
|
if (numAudioWritten > bufSize) {
|
2014-12-26 22:20:50 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (num_vis > numAudioWritten) {
|
|
|
|
num_vis = numAudioWritten;
|
2014-12-26 22:20:50 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
ati_vis->data.assign(resampledOutputData.begin(), resampledOutputData.begin() + num_vis);
|
2014-12-26 22:20:50 -05:00
|
|
|
} else {
|
2015-01-04 13:20:31 -05:00
|
|
|
if (num_vis > bufSize) {
|
|
|
|
num_vis = bufSize;
|
2014-12-26 22:20:50 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
ati_vis->data.assign(demodOutputData.begin(), demodOutputData.begin() + num_vis);
|
2014-12-21 16:08:32 -05:00
|
|
|
}
|
2014-12-26 22:20:50 -05:00
|
|
|
|
|
|
|
// std::cout << "Signal: " << agc_crcf_get_signal_level(agc) << " -- " << agc_crcf_get_rssi(agc) << "dB " << std::endl;
|
2014-12-21 16:08:32 -05:00
|
|
|
}
|
2014-12-21 17:37:41 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
audioVisOutputQueue->push(ati_vis);
|
2014-12-21 16:08:32 -05:00
|
|
|
}
|
|
|
|
if (!threadQueueControl->empty()) {
|
2015-01-01 18:08:54 -05:00
|
|
|
int newDemodType = DEMOD_TYPE_NULL;
|
|
|
|
|
2014-12-21 16:08:32 -05:00
|
|
|
while (!threadQueueControl->empty()) {
|
|
|
|
DemodulatorThreadControlCommand command;
|
|
|
|
threadQueueControl->pop(command);
|
|
|
|
|
|
|
|
switch (command.cmd) {
|
2015-01-03 17:07:39 -05:00
|
|
|
case DemodulatorThreadControlCommand::DEMOD_THREAD_CMD_CTL_SQUELCH_ON:
|
|
|
|
squelchEnabled = true;
|
2014-12-21 16:08:32 -05:00
|
|
|
break;
|
|
|
|
case DemodulatorThreadControlCommand::DEMOD_THREAD_CMD_CTL_SQUELCH_OFF:
|
2015-01-03 17:07:39 -05:00
|
|
|
squelchEnabled = false;
|
2014-12-21 16:08:32 -05:00
|
|
|
break;
|
2015-01-01 18:08:54 -05:00
|
|
|
case DemodulatorThreadControlCommand::DEMOD_THREAD_CMD_CTL_TYPE:
|
|
|
|
newDemodType = command.demodType;
|
|
|
|
break;
|
2014-12-21 16:08:32 -05:00
|
|
|
default:
|
|
|
|
break;
|
|
|
|
}
|
2014-12-18 21:39:32 -05:00
|
|
|
}
|
2015-01-01 18:08:54 -05:00
|
|
|
|
|
|
|
if (newDemodType != DEMOD_TYPE_NULL) {
|
|
|
|
switch (newDemodType) {
|
|
|
|
case DEMOD_TYPE_FM:
|
2015-01-07 21:25:35 -05:00
|
|
|
freqdem_reset(demodFM);
|
2015-01-01 18:08:54 -05:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_LSB:
|
2015-02-12 22:46:45 -05:00
|
|
|
demodAM = demodAM_LSB;
|
2015-01-07 21:25:35 -05:00
|
|
|
ampmodem_reset(demodAM);
|
2015-01-01 18:08:54 -05:00
|
|
|
break;
|
|
|
|
case DEMOD_TYPE_USB:
|
2015-02-12 22:46:45 -05:00
|
|
|
demodAM = demodAM_USB;
|
2015-01-07 21:25:35 -05:00
|
|
|
ampmodem_reset(demodAM);
|
2015-01-01 18:08:54 -05:00
|
|
|
break;
|
2015-01-07 20:23:15 -05:00
|
|
|
case DEMOD_TYPE_DSB:
|
|
|
|
demodAM = demodAM_DSB;
|
2015-01-07 21:25:35 -05:00
|
|
|
ampmodem_reset(demodAM);
|
2015-01-07 20:23:15 -05:00
|
|
|
break;
|
2015-01-01 18:08:54 -05:00
|
|
|
case DEMOD_TYPE_AM:
|
2015-01-03 17:07:39 -05:00
|
|
|
demodAM = demodAM_DSB_CSP;
|
2015-01-07 21:25:35 -05:00
|
|
|
ampmodem_reset(demodAM);
|
2015-01-01 18:08:54 -05:00
|
|
|
break;
|
2015-06-09 11:00:05 -04:00
|
|
|
// case DEMOD_TYPE_QPSK:
|
|
|
|
// std::cout << "reset modem qpsk" << std::endl;
|
|
|
|
// modem_reset(demodQPSK);
|
|
|
|
// break;
|
2015-01-01 18:08:54 -05:00
|
|
|
}
|
|
|
|
demodulatorType = newDemodType;
|
|
|
|
}
|
2014-12-16 21:30:03 -05:00
|
|
|
}
|
2014-12-21 16:08:32 -05:00
|
|
|
|
2014-12-24 01:28:33 -05:00
|
|
|
inp->decRefCount();
|
2014-12-16 21:30:03 -05:00
|
|
|
}
|
2014-12-16 18:27:32 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
if (audioResampler != NULL) {
|
|
|
|
msresamp_rrrf_destroy(audioResampler);
|
2014-12-16 21:30:03 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
if (stereoResampler != NULL) {
|
|
|
|
msresamp_rrrf_destroy(stereoResampler);
|
2014-12-26 21:55:13 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
if (firStereoLeft != NULL) {
|
|
|
|
firfilt_rrrf_destroy(firStereoLeft);
|
2014-12-27 15:04:43 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
if (firStereoRight != NULL) {
|
|
|
|
firfilt_rrrf_destroy(firStereoRight);
|
2014-12-27 15:04:43 -05:00
|
|
|
}
|
2014-12-16 20:33:44 -05:00
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
agc_crcf_destroy(iqAutoGain);
|
|
|
|
firhilbf_destroy(firStereoR2C);
|
|
|
|
firhilbf_destroy(firStereoC2R);
|
|
|
|
nco_crcf_destroy(stereoShifter);
|
2015-02-19 23:51:24 -05:00
|
|
|
resamp2_cccf_destroy(ssbFilt);
|
2015-01-03 17:07:39 -05:00
|
|
|
|
|
|
|
while (!outputBuffers.empty()) {
|
|
|
|
AudioThreadInput *audioDataDel = outputBuffers.front();
|
|
|
|
outputBuffers.pop_front();
|
2014-12-24 01:28:33 -05:00
|
|
|
delete audioDataDel;
|
|
|
|
}
|
|
|
|
|
2015-01-22 23:41:33 -05:00
|
|
|
if (audioVisOutputQueue && !audioVisOutputQueue->empty()) {
|
|
|
|
AudioThreadInput *dummy_vis;
|
|
|
|
audioVisOutputQueue->pop(dummy_vis);
|
|
|
|
}
|
2015-01-22 21:32:32 -05:00
|
|
|
delete ati_vis;
|
|
|
|
|
2014-12-16 21:30:03 -05:00
|
|
|
DemodulatorThreadCommand tCmd(DemodulatorThreadCommand::DEMOD_THREAD_CMD_DEMOD_TERMINATED);
|
|
|
|
tCmd.context = this;
|
|
|
|
threadQueueNotify->push(tCmd);
|
2015-01-11 20:26:51 -05:00
|
|
|
std::cout << "Demodulator thread done." << std::endl;
|
2015-01-20 19:13:49 -05:00
|
|
|
|
|
|
|
#ifdef __APPLE__
|
|
|
|
return this;
|
|
|
|
#endif
|
2014-12-16 18:27:32 -05:00
|
|
|
}
|
|
|
|
|
2015-01-03 17:07:39 -05:00
|
|
|
void DemodulatorThread::setVisualOutputQueue(DemodulatorThreadOutputQueue *tQueue) {
|
|
|
|
audioVisOutputQueue = tQueue;
|
|
|
|
}
|
|
|
|
|
|
|
|
void DemodulatorThread::setAudioOutputQueue(AudioThreadInputQueue *tQueue) {
|
|
|
|
audioOutputQueue = tQueue;
|
|
|
|
}
|
|
|
|
|
2014-12-16 18:27:32 -05:00
|
|
|
void DemodulatorThread::terminate() {
|
2014-12-16 21:30:03 -05:00
|
|
|
terminated = true;
|
2014-12-23 01:12:14 -05:00
|
|
|
DemodulatorThreadPostIQData *inp = new DemodulatorThreadPostIQData; // push dummy to nudge queue
|
2015-01-03 17:07:39 -05:00
|
|
|
iqInputQueue->push(inp);
|
2014-12-16 18:27:32 -05:00
|
|
|
}
|
2014-12-26 20:58:42 -05:00
|
|
|
|
|
|
|
void DemodulatorThread::setStereo(bool state) {
|
|
|
|
stereo = state;
|
2014-12-26 22:20:50 -05:00
|
|
|
std::cout << "Stereo " << (state ? "Enabled" : "Disabled") << std::endl;
|
2014-12-26 20:58:42 -05:00
|
|
|
}
|
|
|
|
|
|
|
|
bool DemodulatorThread::isStereo() {
|
|
|
|
return stereo;
|
|
|
|
}
|
2014-12-31 19:45:01 -05:00
|
|
|
|
|
|
|
float DemodulatorThread::getSignalLevel() {
|
2015-01-03 17:07:39 -05:00
|
|
|
return signalLevel;
|
2014-12-31 19:45:01 -05:00
|
|
|
}
|
|
|
|
|
|
|
|
void DemodulatorThread::setSquelchLevel(float signal_level_in) {
|
2015-01-03 17:07:39 -05:00
|
|
|
if (!squelchEnabled) {
|
|
|
|
squelchEnabled = true;
|
2014-12-31 19:45:01 -05:00
|
|
|
}
|
2015-01-03 17:07:39 -05:00
|
|
|
squelchLevel = signal_level_in;
|
2014-12-31 19:45:01 -05:00
|
|
|
}
|
|
|
|
|
|
|
|
float DemodulatorThread::getSquelchLevel() {
|
2015-01-03 17:07:39 -05:00
|
|
|
return squelchLevel;
|
2014-12-31 19:45:01 -05:00
|
|
|
}
|
2015-01-01 03:48:32 -05:00
|
|
|
|
2015-01-01 18:08:54 -05:00
|
|
|
void DemodulatorThread::setDemodulatorType(int demod_type_in) {
|
2015-01-01 03:48:32 -05:00
|
|
|
demodulatorType = demod_type_in;
|
|
|
|
}
|
|
|
|
|
2015-01-01 18:08:54 -05:00
|
|
|
int DemodulatorThread::getDemodulatorType() {
|
2015-01-01 03:48:32 -05:00
|
|
|
return demodulatorType;
|
|
|
|
}
|
2015-01-01 18:08:54 -05:00
|
|
|
|
2015-06-05 03:51:46 -04:00
|
|
|
void DemodulatorThread::setDemodulatorLock(bool demod_lock_in) {
|
2015-06-09 11:00:05 -04:00
|
|
|
demod_lock_in ? currentDemodLock = true : currentDemodLock = false;
|
2015-06-05 03:51:46 -04:00
|
|
|
}
|
|
|
|
|
|
|
|
int DemodulatorThread::getDemodulatorLock() {
|
|
|
|
return currentDemodLock;
|
|
|
|
}
|
|
|
|
|
2015-06-10 04:50:14 -04:00
|
|
|
void DemodulatorThread::setDemodulatorCons(int demod_cons_in) {
|
|
|
|
std::cout << "Updating constellations" << std::endl;
|
|
|
|
demodulatorCons = demod_cons_in;
|
|
|
|
}
|
|
|
|
|
|
|
|
int DemodulatorThread::getDemodulatorCons() {
|
|
|
|
return demodulatorCons;
|
|
|
|
}
|
|
|
|
|
2015-06-09 11:00:05 -04:00
|
|
|
void DemodulatorThread::updateDemodulatorLock(modem demod, float sensitivity) {
|
|
|
|
modem_get_demodulator_evm(demod) <= sensitivity ? setDemodulatorLock(true) : setDemodulatorLock(false);
|
2015-06-05 03:51:46 -04:00
|
|
|
}
|
|
|
|
|