В качестве альтернативы решению @ dantswain, вы можете перейти непосредственно к линейным индексам, предполагая, что вы всегда выбираете из столбцов:
r = size(A,1);
A( (1:r).' + (b-1) * r)
Это будет быстрее, но не обязательно более четким.
К сожалению, элегантного решения не существует.