function[class1_err_rate,class2_err_rate]=tchfsPe(snr_in_dB)
%噪聲功率
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));
%編碼
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);
   %交織前合成
   C([1:456]+n*456)=in_en_sig;
end;
%交織
u=INTER_WEAVE_8(C);
%%%%%%%%%%%%%%%%%%%%%
%經信道加高斯燥聲
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);
%%%%%%%%%%%%%%%%%%%%%
%反交織
D=COUNTER_WEAVE_8(u);
%錯誤位記錄初始化
err_flag=0;
%解碼
for n=0:(N-1)
	%分離楨
   in_de_sig=D([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;
%驗證
if (length(find(signal~=rev_sig))==0)&(err_flag==0)
   display "good!"
else
   display "error!"
   sig_err_rate=(length(find(signal~=rev_sig)))/length(signal)
end
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)