arcoef.m

来自「谱估计及阵列信号处理算法仿真库」· M 代码 · 共 17 行

function Aq = ARcoef(q) 

% AR model generated signal
% q = model order

rand('state',sum(100*clock));
Vm = 0.8 + 0.1*rand(q/2,1);
rand('state',sum(100*clock));
Va = pi*rand(q/2,1);

j = sqrt(-1);
V=[Vm.*exp(Va*j) ; Vm.*exp(-Va*j)];

%the coefficients of the polynomial
Aq=poly(V);