Further filtering of accelerometer data

2 次查看(过去 30 天)
Hi Star Strider, I find your answers regarding the filters to be quite insightful.
I have a similar problem. I have used the Butterworth filter of 5th order and used 0.7Hz and 50 Hz as cutoff frq for passband filter. I am quite sure the frequncy beyond this frequncy region is not of point of interest. My sampling frequncy : 200Hz.
Still I can see there is noise as it can be seen in the attched figure.
Can you please help me to further denoise the data/filter data? I have attached the .fig file aswell for your reference.
Thanks in advance.acc8 filtered.jpg
  2 个评论
Nikhil Sonavane
Nikhil Sonavane 2019-10-9
It would be great if you can share me your data and filter parameters.
PIYUSH MOHANTY
PIYUSH MOHANTY 2019-10-9
Hi Nikhil,
Please find my data as attached.
Here is my filter parameters:
d=dlmread('acc8.txt');
t=dlmread('Time.txt'); %File having sampling time instances
order = 5;
fcutlow = 0.3;
fcuthigh = 50;
fs=200;
[b,a]=butter(order,[fcutlow fcuthigh]/(fs/2),'bandpass');
X_filt=filtfilt(b,a,d);
figure()
plot(t,d);
hold on
plot(t,X_filt);%2
hold off

请先登录,再进行评论。

回答(1 个)

Nikhil Sonavane
Nikhil Sonavane 2019-10-9
After looking at the Fourier Transform (figure attached) of your input data it can be concluded that it majorly comprises of frequencies from 0Hz to 50Hz. If you use a bandpass filter with cut off frequencies of 0.3Hz and 50Hz, you are allowing most of the frequencies through the filter and hence you are not able to remove the noise. Try changing the cut off frequencies to eliminate the spikes in your output data. fft.jpg
  5 个评论
PIYUSH MOHANTY
PIYUSH MOHANTY 2019-10-11
Further, reducing the frequncy cutoff points is not helping much. What else kind of noise can be there in the signal?
Nikhil Sonavane
Nikhil Sonavane 2019-10-11
Consider reducing the frequency range of the filter to less than 10Hz. Also try using a higher order filter to get a sharp transition.

请先登录,再进行评论。

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by