%   this is the starting point of this function

block_sim = 16;          %   blocks that the program to simulate    

mapping_index = 4;	     %  for modu index:1/2/4/6 bits per symbol

A_tag = diag(gen_tag(m_transmitter, N_point_FFT));
save A_tag A_tag;
load A_tag;

while error_temp < 1000 
    %   generate transmitter symbols
    trans_symbols = transmitter(mapping_index); %   16QAM mapping
    %   get the variation of trans symbols
    signal_level = mean(abs(trans_symbols).^2);
    %   get dB level
    signal_level_dB = 10*log10(signal_level);
    
    trans_symbols_tag = trans_symbols * A_tag;
    signal_level_tag = mean(abs(trans_symbols).^2);
    signal_level_tag_dB = 10*log10(signal_level);
    
    load Channel_H;

   %    function out_bits_H = receiver(in_bits_H, channel_H,index);
    receiver(trans_symbols, Channel_H, index);
end;