Dada una matriz cuadrada de decir tamaño400x400, ¿cómo haría para dividir esto en submatrices constituyentes de20x20 utilizando un for-loop? ¡Ni siquiera puedo pensar por dónde empezar!

Me imagino que quiero algo como:

[x,y] = size(matrix)

for i = 1:20:x
    for j = 1:20:y

pero no estoy seguro de cómo proceder. ¿Pensamientos?