MATLAB_时间序列

MATLAB的时间序列，有文档，有可运行的代码，文档里面有例题，有解答，也有代码。MATLAB的时间序列，有文档，有可运行的代码，文档里面有例题，有解答，也有代码。

%--------------------------------------------------------------------------%
%       仿真例子为：时间序列分析--王振龙著    例题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