ip_9_5.m

二进制通信系统的进行蒙特卡罗仿真程序

echo on
SNRindB1=0:1:12;
SNRindB2=0:0.1:12;
for i=1:length(SNRindB1),
  % 仿真误码率
  smld_err_prb(i)=smldPe54(SNRindB1(i));
end;
for i=1:length(SNRindB2),
    SNR=exp(SNRindB2(i)*log(10)/10);
    % 理论误码率
    theo_err_prb(i)=Qfunct(sqrt(SNR));
end;
% 随后为绘图命令
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);  
pause
A=[0 0.1 1.0 2.0];
for i=1:4,
    sgma=sqrt(A(i));
n0=sgma*randn(1000,1);
n1=sgma*randn(1000,1);
n2=sgma*randn(1000,1);
n3=sgma*randn(1000,1);
x1=1.+n0;
y1=n1;
x2=n2;
y2=1.+n3;
pause
plot(x1,y1,'o',x2,y2,'*')
axis('square')
end;