%閆海停??學(xué)號S11093029
%AR時間序列模型估計

clear
tic
%%%%%%%%%%%%%%%第一步，模擬一個AR模型并繪制ACF，PACF圖%%%%%%%%%%%%%
%s首先設(shè)定AR模型的多項式系數(shù)。AR模型中只有多項式A（q）和C（q），
a1?=?0.6;
a2?=?0.8;
a3?=?0;
a4?=?0;
c1?=?0;
c2?=?0;
c3?=?0;
c4?=?0;?
obv?=?700;????????????????????????%obv是模擬的觀測數(shù)目。?
A?=?[1?a1?a2?a3?a4];
B?=?[];???????????????????????????%因為AR模型沒有輸入，因此多項式B是空的。
C?=?[1?c1?c2?c3?c4];
D?=?[];???????????????????????????%把D也設(shè)為空的。
F?=?[];???????????????????????????%AR模型里的F多項式也是空的。
m?=?idpoly（ABCDF11）????%這樣就生成了AR模型，把它存儲在m中。NoiseVariance被設(shè)定為1，1也是默認值。抽樣間隔Ts設(shè)為1。?
error?=?randn（obv1）;???????????%生成一個obv*1的正態(tài)隨機序列。準(zhǔn)備用作模型的誤差項。
e?=?iddata（[]error1）;????????%用randn函數(shù)生成一個噪聲序列。存儲在e中。抽樣間隔是1秒。???????????????????????。
y?=?sim（me）;?
get（y）????????????????????????????%使用get函數(shù)來查看動態(tài)系統(tǒng)的所有性質(zhì)。
Y=y.OutputData;?%把y.Ou