bianma.m

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

%有内编码的函数
function[sig_err_rate,class1_err_rate,class2_err_rate] =bianma(snr_in_dB,g,k)
%噪声功率
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%桢个数
N=100;
%信息序列d
 d=round(rand(1,N*267));
 for n=0:(N-1)
 d([186:189]+n*267)=0;
end
for n=0:(N-1)
	%内编码
   in_en_sig=T_FS_IN_EN(d([1:267]+n*267));
   %交织前合成
   C([1:456]+n*456)=in_en_sig;
end;
u=C;
%%%%%%%%%%%%%%%%%%%%%
%经信道加高斯燥声
%噪声功率
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
u=2*u-1;
u=sign(u+sgma*randn(size(u))+eps);
u=(u+1)/2;
for n=0:(N-1)
	%分离桢
   in_de_sig=u([1:456]+n*456);
  %内编码解码
	ex_de_sig=T_FS_IN_DE(in_de_sig);
  %桢合成
   rev_sig([1:267]+n*267)=ex_de_sig;
end
   m= rev_sig;
sig_err_rate=(length(find(d~=m)))/length(d);
class1_err_num=0;
class2_err_num=0;
for n=0:(N-1)
   class1_err_num=class1_err_num+length(find(d([1:189]+n*267)~=rev_sig([1:189]+n*267)));
   class2_err_num=class2_err_num+length(find(d([190:267]+n*267)~=rev_sig([190:267]+n*267)));
end
class1_err_rate=class1_err_num/(N*189);
class2_err_rate=class2_err_num/(N*78);