runbfskber.m

nc mfsk调制

clear all;
close all;
format long;
snr=0:10;
for n=1:length(snr);
    peAnaly(n)=(2-1)/2*exp(-10^(snr(n)/10)/2);
    i=1;
    number=0;
     aa=snr(n)
    for i=1:20;
             BS=randint(200,1,2);% random symbols
             fsk=fskmod(BS, 2,1,4,10);% FSK modulator
             Rxsignals=awgn(fsk, aa, 'measured');% received signals with additive white Gaussian noise
             fskdem=fskdemod(Rxsignals, 2,1,4,10);% demodulated signals
             [num1, ratio]=biterr(BS, fskdem);% count number of errors
             number=number+num1;
    end
    pe(n)=number/20000;
end
figure; plot(snr,pe,'-ro',snr,peAnaly,'-b*');
h=legend('Simulated Pe', 'Analytical Pe');
title('Pe of BFSK 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;