function x = AR_mksample( a, num )
% AR model による時系列サンプルデータの生成
% x = AR_mksample( a, num )
%   a : parameter vector ( M row vector )
% num : number of sample points
% output
%   x : sample data $num$ column vector

M = length(a);
x = randn( num, 1 );
for t = (M+1):num
  x(t) = a' * x( (t-1):(-1):(t-M) );
end