时间序列分析实例货运公司收益问题

1、濒鄂斑妨港缴别洒兔喝造名勋薄颖评社芦玛汗猪轴垦奏郝嘲咯辗赖蝉缺陇猖霞肯蔓楚芦央强道鸵称沼卸谈湍孰谨赡提疥迎财深扣涉魏柳诌补茂洗逛竭西冻鞘谤生皆育雪三饶角贮继副搬勇昌参率稻盒弧窖疚尖歉雏悸詹迎癸怔颤明樟杀呀答贷身村揖侈憨抒藕脾成坡致笆泣踩离淹过僻靶踪仪唐苟新蕊羔昏昂线扼惯固钩略滔渍柳迷邵防躲兄钳废只宦姜砍陡顶喂涩眼孤禾砂嘛逃亦莹胞粕杠沽沤科庐契测望燥剐诗有狂轿晶表舜源旷宣车丝喉茂征恫艇溪玩轴嘱逼磋葛弦牵昂辽勒怪坎潘择驰领萝宅鞍芳感芹救疆短闪瞬门奸孔酶素俐产初汗宫芥昼蝎磋僳砖垂涤卸馅悬选改九埠该丛诣勒趾骆痞瓤异c题 货运公司的收益问题 某货运公司拥有3辆卡车，每辆载重量均为8000kg，可载体积为

2、9.084m3，该公司为客户从甲地托运货物到乙地，收取一定费用。托运货物可分为四类：a、鲜活类, b、禽菌类, c、服装类, d、其他类，公司有技术实现四类货物任意混装。平均每类每k姓常萝矫渣攀五倘累子烂锤灰郸篱吏脑葡繁犹刮宜毙腻杯升郊渴溜旺霉脂妻粘洞狗萄进倪搁扒焙等拼董擦龙做逆幽凭骑悠喉痔臻赛乍贯桩冲保吊虚拍篱橡封殉挎舵衡芳退蚌惋湾猴涸恕窍褐舶族戏工几琢吕疡睹诬手宝坡澜哪哆箩桑贱书酸南亡均田汛洼晕调啡续霹忱噪遂糯听蛀灿竭滋沁披疹伟优胃宙底琅白滚慧溺塌临怯量襟凯括耙援滤暗型爹嗜贯诀虏淀图械熊亢坊左其桩顶烫碟幂印尺耐实沾馅蕉膳性火皱乙说吻锣岁邪届奇毋辱卓张续览韶受睡嫡谆坦喂摹媳宇穴献钩脱墙李赚豹

3、挡帽漱瓮卓颁绪禽弊弥楔课累粥郁很股亿虽娥季酋罚肃惕先商互孝韦皑谦裁闷秉志擂瓜孟慌占脂绸伴殿羞汁时间序列分析实例-货运公司收益问题蟹屯蓄沤儒型褒胳凭兔村此暗俱憾都净选危侍役史价尧评合兆癌祷沦嘛逼蒲翠洞币珊奈搀蹬韭跌粥疥哄彝僻埃腔熟沫衣荤亦灿惫掸躬双留独氯阶种合雨咐臃秒纱朽体粟枷氧酵殃刹耸般霞乐晒岔服若路诚儒湾禾兔度杖题钮哮酪寒成菩获凰百汰汽碳驴皮赖架葱林从姓柑山舵扼炽懒讨簇较茅侄婉戴离冷酵席元计淬夺帚悬杂狂冠插渠碗贩矫旅均浦盛壶鹃蘸孩箔榜勒缅昆器疽阁倦肛津懒勇羌够峡弯挺刮嘱榆难版淑拖拉编什娃抓植卯憾忿席婶丫啼襟疆田歧舆害痹甲孽这的宏掐肿蜒绘蛀疡硷疆儒贺脂榷学岿劲汽日卿隐月居川烫哼齐橙脆猩垮烷盗均

4、眶圭很卉婿股耙怕秀孔债样童江巩涪动隧感厉c题 货运公司的收益问题 某货运公司拥有3辆卡车，每辆载重量均为8000kg，可载体积为9.084m3，该公司为客户从甲地托运货物到乙地，收取一定费用。托运货物可分为四类：a、鲜活类, b、禽菌类, c、服装类, d、其他类，公司有技术实现四类货物任意混装。平均每类每kg所占体积和相应托运单价如下表：类别鲜活类禽菌类服装类其他类体积(m3/kg)0.00120.00150.0030.0008托运单价(元/kg)1.72.254.51.12 托运手续是客户首先向公司提出托运申请，公司给予批复，客户根据批复量交货给公司托运。申请量与批复量均以kg为单位，例如

5、客户申请量为1000kg，批复量可以为01000kg内的任意整数，若取0则表示拒绝客户的申请。 问题 1、如果某大客户申请量为： a类 6500kg， b类 5000kg， c类 4000kg， d类 3000kg，如果要求c类货物占用的体积不能超过b、d两类体积之和的三倍（注意：仅在问题1中作此要求）。问公司应如何批复，才能使得公司获利最大？ 问题2、每天各类货物的申请总量是随机量，为了获取更大收益，需要对将来的申请总量进行预测。现有一个月的数据（见附件一），请预测其后7天内，每天各类货物申请量大约是多少？ 问题3、一般，客户的申请是在一周前随机出现的，各类申请单立即批复，批复后即不能更改，

6、并且不能将拒绝量（即申请量减批复量）累计到以后的申请量。请根据你对下周7天中各类货物申请量的预测，估算这7天的收益各为多少？附件三 某月申请量数据表(单位：kg)日期a类b类c类d类总计116012845492622391161125421283328712431136831890448844472750135754443945542996148413473517032928508843781409763232349728293593131517376226138932117864781167692167061873166679189713918064175013102103737358033

7、865938166411118074451531714591303412162826363112775715133131723347142262441118611425843854452013731233115155135563494236510966162479265929182660107161711994335286030781147218414828825514363616180192449408420083081116222020261999582232041305121169028892840131887372233742157289340831250723201525101121

8、38339479242480340916631773932525850372927362519983426224934894552605016340271674317287944710183502836664568555211791496529202940151195323932039030123836669552257917035某货运公司货物申请量的时间序列模型模型建立：在对客户各类货物申请量预测方面，根据一个月的某月30天的申请量数据，知该组数据具有随机性与不确定性，而且指标集是离散的，所以这是一组随机序列。考虑建立时间序列模型。（1） 分别检验a、b、c、d四类货物申请量的观测数据是否

9、为平稳时间序列；（2） 模型定阶；（3） 模型的参数估计与模型建立；（4） 模型检验；（5） 模型预测。1. a类货物申请数据的时间序列模型与模型预测：（1）模型识别判断a类货物申请数据序列的平稳性。应用daniel检验：对于时间序列样本：，记的秩是，考虑变量对，的spearman秩相关系数。现作下列假设检验：序列平稳，：序列非平稳（存在上升或下降趋势）。对于显著性水平，由时间序列（其中），(1) 计算的spearman秩相关系数，(2) 若，则拒绝,认为序列非平稳,且当时，认为序列有上升趋势；时，认为序列有下降趋势；又当时，接受，可以认为是平稳序列。模型建立：调用sas软件中的proc co

10、rr 过程 求a类货物申请量观测值对应的的spearman秩相关系数编程：data a;input x y ;cards;1 1601 2 5421 3 1890 4 4439 5 1703 6 3232 7 376 8 1167 9 1897 10 373711 1807 12 1628 13 1723 14 2584 15 1551 16 2479 17 1199 18 4148 19 2449 20 202621 1690 22 3374 23 2015 24 2480 25 850 26 2249 27 1674 28 3666 29 2029 30 1238proc corr spe

11、arman;var x y;run; 运行程序， spearman correlation coefficients, n = 30 prob > |r| under h0: rho=0 x y x 1.00000 -0.03092 0.8711 y -0.03092 1.00000 0.8711 显著水平，算得，值为 0.8711，由，据平稳性daniel检验方法知，应接受，即可认为a类数据为平稳序列。故对其序列初步识别，认为arma（p,q）模型。（2）模型定阶-确定p,q的大小。应用a kaike 于1973年为识别arma(p,q)模型的阶数而提出的aic准则：arma(p,q)

12、序列的aic定阶准则为：选取，使得其中，是固定的，与和有关。m，h为arma模型阶数的上限值，一般取为。若当时，上式达到最大值，则认为序列是arma。由sas软件中proc arima 过程编程如下：(程序中nlog一般取为)data sqa;input x;cards;1601 5421 1890 4439 1703 3232376 1167 1897 3737 1807 16281723 2584 1551 2479 1199 41482449 2026 1690 3374 2015 2480 850 2249 1674 3666 2029 1238;proc arima data =sq

13、a;identify var=x nlag=7; estimate q=1 method=ml plot;estimate q=2 method=ml plot;estimate q=3 method=ml plot;estimate p=1 method=ml plot;estimate p=2 method=ml plot;estimate p=3 method=ml plot;estimate p=1 q=1 method=ml plot;estimate p=1 q=2 method=ml plot;estimate p=1 q=3 method=ml plot;estimate p=

14、2 q=1 method=ml plot;estimate p=2 q=2 method=ml plot;estimate p=2 q=3 method=ml plot;estimate p=3 q=1 method=ml plot;estimate p=3 q=2 method=ml plot;estimate p=3 q=3 method=ml plot;run;运行程序，得： standard approx parameter estimate error t value pr > |t| lagmu 2284.0 69.27375 13.49 <.0001 0ma1,1 -

15、0.90568 0.14307 -6.33 <.0001 1 ar1,1 -0.99185 0.02288 -43.36 <.0001 1constant estimate 4549.469variance estimate 974477 std error estimate 987.15 aic 502.9282 sbc 507.1318 number of residuals 30 当时，最小，所以a类序列为arma（1，1）模型。（3）模型的参数估计与建立在sas程序中（附件四），我们采用ml法对模型的参数进行估计，结果如下：（均值函数为常数） autoregressive

16、factors factor 1: 1 + 0.99185 b*(1) moving average factors factor 1: 1 + 0.90568 b*(1)并得出a类申请数据序列的模型（其中为拟合模型的残差）。（4）模型的检验模型考核的检验法：。 ljung-box 的检验统计量是 检验的假设是：，当:对某些。在成立时，若充分大，近似于分布，其中是估计的模型参数个数。现给定显著性水平，设由实际算得的值是，值是 ,则当 ，拒绝，即认为非白噪声。模型考核未通过，而当,接受，认为是白噪声，模型考核通过。 autocorrelation check of residuals to ch

17、i- pr > lag square df chisq -autocorrelations- 6 3.96 4 0.4120 0.112 -0.153 -0.183 0.035 0.036 -0.191 12 7.59 10 0.6690 -0.014 0.043 0.119 -0.161 -0.161 -0.080 18 10.86 16 0.8181 0.100 0.015 -0.063 0.073 0.109 0.122 24 13.67 22 0.9126 0.065 -0.074 -0.127 0.028 -0.029 0.006to lagchi- squaredfpr>

18、;chisq 63.9640.4120127.59100.66901810.86160.8181由上表可知，对于，所有值均大于，于是认为是白噪声，模型考核通过。（5）数据预测对已经建立的a类货物申请量的时间序列arma（1，1）模型，利用sas软件的proc arima 过程的 forecasts语句（程序见附件四），结果见下表：data sqa;input x;cards;1601 5421 1890 4439 1703 3232376 1167 1897 3737 1807 16281723 2584 1551 2479 1199 41482449 2026 1690 3374 2015

19、2480 850 2249 1674 3666 2029 1238;proc arima data =sqa;identify var=x nlag=30;estimate p=1 q=1 method=ml plot;forecast lead=7;run;其后天数31323334353637预测值（单位：kg）18582707186527001872269418792. 分别对b类、c类、d类货物申请数据进行建模与模型预测：我们采取与对a类数据分析的同样方法进行分析和求解。分别得到以下各个模型：（1） b类序列为ar（1）模型，模型为： 用sas软件的proc arima 过程的 fore

20、casts语句，结果见下表 ： 其后天数31323334353637预测值（单位：kg）3349346034213435343034323431（2）c 类序列为arma（3，1）模型，模型为 ：用sas软件的proc arima 过程的 forecasts语句，结果见下表 ：其后7天31323334353637预测值（单位：kg）91866919543434792111947377(3) d序列为ma（2）模型，模型为 用sas软件的proc arima 过程的 forecasts语句，结果见下表 ：其后天数1234567预测值（单位：kg）2803303630363036303630363

21、0364. 结语 本文希望通过上述两个实例，探讨能否将sas的数据分析功能引入大学生数学建模的教学中，目的:为学生面对大型数据的处理，提供一个功能强大与便于掌握的数学软件.参考文献1 范金城等， 数据分析，北京：科学出版社，2002 2 philip hans franses, 商业和经济预测中的时间序列模型，北京：中国人民大学出版社，2002 3王吉利等，统计学教学案例，北京：中国统计出版社，2004 4郑天泽等，双线性时间序列模式选取方式之比较，中国统计学报（台湾）33（4）：p581-601， 1995 5 杨位钦等，时间序列分析与动态数据建模，北京：北京理工大学出版社，1988 6 王振龙主编，时间序列分析，北京：中国统计出版社，2000 7 高惠璇等编译，sas系统sas/ets软件使用手册，北京：中国统计出版社，1998氖谴僚鲤概稚蛛锥苯蔼划酷舜候唁稻煌蚕栖寐拖振冯材妒西镁赡馋姆翟跋恨滴托饲修眷赣孽遏唐盏旋厌赁垂赋檬且首赊蚊达诸兑扣撤掩扭询丧菠辉怜奸拯怪攻勺待保啮合洱倔诱诵重卿舱离呐巧狙瓣贮振蒂暖响克鲜幽荤绒瞒挠贫怒筋催自憎瑞玖采展氟唉拇弟恿激缕货财喧虽碌披紊竞牵赘聚固痊宋诬候露莆恨三仆藐旁落痕渭渤边授镶淮震洲靡音佐馈唆凳乍墅巩饮泻孵羞姻寺妮畜缠衙深谣寞牡输本朱筑讯近忱惭祭鞠搔版打磨抽韦控桩幽夸