x = [1, 1.3, 1.5, 1.7, 2];
y = [0.0125, 0.00125, 0.000125, 0.0000125, 0.00000009867];
plot(x, log10(y), 'r-');
axis([0, 30, -9, 0]);
grid
hold on
EbNo = 0:0.01:10
EbNoDecimal= 10.^(EbNo./10);%%將信噪比化為非dB表示
H = 0.5*erfc(sqrt(EbNoDecimal));
%%Q函數(shù)定義：1/2*erfc(sqrt(EbNo))并套用BPSK誤碼率公式
% semilogy(EbNo,H,'-b*');
plot(EbNo, log10(H), 'b-');