24年五一杯数学建模竞赛一等奖附录.pdfVIP

24年五一杯数学建模竞赛一等奖附录

附录1：问题求解代码

代码1：问题1求解代码（SPSS.28）

GETDATA***加载数据

/TYPE=XLSX

/FILE=C:\Users\24769\Desktop\附件1.xlsx

/SHEET=nameSheet1

/CELLRANGE=FULL

/READNAMES=ON

/DATATYPEMINPERCENTAGE=95.0

/HIDDENIGNORE=YES.

EXECUTE.

DATASETNAME数据集2WINDOW=FRONT.

DATEY1970.

****定义时间序列

Thefollowingnewvariablesarebeingcreated:

NameLabel

YEAR_YEAR,notperiodic

DATE_Date.Format:YYYY

*序列图.

TSPLOTVARIABLES=

/ID=时间年

/NOLOG

/FORMATNOFILLREFERENCE.

***SARIMA模型.

PREDICTTHRUYEAR2031.

TSMODEL

/MODELSUMMARYPRINT=[MODELFITRESIDACFRESIDPACF]PLOT=[SRSQUARE]

/MODELSTATISTICSDISPLAY=YESMODELFIT=[SRSQUARERSQUARERMSENORMBIC]

/MODELDETAILSPRINT=[PARAMETERSRESIDACFRESIDPACF]PLOT=[RESIDACF

RESIDPACF]

/SERIESPLOTOBSERVEDFORECASTFIT

/OUTPUTFILTERDISPLAY=ALLMODELS

/SAVEPREDICTED(预测)

/AUXILIARYCILEVEL=95MAXACFLAGS=24

/MISSINGUSERMISSING=EXCLUDE

/MODELDEPENDENT=

PREFIX=模型

/EXPERTMODELERTYPE=[ARIMAEXSMOOTH]

/AUTOOUTLIERDETECT=ONTYPE=[ADDITIVELEVELSHIFTINNOVATIONALTRANSIENT

LOCALTREND

ADDITIVEPATCH].

1

***时间序列模型

PREDICTTHRUYEAR2031.

TSMODEL

/MODELSUMMARYPRINT=[MODELFITRESIDACFRESIDPACF]PLOT=[SRSQUARE]

/MODELSTATISTICSDISPLAY=YESMODELFIT=[SRSQUARERSQUARERMSENORMBIC]

/MODELDETAILSPRINT=[PARAMETERSRESIDACFRESIDPACF]PLOT=[RESIDACF

RESIDPACF]

/SERIESPLOTOBSERVEDFORECASTFIT

/OUTPUTFILTERDISPLAY=ALLMODELS

/SAVEPREDICTED(预测)

/AUXILIARYCILEVEL=95MAXACFLAGS=24

/MISSINGUSERMISSING=EXCLUDE

/MODELDEPENDENT=

PREFIX=模型

/EXPERTMODELERTYPE=[EXSMOOTH]

/AUTOOUTLIERDETECT=OFF.

代码2：问题2预测代码（Stata.16）

clear//清空变量

*(1)导入数据并生成时间序列(*和//都可以用来注释一行)

insheetusingdata1.csv//导入位于和代码同一文件夹下的csv数据文件