Daniel@0: function [coef, C] = SS_to_AR(F, Q, k, diagonal)
Daniel@0: %
Daniel@0: % Extract the parameters of a vector autoregresssive process of order k from the state-space form.
Daniel@0: % [coef, C] = SS_to_AR(F, Q, k, diagonal)
Daniel@0: 
Daniel@0: if nargin<4, diagonal = 0; end
Daniel@0: 
Daniel@0: s = length(Q) / k;
Daniel@0: bs = s*ones(1,k);
Daniel@0: coef = zeros(s,s,k);
Daniel@0: for i=1:k
Daniel@0:   if diagonal
Daniel@0:     coef(:,:,i) = diag(diag(F(block(1,bs), block(i,bs))));
Daniel@0:   else
Daniel@0:     coef(:,:,i) = F(block(1,bs), block(i,bs));
Daniel@0:   end
Daniel@0: end
Daniel@0: C = Q(block(1,bs), block(1,bs));
Daniel@0: if diagonal
Daniel@0:   C = diag(diag(C));
Daniel@0: end
Daniel@0: %C = sqrt(Q(block(1,bs), block(1,bs))); % since cov(1,1) of full vector = C C'