Matlab实现M序列的产生和其自相关序列.docVIP

电子信息工程 专业课程设计任务书 学生姓名 专业班级 电信1班 学号 题 目 M序列产生器的MATLAB设计与实现 课题性质 工程技术研究 课题来源 自拟课题 指导教师 同组姓名 无 主要内容 了解m序列的产生、输出及其自相关序列，观察各种成形信号的波形。 任务要求 利用Matlab实现M序列的产生及其自相关序列。 观察成形信号波形。 参考文献 1．樊昌信，曹丽娜编著，通信原理(第六版)，国防工业出版社，2006 2．吴先用，邹学玉，一种m序列伪码发生器的产生方法[J].测控技术，2003，,22(9) 3．m序列是伪随机序列中最重要的一种，最长线性移位寄存器序列，m序列易于实现，具有优良的自相关特性，在直扩通信系统中用于扩展要传递的信号。n级线性移位寄存器的如图1所示： 图1 n级线性移位寄存器 图中表示反馈线的两种可能连接方式，=1表示连线接通，第n-i级输出加入反馈中；=0表示连接线断开，第n-i级输出未参加反馈。 因此，一般形式的线性反馈逻辑表达式为 将等式左面的移至右面，并将代入上式，则上式可改写为 定义一个与上式相对应的多项式 其中x的幂次表示元素的相应位置。式称为线性反馈移位寄存器的特征多项式，特征多项式与输出序列的周期有密切关系.当F(x)满足下列三个条件时，就一定能产生m序列： (1) F(x)是不可约的，即不能再分解多项式； (2) F(x)可整除,这里; (3) F(x)不能整除，这里qp. 满足上述条件的多项式称为本原多项式.这样产生m序列的充要条件就变成了如何寻找本原多项式. 根据m 序列的特征方程: 并根据其联接多项式编写Matlab 程序. 3 运行环境 硬件环境：Window xp 软件环境：Matlab 6.5 4 开发工具和编程语言 Matlab 6.5 5 详细设计clear all; close all; g=19;%G=10011; state=8;%state=1000 L=1000; %m序列产生 N=15; mq=mgen(g,state,L); %m序列自相关 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序列自相关序列) figure(2) %m序列构成的信号（矩形脉冲） 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(2) %subplot(221); plot(t,mt(1:length(t))); axis([0 63 -0.3 1.2]);title(m序列矩形成形信号) 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); figure(3) %subplot(223) plot(t,rt1(length(st):length(st)+length(t)-1)); axis([0 63 -0.3 1.2]);title(m序列矩形成形信号的自相关);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); figure(4) %subplot(224); t1=-55+dt:dt:Tc*L-dt; plot(t,mt(1:length(t))); plot(t1,rt2(1:length(t1))); axis([0 63 -0.5 1.2]);title(m序列since成形信号的自相关);xlabel(t) 调用的子程序如下： (1)mgen.m： function [out] = mgen(g,state,N) %输入 g:m序列生成多项式（10进制输入） %state:寄存器初始状态（10进制输入） %N:输出序列长度 % test g=11;state=3;N=15; gen =