%-----------------------------------------------------------------------
% Compute the bit error rate for binary signaling over Rician fading
% channels via simulation. Assumes ideal synchronization and slow fading.
% For Gaussian performance set KdB = a large number, e.g. 40 dB
% For Rayleigh performance, delete direct path and set KdB = 0 dB 
%-----------------------------------------------------------------------

BlkSize = 2000;
Nblks   = 100;
KdB     = [6];
EbNodB  = [12:2:20];

Output = [EbNodB'];

for kk= 1:length(KdB)
    %--- assume coherent carrier recover ---
    for k = 1:length(EbNodB)
        errors  = 0;
        K       = 10 .^(-KdB(kk)/10);
        sigma   = 10 .^(-EbNodB(k)/20) /sqrt(2);
        
        for n = 1:Nblks
            Rice   = abs( 1 + sqrt(K/2)* (randn(1,BlkSize) + j*randn(1,BlkSize)));
            Signal = Rice   + sigma * randn(1,BlkSize);
            errors = errors + length(find(Signal<0));
        end
        [kk k]
        BER(k) = errors / (Nblks*BlkSize);
    end
    
    Output = [Output BER']
end