基于MATLAB的m序列产生
更新时间:2024-03-03 22:21:01 阅读量: 综合文库 文档下载
第一章 设计内容及要求
基于MATLAB产生m序列
要求:
1.通过matlab编程产生m序列的产生原理及其产生方法。
2.对特定长度的m序列,分析其性质,及其用来构造其它序列的方法。
1
第二章 m序列设计方案的选择
2.1 方案一
MATLAB编程非常简单,无需进行变量声明,可以很方便的实现m序列。 2.2 方案二
图2.1 Simulink实现m序列
Simulink是MATLAB最重要的组件之一,它提供了一个动态系统建模,仿真和综合分析的集成环境。在此环境中无需大量书写程序,而只需通过简单直观的鼠标操作,就可构造出复杂的系统。Simulink具有适应性广,结构及流程清晰及仿真精细等优点,基于以上优点,Simulink已被广泛的运用到控制理论和数字信号处理的复杂仿真和设计。
通过比较方案一和方案二,发现方案一的有点具有通用性而方案二利用MATLAB的Simulink直接搭建模块,在移位寄存器较少的情况下利用此方法比较简单,可是当移位寄存器的个数增多时,要搭建那么多的模块就显的很繁琐了,缺乏通用性,因此本次实验选择方案一。
2
第三章 m序列的产生及性质
3.1 m序列的产生原理、结构及产生
m序列是最长线性反馈移位寄存器序列的简称,m序列是由带线性反馈的移位寄存器产生的。
由n级串联的移位寄存器和反馈逻辑线路可组成动态移位寄存器,如果反馈逻辑线路只由模2和构成,则称为线性反馈移位寄存器。
带线性反馈逻辑的移位寄存器设定初始状态后,在时钟触发下,每次移位后各级寄存器会发生变化,其中任何一级寄存器的输出,随着时钟节拍的推移都会产生一个序列,该序列称为移位寄存器序列。 n级线性移位寄存器的如图3.1所示:
A ◇
图3.1 n级线性移位寄存器
图中Ci表示反馈线的两种可能连接方式,Ci=1表示连线接通,第n-i级输出加入反馈中;Ci=0表示连线断开,第n-i级输出未参加反馈。
因此,一般形式的线性反馈逻辑表达式为
------表达式3.1
将等式左边的an移至右边,并将an=C0an(C0=1)带入上式,则上式可以写成
-------表达式3.2
定义一个与上式相对应的多项式
--------表达式3.3
其中x的幂次表示元素的相应位置。该式为线性反馈移位寄存器的特征
3
多项式,特征多项式与输出序列的周期有密切关系。当F(x)满足下列三个条件时,就一定能产生m序列:
(1)F(x)是不可约的,即不能再分解多项式; (2)F(x)可整除xn+1,这里p=2n+1; (3)F(x)不能整除xn+1,这里q 满足上述条件的多项式称为本原多项式,这样产生m序列的充要条件就变成了如何寻找本原多项式。 3.2 m序列的基本性质 (1)均衡性。在m序列一个周期中‘1’的个数比‘0’要多1位,这表明 4 序列平均值很小。 (2)m序列与其移位后的序列模2相加,所得的序列还是m序列,只是相位不同而已。例如:1110100与向又移3位的序列1001110相对应模二相加后的序列为0111010,相当于原序列向右移一位后的序列,仍为m序列。 (3)m序列发生器中移位寄存器的各种状态,除全0状态外,其他状态只在m序列中出现一次。 (4)m序列发生器中,并不是任何抽头组合都能产生m序列。理论分析指出,产生的m序列数由下式决定: --------表达式3.4 其中φ(X)为欧拉数。例如5级移位寄存器产生31位m序列只有6个。 (5)m序列具有良好的自相关性,其自相关系数: -------表达式3.5 从m序列的自相关系数可以看出m序列是一个狭义伪随机码。 图3.2 m序列信号的自相关函数 3.3 生成m序列的模块 5 根据m序列的生成原理图,如图3.1所示,由图可知m序列是多级移位寄存器通过线性反馈再进行模二相加产生的,最后一位一位输出观察其波形图。程序中使用的代码如下: N=2^length(reg1)-1; for k=1:N a_n=mod(sum(reg1.*coeff1(1:length(coeff1)-1)),2); reg1=[reg1(2:length(reg1)),a_n]; out1(k)=reg1(1); end 其中N为m序列的长度,值为(2-1)。由程序已定义了移位寄存器的初始状态和抽头系数,在此基础上进行反馈,后进行模二加,所得的结果为输出的第一个值,初始状态向左移移位,而所得的输出值填补上变成新的序列,在进行第二次反馈和模二加,依次循环N次,产生m序列 第四章 m 6 6 序列构造其他序列 Gold序列具有三值互相关函数,其值为: ?这里, 0.5(r?1)?1?2,r为奇数 (式2.11) p?2?1,t(r)??0.5(r?2),r为偶数但不是4的倍数?1?2r1pt(r),?1p,1p[t(r)?2] (式2.10) 当r为奇数时,gold序列中约有50%的码序列归一化相关函数值为-1/p。当r为偶数但又不是4的倍数是,约有75%的码序列归一化互相关函数值为-1/p。 Gold序列是R.Gold于1967年提出来的,它由两个m序列按下述方法演变而来的:把2个码长相同的m序列移位并进行模2加,如果相加的两个m序列是一对优选对,则相加的结果为一个Gold序列。 设有一对周期为N=2r-1的m序列优选对{a},{b},以其中任意一个序列为基准序列,如{a},对另一个序列{b}进行移位i次,得到{b}的移位序列{bi},然后与序列{a}进行模二加得到一个新的周期为N的序列{c},则称新序列{c}为Gold序列,既 {ci}={a}+{b} i=0,1,2,....,N 4.1 Gold序列的产生方框图 产生gold序列的程序代码如下: gold=mod(out1+out2,2); 第五章 程序调试及运行结果 5.1仿真设计流程图 7 否 确认程序无 误 仿真准备 规划程序模块 编写程序 是 否 是否达到要 是 图5.1 实验仿真流程图 结束 求 生成仿真图形 运行,调试及数据分析 设置仿真参数 5.2实验的调试与运行结果 程序中把移位寄存器的初始值定义全为1,抽头系数定义为[1000011]和 8 [1100111],根据公式m序列的长度=2n-1,可知道所得的两个m序列的长度都为63,所利用的移位寄存器为6个。代码如下: reg1=ones(1,6); %寄存器初始状态 coeff1=[1 0 0 0 0 1 1]; %设置系数 reg2=ones(1,7); %寄存器初始状态 coeff2=[1 0 0 1 1 1 1]; %设置系数 程序检测无误后,运行程序,得到图形如下: 图5.2 运行后的两m序列波形仿真图 根据产生Gold序列的原理,运行程序,得到如下Gold序列的仿真图; 9 图5.3 运行后Gold序列的仿真图 自相关性:首先将第一个m序列变成双极性的序列,在与本身进行移位相乘进行积分运算,代码如下: out1=2*out1-1; %变为双极性序列 for j=0:N-1 rho(j+1)=sum(out1.*[out1(1+j:N),out1(1:j)])/N; end j=-N+1:N-1; rho=[fliplr(rho(2:N)),rho]; figure(3) plot(j,rho); axis([-10 10 -0.1 1.2]);title('第一个m序列的自相关函数') 互相关性:第一个m序列的函数与第二个m序列函数的移位相乘进行积分运算。代码如下: 10 for j=0:N-1 R(j+1)=sum(out1.*[out2(1+j:N),out2(1:j)]); 图5.4 m序列自相关性仿真图 11 图5.5 两m序列的互相关性仿真图 5.3 错误排除 实验过程中出了许多错误,特别是在对序列的自相关性上,缺乏了对自相关性的了解,导致出的波形错误,但在查找资料,修改自相关函数,终于得到了正确的仿真图。 12 结论 课程设计是一个十分有价值、有意义的实践活动,把一个课题设计好不是一步到位的,是经过反复修改,不断调试的过程,期间有困难也有乐趣,使对工程实践有了一个初步的认识。 本次课程设计实现了设计要求,利用软件实现m序列的生成,通过这次实验不但加深了对m序列的了解,而且对MATLAB的编程有了很好的掌握,虽然在仿真过程中会出现一些如程序不对或出的仿真图没有达到实验要求,如m序列中‘1’的个数要比‘0’的个数多1,而实际出的仿真图‘1’和‘0’的个数却是相等的,最后在调整了顶层原理图才使得仿真图正确,在不断的程序调整中提高了自己的能力。 m序列可以用软件实现,也可以用硬件实现,但是通过此次实验看出了软件的诸多优点。在课程设计的过程中,查询了大量的资料,通过相关资料的查询,使我对通信领域的有关知识有了一定的了解,扩大了知识面。 13 参考资料 [1]肖国镇,梁传甲.伪随机序列及其应用[M].北京:国防工业出版社.1985 [2]林可祥,汪一飞.伪随机码的原理与应用[M].北京:人民邮电出版社.1998. [3]吴先用,邹学玉.一种m序列伪码发生器的产生方法[J].西安:西安电子科技大学出版社.2003 14 附录 clear; reg1=ones(1,7); %寄存器初始状态 coeff1=[1 0 0 0 0 1 1]; %设置系数 N=2^length(reg1)-1; %产生m序列 for k=1:N a_n=mod(sum(reg1.*coeff1(1:length(coeff1)-1)),2); reg1=[reg1(2:length(reg1)),a_n]; out1(k)=reg1(1); end reg2=ones(1,7); %寄存器初始状态 coeff2=[1 10 0 1 1 1]; %设置系数 N=2^length(reg2)-1; for k=1:N a_n=mod(sum(reg2.*coeff2(1:length(coeff2)-1)),2); %馈 reg2=[reg2(2:length(reg2)),a_n]; % out2(k)=reg2(1); %值输出 end %产生gold序列 gold=mod(out1+out2,2); c=1:N; figure(1) [b1,t1]=stairs(c,out1); subplot(2,1,1);plot(b1,t1); axis([0 130 -0.1 1.1]);title('第一个m序列'); [b2,t2]=stairs(c,out2); subplot(2,1,2);plot(b2,t2); 15 移位,反反馈 取第一个 axis([0 130 -0.1 1.1]);title('第二个m序列'); figure(2) [b3,t3]=stairs(c,gold); plot(b3,t3); axis([0 130 -0.1 1.1]);title('gold序列') out1=2*out1-1; %变为双极性序列 out2=2*out2-1; %自相关函数 for j=0:N-1 rho(j+1)=sum(out1.*[out1(1+j:N),out1(1:j)])/N; end j=-N+1:N-1; rho=[fliplr(rho(2:N)),rho]; figure(3) plot(j,rho); axis([-10 10 -0.1 1.2]);title('第一个m序列的自相关函数') %互相关函数 for j=0:N-1 R(j+1)=sum(out1.*[out2(1+j:N),out2(1:j)]); end j=-N+1:N-1; R=[fliplr(R(2:N)),R]; figure(4) plot(j,R); axis([-N N -20 20]);title('两个m序列的互相关函数'); 16 axis([0 130 -0.1 1.1]);title('第二个m序列'); figure(2) [b3,t3]=stairs(c,gold); plot(b3,t3); axis([0 130 -0.1 1.1]);title('gold序列') out1=2*out1-1; %变为双极性序列 out2=2*out2-1; %自相关函数 for j=0:N-1 rho(j+1)=sum(out1.*[out1(1+j:N),out1(1:j)])/N; end j=-N+1:N-1; rho=[fliplr(rho(2:N)),rho]; figure(3) plot(j,rho); axis([-10 10 -0.1 1.2]);title('第一个m序列的自相关函数') %互相关函数 for j=0:N-1 R(j+1)=sum(out1.*[out2(1+j:N),out2(1:j)]); end j=-N+1:N-1; R=[fliplr(R(2:N)),R]; figure(4) plot(j,R); axis([-N N -20 20]);title('两个m序列的互相关函数'); 16
正在阅读:
基于MATLAB的m序列产生03-03
生物药剂学与药代动力学前六章习题04-19
世界重大劫机事件一览09-05
西城区2011初三二模物理试题 - 图文05-12
新浮桥桥台基础土方开挖支护方案04-22
货币经济学—开放经济的宏观经济学ppt04-21
江苏省泰兴市黄桥镇七年级数学上学期第一次独立作业试题(无答案03-15
表态发言:在省委第X巡视组巡视XX县情况反馈会议上的主持词及表06-15
公司项目工作及工作计划04-29
- 人教新课标必修4 Unit2 Working the land名师导航
- 毕业生“校漂族”大行其道 - 0
- 江苏各市中考作文题出炉 - 0
- 暑期精品班 - 三角形 - 图文
- 情人节送什么礼物好??超强礼物已抵达
- 工程项目管理制度1
- 第四次业务学习 2016
- 会计要素与会计科目
- 欠发达地区小企业会计准则运用问题研究
- 一级锅炉水G4题库
- BBD双进双出筒式磨煤机安装使用说明书 SM-1
- 初一数学有理数教案
- 渝北区房地产评估市场调研报告
- iWebMall 数据字典
- 2018年小学入学教育工作计划
- 计量专业实务与案例分析 - 模拟题三 - 2013年版
- 启示录讲义
- 路基灰土改良(方案)
- 人行反洗钱岗位准入培训测试题集
- 2015电大《学前儿童发展心理学》期末试题及答案
- 序列
- 基于
- 产生
- MATLAB
- 2015年绵阳东辰国际学校小升初招生模拟题
- 毕业论文-溶胶凝胶法制备纳米二氧化钛
- 数学人教B版选修2-1课后训练:3.1.2 空间向量的基本定理 含解析
- 箱梁预制施工技术方案
- 经典宋词名句赏析精选句
- 廉政风险防控机制建设工作培训讲稿
- 现代管理实务案例(10)
- Winhex手动恢复NTFS分区中误删除的文件 - 图文
- 幼儿教师如何做行动研究
- 中山大学物理科学与工程技术学院(2) - 图文
- 政法干警面试历年真题(四)
- 群面经典真题(沙漠求生记 月球求生记 荒岛求生记等)(2)
- 安全生产管理条例(2008.8.18)
- 《西方经济学》院考试题及解答2009
- 中国移动铁通宽带上网操作指导手册 - 图文
- 高中数学知识点《解析几何》《圆锥曲线》《曲线参数方程》精选强
- 明星代言广告的风险及应对策略
- 文化娱乐场所试卷(有答案)
- 浅析FTTx接入网技术应用
- 2013《复变函数论》试题库及答案