test1.m

9.6kbit/s 全速率数据业务信道(TCH/F9.6)仿真

echo off
clear all
close all

snr_in_DB=-4:2:6;
SNR_in_DB=-4:0.1:6;
snr_in_DB2=-4:1:6;
%桢个数
N=10;
%信息序列
signal=round(rand(1,N*260));

for i=1:length(SNR_in_DB),
   SNR=10^(SNR_in_DB(i)/10);
   t_e_r(i)=Qfunct(sqrt(2*SNR)),
end; 

for i=1:length(snr_in_DB),
   s_e_r(i)=PrimePe(snr_in_DB(i),signal),
end;

  

for i=1:length(snr_in_DB2),
   c_e_r(i)=converPe(snr_in_DB2(i),signal),
end;
semilogy(SNR_in_DB,t_e_r,'-',snr_in_DB,s_e_r,'k*-',snr_in_DB2,c_e_r,'go-');
legend('Theory Error rate',' Gaussian channel Error rate',' Juanji Error rate');