ysgspe.m

根据TCH/FS信道的技术标准

 function[class1_err_rate,class2_err_rate,sig_err_rate]=ysgspe(snr_in_dB)
%噪声功率
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%桢个数
N=1000;

%突发间隔
N_Interval=80;
%突发长度
N_Length=10;
%信息序列
signal=round(rand(1,N*260));

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

%验证
rev_sig=u;
sig_err_rate=(length(find(signal~=rev_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)~=rev_sig([1:182]+n*260)));
   class2_err_num=class2_err_num+length(find(signal([183:260]+n*260)~=rev_sig([183:260]+n*260)));
end
class1_err_rate=class1_err_num/(N*182)
class2_err_rate=class2_err_num/(N*78)