matlab - Matlab滤波器与平滑数据
问题描述
我找不到为什么该filter
函数似乎输出与该smoothdata
函数不同的东西的原因。5
他们都应该使用移动平均线来实现平滑的 windowSize 。精通这些功能的人可以解释发生了什么吗?
这是文档的 >链接<,下面的代码改编自那里,添加了smoothdata
(对于其他难题,smooth
函数也不相同):
% from documentation
t = linspace(-pi,pi,100);
rng default %initialize random number generator
x = sin(t) + 0.25*rand(size(t));
windowSize = 5;
b = (1/windowSize)*ones(1,windowSize);
a = 1;
y = filter(b,a,x);
%this is added for this example, you need >Matlab 2017a to run this
y2=smoothdata(x,'movmean',windowSize)
%y3=smooth(x,100); %bug in the code (obsolete)
y3=smooth(x,windowSize);
%now plot data
figure;
plot(t,x)
hold on
plot(t,y)
plot(t,y2)
plot(t,y3)
legend('Input Data','Filtered Data','smoothdata','smooth')
%show obvious parts of plot
xlim([0 3]);
ylim([0 1.25]);
这是我得到的输出:
这是情节的第一部分的一些不一致之处:
%this is added for this example, you need >Matlab 2017a to run this
y2 = smoothdata(x,'movmean',[windowSize-1,0]);
y3 = smoothdata(padarray(x,[0 2]),'movmean',[windowSize-1,0]);
%now plot data
figure(1); clf;
plot(t,x)
hold on
plot(t(1:10),y(1:10))
plot(t(1:10),y2(1:10))
plot(t(1:10),y3(1:10))
legend('Input Data','Filtered Data','smoothdata',['padded ' char(10) 'smoothdata'])
%show obvious parts of plot
xlim([-3.1784 -2.3858]);
解决方案
的window
参数smoothdata
,如果是标量,则产生一个以原点为中心的窗口。filter
实现一个因果过滤器,这意味着它采用windowSize
先前样本的平均值。因此,这两个结果之间的差异是windowSize/2
样本的偏移。您可以在情节中清楚地看到这种转变。使用两个参数来模拟 的结果filter
:
y2 = smoothdata(x,'movmean',[windowSize-1,0])
该smooth
函数还实现了移动平均,第二个参数是窗口大小。您在100
此处使用,而不是windowSize
,其值为5
。因此,这个结果是平均值的 20 倍以上。使用正确的窗口大小复制结果:
y3 = smooth(x,windowSize);
我认为结果y3
会发生变化filter
,就像y2
问题中一样。
推荐阅读
- python - Python子进程终端/空间需要调整为80/24以外的值
- css - Javascript/React:根据状态值更改显示的 SVG 背景
- xml - XSLT 使用属性值排序并创建目录 (TOC)
- javascript - Thunderbird Webextensions 插件获取消息列表
- angular - 更改角度材质多选的选定值 css
- excel - 使用 VBA 更改 XY 散点图上的线条颜色
- java - 在 iframe 中搜索按钮
- javascript - JavaScript 中的渲染何时更新?
- go - 输入自定义电子邮件声明时无法创建 JWT 令牌
- r - 如何绘制从秒到 HH/MM 的时间