%--------------------------------------------------------------------------%
%???????仿真例子為：時間序列分析--王振龍著????例題6.6的源代碼???????????
%???????????????原始數據處理及差分、繪圖、0均值化處理
%--------------------------------------------------------------------------%
D0=[10.93?9.340?11.00?10.98?11.29?11.84?10.62?10.90?12.77?12.15?12.24?12.30?...
????9.910?10.24?10.41?10.47?11.51?12.45?11.32?11.73?12.61?13.04?13.14?14.15?...
????10.85?10.30?12.74?12.73?13.08?14.27?13.18?13.75?14.42?13.95?14.53?14.91?...
????12.94?11.43?14.36?14.57?14.25?15.86?15.18?15.94?16.54?16.90?16.86?18.10?...
????13.70?10.88?15.79?16.36?17.22?17.75?16.62?16.96?17.69?16.40?17.51?19.73?...
????13.73?12.85?15.68?16.79?17.59?18.51?16.80?17.27?20.83?19.18?21.40?23.76?...
????15.73?13.14?17.24?17.93?18.82?19.12?17.70?19.87?21.17?21.44?22.14?22.45?...
????17.88?16.00?20.29?21.03?21.78?22.51?21.55?22.01?22.68?23.02?24.55?24.67?...
????19.61?17.15?22.46?23.19?23.40?26.26?22.91?24.03?23.94?24.12?25.87?28.25];
D00=[20.99?17.04?23.16?24.96?26.2?27.89?24.77?26.32?26.75?26.51?28.57?31.14];
subplot（131）
plot（D0‘r‘）;
title（‘原始工業數據曲線‘）;
grid?on

for?i=13:length（D0）
????D1（i-12）=D0（i）-D0（i-12）;
end
subplot（132）
plot（D1‘r‘）;
title（‘一階差分去周期曲線‘）;
grid?on;

S1=sum（D1）;
V=var（D1）;
S1=S1/length（D1）;
D11=（D1-S1）/sqrt（V）;%去均值化
subplot（133）
plot（D11‘r‘）;
title（‘一階差分去周期和0均值化曲線‘）;
grid?on;

%--------------------------------------------------------------------------%
%????????????????????自相關函數系數計算機繪圖
%--------------------------------------------------------------------------%
R0=D11*D11‘/length（D11）;
for?i=1:10
????r（i）=0;
????for?j=（i+1）:length（D11）
????????r（i）=D11（j）*D11（j-i）/（96-i）+r（i）;
????end
end
P=[R0?r];
P=[R0?r]/R0;????????%得到自相關系數P矩陣

figure
subplot（121）
plot（P‘b‘）;
grid?on
title（‘自相關系數曲線‘）;
axis（[1?9?0?0.5]）;

U=2/sqrt（96）;???????%評價指標曲線，2倍si格馬法則

%--------------------------------------------------------------------------%
%????????????????????偏自相關函數系數計算機繪圖
%--------------------------------------------------------------------------%
KK（1）=1;
PLL1=[P（1）?P（2）];
PLL2=toeplitz（PLL1）;
PRR=[P（2）?P（3）];
kk=inv（PLL2‘*PLL2）*PLL2‘*PRR‘;
KK（2）=kk（2）;????????????????%Y-W方程求解偏相關系數??2

PLL1=[PLL1?P（3）];
PLL2=toeplitz（PLL1）;
PRR=[PRR?P（4）];
kk=inv（PLL2‘*PLL2）*PLL2‘*PRR‘;
KK（3）=kk（3）;????????????????%Y-W方程求解偏相關系數??3

PLL1=[PLL1?P（4）];
PLL2=toeplitz（PLL1）;
PRR=[PRR?P（5）];
kk=inv（PLL2‘*PLL2）*PLL2‘*PRR‘;
KK（4）=kk（4）;????????????????%Y-W方程求解偏相關系數??4

PLL1=[PLL1?P（5）];
PLL2=toeplitz（PLL1）;
PRR=[PRR?P（6）];
kk=inv（PLL2‘*PLL2）*PLL2‘*PRR‘;
KK（5）=kk（5）;????????????????%Y-W方程求解偏相關系數??5

PLL1=[PLL1?P（6）];
PLL2=toeplitz（PLL1）;
PRR=[PRR?P（7）];
kk=inv（PLL2‘*PLL2）*PLL2‘*PRR‘;
KK（6）=kk（6）;????????????????%Y-W方程求解偏相關系數??6

PLL1=[PLL1?P（7）];
PLL2=toeplitz（PLL1）;
PRR=[PRR?P（8）];
kk=inv（PLL2‘*PLL2）*PLL2‘*PRR‘;
KK（7）=kk（7）;????????????????%Y-W方程求解偏相關系數??7

subplot（122）
plot（KK‘b‘）
grid?on
title（‘偏自相關系數曲線‘）;
axis（[1?7?0?1]）;

%-------------

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-14?18:25??MATLAB_時間序列\
?????文件????????4480??2014-08-22?16:15??MATLAB_時間序列\A1.m
?????文件????????2639??2018-04-14?18:04??MATLAB_時間序列\A2.m
?????文件?????7299722??2018-04-14?16:28??MATLAB_時間序列\MATLAB在時間序列分析中的應用.pdf
?????文件??????538528??2018-04-14?16:30??MATLAB_時間序列\MATLAB在時間序列建模預測及程序代碼.pdf
?????文件???????26193??2014-08-22?16:21??MATLAB_時間序列\定階曲線圖.PNG
?????文件??????300574??2018-04-14?16:35??MATLAB_時間序列\時間序列分析模型構建與MATLAB實現_譚贇.pdf
?????文件???????21627??2014-08-22?16:14??MATLAB_時間序列\第一個成果圖.PNG