armacanshuguji.m

利用先后估计法对arma模型进行参数估计-最小二乘法得到的AR部分的估计参数

%---------产生ARMA(2,2)时间序列----------%
ts0=idpoly([1 -1.5 0.7],[],[1 -1.7 0.1]);
ir=sim(ts0,[1;zeros(24,1)]);
e=idinput(200,'rgs');
y0=sim(ts0,e);  
y=y0';
%---------ARMA(2,2)对AR的参数进行估计-------------%
z=fliplr(y);
Ry=conv(y,z);
R=zeros(30,2);%ARMA(2,2)的R
r=zeros(30,1);%ARMA(2,2)的r 
for i=1:30
r(i,1)=Ry(198-i)/200;
for j=1:2
R(i,j)=Ry(198-i+j)/200;
end
end
a=inv(R'*R)*R'*r;%利用最小二乘法得到的AR部分的估计参数
A1=[1,-a']
%--------ARMA(2,2)对MA的参数进行估计-----------%
T=zeros(1,3);
for k=0:2
    for m=1:2
        for n=1:2
        T(1,k+1)=T(1,k+1)+a(m)*a(n)*((Ry(200-abs(k+n-m)))/200);
        end
    end
end
p=T';
r=roots(p);
b1=1/r(1)+1/r(2);
b2=-1/(r(1)*r(2));
B1=[1 -b1 -b2]