c10.2_mcbpskber.m

使用pMatlab改写BPSK和QPSK 的Monte Carlo 仿真程序。在多核PC上实现MC仿真速度翻倍（附原程序）

clear all
EbNodB=0:8;                 %vector of Eb/No(dB) values
z=10.^(EbNodB/10);          %convert to linear scale
delay=54;                    %enter delay value (samples)
BER=zeros(1,length(z));     %initialize BER vector
Errors=zeros(1,length(z));  %initialize Errors vector
BER_T=qfunc(sqrt(2*z));         %theor. (AWGN) BER vector
N=round(1000./BER_T);         %100 errors for ideal (zero ISI) sys
FilterSwitch=1
for k=1:length(z)
    N(k)=max(1000,N(k));    %ensure at least one block processed
    [BER(k),Errors(k)]=c10_MCBPSKrun(N(k),z(k),delay,FilterSwitch)
end;
disp(Errors);
disp(BER);
disp(BER_T);
semilogy(EbNodB,BER,'o',EbNodB,BER_T,'-r')
axis([min(EbNodB),max(EbNodB),1e-4,1]);
xlabel('E_b/N_o-dB');ylabel('Bit Error Rate');grid
legend('System under Studay','AWGN Ref',0)
%End