qpak.m

来自「基于MATLAB实现信号的调制解调源代码的简单事例.」· M 代码 · 共 18 行

SNRindB1=0:2:10;
SNRindB2=0:0.1:10;
for i=1:length(SNRindB1)
    [pb,ps]=cm_sm32(SNRindB1(i));   %仿真比特和误码率
    smld_bit_err_prb(i)=pb;
    smld_symbol_err_prb(i)=ps;
end;
for i=1:length(SNRindB2)
    SNR=exp(SNRindB2(i)*log(10)/10);      %信噪比
    theo_err_prd(i)=Qfunct(sqrt(2*SNR));  %理论误码率
end;
figure
semilogy(SNRindB1,smld_bit_err_prb,'*');
hold
semilogy(SNRindB1,smld_symbol_err_prb,'o');
semilogy(SNRindB2,theo_err_prd);
xlabel('信噪比S/N')
ylabel('误码率')