MATLAB: Faltung einer matrixbewerteten Funktion

Question

Ich habe diesen Code geschrieben, um die 1-d-Faltung einer 2-d-matrixbewerteten Funktion durchzuführen (k ist mein Zeitindex, kend ist in der Größenordnung von 10e3). Gibt es einen schnelleren oder saubereren Weg, dies zu tun, vielleicht mit eingebauten Funktionen?

for k=1:kend
  C(:,:,k)=zeros(3);
  for l=0:k-1
      C(:,:,k)=C(:,:,k)+A(:,:,k-l)*B(:,:,l+1);
  end
end