ch7example22prog1.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 17 行

% 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]);