如何实现一个高/低通滤波器对加速度计数据。

49视图(30天)

显示旧的评论

CAN105 2022年9月28日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1814025-how-to-implement-either-a-high-low-pass-filter-on-accelerometer-data

评论道: 明星黾 2022年9月29日

你好所有的,

我已经收集的数据从一个加速度计和阅读,试图理解过去的帖子现在几个小时但是我不相信我足够了解信号处理妥善implment我需要的过滤器。对于我的加速度计,我将输出数据率设置为50赫兹的带宽将25 hz和采样10赫兹(10次/秒)占奈奎斯特。

我的数据是非常吵,我想提取频率,基于类似的研究在这个领域应该是0.1 - 1 hz范围之间。还从我读过的研究论文,似乎先前的研究要么使用巴特沃斯高通滤波器和低通滤波器。

我实现了一个FFT首先但是我不太清楚如何辨别和accuratley应用过滤信号在所需的范围内通过获得一个过滤数据集。下面我将提供一个形象和代码(不确定是否正确)。

我会很感激,如果有人将我如何实现过滤数据集所以我可以清楚地辨别的峰值频率数据,过滤数据集前进牌汽车。我很业余的信号处理和Matlab !谢谢你,我真的想理解它。

                          X = readmatrix (“30 min.csv”);
                         
                          %的阴谋(X)
                         
                          Y = fft (X);
                         
                          L = 133361;%的长度矩阵
                         
                          Fs = 10;%采样频率(赫兹)
                         
                          P2 = abs (Y / L);
                         
                          P1, P2 (1: L / 2 + 1);
                         
                          P1 (2: end-1) = 2 * P1 (2: end-1);
                         
                          f = f * (0: (L / 2)) / L;
                         
                          情节(f, P1)
                         
                          标题(“单面振幅谱的X (t)”)
                         
                          包含(“f (Hz)”)
                         
                          ylabel (“| P1 (f) |”)

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

答案(1)

明星黾 2022年9月28日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1814025-how-to-implement-either-a-high-low-pass-filter-on-accelerometer-data answer_1062790

” 对于我的加速度计,我将输出数据率设置为50赫兹的带宽将25 hz和采样10赫兹(10次/秒)占奈奎斯特。 ”

让我解释这是困难的。如果采样频率 50 赫兹,奈奎斯特频率 25 赫兹,所以它应该是相对简单的过滤信号的范围 0.1 赫兹, 1 赫兹。

信号似乎宽带噪声。最简单的方法来消除与Savitzky-Golay过滤器( sgolayfilt 函数)因为一个频率选择过滤器将passd噪声以及信号的通频带。使用带通 , 高通滤波或低通滤波器其他过滤功能,并指定 “ImpulseResponse”、“信息检索” ,等待最好的结果。这些功能将设计一个非常有效的椭圆滤波器。

4评论
显示3年长的评论隐藏3年长的评论

CAN105 2022年9月29日

谢谢你的澄清!

我已经尝试了Savitzky-Golar过滤器和得到一些东西但不是entirley确定它是正确的,我不想走太多的兔子洞,所以到目前为止我在这里把我的进步。我说一些一般性的问题,但我当然不希望你回答他们所有人!我主要是想看看我是否在正确的道路上前进。我不确定如何实现过滤数据超过2赫兹和仍然低于0.01……

Ultimatley我需要花10分钟块分析主导frequecy跨越更大的csv文件。

                                   %明星黾建议Savitzky-Golay过滤器首先降低宽带
                                  
                                   %噪音然后紧随其后的是一个带通、高通滤波、低通滤波器函数
                                  
                                   其他过滤…% %
                                  
                                   X = readmatrix (“20220825 _1hr_rollonly.csv”);%建立加速度计数据的矩阵
                                  
                                   % % Savitzky-Golay过滤器适用于光滑的数据
                                  
                                   订单= 3;%,我们如何确定订单?
                                  
                                   framelen = 6001;%为什么framelen需要奇怪吗?
                                  
                                   %在10分钟的数据块
                                  
                                   yg = sgolayfilt (X,秩序,framelen);
                                  
                                   情节(X)
                                  
                                   持有在
                                  
                                   情节(yg)
                                  
                                   传奇(“原始”,“sgolay”)
                                  
                                   % %应用高通滤波器来平滑数据
                                  
                                   fs = 10;%采样率(10 hz)
                                  
                                   本产品=高通滤波(yg、0.01 fs, ImpulseResponse =“信息检索”);
                                  
                                   情节(yh);