伪随机序列matlab仿真
基于Matlab的m序列的仿真与实现
ａ＝［１１０１０１０００１１１０１１１１１００１００１１００００１０］； ｂ＝ａ； Ｎｂ＝ｌｅｎｇｔｈ（ｂ）； Ｎ＝１２０； ｘ＝［１：Ｎ］； ｆｏｒ ｋ＝１：Ｎ ｃ＝ｘｏｒ（ａ，ｂ）； Ｄ＝ｓｕｍ（ｃ）； Ａ＝Ｎｂ－Ｄ； ％ Ｒ（ｋ）＝（Ａ－Ｄ）／（Ａ＋Ｄ）； ｂ＝［ｂ（Ｎｂ），ｂ（１：Ｎｂ－１）］； ｐｌｏｔ（ｋ，Ｒ（ｋ），＇ｒｏ＇）； ｈｏｌｄ ｏｎ ｅｎｄ ｐｌｏｔ（ｘ，Ｒ，＇ｇ＇）； ｇｒｉｄ 运行结果：
－５
图二 ５ 级线性反馈移位寄存器产生 ｍ 序列的仿真系统 （２）利用 Ｓｉｍｕｌｉｎｋｍ观测 ｍ 序列时频域特性
图三 ５ 级线性反馈移位寄存器产生 ｍ 序列的时域波形
图四 ５ 级线性反馈移位寄存器产生 ｍ 序列的频域波形 ＴｏＷｏｒｋｓｐａｃｅ的ｗｒ５就是所得ｍ序列，我们取它一个完整
的周期，其周期是 ３１，序列如下： １１０１０１０００１１１０１１１１１００１００１１００００１０ ２．３ 利用 Ｍａｔｌａｂ 编程环境求 ｍ 序列的自相关特性 程序： ｃｌｅａｒ
Ｆ（ｘ）＝Ｃ ＋Ｃ ｘ１＋Ｃ ｘ２＋…＝
０１
２
（２）
（３）生成函数（亦称序列多项式）可以表示为如下形式
Ｇ（ｘ）＝ａ ＋ａ ｘ１＋ａ ｘ２＋…＝
０１
２
当初始状态为
ａ ＝１，ａ ＝ａ ＝…＝ａ ＝０
－ｍ
－ｍ＋１ －ｒ＋２
－１
联立式（１）、（２）、（３）可得
（３） （４）
国防工业工业出版社，２００６，５． ［４］查光明，熊贤祚．扩频通信［Ｍ］．西安：西安电子科技
大学出版社，２００２． ［５］曾兴雯，刘乃安，孙献璞．扩展频谱通信及其多址技术
［Ｍ］．北京：高等教育出版社，２００４． ［６］樊昌信，张甫翊，吴成柯．现代通信［Ｍ］．北京：高等
matlab论文
本原多项式
初始值
取样时间
3-2 伪随机码产生器模块参数窗口
1
0 .8
0 .6
0 .4
0 .2
幅 度 /v
0
-0 .2
-0 .4
-0 .6
-0 .8
-1
• 图3-3 伪随机码产生器输出的15位序列波形 • （载波频率为 1GH ，码元宽度 T 100 ns ，序列长度15位） • 由上图3-3可以写成15位伪随机序列: • ={ 1 1 1 -1 1 1 -1 -1 1 -1 1 1 1 1 1 } • 3.2 单极性伪码转换成双极性伪码模块 以上伪随机码产生模块形成的伪随机序列是单极性码，不过调制系统需要的是双 极性码，为此必须采用电平转换模块把单极性码变成双极性码。
0.8
幅度 （ V ）
0.6
0.4
0.2
X: 9.83e-005 Y: -0.069 X: 9.85e-005 Y: -0.063 X: 9.98e-005 Y: -0.067 X: 0.0001 Y: -0.068 X: 0.0001013 Y: -0.063 X: 0.0001015 Y: -0.071
• • • • • • • •
clear all; load simout1.mat; load simout2.mat; PN1=simout1; PN2=simout2; R=xcorr_tu (PN1,PN2); plot(R) 以上程序生成的自相关函数波形如下图3-9所示：
1
X: 9.84e-005 Y: 0.987 X: 9.99e-005 Y: 1 X: 0.0001014 Y: 0.987
图3-1 伪码调相过程原理框图 按照原理图可以用MATLAB/Simulink软件包建模实现伪随机码的调相信号框 图上各个组成部分的含义以及形成过程如下：
基于matlab的伪随机序列实现
衡量一个m序列与它的j次移位序列之间的相关 程度，并把它叫做m序列( a1 , a2 , a3 ,, a p)的自相关函 数。记作 R( j ) a a ，当采用二进制数字0和1 A 代表码元的可能取值时，R( j) A D A D D p （其中A表示 a a 0的数目，D表示 ai ai j 1的数目） 由m序列的均衡性和移位相加特性可知： j0 1 R( j ) 1 / p j 1,2,,( p 1)
由上图可以看出，寄存器状态满足： y1(i)=y1(i-1) y4(i-1); y2(i)=y1(i-1); y3(i)=y2(i-1); y4(i)=y3(i-1); 给寄存器组赋予一个非全零的初始状态，y4的输 出即可作为码周期为15的一个m序列。输入 mseq=mseries([1 0 0 1]) 可得 mseq = Columns 1 through 10 0 0 0 1 1 1 1 0 1 0 Columns 11 through 15 1 1 0 0 1
p i 1 i j i
三、伪随机序列的自相关性
i
i j
自相关程序：
a=[ 0 0 0 1 1 1 1 0 1 0 1 1 0 0 1]; b=a; m=length(b); N=100; x=1:N; for k=1:N c=xor(a,b); D=sum(c); A=m-D; R(k)=(A-D)/(A+D); plot(k,R(k),'ro') b=[b(m),b(1:m-1)]; 从图中可以看出在15、30、45、60、75、 hold on 90处得到了最大的相关特性，其他值处自 end 相关值接近于零，此仿真验证了m序列具 plot(x,R,'b') 有良好的自相关特性。 grid on
m序列产生和性能的MATLAB仿真
电厂检修业务的进入战略 就进入新业务战略而言，业务拓展可以分为由
点到线、由线到网的进入策略，存在着购并外部现有 企业、从头开始自创新企业以及与外部业务结盟开 展新业务３种作法，简称购并、自创、结盟。作为电 厂检修业务，浙江火电有着很强的技术基础，人员可 通过短期培训获得，且具备一定的资金实力和融资 能力，自创新企业不存在障碍，但考虑市场资源的获 得，可以考虑与拥有检修市场的单位进行结盟，提高 市场竞争能力，扩大市场规模。因此，最好的进人方 式就是以自主创新为主，以结盟为辅。
尺ｍ）。ｆⅣ
’
【一１
ｍ＝壳Ⅳ＿ｊ｝＝ｏ，±ｌ，±２，…。。
。
其他ｍ
２
ＭＡ吼ＡＢ仿真
为在ＭＡＴＬＡＢ中对ｍ序列发生器和序列性能
进行仿真，可傲絮下定义，将移位寄存器状态设置成 向量ｍｒｅｇｉｓｔｅｒ＝［ａｌ，ａ２，０３，…，ａＮ］，将连接向量设 为向量ｇ＝［ｇ，，９２，９３，…，ｇＮ］，由图ｌ可得模２加法 器的输进等于反馈量ｂａｃｋ，也帮是ａ，的下一个德，
ｍ＝ｚｅｒｏｓ（１，３１）；
ｍｒｅｇｉｓｔｅｒ＝ｍｉｎｉｔｉａｌ； ｎ＝０：ｋ＝０：
万 方数据
水利电力机械
２００７年１２月
５火电施工企业业务组合和选择发展战略
从上述分析中可以看出，浙江火电宜以火电施 工业为基础，大力发展钢结构加工业和电厂检修业 务，退出物流行业。
５．１
警誉芋磊薯君翥望盏霎｛娄馨等嚣喾需嚣嘉兰主
ｎ２＝［０：ｎ一１］；
ｓｕｂｐｌｏｔ（２，１，１） ｓｔｅｍ（ｎｌ，ｍ） ｙｌａｂｅｌ（’ｍ ｓｅｑｕｅｎｃｅ＇） ｓｕｂｐｌｏｔ（２，１，２） ｐｌｏｔ（ｎ２，ｒＥ） ｙｌａｂｅｌ（’Ｒ（ｍ）０
基于matlab的伪随机序列生成及相关函数仿真实验
武夷学院实验报告课程名称：_扩频通信技术及应用_ 项目名称：基于matlab的伪随机序列生成及相关函数仿真实验姓名：翁开传__专业：__通信工程___ 班级：__2__学号：_40__同组成员_无___1注：1、实验预习部分包括实验环境准备和实验所需知识点准备。
2、若是单人单组实验，同组成员填无。
2注：实验过程记录要包含实验目的、实验原理、实验步骤，页码不够可自行添加。
2)函数auto_corr（）计算二进制序列seq的自相关函数，并画出函数曲线。
在函数内调用了生成m序列的函数ms_generator（*）生成的m序列seq，然后在计算。
function seq=ms_generator(registers,connections)registers=[0 0 0 0 1];connections=[1 0 0 1 0 1];n=length(connections);L=2^(n-1)-1;seq(1)=registers(n-1);for i=2:Lsum=0;for m=1:(n-1)sum=mod(sum+registers(m)*connections(m+1),2);endfor k=(n-1):-1:2registers(k)=registers(k-1);endregisters(1)=sum;seq(i)=registers(n-1);endfunction auto_correlation=auto_corr(seq)registers=[1 0 0 0 0 ];%设置初始状态connections=[1 0 1 0 0 1];%设置反馈逻辑seq=ms_generator(registers,connections);seq=-1*(seq*2-1);%负逻辑映射len=length(seq);temp=[seq seq];for i=0:len-1auto_correlation(i+1)=seq*(temp(i+1: i+len))';endauto_correlationplot(0:len-1,auto_correlation);运行的结果如下3)函数cross_corr（）计算二进制序列seq1和seq2的互相关函数并画出曲线。
单片机matlab方式产生伪随机m序列
信息科学与技术学院通信原理课程设计课题名称：伪随机m序列发生器的设计学生姓名：张昕灏2018508087学信息科学与技术学院院：专业年级：电子信息项目2018级指导教师：田敏副教授二O—三年七月十二日完成日期：目录前言1第一章设计内容及要求21.1设计内容21.2设计要求21.3方案选择2第二章m序列的特性分析42.1m序列的原理42.2均衡特性52.3游程分布52.4线性叠加性62.5自相关特性6第三章m序列设计83.1设计流程图83.2特征多项式确定83.3本原多项式确定103.4 m 序列的最终产生＜以五阶移位寄存器举例）11 第四章设计成果分析及总结134.1仿真结果分析134.2设计总结14 心得体会15 参考文献16 附录matlab 程序17 附录51 单片机实现方法18电路图18设计说明18结果验证18C51 代码及与对应matlab 代码20 数模转换输出代码：20 反馈链接状态及波形输出控制代码22 使用器件、八、-前言扩展频谱通信是一种不同于常规通信系统的新调制理论和技术，简称扩频通信[1]。
其设计思想是将待传输的信息信号用特定的扩频码扩展频谱后成为宽带信号进行传输，接收时再采用相应的技术手段将频谱压缩，恢复原来待传信息信号的带宽，从而实现通信。
扩频通信具有两个特点：传输信号的带宽远大于原始信息信号的带宽；传输信号的带宽主要有扩频码决定，此扩频码通常是伪随机码。
伪随机码(pseudo randomcode> 简称PN 码，可以人为产生与复制，具有类似白噪声的性质，相关函数具有尖锐的特性，功率谱占据很宽的频带，易于从其他信号或干扰中分离出来，具有优良的抗干扰特性，其特点是：具有尖锐的自相关函数；互相关函数值应足够小；有足够长的码周期，以确保抗侦破与抗干扰的要求；码的数量足够多，以实现码分多址的要求；平衡性好，以满足抗干扰的要求；项目上易于产生、加工、复制与控制[2]。
扩频通信的优势主要来自于伪随机码具有白噪声的统计特性。
matlab产生伪随机数,改变M序列幅值
1.用下式产生伪随机数：clc;clear;p=35;M=2^p;n=7;A=2^n+1;c=1;length=1000for i=2:length;x(1)=3;x(i)=mod(A*x(i-1)+c,M);endy=x/M;plot([1:length],y,'.')这种产生随机数的方法为混合同余法，当M=2^35，n=21,c=1时，相关系数接近于零。
但当n取值过大时，占用机器时间增加。
Coveyou给出的一组参数M=2^35，n=7,c=1，虽n取值不大，却能有效的减小序列的相关系数，且不占用过多的机器时间。
2.①用给出的等式产生M序列clc;clear;x(1)=0;x(2)=1;x(3)=0;x(4)=1;x(5)=0;x(6)=1;x(7)=0;x(8)=1;x(9)=0;length=500for i=10:500;x(i)=xor(x(i-4),x(i-9));endplot(x)分析:首先对x(1)到x(9)随即赋值，系数都为1，随后的变量由等式计算出来，由图可知，该序列的循环周期很长，0和1出现的概率基本相同。
②以此序列为基础产生逆序列clc;clear;x(1)=0;x(2)=1;x(3)=0;x(4)=1;x(5)=0;x(6)=1;x(7)=0;x(8)=1;x(9)=0;length=500for i=10:500;x(i)=xor(x(i-4),x(i-9));endy(i)=mod(i,2)for i=1:500z(i)=xor(x(i),y(i))endplot(z)逆M序列是M序列和方波序列复合结果，也是周期较长的随机序列，0和1出现概率相同。
③并将逆序M的幅值变为-a和+aclc;clear;x(1)=0;x(2)=1;x(3)=0;x(4)=1;x(5)=0;x(6)=1;x(7)=0;x(8)=1;x(9)=0;length=500for i=10:500;x(i)=xor(x(i-4),x(i-9));endy(i)=mod(i,2)for i=1:500z(i)=xor(x(i),y(i))enda=2for i=1:500b(i)=a*(1-2*z(i))endplot(b)分析：取a=2，幅值变为-2和+2。
伪随机序列的产生与仿真
基于MATLAB 的伪随机序列的产生及相关特性的仿真一、相关概念:平稳随机过程的各态历经性, 随机信号的频谱特性, 自相关函数, 互相关函数二、工程背景与理论基础根据香农的理论，在高斯白噪声干扰情况下，在平均功率受限的信道上，实现有效和可靠通信的最佳信号是具有白噪声统计特性的信号。
扩频通信正是由此而来的，在扩频通信最大的优点就是具有强大的抗噪声性能，使有用信号几乎可以淹没在噪声传播。
故扩频通信对扩频序列一般有如下要求：（1）尖锐的自相关特性（2）尽可能小的互相关值（3）足够多的序列数，具有良好的伪随机性（4）序列均衡性好，0、1等概（5）工程上易于实现伪随机序列具有以上所以有点，故在CDMA 扩频通信系统中，伪随机序列被作为扩频码之一。
下面在理论上阐述下伪随机序列（即m 序列）的产生原理及其所具有的相关数学性质。
然后在用MATLAB 语言实现m 序列的产生，并就其相关特性进行仿真，仿真结果结果表明该方法是可行的。
1、 m 序列简单介绍m 序列是最长线性反馈移位寄存器序列的简称，是由带线性反馈的移位寄存器的周期最长的序列。
它是周期为rN=2-1的伪随机序列，r 是移位寄存器的阶数。
下面是IS-95CDMA 系统中I 信道引导PN 序列的生成多项式和线性反馈移位寄存器的框图。
I 支路生成表达式：15139875()1I P x x x x x x x =++++++ 123456789101112131415输出图1-1 I 路信号产生器m 序列具有以下基本性质：（1）均衡性：在m 序列的一个周期中，“1”的个数之比“0”的个数多一个。
这表明序列平均性很好，即“1”和“0”几乎就是随机出现的，具有较好的随机性。
（2）具有尖锐的自相关特性，相互不同码字之间几乎是完成正交的。
周期函数的自相关函数定义为：/2/201R()()()T s s T T s t s t dt ττ-=+⎰，式中0T 是s()t 的周期。
m序列伪随机码产生及应用仿真matlab
序列伪随机码产生及应用仿真matlab一、概述随机序列是一种具有随机性质的数字序列，可用于信息传输、通信系统、密码学、雷达等领域。
伪随机序列是一种经过数学算法产生的序列，其具有类似与随机序列的统计特性，但实际上是确定性的。
在通信系统中，伪随机序列广泛应用于码分多址技术、扩频通信、混沌通信等领域。
本文将介绍伪随机码的产生原理，并使用matlab进行仿真实现，以加深对该技术的理解。
二、伪随机码的产生原理伪随机码的产生主要包括线性反馈移位寄存器（LFSR）、加法（XOR）运算和乘法（AND）运算等步骤。
下面以LFSR为例，介绍伪随机码的产生原理。
1. LFSR原理LFSR是一种利用移位寄存器和反馈逻辑电路构成的伪随机码产生器。
在LFSR中，存在一个移位寄存器和一个反馈逻辑电路。
移位寄存器中存储了若干比特的信息，经过时钟信号的控制完成信息的移位操作。
而反馈逻辑电路则根据寄存器中的信息产生反馈信号，从而实现伪随机序列的产生。
2. 伪随机码的产生过程（1）初始化LFSR。
将移位寄存器中的初始状态设置为一个非零的值。
（2）循环移位寄存器。
根据时钟信号的控制，移位寄存器中的信息进行移位操作。
（3）根据反馈逻辑电路的输出，生成伪随机序列。
通过反馈逻辑电路生成的反馈信号，即为伪随机序列的一部分。
三、matlab仿真实现matlab是一种功能强大的科学计算软件，具有丰富的绘图和仿真功能。
下面将使用matlab进行伪随机码的产生和仿真实现。
1. 产生伪随机码在matlab中，可以使用shiftregister函数实现LFSR的移位寄存器功能。
结合matlab的位运算功能，可以方便地实现伪随机码的产生。
以下是一个简单的matlab代码示例：```matlab定义LFSR的初始状态state = [1 0 1 0 1];定义LFSR的反馈多项式polynomial = [5 2];产生伪随机码for i = 1:10获取LFSR的输出output = state(end);计算新的状态new_state = mod(sum(state(polynomial)), 2);更新状态state = [new_state, state(1:end-1)];显示输出disp(output);end```上述代码中，通过定义LFSR的初始状态和反馈多项式，使用循环产生了10个伪随机码的输出。
单片机matlab方式产生伪随机m序列
信息科学与技术学院通信原理课程设计课题名称：伪随机m序列发生器的设计学生姓名：张昕灏2018508087学院：信息科学与技术学院专业年级：电子信息项目2018级指导教师：田敏副教授完成日期：二○一三年七月十二日目录前言1第一章设计内容及要求21.1设计内容21.2设计要求21.3方案选择2第二章 m序列的特性分析42.1 m序列的原理42.2均衡特性52.3游程分布52.4线性叠加性62.5自相关特性6第三章 m序列设计83.1设计流程图83.2特征多项式确定83.3本原多项式确定103.4 m序列的最终产生<以五阶移位寄存器举例）11第四章设计成果分析及总结134.1仿真结果分析134.2设计总结14心得体会15参考文献16附录matlab程序17附录51单片机实现方法18电路图18设计说明18结果验证18C51 代码及与对应matlab代码20数模转换输出代码：20反馈链接状态及波形输出控制代码22使用器件23前言扩展频谱通信是一种不同于常规通信系统的新调制理论和技术，简称扩频通信[1]。
其设计思想是将待传输的信息信号用特定的扩频码扩展频谱后成为宽带信号进行传输，接收时再采用相应的技术手段将频谱压缩，恢复原来待传信息信号的带宽，从而实现通信。
扩频通信具有两个特点：传输信号的带宽远大于原始信息信号的带宽；传输信号的带宽主要有扩频码决定，此扩频码通常是伪随机码。
伪随机码(pseudo randomcode>简称PN码，可以人为产生与复制，具有类似白噪声的性质，相关函数具有尖锐的特性，功率谱占据很宽的频带，易于从其他信号或干扰中分离出来，具有优良的抗干扰特性，其特点是：具有尖锐的自相关函数；互相关函数值应足够小；有足够长的码周期，以确保抗侦破与抗干扰的要求；码的数量足够多，以实现码分多址的要求；平衡性好，以满足抗干扰的要求；项目上易于产生、加工、复制与控制[2]。
扩频通信的优势主要来自于伪随机码具有白噪声的统计特性。
基于MATLAB的伪随机序列的产生和分析
收稿日期： !""# $ %! $ %&
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
第 !! 卷
&’(’)*+,-( *(. "(*/01,1 -2 34 5’67’(8’ %*1’. -( !"#$"%
89:; 9<= $ 43>? ，/@ 93>?
（ 8344A?A 3B :4ACDE3>=C <>5 @>B3E1<D=3> ， ;3EDFGAHDAE> 634IDACF>=C<4 J>=KAEH=DI， L=’ <> MF<>N= +%""+!， 8F=><） "%5#9":#： ,-./-0 H3BDG<EA F<H OAA> PHA5 ANDA>H=KA4I => 1<>I A>?=>AAE=>? B=A45H 5PA D3 =DH HDE3>? 3QAE<D=3> B<>CD=3>7 .3 ANQ<>5=>? 3E C31QEAHH=>? DFA H=?><4 HQACDEP1 => HQEA<5 HQACDEP1 HIHDA1，DFA H=?><4 =H ?A>AE<44I 1P4D=Q4A5 OI < HQEA<5 HQACDEP1 HARPA>CA7 .FA CF<E<CDAE 3B HQEA<5 HQACDEP1 HARPA>CA H=?>=B=C<>D4I <BBACDH DFA C311P>=C<D=3> RP<4=DI7 @> <44 6; HAS RPA>CAH，1 $ HARPA>CA <>5 2345 $ HARPA>CA <EA 3BDA> PHA5 <H HQEA<5 HQACDEP1 HARPA>CA7 @> DF=H Q<QAE，DFA OE=AB =>DE35PCS D=3> 3B DFA DFA3EI，QE3QAEDI <>5 C3>HDEPCD=>? 1A<>H 3B DFA DG3 HARPA>CAH <EA ?=KA> B=EHD， <>5 DFA ?A>AE<D=3> <>5 <><4IH=H 3B DFA1 OI QE3?E<11=>? G=DF , 4<>?P<?A => ,-./-0 <EA ?=KA> 4<DAE7 .FA H=1P4<D=3> EAHP4DH HF3G DFA C3EEACD>AHH <>5 BA<H=O=4=S DI 3B DF=H 1ADF357 .FA H=1Q4A <>5 =>DP=D=KA 1ADF35 =H C3>KA>=A>D B3E DFA A>?=>AAE=>? QAEH3>>A47 ;<=>?9@5： 6; HARPA>CA；MARPA>CA；MQEA<5 HQACDEP1 HARPA>CA
基于MATLAB的伪随机序列发生器的设计 (2)
通信工程专业专业综合课程设计题目基于MATLAB的伪随机序列发生器的设计学生姓名学号所在院(系)专业班级指导教师完成地点2013年 12 月 29 日专业综合课程设计任务书院(系)专业班级学生姓名一、专业综合课程设计题目基于MATLAB的伪随机序列发生器的设计二、专业综合课程设计工作自 2013 年 12 月 20 日起至 2014 年 1 月 17 日止三、专业综合课程设计进行地点:四、专业综合课程设计的内容要求：1、要求在MATLAB软件环境下编写程序完成九级m序列反馈系数，本原多项式的求解。
2、要求完成相应的Gold序列发生器的设计;或采用软件编程来实现。
3、要求设置合适的系统参数，输出Gold序列的时域、频域波形。
指导教师系(教研室)通信工程系接受论文 (设计)任务开始执行日期2013年12月20日学生签名基于matlab的伪随机序列发生器的设计[摘要]伪随机序列码越来越受到人们的重视，被广泛用于导弹、卫星、飞船轨道测量和跟踪、雷达、导航、移动通信、保密通信和通信系统性能的测量以及数字信息处理系统中。
本文主是对基于matlab的伪随机序列发生器的设计，及其利用matlab 软件对其进行仿真和利用simulink对其仿真性能的研究，主要阐述了扩频系统中m序列和gold序列的产生。
在第一部分中介绍了课题研究的背景，第二部分中介绍了扩频系统的相关知识，第三部分介绍了m序列和Gold序列产生的原理和方法，第四部分利用matlab和simulink对其进行仿真。
[关键词]伪随机序列 m序列移位寄存器Design of the pseudo-random sequence generatorbased on matlab[Abstract]pseudo-random sequence code more and more get people's attention, is widely used in missiles, satellites, spacecraft orbit measurement and tracking, radar, navigation, mobile communications, and the measurement of the performance of the communication security and communication system of digital information processing system. In this paper, the main is to the design of pseudo-random sequence generator based on matlab, and the use of matlab simulation and the use of simulink software to the study of the simulation performance, mainly expounds the m sequence and gold sequence in spread spectrum system. In the first part introduces the research background, the second part introduces the related knowledge of spread spectrum system, in the third part introduces the m sequence and Gold sequence principle and method of the fourth part carries on the simulation using matlab and simulink.[key words]pseudo random sequence m sequence shift register目录1．引言 (4)1.1研究的背景及意义 (4)2．扩频通信系统简介 (5)2.1扩频通信的基本概念及相关模型 (5)2.1.1基本概念 (5)2.1.2 数学模型 (8)2.2扩频通信系统的主要特点 (10)2.3扩频通信系统分类 (11)2.4伪随机序列在扩频通信中的应用 (11)3.m序列 (11)3.1m序列的定义 (11)3.2m序列的原理 (12)3.3m序列的性质 (15)4.Gold序列 (17)4.1 Gold序列的产生原理 (17)4.2 Gold序列的性质 (18)5.MATLAB仿真实现 (18)5.1 MATLAB软件介绍 (18)5.2 m序列的仿真及分析 (19)5.2.1程序分析 (19)5.2.2 simulink分析 (20)5.3 Gold的仿真及分析 (21)5.3.1程序分析 (21)5.3.2 simulink分析 (22)6.致谢 (23)7.参考文献 (24)附录 (25)附录A (25)附录B (27)1．引言1.1研究的背景及意义移动通信由于具有时实性、机动性、具有不受时空限制等特点，己经成为一种深受人们欢迎的通信方式，并融入了现代生活当中。
常用伪随机码序列的相关性分析与MATLAB仿真
伪随机序列仿真实验报告
伪随机序列仿真仿真实验报告报告时间：2013/5/13姓名：刘敬亚学号：2010101026电子邮件：403555418@一、实验目的1、熟悉MATLAB 仿真的应用；2、掌握伪随机序列的原理、软件产生，并仿真分析其相关特性。
二、实验内容1、MATLAB中产生127位gold序列，并分析其自相关和互相关特性。
2、采用MATLAB的伪随机函数，生成127位二进制伪随机序列，并分析其自相关和互相关特性。
3、收集学习技术资料，列举五种以上通信系统中采用的伪随机序列。
三、实验过程记录根据搜集的资料在matlab中编写代码，然后运行，得到相应的特性图。
1、产生127位的gold序列2、gold序列的自相关性close all;clear all;clc;n=7;%移位寄存器级数P=2^n-1%m列的周期%产生m1序列%c1=[1 0 1 0 0 1 1];%反馈系数==211reg1=[zeros(1,n-1) 1];%初始化移位寄存器m1(1)=reg1(n);for i=2:Pnew_reg1(1)=mod(sum(c1.*reg1),2);for j=2:nnew_reg1(j)=reg1(j-1);endreg1=new_reg1;m1(i)=reg1(n);endsubplot(311);stairs(m1);xlabel('k');title('m1序列');%产生m2序列%c2=[1 0 1 1 0 0 1];%反馈系数==217reg2=[zeros(1,n-1) 1];%初始化移位寄存器m2(1)=reg2(n);for i=2:Pnew_reg2(1)=mod(sum(c2.*reg2),2); for j=2:nnew_reg2(j)=reg2(j-1);endreg2=new_reg2;m2(i)=reg2(n);endsubplot(312);stairs(m2);xlabel('k');title('m2序列');%产生gold序列%for i=0:P-1shift_m2=[m2(i+1:P) m2(1:i)];%m2g(i+1,:)=mod(m1+shift_m2,2);%模2加endsubplot(313)stairs(g(1,:));xlabel('k');title('gold序列的第1行');%gold序列的自相关特性和互相关特性%g1=g(1,:);g2=g(2,:);%取gold序列的第2行% for i=-P+1:-1% shift=[g1(i+P+1:P) g1(1:i+P)];% r1(P+i)=sum(g1(P+i).*shift(P+i)); % end% for i=0:P-1% shift=[g1(i+1:P) g1(1:i)];% r1(P+i)=sum(g1(i+1).*shift(i+1)); % endr1=conv(g1,g1);figure(2)subplot(211)stem(r1/P);xlabel('k');title('g1的自相关特性');% for i=-P+1:-1% shift=[g1(i+P+1:P) g1(1:i+P)]; % r2(P+i)=g2(P+i).*shift(P+i); % end% for i=0:P-1% shift=[g1(i+1:P) g1(1:i)];% r2(P+i)=g2(i+1).*shift(i+1); % endr2=conv(g1,g2);subplot(212);stem(r2/P);xlabel('k');title('g1与g2的互相关特性');产生的图形：四、分析及结论从实验结果可以得出，Gold序列间不仅有良好的互相关性能，而且个数远比m序列多。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB实验报告
姓名：***
学号：********
班级：14184111
实验七伪随机序列的产生与相关特性分析
一、实验目的
1、了解伪随机序列的相关知识。
2、了解m 序列的相关知识，了解其相关性质。
3、学会用matlab 实现方框图描述的系统，并由此产生m 序列。
二、相关理论知识
（1）伪随机序列
伪随机序列, 又称伪随机码, 是一种可以预先确定并可以重复地产生和复制, 又具有随机统计特性的二进制码序列。
在现代工程实践中, 伪随机信号在移动通信、导航、雷达和保密通信、通信系统性能的测量等领域中有着广泛的应用。
例如，在连续波雷达中可用作测距信号, 在遥控系统中可用作遥控信号, 在多址通信中可用作地址信号, 在数字通信中可用作群同步信号, 还可用作噪声源以及在保密通信中的加密作用等。
伪随机发生器在测距、通信等领域的应用日益受到人们重视。
伪随机信号与随机信号的区别在于: 随机信号是不可预测的, 它在将来时刻的取值只能从统计意义上去描述；伪随机序列实质上不是随机的, 而是收发双方都知道的确定性周期信号。
之所以称其为伪随机序列, 是因为它表现出白噪声采样序列的统计特性, 在不知其生成方法的侦听者看来像真的随机序列一样。
m 序列作为一种基本的PN 序列,具有很强的系统性、规律性和相关性。
（2）m 序列的产生
①线性反馈移位寄存器
m 序列发生器的系统框图如图。
其中加法器为“模2 相加”运算，寄存器与反馈的每一位只有1、0 两种状态。
由于带有反馈，因此在移位脉冲作用下，移位寄存器各级的状态将不断变化，通常移位寄存器的最后一级做输出，输出序列为
输出序列是一个周期序列。
其特性由移位寄存器的级数、初始状态、反馈逻辑以及时钟速率(决定着输出码元的宽度)所决定。
当移位寄存器的级数及时钟一定时，输出序列就由移位寄存器的初始状态及反馈逻辑完全确定。
当初始状态为全零状态时，移位寄存器输出全0序列。
为了避免这种情况，需设置全0 排除电路。
三、实验任务
编写程序利用5级移位寄存器产生m序列，初始状态全1，抽头系数[1 0 0 1 1 1]。
四.我自己的实验代码
clear all
close all
solve=0;
zk=0
a=[1,1,1,1,1]%an
c=[1,1,1,0,1,1]
lena=length(a)
lenb=length(c)
newc=fliplr(c)
shuchu=zeros(1,lena)
for zk=1:2^(lena+1)
for k=1:lena
zhong(k)=a(k)*newc(k);
lenz=length(zhong)
%solve=0;
end%%%%得到每个被加数
for t=1:lenz
solve=solve+zhong(lenz-t+1)
if solve>=2
solve=0
end
end%%%%得到最后一个加数
shuchu(zk)=a(1)
for i=2:lena
a(i-1)=a(i)
end
a(lena)=solve
%a(lena)=solve
solve=0;
end
stem(solve)
五.实验心得：
通过此次实验我学到了如何使用matlab产生伪随机序列，深刻体会了产生伪随机序列的原理，此次实验较之前的小实验较为综合复杂，同时也是对之前所学内容的回顾与巩固，在编写程序时也发现了自己之前一些表达方式掌握不够，将在接下来的实验中进一步巩固。