runbpskber.m

bpsk解调

clear all;
format long;
snr=0:10;
for n=1:length(snr);
    peAnaly(n)=erfc(sqrt(10^(snr(n)/10)))/2;
    i=1;
    number=0;
    for i=1:100000;
        number=number+BPSKsimple(snr(n));
    end
    pe(n)=number/1000000;
end
figure(1); plot(snr,pe,'-ro',snr,peAnaly,'-b*');
h=legend('Simulated Pe', 'Analytical Pe');
title('Pe of BPSK with Additive White Gaussian Noise');
axis([0 10 1.0e-006 1]);
set(gca, 'YScale', 'log');
xlabel('Eb/No (dB)');
ylabel('probability of err Pe');
grid on;