gbconv.m

来自「依据GSM 05.03建议」· M 代码 · 共 49 行

function sig_err_rate=gbconv(snr_in_dB)
%噪声功率
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%桢个数
N=50;
%突发间隔
N_Interval=80;
%突发长度
N_Length=6;
%信息序列
signal=round(rand(1,N*260));
%编码
for n=0:(N-1)
	%外编码
	ex_en_sig=T_FS_EX_EN(signal([1:260]+n*260));
	%内编码
   in_en_sig=T_FS_IN_EN(ex_en_sig);
   %交织前合成
   u([1:456]+n*456)=in_en_sig;
end;
%%%%%%%%%%%%%%%%%%%%%
%经信道加高斯燥声
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);
%%%%%%%%%%%%%%%%%%%%%
%错误位记录初始化
err_flag=0;
%解码
for n=0:(N-1)
	%分离桢
   in_de_sig=u([1:456]+n*456);
   %内编码解码
	ex_de_sig=T_FS_IN_DE(in_de_sig);
	%外编码解码
   [newsig,flag]=T_FS_EX_DE(ex_de_sig);
   %错误位记录，如果有错，err_flag会在相应的位置记1
   err_flag=err_flag*10+mod(flag,2);
   %桢合成
   rev_sig([1:260]+n*260)=newsig;
end;
%验证
   sig_err_rate=(length(find(signal~=rev_sig)))/length(signal);