www.pudn.com > mentos_fs.rar > PrimePe_iw.m

%计算经过卷积编码和交织之后的误码率，检测内编码和交织的抗早性能,加高斯噪声及突发干扰 
function [sig_err_rate2]=PrimePe_iw(snr_in_dB) 
% snr_in_dB=0; 
%噪声功率 
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*267));      %内编码器输入序列，一级189bit,二级78bit 
 
in_en_sig=zeros(1,N*456);         %初始化卷积编码器输出矩阵 
in_de_sig=zeros(1,N*267);         %初始化卷积译码器输出矩阵 
 
inter=zeros(1,N*456+456);         %初始化交织器输出矩阵 
counter=zeros(1,N*456+456);       %初始化去交织输出矩阵 
 
for n=0:N-1   %对每桢进行编码 
  %内编码 
  in_en_sig(n*456+1:n*456+456)=T_FS_IN_EN(signal(n*267+1:n*267+267)); 
end 
 
%交织 
inter=INTER_WEAVE_8(in_en_sig); 
 
 %加高斯燥声 
 inter=2*inter-1; 
 inter=sign(inter+sgma*randn(size(inter))+eps); 
 inter=(inter+1)/2; 
 %加突发噪声 
 Noise=BurstNoise(N_Interval,N_Length,length(inter)); 
 inter=mod(inter+Noise,2); 
 
% %加高斯燥声 
% in_en_sig=2*in_en_sig-1; 
% in_en_sig=sign(in_en_sig+sgma*randn(size(in_en_sig))+eps); 
% in_en_sig=(in_en_sig+1)/2; 
    
%去交织 
counter=COUNTER_WEAVE_8(inter); 
 
for n=0:N-1  
  %译码 
  %内译码 
in_de_sig(n*267+1:n*267+267)=T_FS_IN_DE(counter(n*456+1:n*456+456)); 
% in_de_sig(n*267+1:n*267+267)=T_FS_IN_DE(in_en_sig(n*456+1:n*456+456)); 
end 
%验证 
sig_err_rate2=(length(find(signal~=in_de_sig)))/length(signal)   %误码率