How to create a moving average filter to a specific signal
231 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Hi
I wanna apply a moving average filter to my temperature dataset, it is a 92-day daily temperature dataset. How can I create the filter?
Thanks a lot!
0 comentarios
Respuestas (2)
Chunru
el 28 de Nov. de 2021
Editada: Chunru
el 28 de Nov. de 2021
x = readmatrix("2013.txt");
n = 7; % filter order
y = filter(ones(n, 1)/n, 1, x);
plot(x); hold on;
plot(y);
legend('Original', 'MovAve')
To get rid of delay:
figure
y = filtfilt(ones(n, 1)/n, 1, x);
plot(x); hold on;
plot(y);
legend('Original', 'MovAve')
3 comentarios
Image Analyst
el 28 de Nov. de 2021
Look at the movmean() function.
2 comentarios
Image Analyst
el 28 de Nov. de 2021
@Yian Chen I don't see any zeros, and you might want to use more than 3 elements:
y = readmatrix('2013.txt')
plot(y, 'b.-', 'LineWidth', 2, 'MarkerSize', 20);
grid on;
xlabel('Index', 'FontSize',fontSize)
% Filter it
windowWidth = 9
smoothy = movmean(y, windowWidth);
hold on;
plot(smoothy, 'r.-', 'LineWidth', 2, 'MarkerSize', 20);
legend('Original', 'Smoothed')
You could also try sgolayfilt() to smooth it. You can choose an order like 2 or 3 to follow the original data better than movmean() which uses a linear fit. An order of 2 or 3 will fit a quadratic or cubic to the data in each window position.
Ver también
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!