SS_to_AR.m
上传用户:mozhenmi
上传日期:2008-02-18
资源大小:13k
文件大小:1k
源码类别:

其他小程序

开发平台:

Matlab

  1. function [coef, C] = SS_to_AR(F, Q, k, diagonal)
  2. %
  3. % Extract the parameters of a vector autoregresssive process of order k from the state-space form.
  4. % [coef, C] = SS_to_AR(F, Q, k, diagonal)
  5. if nargin<4, diagonal = 0; end
  6. s = length(Q) / k;
  7. bs = s*ones(1,k);
  8. coef = zeros(s,s,k);
  9. for i=1:k
  10.   if diagonal
  11.     coef(:,:,i) = diag(diag(F(block(1,bs), block(i,bs))));
  12.   else
  13.     coef(:,:,i) = F(block(1,bs), block(i,bs));
  14.   end
  15. end
  16. C = Q(block(1,bs), block(1,bs));
  17. if diagonal
  18.   C = diag(diag(C));
  19. end
  20. %C = sqrt(Q(block(1,bs), block(1,bs))); % since cov(1,1) of full vector = C C'