%matlab程序:生成m序列和相关性分析
clc;
N=5;
connections=gfprimfd(N,'all');%生成级数为5时所有本原多项式系数序列矩阵f1=connections(2,:); %取一组本原多项式序列,此系数为45(100101)f2=connections(3,:); %取另一组本原多项式序列,此系数为75(111101)registers1=[0 0 0 0 1];%给定寄存器的初始状态
registers2=[0 0 0 0 1];%取相同的初始状态
L=2^N-1; %周期长度
sum2=0;
sum1=0;
for k=1:L
seq1(k)=registers1(N); %第一组m序列
seq2(k)=registers2(N); %第二组序列
for j=1:N %进行模2加
sum1=sum1+f1(j+1)*registers1(j); %各级寄存器送参与模2加的值
sum1=mod(sum1,2);
sum2=sum2+f2(j+1)*registers2(j); %各级寄存器送参与模2加的值
sum2=mod(sum2,2);
end
for t=N:-1:2 %寄存器移位
registers1(t)=registers1(t-1);
registers2(t)=registers2(t-1);
end
registers1(1)=sum1; registers2(1)=sum2;
sum1=0; sum2=0;
end
disp(f1); %显示反馈系数序列
disp(seq1); %显示第一组m序列
m1=xcorr(seq1,’unbiased’);%计算自相关函数
figure;stem(m1/max(m1));
title('m1的自相关函数'); %画出m1序列的自相关函数
disp(f2);
disp(seq2);
m2=xcorr(seq2,’unbiased’);
figure;stem(m2/max(m2));
title('m2的自相关函数');%画出m2序列的自相关函数
R12= xcorr(seq1,seq2,’unbiased’)
figure;plot(R12/max(R12));
title('m2,m1的互相关函数');%画出m1,m2序列的互相关特性
s=fftshift(abs(fft(seq1,2*L)).^2); s=s/max(s);
figure;plot(s);
如果觉得《matlab中m序列的反馈系数 m序列的生成及其相关性的matlab分析》对你有帮助,请点赞、收藏,并留下你的观点哦!