matlab時(shí)間序列分析在測繪中的應(yīng)用

《matlab時(shí)間序列分析在測繪中的應(yīng)用》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。

1、《現(xiàn)代測繪數(shù)據(jù)處理方法》課程課間實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)項(xiàng)目：MATLAB時(shí)間序列分析在測繪中的應(yīng)用班級(jí)：測繪工程專業(yè)指導(dǎo)教師：一、實(shí)驗(yàn)?zāi)康募八密浖姹?、實(shí)驗(yàn)?zāi)康蘑倭私釳ATLAB時(shí)間序列分析的基本原理及應(yīng)用②學(xué)會(huì)用MATLAB時(shí)間序列的分析方法解決測繪工程中的實(shí)際問題2、實(shí)驗(yàn)軟件所用版本MATLAB2011bWindows2007二、實(shí)驗(yàn)內(nèi)容及問題背景1、實(shí)驗(yàn)內(nèi)容為了驗(yàn)證MATLAB回歸分析的合理性與正確性,為了考察變形量在測量時(shí)隨時(shí)間變化x的影響,利用自回歸移動(dòng)平均模型ARMA，選擇一點(diǎn)在相同時(shí)間段內(nèi)，觀察并記錄其變形量h。待

2、解決實(shí)例如下：2、實(shí)驗(yàn)內(nèi)容所涉及的問題背景現(xiàn)代變形監(jiān)測對生態(tài)民生有著極其重要的作用，在一些地區(qū)塌陷地從無到有,并且增長速度很大,如何治理塌陷地是地區(qū)下一步土地復(fù)墾工作的重點(diǎn)。林地面積的減少,對保持區(qū)域的生態(tài)平衡提出了嚴(yán)峻的挑戰(zhàn)。因此通過變形監(jiān)測可以預(yù)測其變化趨勢。三、實(shí)驗(yàn)原理與數(shù)學(xué)模型1.自回歸移動(dòng)平均模型ARMA設(shè)K是零均值的序列值，則：稱為序列的自回歸移動(dòng)平均模型，簡記ARMA(p，q)。其中p，q為正整數(shù)，分別為自回歸、移動(dòng)平均階次(P(1～p)、廿(1～q)分別為模型的自回歸系數(shù)和移動(dòng)平均系數(shù)為白噪聲序列，也稱殘差

3、序列，為獨(dú)立正態(tài)同分布序列口，2.ARMA線性模型的概率特性對于零均值化序列，自協(xié)方差函數(shù)，自相關(guān)函數(shù)。對于一個(gè)有限樣本數(shù)據(jù)，利用樣本數(shù)據(jù)計(jì)算自協(xié)方差函數(shù)為，其中k=1，2，3，．．．，N一1。對于非零均值序列，由于理論均值往往未知，通常用樣本均值近似代替，這時(shí)的樣本自協(xié)方差用下式計(jì)算：利用自協(xié)方差函數(shù)的計(jì)算公式可得MA(q)模型的自協(xié)方差函數(shù)。當(dāng)k>q時(shí)，，那么，MA(q)模型的自相關(guān)函數(shù)，在k>q后全為o，表現(xiàn)為截尾性，“截尾”處的k值就是MA模型的階。但對于AR、ARMA模型沒有此性。四、實(shí)驗(yàn)過程記錄（含基本步驟、程

4、序清單及異常情況記錄等）某一變形點(diǎn)的變形監(jiān)測數(shù)據(jù)表一周期變形量周期變形量周期變形量周期變形量周期變形量周期變形量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將數(shù)據(jù)導(dǎo)入MATLAB>>x=A(:,1);y=A(:,2);>>plot(x,y)則其變形的線形圖如圖2.ARMA模型類型的判定以及階數(shù)的確定>>y=A(1:44,2);>>[ACF,lags,bounds]=autocorr(y);>>[ACF,lags,bounds]=parcorr(y);>>plot(autocorr(y))>>plo

6、t(parcorr(y))自相關(guān)函數(shù)和偏相關(guān)函數(shù)均具有拖尾性，所以可以判定是符合ARMA模型的。并且可以初步判定其階數(shù)為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模型進(jìn)行預(yù)測分析>>predict(A(45:52,2),m)>>y2=predict(A(45:52,2),m)y2=[8x1double]>>g=y2{1,1}周期4546474849505152實(shí)測值0.63000.69000.76000.64000.40000.66000.74000.8200預(yù)測值0.63000.59940.68550.74090.54870.70840.68120.7500五、實(shí)驗(yàn)結(jié)果實(shí)驗(yàn)結(jié)果見上圖，（藍(lán)色為預(yù)測值曲線，綠色為實(shí)測值曲線）可看出

8、預(yù)測曲線與實(shí)測曲線相近，同時(shí)也可以看到隨著預(yù)測次數(shù)的增加精度也有所減小。六、深入探索與思考對于觀測數(shù)據(jù)的預(yù)測應(yīng)該有一定的次數(shù)限制，同時(shí)，對于精度要求高的監(jiān)測項(xiàng)目MATLAB所進(jìn)行的預(yù)測只可作為監(jiān)測數(shù)據(jù)的參考，不可應(yīng)用于實(shí)際記錄。七、實(shí)驗(yàn)總結(jié)實(shí)驗(yàn)中，有關(guān)的參量應(yīng)區(qū)別對應(yīng)，同時(shí)應(yīng)該注意字符的屬性。在最后的成