资源描述:
《将simulink的scope波形数据保存到workspace》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、在用Simulink做仿真时,我们经常会用到示波器Scope来观察波形,它可以对波形进行局部放大、按横、纵座标放大,非常方便,但是如果我们要保存波形时,就最好别直接拷贝Scope波形了,因为它的背景是黑的,而且不能进行线形修改和标注,不适合作为文档用图。一般的做法是将数据输出到工作空间,然后用画图指令Plot画图。输出到工作空间的方法一般有这么几种:1.添加ToWorkspace模块;2.添加out模块;3.直接用Scope输出。本人比较懒,一般不再添加其他输出模块,直接选用方法3。当然不是说放一个Scope就能数出数
2、据的,需要对Scope进行设置。设置界面如下:在仿真结束后,打开scope,点击第二个图标,叫做parameter,选择datahistory,再在savedatatoworkspace前勾上,再仿真一次,数据就保存在workspace里了。这里最好把Limitdatapointstolast勾掉,因为很有可能你的数据会超过5000个。勾选SavedatatoWorkspace,变量类型可以选结构体,结构体带时间,以及向量(后面我们会分别介绍这几种变量类型的画图方法)。运行Simulink,输出完数据,你就可以利用Ma
3、tlab的画图工具随心所欲的画图了。下面以一个例子分别介绍三种变量类型的画图方法。1.输出类型为向量形式。从图上看到,输出了两维时间序列,而实际输出到工作空间的变量ScopeData为三维序列,其中第一列为时间,这正好为我们画图提供了方便。我们可以采用画图命令如下:figure;plot(ScopeData(:,1),ScopeData(:,2),'LineWidth',1.5);holdon;plot(ScopeData(:,1),ScopeData(:,3),'r:','LineWidth',1.5);legend
4、('正弦波','锯齿波');holdoff;当然你还可以采用其他绘图方式,如采用Subplot方式。2.输出类型为StructurewithTime。即结构体带时间。我们可以看一下这个结构体包含哪些东西。在CommandWindow里直接输入变量名。ScopeData=time:[51x1double]signals:[1x1struct]blockName:'untitled/Scope'可见,该结构体包含了时间序列,信号结构体,以及我的框图名。实际上我们的输出信号都包含在signals这个结构体里了,我们接着可以再
5、看看signals结构体的组成部分。我们输入ScopeData.signals(这点与C语言是类似的)。ans=values:[51x2double]dimensions:2label:''title:''plotStyle:[00]可以看到,values是一个51x2的double型矩阵,它正好是我们输出的数据。我们采用这样的画图命令即可完成画图:figure;plot(ScopeData.time,ScopeData.signals.values(:,1),'LineWidth',1.5);holdon;plot(
6、ScopeData.time,ScopeData.signals.values(:,2),'r:','LineWidth',1.5);legend('正弦波','锯齿波');holdoff;结果同上。3.对于Structure类型,正好是Structurewithtime的精简版,因为它的时间为空,因此你必须用其他方式获得时间,这里就不介绍了。上述三种输出方式,1最简单,但有时候你又不得不用结构体形式画图,例如下面这种情况,这时候你就只能用2和3了,当然最好还是2。可以试试下面这段程序:f=[1,2,3,4,5];%示
7、波器对应坐标系的序号,比如一个示波器有8个量显示,你想选择第1,2,3,4,5个量另外输出到新的figure窗口。a=1;%取起始点,这里从第一个点开始,当然你也可以从其他比如第一100个点开始作为起始点。b=size(ScopeData.time,1);%时间采样点总个数,即从0开始,以系统默认步长为间隔到末时间点的总个数。d=fix(5/5*b);%选择数据范围,这里5/5表示整个数据范围,3/5表示取原范围的3/5,c=1;%取点步长,1表示每个点都取,如c=5则表示每隔5个点取数据。H=figure(1);%新
8、建图形窗口fori=1:size(f,2)h(i)=subplot(size(f,2),1,i);%将几个曲线图画在一个figure里,并竖排,并且每个坐标图都有一个句柄值,以便后续处理需要,比如下面的xlabeltitle之类的。plot(ScopeData.time(a:c:d),ScopeData.signals(f(i))