matlab时间序列分析在测绘中的应用

《matlab时间序列分析在测绘中的应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、《现代测绘数据处理方法》课程课间实验报告实验项目：MATLAB时间序列分析在测绘中的应用班级：测绘工程专业指导教师：一、实验目的及所用软件版本1、实验目的①了解MATLAB时间序列分析的基本原理及应用②学会用MATLAB时间序列的分析方法解决测绘工程中的实际问题2、实验软件所用版本MATLAB2011bWindows2007二、实验内容及问题背景1、实验内容为了验证MATLAB回归分析的合理性与正确性,为了考察变形量在测量时随时间变化x的影响,利用自回归移动平均模型ARMA，选择一点在相同时间段内，观察并记录其变形量h。待

2、解决实例如下：2、实验内容所涉及的问题背景现代变形监测对生态民生有着极其重要的作用，在一些地区塌陷地从无到有,并且增长速度很大,如何治理塌陷地是地区下一步土地复垦工作的重点。林地面积的减少,对保持区域的生态平衡提出了严峻的挑战。因此通过变形监测可以预测其变化趋势。三、实验原理与数学模型1.自回归移动平均模型ARMA设K是零均值的序列值，则：称为序列的自回归移动平均模型，简记ARMA(p，q)。其中p，q为正整数，分别为自回归、移动平均阶次(P(1～p)、廿(1～q)分别为模型的自回归系数和移动平均系数为白噪声序列，也称残差

3、序列，为独立正态同分布序列口，2.ARMA线性模型的概率特性对于零均值化序列，自协方差函数，自相关函数。对于一个有限样本数据，利用样本数据计算自协方差函数为，其中k=1，2，3，．．．，N一1。对于非零均值序列，由于理论均值往往未知，通常用样本均值近似代替，这时的样本自协方差用下式计算：利用自协方差函数的计算公式可得MA(q)模型的自协方差函数。当k>q时，，那么，MA(q)模型的自相关函数，在k>q后全为o，表现为截尾性，“截尾”处的k值就是MA模型的阶。但对于AR、ARMA模型没有此性。四、实验过程记录（含基本步骤、程

4、序清单及异常情况记录等）某一变形点的变形监测数据表一周期变形量周期变形量周期变形量周期变形量周期变形量周期变形量10.19110.53210.36310.71410.85510.7420.27120.53220.54320.65420.66520.8230.25130.47230.6330.44430.4740.18140.54240.37340.53440.3550.32150.29250.49350.75450.6360.34160.41260.35360.42460.6970.28170.56270.37370.63

5、470.7680.28180.28280.5380.71480.6490.25190.38290.5390.46490.4100.44200.55300.59400.75500.66将数据导入MATLAB>>x=A(:,1);y=A(:,2);>>plot(x,y)则其变形的线形图如图2.ARMA模型类型的判定以及阶数的确定>>y=A(1:44,2);>>[ACF,lags,bounds]=autocorr(y);>>[ACF,lags,bounds]=parcorr(y);>>plot(autocorr(y))>>plo

6、t(parcorr(y))自相关函数和偏相关函数均具有拖尾性，所以可以判定是符合ARMA模型的。并且可以初步判定其阶数为p=4q=3。>>m=armax(y,[43])Discrete-timeIDPOLYmodel:A(q)y(t)=C(q)e(t)A(q)=1+0.3373q^-1-0.7136q^-2-0.5233q^-3-0.1489q^-4C(q)=1+0.8961q^-1-0.8413q^-2-0.878q^-3EstimatedusingARMAXondatasetyLossfunction0.andFPE0

7、.Samplinginterval:13、用ARMA模型进行预测分析>>predict(A(45:52,2),m)>>y2=predict(A(45:52,2),m)y2=[8x1double]>>g=y2{1,1}周期4546474849505152实测值0.63000.69000.76000.64000.40000.66000.74000.8200预测值0.63000.59940.68550.74090.54870.70840.68120.7500五、实验结果实验结果见上图，（蓝色为预测值曲线，绿色为实测值曲线）可看出

8、预测曲线与实测曲线相近，同时也可以看到随着预测次数的增加精度也有所减小。六、深入探索与思考对于观测数据的预测应该有一定的次数限制，同时，对于精度要求高的监测项目MATLAB所进行的预测只可作为监测数据的参考，不可应用于实际记录。七、实验总结实验中，有关的参量应区别对应，同时应该注意字符的属性。在最后的成