primepe.m

来自「9.6kbit/s 全速率数据业务信道(TCH/F9.6)仿真」· M 代码 · 共 22 行

function[sig_err_rate]=PrimePe(snr_in_dB,signal)
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%桢个数
N=100;
%突发间隔
N_Interval=80;
%突发长度
N_Length=5;
%%%%%%%%%%%%%%%%%%%%%
%经信道加高斯燥声
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_rate=(length(find(signal~=u)))/length(signal);