BER output problem for Radio over Fiber system/Optical Comm

%*****************************Defining constant parameter****************
%Rb = 1e+6;                  %bit rate
%Tb = 1/Rb;                  %bit duration / period
Num_bits = 1e4;             %number of bits

nsamp = 20;                 %samples per symbols
%Tsamp = Tb/nsamp;           %sampling rate, fsamp = 20MHz.
a = 0.2;                    %channel attenuation, dB/km
L = 20;                     % Length of fiber
Gain = 10;                  %Gain in dB

%*****************************Filter definations************************

snr_dB=0:20;                            %SNR in dB from 0 to 20

SNR=10.^(snr_dB./10);                   %SNR in linear scale

Tx_filter = ones(1,nsamp);              %transmitter filter, in terms of 
                                        % h(t),match filter
Rx_filter = fliplr(Tx_filter/nsamp);    % receiver filter (matched filter),
                                        % flip from left to right to make
                                        % it h(T-t), match filter

Tx_OOK = randint(1,Num_bits);           % randomly generating 1 and 0 in a 
                                        % row of matrix for 1000
                                        % column,basically ON-OFF keying
                                        % with random bits
                                        
                                        
%*****************************Channel input*****************************                                        
Tx_OOKin = 0+Tx_OOK;                    % Tx+with power, 0dBm,ideal.

%*****************************Channel properties************************
                                        
 Txch = Tx_OOKin -(a*L);                %Signal in channel, deducting with 
                                        %channel attenuation=4dB

Txout = 10.^(Txch./10);                 %Changing back the signal to linear
                                        %scale

Tx_signal = rectpulse(Txout,nsamp);     % generate rectpulse for nsamp time
                                        % per symbol
                                        
                                      
%*****************************Receiver**********************************                                        
for i=1:length(snr_dB)                                        
[b]Rx_signal = awgn(Tx_signal,snr_dB(i)+3-10*log10(nsamp),'measured'); [/b]
                                   


%*****************************Match filtering***************************
MF_out = conv(Rx_signal,Rx_filter);          %equivalent to y(t)=x(t)*h(t)

MF_out_downsamp = MF_out(nsamp:nsamp:end);   %sampling and hold



%try without truncation

%*****************************Thresholding*******************************

Rxth = zeros(1,Num_bits);               %Reset receiver before thresholding

Rxth(find(MF_out_downsamp>0.45))=1;     %value above 0.45 equivalent to one
                                        %at the receiver, thresholding.
                                        
%*****************************BER CALC**********************************
ber(i) = biterr(Tx_OOK,Rxth);           %Matlab function for ber
%******************Theoritical BER CALC**********************************
ber_pr(i) = Qfunct(sqrt(SNR(i)));       %theoritical for ber
end

figure;
semilogy(snr_dB,ber,'b');               %Simulation 
hold on;
semilogy(snr_dB,ber_pr,'r');
xlabel('SNR(dB)');
ylabel('BER');