primepe.m

GSM TCH全速话音信道仿真（基带）

%分别计算一级二级和系统的误码率，在加高斯噪声和突发干扰的情况下
function [class1_err_rate,class2_err_rate,sig_err_rate]=PrimePe(snr_in_dB)
echo off
close all
% snr_in_dB=0;
%噪声功率
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%桢个数
N=10;
%突发间隔
N_Interval=80;
%突发长度
N_Length=10;
%信息序列
signal=round(rand(1,N*260));

%初始化
ex_en_sig=zeros(1,N*267);
in_en_sig=zeros(1,N*456);
in_de_sig=zeros(1,N*267);
ex_de_sig=zeros(1,N*260);

inter=zeros(1,N*456+456);
counter=zeros(1,N*456+456);

for n=0:N-1
  %编码
  %外编码
  ex_en_sig(n*267+1:n*267+189)=T_FS_EX_EN(signal(n*260+1:n*260+182));
  ex_en_sig(n*267+190:n*267+267)=signal(n*260+183:n*260+260);
  %内编码
  in_en_sig(n*456+1:n*456+456)=T_FS_IN_EN(ex_en_sig(n*267+1:n*267+267));
end

%交织
inter=INTER_WEAVE_8(in_en_sig);

 %加高斯燥声
 inter=2*inter-1;
 inter=sign(inter+sgma*randn(size(inter))+eps);
 inter=(inter+1)/2;
 %加突发噪声
 Noise=BurstNoise(N_Interval,N_Length,length(inter));
 inter=mod(inter+Noise,2);

% % 加高斯燥声
% in_en_sig=2*in_en_sig-1;
% in_en_sig=sign(in_en_sig+sgma*randn(size(in_en_sig))+eps);
% in_en_sig=(in_en_sig+1)/2;
% %加突发噪声
% Noise=BurstNoise(N_Interval,N_Length,length(in_en_sig));
% in_en_sig=mod(in_en_sig+Noise,2);
   
%去交织
counter=COUNTER_WEAVE_8(inter);

for n=0:N-1 
  %译码
  %内译码
  in_de_sig(n*267+1:n*267+267)=T_FS_IN_DE(counter(n*456+1:n*456+456));
%   in_de_sig(n*267+1:n*267+267)=T_FS_IN_DE(in_en_sig(n*456+1:n*456+456));
  %外译码
  ex_de_sig(n*260+1:n*260+182)=T_FS_EX_DE(in_de_sig(n*267+1:n*267+189));
  ex_de_sig(n*260+183:n*260+260)=in_de_sig(n*267+190:n*267+267);
end
%验证
sig_err_rate=(length(find(signal~=ex_de_sig)))/length(signal)  %系统误码率
class1_err_num=0;
class2_err_num=0;
for n=0:(N-1)
   class1_err_num=class1_err_num+length(find(signal([1:182]+n*260)~=ex_de_sig([1:182]+n*260)));
   class2_err_num=class2_err_num+length(find(signal([183:260]+n*260)~=ex_de_sig([183:260]+n*260)));
end
class1_err_rate=class1_err_num/(N*182)  %一级误码率
class2_err_rate=class2_err_num/(N*78)   %二级误码率