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

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