【仿真总结】使用FFT对仿真数据滤波

问题描述

在simulink仿真的时候希望对一组数据进行FFT分析，并且用plot绘制出图像，主要问题在以下两部分：

1. 仿真模型输出的是1，1，n的数组，而plot不能处理三维数组
2. fft分析的程序

解决方案：

对于第一个问题，
首先在模块参数处按照如下方式选择，将out前缀去掉，不然记录的所有数据都会保存在这个结构体中，不方便后续的处理
之后通过scope记录数据，之后再设置-配置属性-记录中勾选记录数据到工作区，同时编辑相应的数据名称，将其保存为“带时间的结构体”格式

到此可以实现的是将仿真中的数据保存至matlab工作区，但是由于之后想通过plot将数据通过图标的格式表示出来，但是plot在上图中的“保存格式”不管选择什么格式最终都将数据保存为[1,1,x]的三维数组，而plot无法处理三维数组的成像。
其次通过reshape函数将三维数组变为二维，可以看到我仿真得到的数据保存至工作区如下：

具体数据保存在结构体signals的values中，

通过如下函数可将三维数据转为二维：

load PeData.mat
relPe = reshape(PeData.signals.values(1,1,:),[230366,1]);
time = PeData.time;
figure(1)
plot(time, relPe)

对于第二个问题
主要参考的是help文档中fft相关的实例

%计算双侧频谱与单侧频谱
y = fft(relPe);
P2 = abs(y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
%绘制频域图像
f = Fs*(0:(L/2))/L;
figure(2)
plot(f,P1) 

这里的relPe就是之前第一部处理好的数据，我的问题就解决好了，能够看到相应的频域图像。
PS：本身打算通过fft将扰动通过频域图的形式观测出来，之后简单的做一个滤波，之后发现可能我的仿真并不是出现扰动hh又是白干的一天。