實(shí)驗(yàn)三 m序列的仿真驗(yàn)證.doc

實(shí)驗(yàn)三 m序列的仿真驗(yàn)證1、 實(shí)驗(yàn)?zāi)康模豪胢atlab驗(yàn)證m序列的產(chǎn)生方法及其自相關(guān)特性。2、 實(shí)驗(yàn)要求：設(shè)m序列的生成多項(xiàng)式為， 求（1）m序列的輸出及其自相關(guān)序列；（2）設(shè)脈沖成形為P（t）=畫出其m序列信號(hào)的自相關(guān)函數(shù)；（3）設(shè)脈沖波形為升余弦成形（a=0），畫出其m序列信號(hào)的自相關(guān)函數(shù)。3、 實(shí)驗(yàn)原理：M序列即為maximal length linear feedback shift register sequence 它具有類似與隨機(jī)噪聲的某些統(tǒng)計(jì)特性，有時(shí)可以重復(fù)產(chǎn)生的。四、實(shí)驗(yàn)源碼%m序列發(fā)生器及其自相關(guān) mseq.mclear all;close all;g=19;%G=10011;state=8;%state=1000L=1000;%m序列產(chǎn)生N=15;mq=mgen(g,state,L);%m序列自相關(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.3 1.2);title(m序列自相關(guān)序列)%m序列構(gòu)成的信號(hào)（矩形脈沖）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);figure(1)subplot(221);plot(t,mt(1:length(t);axis(0 63 -0.3 1.2);title(m序列矩形成型信號(hào))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.3 1.2);title(m序列矩形成型信號(hào)的自相關(guān));xlabel(t);Tc=1;dt = Tc/N_sample;t=-20:dt:20;gt=sinc(t/Tc);mt = sigexpand(1-2*mq, N_sample);mt = conv(mt,gt);st2=sigexpand(1-2*mq(1:15),N_sample);s2=conv(st2,gt);st2=s2;rt2=conv(mt,st2(end:-1:1)/(N*N_sample);subplot(224);t1=-55+dt:dt:Tc*L-dt;%plot(t,mt(1:length(t);plot(t1,rt2(1:leng