 

用户972891

文章：75 阅读：208804 评论：14 赞：745

 好友  私信个人主页

文章 75
原创 0
阅读 208804
评论 14
赞 745

原创 [原创]simulink快捷入门--7阶PN序列

 2009-3-23 23:12  5363 8 7 分类: 通信

晚上还有点儿时间，把matlab拿出来倒腾了一下，今天接触了下simulink，在这里把过程记录下来！

其实本不打算这么早折腾simulink的，一切源于昨天看到的一篇文章《基于Matlab的m序列的仿真与实现》，在这篇文章的第二页有个图，如下：

图1

看上去有点儿意思，于是我就想就把这个图实现了看看输出的是啥玩意儿！

但是那个论文上只是给出了拓扑，并没有手把手的教（人家肯定不会在论文上写这个了^_^），关键问题是我连怎么启动simulink都还是未解之谜...

这不是网络世界么，于是上网搜了个simulink的入门文档《simulink快速入门》，其中P14有个例子，那个例子是这样滴：

图2

从图2对比来看，要实现图1的结构只要找到图1中用到的三个原件就ok了！

不着急，首先在simulink中使用search功能搜了一把，都齐了！又不着急的翻了翻所有的simulink库。基本上知道这个玩意儿是咋回事儿了，hoho，东西看来都基本相通，看起来和Pspice一个道理的东东！

下面进入操作部分，就是实现图1的拓扑的操作步骤了（稍有不一样，我们这里使用的PN码本源多项式为X6+X5+1）：

1.启动simulink；

两种方式：a.在matlab的command中输入simulink；

b.点击matlab主窗口中的这个图标（下图中红框中的）

图3

这样，我们就进入了simulink的世界，如下：

图4

2.新建simulink文件，作出如图1的仿真拓扑。

a.在图4中点击新建文件图标（显然是那个纯净的白图标）,弹出如下窗口：

图5

b.添加仿真模块；

b.1在图4所示的窗口中找到Commonly Used Blocks（显然是simulink下拉的第一个），找到下面红框标注的东东：

图6

使用鼠标左键把它们一个个拖到图5的窗口中（当然unit delay要拖7次，如果你知道拖一次以后可以使用经典的两键式操作，那我想花时间看这篇文档有点儿浪费了）；

b.2修改参数。

b.2.1双击标示为logical operator的模块，配置参数为2输入XOR门如下：

图7

b.2.2双击标示为unit delay的模块，配置初始化值为1，如下：

图8

b.2.3修改仿真时间为127，如下：

图9

c.连线（鼠标左键先放在输入/输出，按下左键不放，拖到欲连接模块的输出/输入端口，注意必须入/出配对，否则连不上）。连好的图形如下:

图10

d.保存文件为m_seq_ex.mdl；

图11

e.运行仿真（simulation--start或者点击图11中标示的按钮）；

f.查看结果（双击scope模块），如下:

图12

看看结果，是不是7阶PN码？

我们可以用下面的程序比较下:

%m伪随机序列Matlab源代码
%5阶m序列
% 在MATLAB命令窗口输入以下：
% fbconnection=[0 1 0 0 1];
% mseq="m"_sequence(fbconnection);
% mseq

function mseq="m"_sequence(fbconnection)
n=length(fbconnection);
N=2^n-1;
register=[ones(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

保存为m_sequence.m;

下面的程序，

fbconnection =[0 0 0 0 0 1 1 ];
mseq = m_sequence(fbconnection);

保存为m_run.m

在m_run.m中点击运行，在matlab的command中输入

>> mseq

结果如下：

mseq =

Columns 1 through 18

1 1 1 1 1 1 1 0 0 0 0 0 0 1 0 0 0 0

Columns 19 through 36

0 1 1 0 0 0 0 1 0 1 0 0 0 1 1 1 1 0