L=1000;?
M=4;%每個符號的比特數
N=2^M-1;%編碼后碼字長度
K=N-4;%信息長度
MSG=randint（L1）;%隨機產生L比特信號
TP=gftuple（[-1:N-1]‘M）;%產生加羅華域元素
PG=rspoly（NK）;%產生生成式
[CODEADDED]=rsenco（MSGTPK）;%編碼
NOI=rand（length（CODE）/M1）<.03;%加入3%的噪聲
NOI=（NOI*ones（1M））‘;%產生突發錯誤
NOI=NOI（:）;
CODE_NOI=rem（CODE+NOI2）;%噪聲加入信號
[DECERRCCODEERR_C]=rsdeco（CODE_NOITPK）;%譯碼
MSG=[MSG;zeros（ADDED1）];%調整長度
max（abs（DEC-MSG））;%比較
X=[1:length（NOI）];
Z=[1:M*N:length（NOI）];
Y=zeros（1length（Z））;
Z=[Z;Z];
Y=[Y+min（ERR_C）;Y+max（ERR_C）];
subplot（211）;
plot（XNOI‘yo‘XER