clear;
clc;
%產生一m序列，其特征多相式:1+x^3+x^5
Np=62;%初始化m序列的長度
a（1）=1;
a（2）=0;
a（3）=0;
a（4）=0;
a（5）=0;
for?i=6:Np;%根據m序列特征式計算m序列的其余項
????a（i）=xor（a（i-3）a（i-5））;
end
%產生一加擾序列c（k）;滿足c（k）=a（k）+c（k-3）+c（k-5）
Nq=Np;%初始化加擾序列的長度
c（1）=1;
c（2）=0;
c（3）=0;
c（4）=1;
c（5）=0;
for?n=6:Nq;%計算加擾序列的其余項
????c（n）=xor（c（n-3）c（n-5））;
????c（n）=xor（c（n）a（n））;
end
%產生解擾序列d（k）恢復出a（k）滿足d（k）=c（k）+c（k-3）+c（k-5）
Nr=Np;%初始化解擾序列的長度
d