c10.5_mcqpskber.m

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

clear all
t=cputime;
Eb=22:0.5:26;No=-50                 %Eb(dBm) and No(dBm/Hz)
ChannelAttenuation=70               %Channel attenuation in dB
EbNodB=(Eb-ChannelAttenuation)-No;  %Eb/No in dB
EbNo=10.^(EbNodB./10);              %Eb/No in linear units
BER_T=0.5*erfc(sqrt(EbNo));         %BER (theor)
N=round(150./BER_T);                %Symbols to transmit
BER_MC=zeros(size(Eb));             
for k=1:length(Eb)                  %Main Loop
    BER_MC(k)=c10_MCQPSKrun(N(k),Eb(k),No,ChannelAttenuation,0,0,0,0);
    disp(['Simulation',num2str(k*100/length(Eb)),'%Complete']);
end;
disp(BER_MC);
semilogy(EbNodB,BER_MC,'o',EbNodB,2*BER_T,'-r')
xlabel('Eb/No(dB)');ylabel('Bit Error Rate');
legend('MC BER Estimate','Theoretical BER');grid;
e=cputime-t;
disp(e);
%End