if ~exist('dvbs2xLDPCParityMatrices.mat','file')
    if ~exist('s2xLDPCParityMatrices.zip', 'file')
        url = 'https://ssd.mathworks.com/supportfiles/spc/satcom/DVB/s2xLDPCParityMatrices.zip';
        websave('s2xLDPCParityMatrices.zip',url);
        unzip('s2xLDPCParityMatrices.zip');
    end
addpath('s2xLDPCParityMatrices');
end

numFrames = 1;

s2WaveGen = dvbs2WaveformGenerator;
s2WaveGen.NumInputStreams = 2;
s2WaveGen.MODCOD = [21 16];
s2WaveGen.DFL   = 47008;
s2WaveGen.ISSYI = true;
s2WaveGen.SamplesPerSymbol = 2;
disp(s2WaveGen)

info(s2WaveGen)

syncBits = [0 1 0 0 0 1 1 1]';       % Sync byte for TS packet is 47 Hex
pktLen = 1496;                       % UP length without sync bits is 1496
data =  cell(1,s2WaveGen.NumInputStreams);
for i = 1:s2WaveGen.NumInputStreams
    numPkts = s2WaveGen.MinNumPackets(i)*numFrames;
    txRawPkts = randi([0 1],pktLen,numPkts);
    ISSY = randi([0 1],16,numPkts);   % ISCRFormat is 'short' by default
                                      % 'short' implies the default length of ISSY as 2 bytes
    txPkts = [repmat(syncBits,1,numPkts);txRawPkts;ISSY];    % ISSY is appended at the end of UP
    data{i} = txPkts(:);
end

txWaveform = [s2WaveGen(data)];

flushFilter(s2WaveGen)

if ~exist('dvbs2xLDPCParityMatrices.mat','file')
    if ~exist('s2xLDPCParityMatrices.zip','file')
        url = 'https://ssd.mathworks.com/supportfiles/spc/satcom/DVB/s2xLDPCParityMatrices.zip';
        websave('s2xLDPCParityMatrices.zip',url);
        unzip('s2xLDPCParityMatrices.zip');
    end
addpath('s2xLDPCParityMatrices');
end

numFrames = 2;

s2xWaveGen = dvbs2xWaveformGenerator();
s2xWaveGen.HasTimeSlicing = true;
s2xWaveGen.NumInputStreams = 2;
s2xWaveGen.PLSDecimalCode = [135 145];   % QPSK 9/20 and 8PSK 25/36
s2xWaveGen.DFL = [18048 44656];
s2xWaveGen.PLScramblingIndex = [0 1];
disp(s2xWaveGen)

info(s2xWaveGen)

syncBits = [0 1 0 0 0 1 1 1]';       % Sync byte for TS packet is 47 Hex
pktLen = 1496;                       % UP length without sync bits is 1496
data =  cell(1, s2xWaveGen.NumInputStreams);
for i = 1:s2xWaveGen.NumInputStreams
    numPkts = s2xWaveGen.MinNumPackets(i)*numFrames;
    txRawPkts = randi([0 1], pktLen, numPkts);
    txPkts = [repmat(syncBits, 1, numPkts); txRawPkts];
    data{i} = txPkts(:);
end

txWaveform = s2xWaveGen(data);

flushFilter(s2xWaveGen)

tmWaveGen = ccsdsTMWaveformGenerator;
tmWaveGen.WaveformSource = "synchronization and channel coding";
tmWaveGen.ChannelCoding = "LDPC";
tmWaveGen.NumBitsInInformationBlock = 1024;
tmWaveGen.Modulation = "QPSK";
tmWaveGen.CodeRate = "1/2";
disp(tmWaveGen)

numTF = 20;

info(tmWaveGen)

bits = randi([0 1], tmWaveGen.NumInputBits*numTF,1);
waveform = tmWaveGen(bits);

flushFilter(tmWaveGen)

chan = etsiRicianChannel;
chan.SampleRate = 2e5;     
chan.KFactor = 10;
chan.MaximumDopplerShift = 20;
chan.NumSinusoids = 58;
disp(chan)

txWaveform = randi([0 1],500,1);
rxWaveform = chan(txWaveform);

info(chan)