%  2 transmit,1 receive(1bit/s/Hz)
%  QPSK  2bit/s/Hz orthogonal coded
clear all
clc
%SNR upto 20 dBs
EbNo=-5:2;

%N, M: number of transmit and receive antennas
N=2;
M=1

Num=200;% number of frame
K=100;  % numver of symbols in a frame

% initialize count
idx=1;
BIT=2;
%  BIT=2 for QPSK
%  BIT=3 for 8PSK
%  BIT=4 for 16QAM

%  the total bit is K*BIT
for en=1:length(EbNo)
       EbNo(en)
       sigma=0.5/(10^(EbNo(en)/10));  

        % Num -> number of packets
       for packet_count=1:Num
           packet_count
           data=randint(1,K*BIT);%%%%%%%%%%%%%%%%%%%%%%
           temp=EN_QPSK(data);
           for i=1:K/2
               x_time1(:,i)=[temp(1,(2*i-1));temp(1,2*i)];
               x_time2(:,i)=[-conj(temp(1,2*i));conj(temp(1,(2*i-1)))];
           end
           %----------------------------------------------------------
           % Rayleigh
           H=(randn(1,N)+j*randn(1,N))/sqrt(2);  
           %---------------------------------------
           %  Rice channel
           %KK=1;
           %H=sqrt(KK/(1+KK))*ones(M,N)+sqrt(1/(1+KK))*(randn(M,N)+j*randn(M,N))/sqrt(2);
           %---------------------------------------
           % Nakagami-m
           m=2;
           H=nakagamirnd_gamma(m,1,N);
           %----------------------------------------------------------
           %  correlation channels
           Rt=[1 0.1;0.1 1];
           Rt=chol(Rt);
           Rr=[1 0;0 1];
           Rr=chol(Rr);
           H=H*Rt;
          %________________________________________________________________
          %________________________________________________________________
           
           
           
           r_time1=H*x_time1+sqrt(sigma*N)*(randn(1,K/2)+j*randn(1,K/2));
           r_time2=H*x_time2+sqrt(sigma*N)*(randn(1,K/2)+j*randn(1,K/2));
           

           x1_est=conj(H(1,1))*r_time1+H(1,2)*conj(r_time2);
           x2_est=conj(H(1,2))*r_time1-H(1,1)*conj(r_time2);

           for i=1:K/2
               x_est(1,(2*i-1))=x1_est(1,i);
               x_est(1,2*i)=x2_est(1,i);
           end
           [x_est_bit x_est_symbol]=DE_QPSK(x_est);
           
           error_bit(packet_count)=sum(x_est_bit~=data);
           error_symbol(packet_count)=sum(x_est_symbol~=temp);

       end
       BER(idx)=sum(error_bit)/(BIT*Num*K);
       SER(idx)=sum(error_symbol)/(Num*K);
       idx=idx+1
   end
   
% ***************************************************
%figure(1);
semilogy(EbNo,SER,'r-');
grid on;
xlabel('SNR[dB]')
ylabel('SER')

%figure(2)
%semilogy(EbNo,BER,'bo');
%grid on;
%xlabel('SNR[dB]')
%ylabel('BER')