MATLAB For 循环优化

Question

我是 MATLAB 的新手，还不知道它的功能和编程方式。

我有一个 3D 矩阵：19*1000*134，其中包含我想用巴特沃思过滤的数据。 我做了一个双循环的代码，效果很好：

%% Filtering
w1 = 1/(fs/2);
w2 = 45/(fs/2);
[b, a] = butter(6, [w1 w2], 'bandpass');

for i=1:channel
    for j=1:nb_epoch
        Data_filtered(i,:,j) = filtfilt(b, a, s_EEG.data(i,:,j));
    end
end

有没有其他方法可以更好地做同样的事情？由于我不了解 MATLAB，因此我不确定这是实现此带通滤波器的最佳方式。

谢谢！

Answer 1

我能想到的唯一优化是使用 parfors：

w1 = 1/(fs/2);
w2 = 45/(fs/2);
[b, a] = butter(6, [w1 w2], 'bandpass');

parfor i=1:channel
    parfor j=1:nb_epoch
        Data_filtered(i,:,j) = filtfilt(b, a, s_EEG.data(i,:,j));
    end
end

请记住，第一次运行 parfors 时，由于 worker 池是空的，因此脚本会慢得多，因为 worker 必须被实例化和初始化。多线程应该会有所帮助。