BERTool_MPSK_Equalizer_T1.m

function [BER, numBits] = BERTool_MPSK_Equalizer_T1(EbNo, maxNumErrs, maxNumBits)

    persistent FullOperatingTime

    % Display Line on the Start of Imitation Modeling
    disp('======================================');
    % Start time
    tStart = clock;
    % Total duration of Imitation Modeling
    % Saving for each trials. To restart need 'clear all' command.
    if isempty(FullOperatingTime)
        FullOperatingTime = 0;
    end
    
    
    %%%%% Initial Information Source %%%%%
    
    % Symbol Rate
    Rs = 1e3;
    
    % Number of Symbols per Frame
    numSymbols = 1e4;
    
    
    %%%%% M-PSK Modulation %%%%%
    
    % Modulation Order
    M = 2;
    % Number of Bits in Symbol
    k = log2(M);
    
    % M-PSK Modulator Object
    MPSKModulator = comm.PSKModulator( ...
        'ModulationOrder', M, ...
        'BitInput', true, ...
        'PhaseOffset', 0, ...
        'SymbolMapping', 'Gray' ...
    );

    % M-PSK Demodulator Object
    MPSKDemodulator = comm.PSKDemodulator( ...
        'ModulationOrder', MPSKModulator.ModulationOrder, ...
        'BitOutput', MPSKModulator.BitInput, ...
        'PhaseOffset', MPSKModulator.PhaseOffset, ...
        'SymbolMapping', MPSKModulator.SymbolMapping, ...
        'DecisionMethod', 'Hard decision' ...
    );


    %%%%% HF Communication Channel with Multipath and Signal Fading %%%%%
    
    % Rayleigh Channel Object
    RayleighChannel = stdchan( ...
        'iturHFMQ', ...
        Rs, ...
        1 ...
    );
    % Delay in Channel Object
    ChanDelay = info(RayleighChannel).ChannelFilterDelay;

    % AWGN Channel Object
    AWGNChannel = comm.AWGNChannel( ...
        'NoiseMethod', 'Signal to noise ratio (Eb/No)', ...
        'EbNo', EbNo, ...
        'SignalPower', 1, ...
        'BitsPerSymbol', k ...
    );


    %%%%% Adaptive Equalizer %%%%%
    
    % Decision Feedback Equalizer Object
    DFEqualizer = comm.DecisionFeedbackEqualizer( ...
        'Algorithm', 'RLS', ...
        'NumForwardTaps', 5, ...
        'NumFeedbackTaps', 3, ...
        'ReferenceTap', 3, ...
        'ForgettingFactor', 0.9, ...
        'InitialInverseCorrelationMatrix', 0.1, ...  
        'Constellation', constellation(MPSKModulator), ...
        'InputDelay', ChanDelay ...
    );

    % Delay in Equalizer Object
    EqualizerDelay = info(DFEqualizer).Latency;

    % Number of Training Symbols per Frame
    numTrainingSymbols = numSymbols/5;
    % Number of Data Symbols per Frame
    numDataSymbols = numSymbols - numTrainingSymbols;  
    
    
    %%%%% Imitation Modeling %%%%%
    
    % Import Java class for BERTool
    import com.mathworks.toolbox.comm.BERTool;
    
    % Total Delay in Channel & Equalizer Objects
    TotalDelay = ChanDelay + EqualizerDelay;
    
    % BER Calculator Object
    BERCalculater = comm.ErrorRate;
    % BER Intermediate Variable
    BERIm = zeros(3,1);
    
    % Generation of Trainig Bits
    TrainingBits = randi([0 1], k*numTrainingSymbols, 1);
    % M-PSK modulation of Trainig Bits
    TrainigSignal = MPSKModulator(TrainingBits);
    
    % Imitation Modeling Loop
    tLoop1 = clock;
    while BERIm(2) < maxNumErrs && BERIm(3) < maxNumBits
        
        % Check of User push Stop
        if BERTool.getSimulationStop
            break;
        end
        
        % >>> Transmitter >>>
        
        % Generation of Data Bits
        DataBitsTx = randi([0 1], k*numDataSymbols, 1);
        % M-PSK modulation of Data Bits
        DataSignal = MPSKModulator(DataBitsTx);
        % Merging the Trainig & Data Signals in Transmitted Signal
        SignalTx = [TrainigSignal; DataSignal];
        
        
        % >>> HF Communication Channel with Multipath and Signal Fading  >>>
        
        % Rayleigh Channel
        SignalChan1 = RayleighChannel(SignalTx);
        % AWGN Channel
        SignalChan2 = AWGNChannel(SignalChan1);
        
        
        % >>> Receiver >>>
        
        % Equalization
        SignalRx1 = DFEqualizer(SignalChan2, TrainigSignal);
        % M-PSK Demodulation
        BitsRx = MPSKDemodulator(SignalRx1(TotalDelay + 1 : end));
        
        % Selection of Data Bits
        DataBitsRx = BitsRx(k*numTrainingSymbols + 1 : end);
        % BER Calculation
        BERIm = BERCalculater( ...
            DataBitsTx(1 : end - k*TotalDelay), ...
            DataBitsRx ...
        );
        
    end
    tLoop2 = clock;
    
    % BER Results
    BER = BERIm(1);
    numBits = BERIm(3);
    disp(['BER = ', num2str(BERIm(1), '%.5g'), ' at Eb/No = ', num2str(EbNo), ' dB']);
    disp(['Number of bits = ', num2str(BERIm(3))]);
    disp(['Number of errors = ', num2str(BERIm(2))]);
    
    
    % Performance of Imitation Modeling
    Performance = BERIm(3) / etime(tLoop2, tLoop1);
    disp(['Performance = ', num2str(Performance), ' bit/sec']);    
    
    % Duration of this Imitation Modeling
    Duration = etime(clock, tStart);
    disp(['Operating time = ', num2str(Duration), ' sec']);
    
    % Total duration of Imitation Modeling
    FullOperatingTime = FullOperatingTime + Duration;
    assignin('base', 'FullOperatingTime', FullOperatingTime);

end