用MATLAB实现m序列的加扰与解扰

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