You have a problem in that the structure of your plot (hardcoded as 7x8) is not linked to the structure of your data (unique_rotation and unique_SF). But the general approach I would take is this:
if i == 1
ylabel(sprintf('Orientation = %d', unique_rotation(j))
end
if j == 1
title('construct your string here')
end
if i == length(unique_SF)
xlabel('construct your string here')
end