primeerr.m

来自「GSM信道编码设计」· M 代码 · 共 21 行

echo off
clear all
close all
%snr_in_DB
snr_in_DB=-4:2:6;
SNR_in_DB=-4:0.1:6;
for i=1:length(snr_in_DB),
   s_e_r(i)=PrimePe(snr_in_DB(i));
end;
for i=1:length(SNR_in_DB),
   SNR=10^(SNR_in_DB(i)/10);
   t_e_r(i)=Qfunct(sqrt(2*SNR));
end;   
semilogy(snr_in_DB,s_e_r,'r*');
hold
semilogy(SNR_in_DB,t_e_r);