Построение переходной матрицы цепи Маркова нескольких порядков в Matlab

Переходная матрица первого порядка из 6 состояний может быть построена очень элегантно, как следует за

 x = [1 6 1 6 4 4 4 3 1 2 2 3 4 5 4 5 2 6 2 6 2 6]; % the Markov chain
 tm = full(sparse(x(1:end-1),x(2:end),1)) % the transition matrix.

. Итак, вот моя проблема: как элегантно построить матрицу перехода второго порядка? Решение, к которому я пришел, выглядит следующим образом

 [si sj] = ndgrid(1:6);
 s2 = [si(:) sj(:)]; % combinations for 2 contiguous states
 tm2 = zeros([numel(si),6]); % initialize transition matrix
 for i = 3:numel(x) % construct transition matrix
   tm2(strmatch(num2str(x(i-2:i-1)),num2str(s2)),x(i))=...
   tm2(strmatch(num2str(x(i-2:i-1)),num2str(s2)),x(i))+1;
 end

Существует ли одно- или двухстрочная альтернатива без циклов?

--

Редактировать: Я попытался сравнить свое решение с решением Амро с помощью "x=round(5*rand([1,1000])+1);"

 % ted teng's solution
 Elapsed time is 2.225573 seconds.
 % Amro's solution
 Elapsed time is 0.042369 seconds. 

Какая разница! К вашему сведению, grp2idxдоступен онлайн.