How to compute the moving average of a signal with 6-hour window length and shifted 15 min at a time?

2 次查看（过去 30 天）
Mohanad Alkhodari
Mohanad Alkhodari 2019-9-22
评论： Mohanad Alkhodari 2019-9-27
Hello,
I am working on ECG signal and I want to determine the HR.
How to compute the moving average of a signal with 6-hour window length and shifted 15 min at a time?
Sampling frequency is 200.
Thank you

请先登录，再回答此问题。

采纳的回答

Dimitris Kalogiros
Dimitris Kalogiros 2019-9-22
Hi Mohanad
Bellow, you will find an example, where a moving average corresponding to a 6h window, is applied upon a signal taken within an interval of 24h with Fsampling=200Hz
clearvars; close all; clc;
%% input data
Fs = 200; % 200 Hz sampling frequency
time = 0:(1/Fs):(24*3600-1); % data collected within 24 hours
data = sin(2*pi*(1/2E5).*time)+0.2*randn(size(time)); % an example of input data: a sine with some noise
%% averanging
Kstep = (15*60)*Fs; % step of 15 minutes
Kwindow = (6*3600)*Fs; % aneranging window is 6 hours
meanData=[]; meanTime=[];
for k=1:Kstep:length(data)-Kwindow
dataWindow=data(1:k+Kwindow);
meanData=[meanData mean(dataWindow)];
meanTime=[meanTime time(k+floor(Kwindow/2))];
end
figure;
%mean values are depocted at the middle of their corresponding 6h window
plot(time, data); hold on;
plot(meanTime, meanData, '-*'); zoom on; grid on;
xlabel('time'); ylabel('data values');
  1 个评论
Mohanad Alkhodari
Mohanad Alkhodari 2019-9-27
I tweaked this code and it worked perfectly.
Fs = 200; % 200 Hz sampling frequency
time = ((0:17280000)/Fs)'; %24 hours in seconds
Kstep = (15*60)*Fs; % step of 15 minutes
Kwindow = (6*3600)*Fs; % averaging window is 6 hours
meanData=[];
meanTime=[];
k_old = 1;
for k=1:Kstep:length(data)-1
dataWindow=data(k_old:k+Kstep);
meanData=[meanData mean(dataWindow)];
meanTime=[meanTime time(k+floor(Kstep/2))];
k_old = k+Kstep;
end
figure;
%mean values are depocted at the middle of their corresponding 6h window
plot(time, data); hold on;
plot(meanTime, meanData, '-*'); zoom on; grid on;
xlabel('time'); ylabel('data values');

请先登录，再进行评论。

更多回答（1 个）

Steven Lord
Steven Lord 2019-9-22
Use movmean with a datetime or duration SamplePoints input. Specify the window using a duration array.
  1 个评论
Mohanad Alkhodari
Mohanad Alkhodari 2019-9-27
I did not get this.
The signal is in samples, I use movemean but it only provides the window length as inputs.
How to use SamplePoints for steps of 15 minutes each?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Signal Generation and Preprocessing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by