Лучшие решения перечислены Лорен . Также возможно создать эти матрицы, используя SPDIAGS :
vec = 101:105;
A = full(spdiags(repmat(vec,5,1),0:4,5,5)); % The second matrix
B = fliplr(full(spdiags(repmat(fliplr(vec),5,1),0:4,5,5))); % The first matrix
Я помню, как создавал подобные матрицы, прежде чем узнал о некоторых встроенных функциях, упомянутых Лорен. Это не так просто и чисто, как их использование, но это сработало. =)