c10_mcqpsksymjitter.m

对于QPSK通信系统的蒙特卡洛仿真

% File: c10 _MCQPSKSymJitter.m
% 
%
SymJitter = 0:0.02:0.2;
Eb = 24; No = -50;              	% Eb (dBm) and No (dBm/Hz)
ChannelAttenuation = 70; 			% channel attenuation in dB
EbNo = 10.^((Eb-ChannelAttenuation-No)/10); 
BER_T = 0.5*erfc(sqrt(EbNo)*ones(size(SymJitter)));
N=round(100./BER_T); 
BER_MC = zeros(size(SymJitter)); 
for k=1:length(SymJitter)
  BER_MC(k) = c10_MCQPSKrun(N(k),Eb,No,ChannelAttenuation,0,SymJitter(k),0,0);
  disp(['Simulation ',num2str(k*100/length(SymJitter)),'% Complete']);
end
semilogy(SymJitter,BER_MC,'o',SymJitter,2*BER_T,'-')
xlabel('Symbol Timing Error Std. Dev. (Symbols)'); 
ylabel('Bit Error Rate'); 
legend('MC BER Estimate','Theoretical BER'); grid;
% End of script file.