/實(shí)驗(yàn)報(bào)告課程名稱：3G移動(dòng)通信實(shí)驗(yàn)實(shí)驗(yàn)名稱:偽隨機(jī)序列生成編程(２）學(xué)生姓名：章仕波學(xué)生學(xué)號(hào)：11081522實(shí)驗(yàn)日期:5月15日實(shí)驗(yàn)六偽隨機(jī)序列生成編程(二)【實(shí)驗(yàn)?zāi)康摹考由顚?duì)偽隨機(jī)序列性能的理解;能夠使用Maｔlab(或者C語(yǔ)言)生成偽隨機(jī)序列;比較這些偽隨機(jī)序列的自相關(guān)函數(shù)和互相關(guān)函數(shù)特性.【實(shí)驗(yàn)內(nèi)容】畫出m序列的自相關(guān)函數(shù)和互相關(guān)函數(shù)并進(jìn)行分析；【實(shí)驗(yàn)步驟】畫出兩個(gè)m序列各自的自相關(guān)函數(shù)，以及這兩個(gè)序列的互相關(guān)函數(shù).（本原多項(xiàng)式為（4５）8，(75）8。）【實(shí)驗(yàn)報(bào)告】%編程產(chǎn)生兩個(gè)周期長(zhǎng)度為31的m序列；（本原多項(xiàng)式為（45)8，(７5）８.％(4５)8=100１０1，(75)8=１1１1０１clearａｌｌ；cｌc;ｒ=5;N=2^r-1；%周期Ｔ=30;Tmｉn＝1;a1=oｎes(1，r);％移位寄存器ｍ1＝zeroｓ(1，2^ｒ-１)；％生成的ｍ序列fｏri=1:(2＾r(nóng)-1）ｔemp=ｍod((ａ１（5）+a1(２)），2)；foｒj＝ｒ:－1：2a1(ｊ）=a1(j—1）；eｎdａ１(1)=ｔｅmp；m1（ｉ)＝a1(ｒ);endsｕbplｏt（2,1,1）steｍ（m１）;％序列1ａｘis([0,2^r，0，2]）；tｉtｌｅ('本原多項(xiàng)式m１為(45）8'）;a2＝oｎes(１，r)；％移位寄存器ｍ2＝zeroｓ(1,2^r—1）；％生成的m序列fｏｒi=1:（2^ｒ-１)tｅｍp=mｏｄ(（a2（5)+ａ2（４）+a2（3）+a2(2））,２）;foｒj=r:—１:2ａ2(ｊ）=a2（ｊ－1);ｅnｄａ2(１)=ｔemｐ;m２(ｉ）=a2（r）;％序列2endsｕbplot(2，１，2）steｍ(m2）；ａxiｓ(［０,2＾r(nóng)，0，２］)；tiｔlｅ（’本原多項(xiàng)式m2為（75)8’)；ｍ１=m1＊２—1；m2＝m2*2—1;％*＊*＊＊＊＊*＊＊＊**＊＊*＊＊＊*％ｆigｕｒｅRａa＝zeros（１，Ｔ+1）;fori=0:TRaa（i+１）＝0;fｏｒｊ＝1:Ｎ;k＝i+j；wｈiｌe（k>N）k=k－N;endRaa（i＋1）=Rａa（ｉ+１）+m1(ｊ）＊m1(k）；endＲaa(ｉ+1)＝Ｒａa(i+1)／N;eｎdsｕbpｌｏt(３,１,１)stｅm（0:T,Ｒａａ）；tiｔle('ｍ１自相關(guān)函數(shù)＇）；%＊**＊＊＊*＊***＊*＊＊***＊＊***＊Rbb=zｅｒｏs(1，T+1)；fｏrｉ=0：ＴＲｂb(i+1）=0;forj=1:N；k=i＋j；ｗhｉle(k＞Ｎ）k=k—N;endＲbb（ｉ+１）＝Rｂb(i＋１）＋m2（j)＊m２（ｋ);ｅｎdRbｂ（i+1)=Rbb（ｉ+１)/Ｎ;ｅndsｕbｐlot(3，1，２）stem(０：T,Rｂb)；title（'m２自相關(guān)函數(shù)＇);％＊＊**＊＊***＊＊＊＊*＊＊＊＊＊＊*＊＊*＊＊＊＊＊＊**＊*Ｒab=ｚeros(1，T+1)；foｒi=0：TRａb（i＋１）=0;forj＝1:Ｎ;k=i+j；whiｌe(k＞N)ｋ＝ｋ－N;eｎdRab(i＋1)=Rab（ｉ+１）＋ｍ１(j）*m2(ｋ）;endＲab（i+1）=Rab(i+１）/N；endｓubpｌｏt（3,1,3）steｍ(0:T,Ｒab）;ｔiｔlｅ（＇m1,ｍ2互相關(guān)函數(shù)’);【背景知識(shí)】【偽代碼】a＝ｚeｒｏs(１，r）；seq＿m=zeｒｏs（1，2^r-１);ｆｏrj=１：（