Add DVBS2 QPSK implementation from G4EWJ Brian

v2beta
F5OEO 2018-09-03 19:41:33 +00:00
rodzic 66f1178486
commit f426481e8f
3 zmienionych plików z 2874 dodań i 66 usunięć

Wyświetl plik

@ -38,10 +38,11 @@ LDFLAGS_Pissb = librpitx/src/librpitx.a -lm -lrt -lpthread -lsndfile -lliquid
../freedv : freedv/freedv.cpp
$(CCP) $(CFLAGS) -o ../freedv freedv/freedv.cpp $(LDFLAGS)
../dvbrf : dvb/dvbrf.cpp dvb/dvbsenco8.s dvb/fec100.c librpitx/src/librpitx.a
../dvbrf : dvb/dvbrf.cpp dvb/dvbsenco8.s dvb/fec100.c dvb/dvbs2arm_1v30.s librpitx/src/librpitx.a
$(CC) $(CFLAGS) -c -o dvb/dvbsenco8.o dvb/dvbsenco8.s
$(CC) $(CFLAGS) -c -o dvb/dvbs2arm_1v30.o dvb/dvbs2arm_1v30.s
$(CC) $(CFLAGS) -c -o dvb/fec100.o dvb/fec100.c
$(CCP) $(CFLAGS) -o ../dvbrf dvb/dvbrf.cpp dvb/dvbsenco8.o dvb/fec100.o $(LDFLAGS)
$(CCP) $(CFLAGS) -o ../dvbrf dvb/dvbrf.cpp dvb/dvbsenco8.o dvb/fec100.o dvb/dvbs2arm_1v30.o $(LDFLAGS)
CFLAGS_Pifm = -Wall -g -O2 -Wno-unused-variable
LDFLAGS_Pifm = librpitx/src/librpitx.a -lm -lrt -lpthread -lsndfile

Wyświetl plik

@ -3,18 +3,22 @@
#include "stdio.h"
#include <cstring>
#include <signal.h>
// DVBS ENCODER
extern "C"
{
#include "fec100.h"
}
extern "C" void dvbsenco_init (void) ;
extern "C" uchar* dvbsenco (uchar*) ;
extern "C" void energy (uchar* input,uchar *output) ;
extern "C" void reed (uchar *input188) ;
extern "C" uchar* interleave (uchar* packetin) ;
// DVBS2 ENCODER
extern "C" uint32_t _dvbs2arm_control (uint32_t command, uint32_t param1) ;
extern "C" uchar* _dvbs2arm_process_packet (uchar *packet188) ;
bool running=true;
#define PROGRAM_VERSION "0.1"
@ -35,7 +39,7 @@ Usage:\ndvbrf [-i File Input][-s Samplerate][-l] [-f Frequency] [-h Harmonic num
-c Fec : {1/2,3/4,5/6,7/8} \n\
-f float central frequency Hz(50 kHz to 1500 MHz),\n\
-l loop mode for file input\n\
-h Use harmonic number n\n\
-n Use harmonic number n\n\
-? help (this help).\n\
\n",\
PROGRAM_VERSION);
@ -67,9 +71,11 @@ int main(int argc, char* argv[])
int InputType=typeiq_i16;
int Decimation=1;
int FEC=-1;
enum {DVBS,DVBS2};
int DVB_Mode=DVBS;//
while(1)
{
a = getopt(argc, argv, "i:f:s:h:lt:c:");
a = getopt(argc, argv, "i:f:s:n:lt:c:m:");
if(a == -1)
{
@ -106,7 +112,11 @@ int main(int argc, char* argv[])
if(strcmp("carrier",optarg)==0) {fprintf(stderr,"Carrier mode\n");FEC=0;}//CARRIER MODE
break;
case 'h': // help
case 'm': // FEC
if(strcmp("dvbs",optarg)==0) DVB_Mode=DVBS;
if(strcmp("dvbs2",optarg)==0) DVB_Mode=DVBS2;
break;
case 'n': // help
Harmonic=atoi(optarg);
break;
case 'l': // loop mode
@ -155,67 +165,134 @@ int main(int argc, char* argv[])
exit(0);
}
int FifoSize=204*100*4;
int NumberofPhase=4;
phasedmasync dvbsmodul(SetFrequency/Harmonic,SampleRate,NumberofPhase,14,FifoSize);
padgpio pad;
pad.setlevel(7);
#define BURST_PACKET 10
#define BURST_SIZE 204*BURST_PACKET
int PhaseBuffer[BURST_SIZE*4*2]; //4 dibit by byte
static uint32_t BuffAligned[256*2];
static uchar *TsBuffer=(uchar *)BuffAligned ;
static unsigned char BuffIQ[204*2]; //QPSK 2bit per sample : 4 dibits by byte
static uchar *pRS204=NULL;
dvbsenco_init() ;
viterbi_init(FEC);
int MapDVBS[4]={0,3,1,2};//0,-pi/2,pi/2,pi
while(running)
if(DVB_Mode==DVBS)
{
int nbread=fread(TsBuffer,sizeof(uchar),BURST_PACKET*188,tsfile);
if(nbread>0)
{
if(TsBuffer[0]!=0x47) fprintf(stderr,"SYNC ERROR \n");
for(int p=0;p<BURST_PACKET;p++)
{
pRS204 = dvbsenco (TsBuffer+p*188) ;
int NbIQOutput=viterbi (pRS204,BuffIQ);
int dibit=0;
for(int k=0;k<NbIQOutput;k++)
{
for(int i=3;i>=0;i--)
{
PhaseBuffer[dibit++]=MapDVBS[(BuffIQ[k]>>(i*2))&0x3];
}
}
dvbsmodul.SetPhaseSamples(PhaseBuffer,NbIQOutput*4);
}
}
else
{
printf("End of file\n");
if(loop_mode_flag)
fseek ( tsfile , 0 , SEEK_SET );
else
running=false;
}
}
int FifoSize=204*100*4;
int NumberofPhase=4;
phasedmasync dvbsmodul(SetFrequency/Harmonic,SampleRate,NumberofPhase,14,FifoSize);
padgpio pad;
pad.setlevel(7);
dvbsmodul.stop();
#define BURST_PACKET 10
#define BURST_SIZE 204*BURST_PACKET
int PhaseBuffer[BURST_SIZE*4*2]; //4 dibit by byte
static uint32_t BuffAligned[256*2];
static uchar *TsBuffer=(uchar *)BuffAligned ;
static unsigned char BuffIQ[204*2]; //QPSK 2bit per sample : 4 dibits by byte
static uchar *pRS204=NULL;
dvbsenco_init() ;
viterbi_init(FEC);
int MapDVBS[4]={0,3,1,2};//0,-pi/2,pi/2,pi
while(running)
{
int nbread=fread(TsBuffer,sizeof(uchar),BURST_PACKET*188,tsfile);
if(nbread>0)
{
if(TsBuffer[0]!=0x47) fprintf(stderr,"SYNC ERROR \n");
for(int p=0;p<BURST_PACKET;p++)
{
pRS204 = dvbsenco (TsBuffer+p*188) ;
int NbIQOutput=viterbi (pRS204,BuffIQ);
int dibit=0;
for(int k=0;k<NbIQOutput;k++)
{
for(int i=3;i>=0;i--)
{
PhaseBuffer[dibit++]=MapDVBS[(BuffIQ[k]>>(i*2))&0x3];
}
}
dvbsmodul.SetPhaseSamples(PhaseBuffer,NbIQOutput*4);
}
}
else
{
printf("End of file\n");
if(loop_mode_flag)
fseek ( tsfile , 0 , SEEK_SET );
else
running=false;
}
}
dvbsmodul.stop();
}
else
{
int FifoSize=546*30*4;
int NumberofPhase=4;
phasedmasync dvbs2modul(SetFrequency/Harmonic,SampleRate,NumberofPhase,14,FifoSize);
padgpio pad;
pad.setlevel(7);
uint32_t BuffAligned[8832];
uchar *IqBuffer=(uchar *)BuffAligned ;
unsigned char TsBuffer[188];
uint32_t Result=-1;
Result=_dvbs2arm_control (1,(uint32_t)IqBuffer) ; // set DVB-S2 mode by passing buffer address
Result=_dvbs2arm_control (2,0x34) ; // set FEC 3/4
Result=_dvbs2arm_control (3,0x35) ; // set rolloff 0.35
Result=_dvbs2arm_control (4,0) ; // Get Efficiency
printf("Net Bitrate should be %f\n",Result*1e-6*SampleRate);
int PhaseBuffer[546*30]; //4 dibit by byte
unsigned char *IQTempFrame=NULL;
int MapDVBS[4]={0,3,1,2};//0,-pi/2,pi/2,pi
while(running)
{
int nbread=fread(TsBuffer,sizeof(uchar),188,tsfile);
if(nbread>0)
{
IQTempFrame =_dvbs2arm_process_packet (TsBuffer);
if(IQTempFrame!=NULL)
{
uint32_t *IQFrame=(uint32_t*)IQTempFrame;
int symbol=0;
for(int i=0;i<546/2;i++)
{
for(int j=31;j>1;j--)
{
PhaseBuffer[symbol++]=MapDVBS[(((IQFrame[2*i]>>j)&1)<<1)|(((IQFrame[2*i+1]>>j)&1))];
//printf("Symbol=%d\n",PhaseBuffer[symbol]);
}
}
//printf("Symbol=%d\n",symbol);
dvbs2modul.SetPhaseSamples(PhaseBuffer,symbol);
}
else
{
//printf("more\n");
}
}
else
{
printf("End of file\n");
if(loop_mode_flag)
fseek ( tsfile , 0 , SEEK_SET );
else
running=false;
}
}
dvbs2modul.stop();
}
}

Plik diff jest za duży Load Diff