converpe2.m

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

  function[err_rate,err_rate1,err_rate2]=converPe2(N_Length)

%桢个数
N=50;
rate=0.05;
%信息序列
signal=round(rand(1,N*267));
q=signal;%用于不做任何编码的测试

%编码
for n=0:(N-1)
	%外编码
	%ex_en_sig=T_FS_EX_EN(signal([1:260]+n*260));
	%内编码
   in_en_sig=T_FS_IN_EN(signal([1:267]+n*267));
   %交织前合成
   C([1:456]+n*456)=in_en_sig;
end;


%交织
u=INTER_WEAVE_8(C);

%加突发噪声
Noise=BurstNoise1(rate,N_Length,length(u));
u=mod(u+Noise,2);
Noise=BurstNoise1(rate,N_Length,length(C));
C=mod(C+Noise,2);
Noise=BurstNoise1(rate,N_Length,length(q));
q=mod(q+Noise,2);

%反交织
D=COUNTER_WEAVE_8(u);

%解码
for n=0:(N-1)
	%分离桢
   in_de_sig=D([1:456]+n*456);
   %内编码解码
	ex_de_sig=T_FS_IN_DE(in_de_sig);
    in_de_sig1=C([1:456]+n*456);
    ex_de_sig1=T_FS_IN_DE(in_de_sig1);
    
	
   %桢合成
   rev_sig([1:267]+n*267)=ex_de_sig;
   rev_sig1([1:267]+n*267)=ex_de_sig1;
end;

%误比特率计算
    err_rate=(length(find(signal~=rev_sig)))/length(signal)%内编码+交织编码后的误码率
    err_rate1=(length(find(signal~=q)))/length(signal)%未经编码的信道误码率
    err_rate2=(length(find(signal~=rev_sig1)))/length(signal)%内编码后的误码率