echo on
SNRindB1=0:0.5:15;
SNRindB2=0:0.1:15;
M=16;
k=log2(M);
for i=1:length(SNRindB1),
    smld_err_prb(i)=qammoto(SNRindB1(i));
end;

for i=1:length(SNRindB2),
    SNR=exp(SNRindB2(i)*log(10)/10);
    theo_err_prb(i)=4*Qfunct(sqrt(3*k*SNR/(M-1)));
end;

semilogy(SNRindB1,smld_err_prb,'mo');   % 實(shí)際的信噪比—誤碼率曲線
hold
semilogy(SNRindB2,theo_err_prb,'r-');   % 理論的曲線