%計算在高斯噪聲及突發錯誤干擾下未編碼系統誤碼率
function [sig_err_rate1]=PrimePe(snr_in_dB)
echo off
close all
%噪聲功率
E=1;
% snr_in_dB=0;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%楨個數
N=30;
%突發間隔
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)   %系統誤碼率