Apply a filter to a dataset of ecg signals

Question

Sabaudian 2022-6-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1742355-apply-a-filter-to-a-dataset-of-ecg-signals

评论： Sabaudian 2022-6-17

dataset.zip

I have a dataset of 48 ECG signals (attached here as a zip file, but I use it unzipped), placed in a directory. I have use a for-cicle to load all the signal in the workspace, but now I am a bit confunsed.

I want to apply to the signals a Band-Pass Butterworth filter, so I use butter to design the filter and filtfilt to apply it. But I have a problem with size in "y(i)". This is the error I get:

"Unable to perform assignment because the indices on the left side are not compatible

with the size of the right side."

this is my code:

Fs = 360; % sample frequency of signals [Hz]
files = dir(fullfile("dataset/","*.mat"));
for i = 1:numel(files)
    data(i) = load(files(i).name); % load signals
    [b,a] = butter(3,[1 30]/(Fs/2),"bandpass"); % bandpass butterworth filter of third order with 1-30 Hz 
    y(i) = filtfilt(b,a,data(i).val); % apply the filter to all signals 
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Karim 2022-6-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1742355-apply-a-filter-to-a-dataset-of-ecg-signals#answer_987440

在 MATLAB Online 中打开

Hello, the problem is that you didn't initalize the variable "y". Hence you are trying to fit an array of 1x3600 into a single position (i.e. y(i) ).

See below for a quick fix:

files = dir(fullfile("dataset/","*.mat"));
Fs = 360; % sample frequency of signals [Hz]
[b,a] = butter(3,[1 30]/(Fs/2),"bandpass"); % bandpass butterworth filter of third order with 1-30 Hz 
y = zeros(numel(files), 3600);
for i = 1:numel(files)
    data(i) = load(fullfile("dataset/",files(i).name)); % load signals
    y(i,:) = filtfilt(b,a,data(i).val); % apply the filter to all signals 
end