Как я могу ускорить этот цикл?

Question

У меня огромная матрица (5000x5000x100), и я пытаюсь сгладить каждый индекс по третьему измерению, но это занимает ЧАСЫ. Я делаю что-то неэффективно.

new_mat=zeros(size(my_mat));
for i = 1:length(mymat)
   for j = 1:length(mymat)
      new_mat(i,j,:) = wdenoise(squeeze(mymat(i,j,:)));
   end
end

Я знаю, что массивы и индексация могут помочь, но я не уверен, как их применять здесь. Спасибо за любую помощь.

Answer 1

У меня нет доступа к более новой Wavelet Toolbox, необходимой для wdenoise, но так как функция будет работать по столбцам, если вы предоставите матрицу, вы сможете удалить внутренний цикл, который может немного увеличить скорость:

new_mat=zeros(my_mat)
for i = 1:length(mymat)
    new_mat(i,:,:) = wdenoise(squeeze(mymat(i,:,:)));
end