第1页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第1页共8页2011数学建模竞赛c题程序x=[566632745755769789985111013131428178219202150229226013149433851455809624168547656877210007113741256714332166141922822844264042968832074];>>y=diff(x);>>plot(y)>>a=lpc(y,3)a=1.0000-1.23370.23020.1088>>estx=filter([0-a(2:end)],1,y);>>e=y-estx;>>subplot(122);>>plot(lags,acs)2011数学建模竞赛c题（1）数据初步处理：对一般的时间序列而言为得到较为平稳的序列都会对原来的序列进行取对和差分处理并010203040-50005001000150020002500300035004000-40-2002040-0.200.20.40.60.811.2预测误差的自相关函数第2页共8页第1页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第2页共8页对新的数列进行白噪声检验cdata=[566,632,745,755,769,789,985,1110,1313,1428,1782,1920,2150,2292,2601,3149,4338,5145,5809,6241,6854,7656,8772,10007,11374,12567,14332,16614,19228,22844,26404,29688,32074];>>x=log(cdata);>>plot(x)>>x1=diff(x);>>subplot(221)>>plot(x)>>subplot(222)>>plot(x1)第3页共8页第2页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第3页共8页a=ar(x2,2)Discrete-timeIDPOLYmodel:A(q)y(t)=e(t)A(q)=1-0.2034q^-1-0.02065q^-2y2(1)=-0.0159;>>y2(2)=0.0383;>>fori=3:32y2(i)=0.2034*y2(i-1)+0.02065*y2(i-2);end>>plot(x2(1:32))>>holdon>>plot(y2)（2）真实值与可能值比较01020304000.511.522.533.5x1040102030406789101101020304000.050.10.150.20.250.30.35010203040-0.15-0.1-0.0500.050.10.150.205101520-0.500.51LagSampleAutocorrelationSampleAutocorrelationFunction(ACF)05101520-1-0.500.511.5LagSamplePartialAutocorrelationsSamplePartialAutocorrelationFunction第4页共8页第3页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第4页共8页、05101520253035-0.500.511.522.533.5x104第5页共8页第4页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第5页共8页0102030405060-505101520x104真实值估计值第6页共8页第5页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第6页共8页程序如下0510152056789101112-20-1001020-0.200.20.40.60.811.2-40-30-20-1001020304000.10.20.30.40.50.60.70.80.91预测误差自相关函数第7页共8页第6页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第7页共8页程序如下x=[566,632,745,755,769,789,985,1110,1313,1428,1782,1920,2150,2292,2601,3149,4338,5145,5809,6241,6854,7656,8772,10007,11374,12567,14332,16614,19228,22844,26404,29688,32074];>>y=log(x);x1=diff(y);a=1.0000-0.92600.03440.0157>>estx=filter([0-a(2:end)],1,x1);e=x1-estx;>>[acs,lags]=xcorr(e,'coeff');>>figure>>plot(x1,'-*')>>holdon>>plot(estx,'-r*')第8页共8页第7页共8页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第8页共8页0510152025303566.577.588.599.51010.5
