arma.m

来自「arma模型的功率谱估计代码。还包括系统的极零点分析等完整软件包。」· M 代码 · 共 9 行

function [a, b, e] = ARMA(x, p ,q)
% 估计X(n)的ARMA模型参数,返回e*B/A中的A B e,e是sigma^2
N = length(x);
a = YuleWalker(x, p, q);
a = [1; a];
y = filter(a, 1, x);
y = y(p+1:N);
[b, e] = MA(y, q);
b = [1; b];