基于matlab的伪随机序列生成及相关函数仿真试验

1、武夷学院实验报告课程名称：扩频通信技术及应用项目名称：基于 matlab 的伪随机序列生成及相关函数仿真实验姓名：翁开传专业:通信工程班级：2 学号：40 同组成员二一、实验预习部分1：1,实验环境准备实验室的计算机和 matlab 软件和书本上的程序实验过程记录2:1.实验目的1) 了解伪随机序列的相关知识。2) 了解 m 序列的相关知识，了解其相关特性。3)熟悉 MATLA 的真的应用。4)掌握伪随机码的原理、软件产生，仿真分析其相关特性。2 .实验步骤1)函数 ms_generator(registers,connections)是 m 序歹!J 的生成函数，其中参数 registers

2、 给出了移位寄存器的初始状态，connections 给出了 m 序列的发生器。functionseq=ms_generator(registers,connections)registers=00001;connections=100101;n=length(connections);L=2A(n-1)-1;seq(1)=registers(n-1);注：1、实验预习部分包括实验环境准备和实验所需知识点准备。2、若是单人单组实验，同组成员填无。注：实验过程记录要包含实验目的、实验原理、实验步骤，页码不够可自行添加fori=2:Lsum=0;form=1:(n-1)sum=mod(sum+re

3、gisters(m)*connections(m+1),2);endfork=(n-1):-1:2registers(k)=registers(k-1);endregisters(1)=sum;seq(i)=registers(n-1);end运行输出的结果如下=C4ltiiuif11.9jooooioo1oijooi1iLColunns20iTiroulh31IOOQJIOLllOlO2)函数 auto_corr()计算二进制序列 seq 的自相关函数，并画出函数曲线。在函数内调用了生成 m 序列的函数 ms_generator(*)生成的 m 序列 seq,然后在计算。functionse

4、q=ms_generator(registers,connections)registers=00001;connections=100101;n=length(connections);L=2A(n-1)-1;seq(1)=registers(n-1);fori=2:Lsum=0;form=1:(n-1)sum=mod(sum+registers(m)*connections(m+1),2);endfork=(n-1):-1:2registers(k)=registers(k-1);endregisters(1)=sum;seq(i)=registers(n-1);endfunctionau

5、to_correlation=auto_corr(seq)registers=10000;%S 置初始状态connections=101001;%S 置反馈逻辑seq=ms_generator(registers,connections);seq=-1*(seq*2-1);%S 逻辑映射len=length(seq);temp=seqseq;fori=0:len-1auto_correlation(i+1)=seq*(temp(i+1:i+len);endauto_correlationplot(0:len-1,auto_correlation);运行的结果如下autQ.correlatiQn

6、=Colwui1through1931-1T-1-1-1-1-1-1-i-I-1T-1-1-1-1-1-1Colwws2Q31I-17T-1-1-177-1-I-IT3)函数 cross_corr()计算二进制序列 seq1 和 seq2 的互相关函数并画出曲线。在函数内调用 ms_generator(*)函数分别生成等长的M 序列 seq1 和 seq2 后在计算。程序也画出了自相关函数。functionseq=ms_generator(registers,connections)%registers=00001;%connections=100101;n=length(connection

7、s);L=2A(n-1)-1;seq(1)=registers(n-1);fori=2:Lsum=0;form=1:(n-1)sum=mod(sum+registers(m)*connections(m+1),2);endfork=(n-1):-1:2registers(k)=registers(k-1);endregisters(1)=sum;seq(i)=registers(n-1);endfunctionauto_correlation=auto_corr(seq)registers=10000;%iS 置初始状态connections=101001;%iS 置反馈逻辑seq=ms_ge

8、nerator(registers,connections);seq=-1*(seq*2-1);%5i 逻辑映射len=length(seq);temp=seqseq;fori=0:len-1auto_correlation(i+1)=seq*(temp(i+1:i+len);endauto_correlationplot(0:len-1,auto_correlation);functioncross_correlation=cross_corr(seq1,seq2)registers1=00001;%iS 置初始状态connections1=100101;%iS 置反馈逻辑seq1=ms_g

9、enerator(registers1,connections1);seq1=-1*(seq1*2-1);%S 逻辑映射registers2=11111;%iS 置初始状态connections2=111101;%iS 置反馈逻辑seq2=ms_generator(registers2,connections2);seq2=-1*(seq2*2-1);%逻辑映射temp=seq2seq2;len=length(seq2);fori=0:len-1cross_correlation(i+1)=seq1*(temp(i+1:i+len);end4)函数 gold()生成一族 gold 序列。程序中

10、首先生成一对 m 序列优先对 seql 和 seq2,然后按照并结构生成了一族 gold 序列，最后画出一个 gold 序列的自相关函数以及族内一对 gold 序列的互相关函数functionseq=ms_generator(registers,connections)%registers=00001;%connections=100101;auto_correlation=auto_corr(seq1);plot(0:len-1,auto_correlation,legend(自相关函数,互相关函数);xlabel(j);ylabel(相关函数);运行的结果如下,0:len-1,cross_

11、correlation);Colunms1through19-9-17-1-9-97-1Colunn20through31TT7一17TT7-1-11-97-977-1-1菰图林哭自相关函数互相关函数n=length(connections);L=2A(n-1)-1;seq(1)=registers(n-1);fori=2:Lsum=0;form=1:(n-1)sum=mod(sum+registers(m)*connections(m+1),2);endfork=(n-1):-1:2registers(k)=registers(k-1);endregisters(1)=sum;seq(i)=

12、registers(n-1);endfunctionauto_correlation=auto_corr(seq)registers=10000;%S 置初始状态connections=101001;%S 置反馈逻辑seq=ms_generator(registers,connections);seq=-1*(seq*2-1);%S 逻辑映射len=length(seq);temp=seqseq;fori=0:len-1auto_correlation(i+1)=seq*(temp(i+1:i+len),;endauto_correlationplot(0:len-1,auto_correla

13、tion);functioncross_correlation=cross_corr(seq1,seq2)registers1=00001;%S 置初始状态connections1=100101;%S 置反馈逻辑seq1=ms_generator(registers1,connections1);seq1=-1*(seq1*2-1);%S 逻辑映射registers2=11111;%S 置初始状态connections2=111101;%S 置反馈逻辑seq2=msgenerator(registers2,connections2);seq2=-1*(seq2*2-1);%逻辑映射temp=s

14、eq2seq2;len=length(seq2);fori=0:len-1cross_correlation(i+1)=seq1*(temp(i+1:i+len);endauto_correlation=auto_corr(seq1);plot(0:len-1,auto_correlation,-,0:len-1,cross_correlation);legend(自相关函数,互相关函数);xlabel(j);ylabel(相关函数);functiongold=gold()registers1=100000;connections1=1100001;seq1=ms_generator(regi

15、sters1,connections1);registers2=100000;connections2=1110011;seq2=ms_generator(registers2,connections2);temp=seq2seq2;len=length(seq2);fori=1:lengold(i,:)=mod(seq1+temp(i:i+len-1),2);endgold(len+1,:)=seq1;gold(len+2,:)=seq2;plot(0:len-1,auto_corr(gold(1,:),-,0:len-1,cross_corr(gold(1,:),gold(2,:),-.);legend(自相关函数,互相关函数);xlabel(j);ylabel(相关函数);运行的结果如下实验结果与讨论：3实验心得体会注：1、实验小结应包含实验所需知识点和实验方法的总结，实验心得