clc
clear

num=100000;
snr_db=-100.116;
snr=10.^(snr_db10);
snr_th_db=-5;
snr_th=10^(snr_th_db10);
p_mrc_e=(1-exp(-snr_th.snr)).^2;
p_sc_e=1-exp(-snr_th.snr).(1+snr_th.(snr));
semilogy(snr_db,p_mrc_e);
hold on
semilogy(snr_db,p_sc_e);

p_mrc_sim=zeros(1,length(snr));
p_sc_sim=zeros(1,length(snr));
for i=11length(snr)
    h1=randn(1,num)+1irandn(1,num);
    h2=randn(1,num)+1irandn(1,num);
    for j=11num
        if (abs(h1(j))^2+abs(h2(j)^2))snr(i)2snr_th
            p_mrc_sim(i)=p_mrc_sim(i)+1;
        end
        if ((abs(h1(j))^2)(0.5snr(i))snr_th)&&((abs(h2(j))^2)(0.5snr(i))snr_th)
            p_sc_sim(i)=p_sc_sim(i)+1;
        end
    end
    p_mrc_sim(i)=p_mrc_sim(i)num;
    p_sc_sim(i)=p_sc_sim(i)num;
end
semilogy(snr_db,p_mrc_sim,'o');
semilogy(snr_db,p_sc_sim,'o');
legend('MRC-exact','SC-exact','MRC-sim','SC-sim')