%m序列發(fā)生器及器自相關(guān) mseq.m
clear all;
close all;
g=19;        %G=10011
state=8;     %state=1000
L=1000;

%m序列發(fā)生器
N=15;
mq=mgen(g,state,L);

%求序列自相關(guān)
ms=conv(1-2*mq,1-2*mq(15:-1:1))/N;

figure(1)
subplot(222)
stem(ms(15:end));
axis([0 63 -0.2 1.2]);title('m序列自相關(guān)序列')

%m序列構(gòu)成的信號（矩形脈沖）
N_sample=8;
Tc=1;
dt=Tc/N_sample;
t=0:dt:Tc*L-dt;

gt=ones(1,N_sample);
mt=sigexpand(1-2*mq,N_sample);
mt=conv(mt,gt);
mt1=sigexpand(mq,N_sample);
mt1=conv(mt1,gt);
figure(1)
subplot(221)
plot(t,mt1(1:length(t)));
axis([0 63 -0.2 1.2]);title('m序列矩形形成信號')


st=sigexpand(1-2*mq(1:15),N_sample);
s=conv(st,gt);
st=s(1:length(st));

rt1=conv(mt,st(end:-1:1))/(N*N_sample);

subplot(223)
plot(t,rt1(length(st):length(st)+length(t)-1));
axis([0 63 -0.2 1.2]);title('m序列矩形形成信號的自相關(guān)');xlabel('t');


subplot(224);
[f,mf1]=T2F(t,mt1);
plot(f,10*log10(abs(mf1).^2));
axis([-2 2 0 50]);title('輸入信號功率譜密度（dB/Hz）');