xx.m

arma模型的功率谱估计代码。还包括系统的极零点分析等完整软件包。

function x = xx(a, b, N1, N2)
% 由白噪声通过系统函数H(z)=B/A产生N1+N2个点的数据，前N1个点丢弃，输出N2个点

%load xdata; %use the same data
%return;

if nargin==0 
    N1 = 200;
    N2 = 256;
    a = [1, -1.3817, 1.5632, -0.8843, 0.4096];
    b = [1, 0.3544, 0.3508, 0.1736, 0.2401];
end
x0 = wgn(N1+N2, 1, 1);
x1 = x0(N1+1:N1+N2);
x = filter(b, a, x1);
save xdata x;
%test
%xx