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

ID:22859780

大?。?98.50 KB

頁(yè)數(shù):10頁(yè)

時(shí)間:2018-11-01

matlab時(shí)間序列分析在測(cè)繪中的應(yīng)用_第1頁(yè)
matlab時(shí)間序列分析在測(cè)繪中的應(yīng)用_第2頁(yè)
matlab時(shí)間序列分析在測(cè)繪中的應(yīng)用_第3頁(yè)
matlab時(shí)間序列分析在測(cè)繪中的應(yīng)用_第4頁(yè)
matlab時(shí)間序列分析在測(cè)繪中的應(yīng)用_第5頁(yè)
資源描述:

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

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

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

3、模型的自回歸系數(shù)和移動(dòng)平均系數(shù)為白噪聲序列,也稱殘差序列,為獨(dú)立正態(tài)同分布序列口,2.ARMA線性模型的概率特性對(duì)于零均值化序列,自協(xié)方差函數(shù),自相關(guān)函數(shù)。對(duì)于一個(gè)有限樣本數(shù)據(jù),利用樣本數(shù)據(jù)計(jì)算自協(xié)方差函數(shù)為,其中k=1,2,3,...,N一1。對(duì)于非零均值序列,由于理論均值往往未知,通常用樣本均值近似代替,這時(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值就是M

4、A模型的階。但對(duì)于AR、ARMA模型沒(méi)有此性。四、實(shí)驗(yàn)過(guò)程記錄(含基本步驟、程序清單及異常情況記錄等)某一變形點(diǎn)的變形監(jiān)測(cè)數(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.63

5、60.34160.41260.35360.42460.6970.28170.56270.37370.63470.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]=au

6、tocorr(y);>>[ACF,lags,bounds]=parcorr(y);>>plot(autocorr(y))>>plot(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.8

7、961q^-1-0.8413q^-2-0.878q^-3EstimatedusingARMAXondatasetyLossfunction0.0109749andFPE0.0148161Samplinginterval:13、用ARMA模型進(jìn)行預(yù)測(cè)分析>>predict(A(45:52,2),m)>>y2=predict(A(45:52,2),m)y2=[8x1double]>>g=y2{1,1}周期4546474849505152實(shí)測(cè)值0.63000.69000.76000.64000.4000

8、0.66000.74000.8200預(yù)測(cè)值0.63000.59940.68550.74090.54870.70840.68120.7500五、實(shí)驗(yàn)結(jié)果實(shí)驗(yàn)結(jié)果見(jiàn)上圖,(藍(lán)色為預(yù)測(cè)值曲線,綠色為實(shí)測(cè)值曲線)可看出預(yù)測(cè)曲線與實(shí)測(cè)曲線相近,同時(shí)也可以看到隨著預(yù)測(cè)次數(shù)的增加精度也有所減小。六、深入探索與思考對(duì)于觀測(cè)數(shù)據(jù)的預(yù)測(cè)應(yīng)該有一定的次數(shù)限制,同時(shí),對(duì)于精度要求高的監(jiān)測(cè)項(xiàng)目MATLAB所進(jìn)行的預(yù)測(cè)只可作為監(jiān)測(cè)數(shù)據(jù)的參考,不可應(yīng)用于實(shí)際記錄。七、實(shí)驗(yàn)總結(jié)實(shí)驗(yàn)中,有關(guān)的參量應(yīng)區(qū)別對(duì)應(yīng),同時(shí)應(yīng)

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫的文件,查看預(yù)覽時(shí)可能會(huì)顯示錯(cuò)亂或異常,文件下載后無(wú)此問(wèn)題,請(qǐng)放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫(kù)負(fù)責(zé)整理代發(fā)布。如果您對(duì)本文檔版權(quán)有爭(zhēng)議請(qǐng)及時(shí)聯(lián)系客服。
3. 下載前請(qǐng)仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無(wú)法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請(qǐng)聯(lián)系客服處理。
关闭