第三个到第八十九个数
dataf1expdata41反对数
rmse1sumdata13391dataf132891计算均方根误差
data4Ycdata2p1e
d1反差分第三个到第八十九个数datafexpdata4反对数
Udataf1Tdata1p1TpsigmaU1TUUaiclogdetsigmaU2pKKTifaicrecordaic
recordaicaicmpYc_bYcdataf_bdataf
e
de
dsubplot211plot1921mYc_b3b1921mdata33m191r差分数据画图subplot212plot1921mdata13m191b1921mdataf_b3r实际数据画图rmsesumdata13m191dataf3291m1计算均方根误差
运行结果如下：
f实验分析：效果图如下图所示上面的子图是对数差分数据预测效果，下面的子图是反差分反对数的预测效果图，预测均方根误差为可以看出，预测效果很好。
例题2：课件ARMA_AR_MA_Represe
tGra
gerA
alysisP34的例题
fARp模型转换成MA∞模型。代码如下：fu
ctio
psiARcha
geMAAR_fi
um_psi输入AR_fi和
um_psi，求MA模型中的psipsipsizeAR_fi2AR_fi的列数ifp0
errordlgmodelis
otcorrectede
difp0
ksizeAR_fi1AR_fi的行数ppk求pDeyekIeyekOzeroskkp1JIO求Jifp1
xiAR_fiJIelsexiAR_fieyekp1zeroskk求xixi为kp×kp阶矩阵e
d判断是否平稳lamtaeigxifori1pkifabslamtai1判断此模型成立条件lamta为xi的特征值，lamta1z
retur
e
de
d求psipsieyekfori1
um_psi1psipsiJxiiJe
de
d运行结果：
f例题3：卡尔曼滤波
仿真中所选的模型的状态方程是AR2模型，基于老师所给数据。
代码如下：clcclearallHdiag111000Gzeros66vphai_1phai_2Bzeros66Fdiago
es16B1313phai_1B1346phai_2B4613diago
es13Origi
al_datazeros923Origi
al_data1180179185192211202207214231229234237206250259263264280282292286302304307317314306304292275273301280289303322315339364371375432453460475496494498526519516531573551538532558524525519526510519538549570559584611597603619635658675700692759782816844830853852833860870830801824831830Origi
al_data24514654854935095205215405485585745835915996106276426536606947097347517637667798087857947997998128378538768979229499799881025106311041131113711781211125612901314134613
f851416143614621493151615571613164216901759175617801807183118731897191019431976201820402070212121322199225322762318236924232457247025212545258026202639261826282651Origi
al_data341542143444845945847948749751051652552953854655557457458660261763965366r