实验环境为matlab2013b

1、首先编写一个mseq.m文件,内容为:

function[mseq]=m_sequence(fbconnection) 

n=length(fbconnection); 

N=2^n-1;  

register=[zeros(1,n-1) 1];  %移位寄存器的初始状态 

mseq(1)=register(n);        %m序列的第一个输出码元 

for i=2:N      

newregister(1)=mod(sum(fbconnection.*register),2);     

for j=2:n          

newregister(j)=register(j-1);     

end;      

register=newregister;     

mseq(i)=register(n); 

end

2、在matlab中打开这个文件

3、在matlab的命令行窗体中执行

primpoly(7,'all')

这个命令的目的是得到7阶M序列的全部本原多项式。假设想得到其它阶的m序列,自行改动数字7就可以。

4、得到的结果为

Primitive polynomial(s) = 

 

D^7+D^1+1

D^7+D^3+1

D^7+D^3+D^2+D^1+1

D^7+D^4+1

D^7+D^4+D^3+D^2+1

D^7+D^5+D^2+D^1+1

D^7+D^5+D^3+D^1+1

D^7+D^5+D^4+D^3+1

D^7+D^5+D^4+D^3+D^2+D^1+1

D^7+D^6+1

D^7+D^6+D^3+D^1+1

D^7+D^6+D^4+D^1+1

D^7+D^6+D^4+D^2+1

D^7+D^6+D^5+D^2+1

D^7+D^6+D^5+D^3+D^2+D^1+1

D^7+D^6+D^5+D^4+1

D^7+D^6+D^5+D^4+D^2+D^1+1

D^7+D^6+D^5+D^4+D^3+D^2+1





ans =





   131

   137

   143

   145

   157

   167

   171

   185

   191

   193

   203

   211

   213

   229

   239

   241

   247

   253

5、输入mseq([1 0 0 0 0 0 1])

6、便可得到一个127位的m序列,整理之后为

说明。

1 0 0 0   0 0 1这个输入參数的物理意义为各寄存器的初始状态,其数学意义是所选用的本原多项式D^7+D^1+1中D^1到D^7的各个系数。

读者能够试试自行选用其它本原多项式。得到不同的m序列。

比如假设选用D^7+D^3+1,则仅仅需调用

mseq([0 0 1 0 0 0 1]),得到的m序列为

在matlab中生成m序列的更多相关文章

  1. matlab中生成随机数的相关知识

    randperm()函数: 功能:用于生成从1到N的随机整数,并且没有重复,它本质上是一个随机排序的函数: 用法:1.  randperm(n)     随机生成从1到n的不重复的整数: 2. ran ...

  2. Matlab生成M序列的伪随机码

    伪随机编码中较常用的是m序列,它是线性反馈移位寄存器序列的一种,其特点是在相同寄存器级数的情况下输出序列周期最长.线性反馈移位寄存器的工作原理是,给定所有寄存器一个初始值,当移位脉冲到来时,将最后一级 ...

  3. matlab中的卷积——filter,conv之间的区别

    %Matlab提供了计算线性卷积和两个多项式相乘的函数conv,语法格式w=conv(u,v),其中u和v分别是有限长度序列向量,w是u和v的卷积结果序列向量. %如果向量u和v的长度分别为N和M,则 ...

  4. MATLAB中白噪声的WGN和AWGN函数的使用

    MATLAB中白噪声的WGN和AWGN函数的使用如下: MATLAB中产生高斯白噪声非常方便,可以直接应用两个函数,一个是WGN,另一个是AWGN.WGN用于产生高斯白噪声,AWGN则用于在某一 信号 ...

  5. Matlab中rand('state',sum(clock))解析

    一.问题来源 来自于一份PSO代码,PSO中需要初始化粒子位置和速度. 二.问题探究 众所周知,Matlab中的rand()函数产生的是伪随机数,但一般用来也可以接受.但是,如果我们知道伪随机数的初始 ...

  6. matlab中矩阵的表示与简单操作

    原文地址为:matlab矩阵的表示和简单操作 一.矩阵的表示在MATLAB中创建矩阵有以下规则: a.矩阵元素必须在”[ ]”内: b.矩阵的同行元素之间用空格(或”,”)隔开: c.矩阵的行与行之间 ...

  7. matlab中fspecial Create predefined 2-D filter以及中值滤波均值滤波以及高斯滤波

    来源: 1.https://ww2.mathworks.cn/help/images/ref/fspecial.html?searchHighlight=fspecial&s_tid=doc_ ...

  8. matlab中fopen 打开文件或获得有关打开文件的信息

    参考:https://ww2.mathworks.cn/help/matlab/ref/fopen.html?searchHighlight=fopen&s_tid=doc_srchtitle ...

  9. 如何在高并发分布式系统中生成全局唯一Id

    月整理出来,有兴趣的园友可以关注下我的博客. 分享原由,最近公司用到,并且在找最合适的方案,希望大家多参与讨论和提出新方案.我和我的小伙伴们也讨论了这个主题,我受益匪浅啊…… 博文示例: 1.     ...

随机推荐

  1. 2018美赛准备之路——Matlab基础——基本运算符号表示

    π pi ln(x) log(x)   lg(x) log10(x) log2(x) log2(x) 根号 sqrt(x) x的y次方 x^y e的y次方 exp(y)    

  2. JS打包与代码分割

    参考来源:https://github.com/ruanyf/webpack-demos#demo01-entry-file-source 后面的代码:https://github.com/94713 ...

  3. 如何禁用python警告

    有-W选项. python -W ignore foo.py 所属网站分类: python基础 > 综合&其它 作者:jiem 链接:http://www.pythonheidong.c ...

  4. ServletContext作用功能详解

    ServletContext,是一个全局的储存信息的空间,服务器开始,其就存在,服务器关闭,其才释放.request,一个用户可有多个:session,一个用户一个:而servletContext,所 ...

  5. .net中的协变和逆变

    百度:委托中的协变和逆变. 百度:.net中的协变和逆变. 协变是从子类转为父类. 逆变是从父类到子类. 这样理解不一定严谨或者正确.需要具体看代码研究.

  6. 修改K/3 Cloud管理中心端口

    有时候可能会应为端口号被占用或者数据隔离等等,不会使用K/3 Cloud默认的8000端口,这时候就设计到要修改端口号了,具体步骤如下: 1. 2. 打开{安装目录}\ManageSite\App_D ...

  7. 背包!背包!HDU 2602 Bone Collector + HDU 1114 Piggy-Bank + HDU 2191 512

    http://acm.hdu.edu.cn/showproblem.php?pid=2602 第一题 01背包问题 http://acm.hdu.edu.cn/showproblem.php?pid= ...

  8. ie6 js报错汇总

    假如ie6是一个非美女的人,那不知道为多少人乱刀砍多少次,假如不幸它是一个美女,那不知道要被多少人轮奸完之后还是乱刀砍死. ie6的js的错误说明总是很隐晦的,它喜欢和捉迷藏.报的行数和错误,常常不知 ...

  9. msp430项目编程05

    msp430中项目---TFT彩屏显示(续) 1.TFT彩屏工作原理 2.电路原理说明 3.代码(静态显示) 4.代码(动态显示) 5.项目总结 msp430项目编程 msp430入门学习

  10. [Android] 随时拍图像处理部分总结及源码分享

    http://blog.csdn.net/eastmount/article/details/45492065#comments [Android] 图像各种处理系列文章合集 http://blog. ...