mseq.m

通信中一些基本实验的程序AM、FM、ASK、FSK、PSK等等

%m序列发生器及器自相关 mseq.m
clear all;
close all;
g=19;        %G=10011
state=8;     %state=1000
L=1000;

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

%求序列自相关
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序列自相关序列')

%m序列构成的信号（矩形脉冲）
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序列矩形形成信号的自相关');xlabel('t');


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