货运数学建模论文分解.doc

托运货物获利的研究摘要 本文运用线性规划，一元线性以及三次回归预测等方法，综合运用Excel,SPSS等数学软件，解决了托运货物最大获益，未来客户申请量的预测问题以及未来每天最大收益的问题。 对于问题一，首先根据客户申请量，考虑到车辆自身载重、体积的约束和问题提出的约束条件，建立线性规划模型，运用Excel软件求解，最终结果为公司全部批复客户各类物品的申请量， 并得到最大获利为44275元，同时给出了合理的车辆承载方式。问题二，首先用Excel绘制出一个月客户申请量的折线图，观察其变化趋势，发现其波动性较大不便于求出其符合的函数关系，因此考虑到用此日之前所有申请量之和作为当天数据，绘图并通过观察发现可以分别用线性、二次和三次拟合，得到其拟合系数，并由此求出未来七天的申请量的预测值。问题三中，要求估算这7天的收益。根据第二问中预测的7天的申请量，并运用第一问中建立的模型，求出7天的批复量，进一步求出公司的获益情况。 在文章最后对模型进行评价及其适用范围做出推广，在实际应用中具有较大的参考价值。关键字：线性规划 回归预测法 Excel 规划求解 SPSS时间序列 ARIMA预测一、问题重述。某公司拥有3辆汽车，每辆载重量均为9000kg，可载体积为10.3，该公司为客户从甲地托运货物到乙地，收取一定费用。托运货物可分为四类：A、鲜活类 B、禽苗类 C、服装类 D、其他类，四类货物可以实现任意混装。托运手续是客户首先向公司提出托运申请，公司给予批复，客户根据批复量交货给公司托运。申请量与批复量均以公斤为单位，例如客户申请量为1000kg，公司根据平均每类每kg所占体积和相应托运单价，批复量可以为01000kg内的任意整数，若取0则表示拒绝客户的申请。问题1、如果某天客户申请量为：A 类 6200kg，B类 5500kg，C 类 4000kg，D 类3000kg，如果要求C类货物占用的体积不能超过B、D两类体积之和的三倍 (注意：仅在问题1中作此要求)。问公司应如何批复，才能使得公司获利最大？问题2、每天各类货物的申请总量是随机量，为了获取更大收益，需要对将来的申请总量进行预测。现有一个月的数据（见附件一），请预测其后7天内，每天各类货物申请量大约是多少？问题3、一般，客户的申请是在一周前随机出现的，各类申请单立即批复，批复后即不能更改，并且不能将拒绝量（即申请量减批复量）累计到以后的申请量。请根据你对下周7天中各类货物申请量的预测，估算这7天的收益各为多少？二、问题分析2.1问题一的分析 问题一要求批复的方法能使公司获利最大，这是一个优化问题。因此可以根据题目所给数据及其约束条件，建立合理的线性规划模型。2.2问题二的分析问题二中为了使公司获取更大利益，要求利用现有一个月的数据，对将来7天每天的申请总量进行预测。根据题意，每天各类货物的申请总量是随机量，而且需要预测的数据是未来短期内的，可以根据时间与货物申请量之间的发展趋势通过数据分析，从而建立合适的预测模型。2.3问题三的分析 问题三在问题二预测的基础上再对其进行最优化求解，在已知七天的预测值的情况下，重复运用问题一的方法，改变约束条件，便可求解到七天的最优解，使得收益最大。三、模型假设（1） 假设各物品混装后总体积等于混装前体积。（2） 无论如何混装均可保证货物的完整性。（3） 每天每辆车只运送一次且当日往返。（4） 不存在油费及路费问题（即只要运量不变，无论用几辆车其成本不变）。（5） 所有数据真实可信。四、符号说明序号符号符号说明12345678第i种货物每千克的体积第i种货物每千克的托运单价第t天的申请量前t天的申请总量A类产品前t天的申请总量B类产品前t天的申请总量C类产品前t天的申请总量D类产品前t天的申请总量五、模型的建立及求解5.1问题一的模型及求解考虑到第一问问题中的约束条件，C类货物占用的体积不能超过B、D两类体积之和的三倍,车的最大载重量及其最大容积，建立线性规划模型如下： 表示第i种货物被j辆货车运输，表示第i种货物的申请量。应用excel软件进行模型求解1，结果如下： 列1A类B类C类D类总计最优批复量（kg）620055004000300018700第一辆的承载量(kg)60000030009000第二辆的承载量(kg)200550078606486第三辆的承载量(kg)00321403214第一辆的承载体积(m3)10.2第二辆的承载体积(m3)10.298第三辆的承载体积(m3)9.642公司获得最大利润（元）44275结果表明，要使公司获得最大利润，则把申请的货物全部批复，并且第一辆车装A类产品6000kg,B,C产品不装，D产品装3000kg;第二辆车A、B、C、D类产品分别装200kg,5500kg,786kg和0kg;第三辆车只装C产品3214kg。共获益44275元。进一步分析可得，在保持最大收益的同时,第一和第二辆车均接近满车, ,第三辆车是没有装满的,是有剩余空间的,所以四类产品比不一定要按照表格的安排来分配货物，正常情况下，如果一个公司同时将三辆车安排出去送货，为了减少风险，通常将货物同时均分给三辆车去运送。所以，公司应该根据实际情况将货物分配给三辆车。5.2问题二的模型及求解 对于问题二，首先检查所给数据，确定无异常后，利用excel画出其变化的趋势图（如图一）。 图一05000100001500020000250001357911131517192123252729A类B类C类D类总计sin 根据图一，并不能看出其发展趋势存在特定的函数关系。因为每天各类货物的申请总量都是随机量，因此考虑将某一天之前的审批量都相加，作为这一天的数据，再利用excel绘制图像（如图二）。 根据曲线趋势判断，可以利用拟合直线或二次,三次方程法建立预测模型2。（1）设一次线性模型为：相关概念： 离差：离差和：离差平方和：因为，且有则（2）三次多项式曲线预测模型为：设统计的数据为，令即解此四元一次方程可求出系数。 假设各类产品各自的第天前申请的总量为,则其第t天的申请量 。利用SPSS软件对各类货物曲线分别进行分析：70000.0060000.0050000.0040000.0030000.0020000.0010000.000.00302520151050日期LinearObservedA类货物申请总量由其散点图可以看出,有比较明显的线性关系,所以用SPSS对其进行线性回归得到上图,并且得到相关分析:由Model Summary表格里的相关系数R=0.998,而调整决定系数R2=0.997,且回归估计的标准差S=1092.824可知模型的拟合效果很理想。从输出的第二个表格:ANOVA方差分析表可以看出回归平方和=33439429.回归方程的显著性检验中,统计量F=9179.320,对应的置信水平为0.000,远低于置信水平0.05,所以可认为显著性水平极高3。由第三个表格:Cofficients回归系数分析表,可以得出回归方程的常数项,回归系数。由此可知回归方程:则预测七天的A类货物申请量(取整处理)为：A类拟合1234567拟合结果71087732967550477712799218213084339预测值2209220922082208220922092209对每天总量取整后得到每天的预测值,因为回归系数不为整,所以最后经过总量取整的预测值并不总是等于其增长幅度。B类货物申请总量 对B类拟合的结果如图所示,根据SPSS输出的表格可以分析得(详见附录),用三次(Cubic)拟合的效果最好,并得其拟合方程为： 并预测七天B类货物申请量为：B类拟合1234567拟合结果105520109219113004116882120858124939129132预测值3620369937853878397640814193 140000.00120000.00100000.0080000.0060000.0040000.0020000.000.00302520151050日期CubicObservedC类货物申请总量对C类产品的总量用SPSS分析后,可得如图所示三次回归,观察图形以及由SPSS输出的表格可得其拟合方程为： 并得预测七天C类货物申请量为：C类拟合1234567拟合结果135970142843150155157932166203174993184329预测值6462687373127777827187909336100000.0080000.0060000.0040000.0020000.000.00302520151050日期LinearObservedD类货物的申请总量对D类产品的总量用SPSS分析后,可得如图所示三次回归,观察图形以及由SPSS输出的表格可得其拟合方程为： 并且预测七天D类货物申请量为：D类拟合1234567拟合结果9405397132100211103290106369109448112526预测值30793079307930793079307930785.3问题三的模型及求解 根据第一问已经建立的模型,第三问的求解只是将第一问的约束条件更改了一部分,将每天的申请量与第一问的特殊条件修改即可。得到新的约束条件: 通过Excel线性规划求解可得七天内最大的收益如下表所示:时间车辆分配及批复情况A类B类C类D类获得收益（元）第一天装载分配批复量220936206462307944427.78车1装载220936208600车2装载0034330车3装载0021693079第二天装载分配批复量2209369968732451.645752.34车1装载220936998230车2装载0034330车3装载0026172451.6第三天装载分配批复量220837857311101646307.27车1装载220837854451016车2装载0034330车3装载0034330第四天装载分配批复量178338787776046748.6车1装载178338789100车2装载0034330车3装载0034330第五天装载分配批复量43439768270046898.8车1装载434397614040车2装载0034330车3装载0034330第六天装载分配批复量040828394046957.5车1装载0408215280车2装载0034330车3装载0034330第七天装载分配批复量041948342046975.5车1装载0419414760车2装载0034330车3装载0034330由表格可知,从第四天起,随着一些产品申请量的增加,公司已经开始拒绝一些请求了。同时随着申请量的进一步增大,第六天到第七天获得的收益已经相同,这个时候,公司的运输水平决定了公司的收益。因此建议公司准备购买或租赁一批运输用车辆,增强公司的运营能力。观察表格,可以发现,四类产品的批复优先度不一样,经过计算其单位体积获得收益,可得到如下逻辑关系式: 因此，在受运输力的影响的情况下，优先满足B类产品的需求，优先拒绝D类产品的申请，再尽可能地满足C类及A类产品的需求，以获得最大的收益。六、模型改进在第二问中,采取的线性模型使得预测值每天都相等,而二次或三次模型可以更加精确地拟合，同时也避免了造成每天的申请量相等和数据不具有很强的波动性的问题。但本论文所建立的方程模型,都对预测值存在递增的趋势,与实际情况的反复落差可能会有出入。查阅资料可知,在相关问题的求解中,有一种方法是根据时间序列中的ARIMA模型或指数平滑等模型来建立。有SPSS软件可直接求得其结果，具体预测参看附录。其中ARIMA模型解决了销量随时间的随机性,具有一定的可参考性。七、模型评价与推广虽然本文中所建的线性模型在短时间内对每天具体的预测可能并不是十分精准,但是如果题目是要求对长时间的预测,例如求解一年的总的各类产品申请量,则精确度会大大提高。同时,对B,C二类产品的三次模型随不能长期预测,但对于短期预测特别适用,能够体现数据的变动。本文对A,B,C,D类产品的申请量根据其实际情况采取了不同的模型对其进行拟合,与现实生活中每类产品有自己不同的需求情况比较契合。另外，此类模型可以进一步推广至股市,在短期类,股市是波动比较大的一组数据,但在一定时间内,例如经济平稳发展的时候,是具有一定的线性规律。也可推广至其他方面，如蔬菜价格的预测等。八、参考文献1 郑丽敏，Excel数据处理与分析，北京：人民邮电出版社，2011.262-2742 百度文库，数学建模统计预测方法及预测模型，3 丁正生，概率论与数理统计简明教程，北京：高等教育出版社，2005.139-156九、附录B类产品SPSS拟合输出表格:线性分析: Model SummaryRR SquareAdjusted R SquareStd. Error of the Estimate.999.998.9981326.825二次分析:三次分析:C类分析结果:D类分析结果:线性分析:二次分析:三次分析:根据表格给出的相关系数R,决定系数R2,回归估计的标准差S,以及置信水平综合考虑可得出A,B,C,D类更适合的模型。Excel处理线性规划如图所示:SPSS时间序列ARIMA得到的第二问数据:时间A类预测B类预测C类预测D类预测12283.443431.4194897.8433018.87322535.573638.3054912.2243164.26531251.3333669.3244031.1033740.20343169.9343068.7034431.9133705.85151173.9142987.3873948.4874012.30463301.2823572.7244675.0183544.44971283.033420.9653955.8793334.98283392.7293857.3364174.3353563.90291447.5172199.1785428.1533854.096103062.7554058.2296399.8154169.319111438.3263436.9614740.6253204.134123084.6033060.0535024.1273671.796131631.6513686.1574186.3882155.561142921.5773441.2594388.0512387.678151681.6213278.224576.1482945.939162892.6643371.3774204.2673159.585171716.0373707.4383843.1573268.681182892.4533128.1423703.5343237.604191568.8313602.9174776.1813025.425202916.1243209.7333642.5593004.66211732.3323932.774886.3792943.9222834.0173677.0484027.9443513.285231690.9723913.113776.383138.306242843.9243813.4162950.2042860.457251757.9743476.0022914.7753286.123262881.523327.1253355.1473445.925271743.0913400.4774260.0822402.893282825.513523.0946334.4461847.353291675.0123025.5865631.542512.503302857.6453179.8798101.7642742.292311848.983321.8827879.1692901.49322714.623460.9347103.982945.277331855.5153423.4666530.3522972.73342708.1343433.5626105.8752989.943351861.9523430.8425791.7693000.735362701.7463431.5755559.3353007.501371868.2923431.3775387.3373011.743