%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%                                                                       %
%%         Name: wimax.m                                                 %
%%                                                                       %
%%     In this file, a menu is defined which is able to call the         %
%%     different functions that carry out the simulations of the         %
%%     system by varying different parameters such as the channel        %
%%     model, the modulation technique, the size of the cyclic           %
%%     prefix and also the channel bandwidth.                            %
%%                                                                       %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
close all;
clear all;

% The variables used in the simulation of the system are defined here
SUI = [];
G = [];
n_mod_type = [];

run = 'y';
while run=='y'
    clc;
    disp('   ********************************************************************************************');
    disp('   *                                                                                          *');
    disp('   *                                   FINAL YEAR PROJECT                                     *');
    disp('   *                                                                                          *');
    disp('   *                                                                                          *');    
    disp('   *              Simulation and Performance Analysis of The 802.16 Physical Layer            *');
    disp('   *                                                                                          *');
    disp('   *                         Study Done by Carlos Batll閟 Ferrer                              *');
    disp('   *                                                                                          *');
    disp('   *       Next, all the possible variants of the simulation shall appear.                    *');
    disp('   *       In every option the user will be asked to define some parameters correctly         *');
    disp('   *       for the simulation.                                                                *');
    disp('   *                                                                                          *');
    disp('   ********************************************************************************************');
    disp('  ');
    disp('      --> The different tests:');
    disp('  ');
    disp('         1) Simulation in which all the modulations are used (BPSK,QPSK,16QAM and 64QAM).'); 
    disp('         2) When we change the size of the "cyclic prefix" (1/4 1/8 1/16 1/32).');
    disp('         3) We realise the simulation WITH and WITHOUT encoding of the bits and study the difference.');
    disp('         4) We carry out the simulation through different SUI channels (1 al 6).');
    disp('         5) We do a simulation with different values of the nominal BW of the system.');
    disp('         6) To exit the program.');
    disp('  ');
    option = input('     Please enter your choice: ');



    switch option
        case 1
            disp('  ');
            G = input('     Please enter the value of G (Cyclic Prefix) [1/4 1/8 1/16 1/32]: ');
            disp('  ');
            SUI = input ('     Please enter which channel you wish to simulate (1 al 6) [AWGN = 0]: ');
            disp('  ');
            disp('     Please enter the nominal BandWidth of the system (BW)');
            BW = input('     Possible Values:28,24,20,15,14,12,11,10,7,6,5.50,5,3.50,3,2.50,1.75,1.5,1.25 [MHz]:');
            disp('  ');
            figur = input('     From which number you wish to start calling the resultant figures:');
            disp('  ');
            samples = input('     Finally enter the number of OFDM symbols to simulate (total bits = 20*symbols):');
            disp('  ');
            disp('     Realizing the simulation.....Please wait a while');
            tic
            TestMods(G,SUI,samples,BW,figur);
        case 2
            disp('  ');
            n_mod_type = input('     Please enter the modulation to use (1->BPSK, 2->QPSK, 3->16QAM