ch7example22prog1.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

% ch7example22prog1.m
clear;
reg=ones(1,9);     % 寄存器初始状态：全1，寄存器级数为9
coeff=[1,0,0,1,0,1,1,0,0,1]; % 抽头系数a0 a1 ...ar，取决于特征多项式
N=2^length(reg)-1; % 周期
for k=1:N          % 计算一个周期的m序列输出
    a_n=mod(sum(reg.*coeff(1:length(coeff)-1)),2);% 反馈系数
    reg=[reg(2:length(reg)),a_n];                 % 寄存器移位,反馈
    out(k)=reg(1);                                % 寄存器最低位输出
end
out=2*out-1; % 转换为双极性序列
for j=0:N-1
    rho(j+1)=sum(out.*[out(1+j : N),out(1:j)])/N;
end
j=-N+1:N-1;
rho=[fliplr(rho(2:N)),rho];
plot(j, rho);axis([-10 10 -0.1 1.2]);