时间序列分析报告实验的报告材料(3)

《时间序列分析报告实验的报告材料(3)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实用标准文案《时间序列分析》课程实验报告项目名称：非平稳序列确定性分析组员姓名：李菲指导教师：牛宪华完成日期：2013年4月20日精彩文档实用标准文案一、上机练习（P124）1.拟合线性趋势12.7914.0212.9218.2721.2218.8125.7326.2726.7528.7331.7133.95程序：dataxiti1;inputx@@;t=_n_;cards;12.7914.0212.9218.2721.2218.8125.7326.2726.7528.7331.7133.95;procgplotdata=xiti1;plot

2、x*t;symbolc=redv=stari=join;run;procautoregdata=xiti1;modelx=t;outputpredicted=xhatout=out;run;procgplotdata=out;plotx*t=1xhat*t=2/overlay;symbol2c=greenv=stari=join;run;运行结果：精彩文档实用标准文案分析：上图为该序列的时序图，可以看出其具有明显的线性递增趋势，故使用线性模型进行拟合：xt=a+bt+It,t=1,2,3,…,12分析：上图为拟合模型的参数估计值，其中a=9.

3、7086,b=1.9829，它们的检验P值均小于0.0001，即小于显著性水平0.05，拒绝原假设，故其参数均显著。从而所拟合模型为：xt=9.7086+1.9829t.精彩文档实用标准文案分析：上图中绿色的线段为线性趋势拟合线，可以看出其与原数据基本吻合。2.拟合非线性趋势1.857.4814.2923.0237.4274.27140.72265.81528.231040.272064.254113.738212.2116405.95程序：dataxiti2;inputx@@;t=_n_;cards;1.857.4814.2923.0237

4、.4274.27140.72265.81528.231040.272064.254113.738212.2116405.95;procgplotdata=xiti2;plotx*t;symbolc=redv=stari=none;run;procnlinmethod=gauss;modelx=a*b**t;parametersa=0.1b=1.1;der.a=b**t;der.b=a*t*b**(t-1);outputpredicted=xhout=out;run;procgplotdata=out;plotx*t=1xh*t=2/overl

5、ay;精彩文档实用标准文案symbol2c=greenv=nonei=join;run;运行结果：分析：上图为该时间序列的时序图，可以很明显的看出其基本是呈指数函数趋势慢慢递增的，故我们可以选择指数型模型进行非线性拟合：xt=abt+It,t=1,2,3,…,12精彩文档实用标准文案分析：由上图可得该拟合模型为：xt=1.0309*1.9958t+It分析：图中的红色星号为原序列值，绿色的曲线为拟合后的拟合曲线，可以看出原序列值与拟合值基本上是重合的，故该拟合效果是很好的。3.X—11过程407774177843160458974194744

6、0614437847237433154339644843468354283343548446374710742552435264503947940437404500746667493254487846234470555031846354472604888352605485275023751592551525045152294546335880253990554775785061978程序：dataxiti3;inputx@@;t=intnx('quarter','1jan1978'd,_n_-1);formattyyq4.;cards;407

7、774177843160458974194744061443784723743315433964484346835428334354844637471074255243526450394794043740450074666749325精彩文档实用标准文案4487846234470555031846354472604888352605485275023751592551525045152294546335880253990554775785061978;procgplotdata=xiti3;plotx*t;symbolc=redv=stari

8、=join;run;procx11data=xiti3;quarterlydate=t;varx;outputout=outb1=xd10=seasond11=ad