% 二進制信道容量 2007-6-4
% gamma=E/N
% 最佳檢測的BPSK的差錯概率是： p=Q(sqrt(2r)) ,Q()=(1/2)*erfc(x/sqrt(2))自寫函數 
% C=1-H(p)
% entropy2熵函數自寫函數 
echo on
gamma_db=[-20:0.1:20];
gamma=10.^(gamma_db./10);
p_error=q(sqrt(2.*gamma));
capacity=1.-entropy2(p_error);
pause % 看看錯誤概率 p vs S/N
clf
semilogx(gamma,p_error)
xlabel('SNR/bit')
title('Error probability versus SNR/bit')
ylabel('Error Prob')
pause % 查看信道容量VS信噪比
clf
semilogx(gamma,capacity,'xr')
%plot(gamma,capacity)
%loglog(gamma,capacity)
xlabel('SNR/bit')
title('Channel capacity versus SNR/bit')
ylabel('Channel capacity')