ch7example23prog1.m

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

% ch7example23prog1.m
clear;
reg=ones(1,6);     % 寄存器初始状态：全1，寄存器级数为6
coeff=[1,0,0,0,0,1,1]; % 抽头系数cr  ...c1 c0，取决于特征多项式
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];                 % 寄存器移位,反馈
    out1(k)=2*reg(1)-1;                           % 寄存器最低位输出,转换为双极性序列
end
reg=ones(1,6);
coeff=[1,1,0,0,1,1,1];% 抽头系数
for k=1:N             % 计算一个周期的m序列输出
    a_n=mod(sum(reg.*coeff(1:length(coeff)-1)),2);% 反馈
    reg=[reg(2:length(reg)),a_n];                 % 寄存器移位,反馈
    out2(k)=2*reg(1)-1;                           % 寄存器最低位输出,转换为双极性序列
end
% 得出两个双极性电平的m序列
for j=0:N-1
    R(j+1)=sum(out1.*[out2(1+j : N),out2(1:j)]);  % 计算相关函数
end
j=-N+1:N-1;% 相关函数自变量
R=[fliplr(R(2:N)),R];% 利用相关函数的偶函数特性计算j为负值的情况
plot(j, R);axis([-N N -20 20]);xlabel('j');ylabel('R(j)');% 作图
max(abs(R)) % 计算相关函数绝对值的最大值