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

%计算不加任何编码的原始误码率，加高斯噪声及突发干扰 
function [sig_err_rate1]=PrimePe(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)); 
%加高斯燥声 
u=2*signal-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); 
%验证 
sig_err_rate1=(length(find(signal~=u)))/length(signal)   %误码率