tufa_wujiaozhi_pe.m

来自「matlab信道编码程序」· M 代码 · 共 45 行

function[class1_err_rate]=tufa_wujiaozhi_pe(snr_in_dB,in)
%噪声功率
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%桢个数
N=100;
%突发间隔
N_Interval=80;
%突发长度
N_Length=5;
signal=reshape(in',1,267*N);
%编码
for n=0:(N-1)
    %卷积编码
    neibian_sig(1:456)=nei_bianma(in(n+1,:));
   %交织前合成
   C([1:456]+n*456)=neibian_sig(1:456);
end;
u=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);
%%%%%%%%%%%%%%%%%%%%%
%解码
for n=0:(N-1)
	%分离桢
   neiyi_sig=u([1:456]+n*456);
   %内编码解码
	new_sig([1:267]+n*267)=nei_yima(neiyi_sig);
end;

class1_err_num=0;

for n=0:(N-1)
   class1_err_num=class1_err_num+length(find(signal([1:182]+n*267)~=new_sig([1:182]+n*267)));
 end
class1_err_rate=class1_err_num/(N*182);