数学建模：高压锅的销售量

1、数学建模辽宁工程技术大学数 学 建 模 课 程 成 绩 评 定 表学 期2014-2015学年1学期姓 名 专 业班 级课程名称 数学建模论文题目 高压锅的销售量评定标准评定指标分值得分知识创新性20理论正确性20内容难易性15结合实际性10知识掌握程度15书写规范性10工作量10总成绩100评语:任课教师时 间2014年11月17 日备 注高压锅的销售量摘要Logistic增长曲线模型和Gompertz增长曲线模型是计量经济学等学科中的两个常用模型，可以用来拟合销售量的增长趋势。本文运用以上两种曲线模型研究了某地区高压锅的销售量的变化规律，并根据给定的1981年到1993年高压锅销售量的数据

2、，运用趋势分析法，分别建立指数增长模型、Logistic增长曲线模型和Gompertz增长曲线模型来对高压锅的销售量进行预测分析，并对各模型进行了比较分析。关键字:曲线 增长率 趋势 MATLAB1 模型的背景问题描述趋势分析法又叫比较分析法、水平分析法，它是通过对财务报表中各类相关数字资料，将两期或多期连续的相同指标或比率进行定基对比和环比对比，得出它们的增减变动方向、数额和幅度，以揭示企业财务状况、经营情况和现金流量变化趋势的一种分析方法。趋势分析法在定量预测。趋势分析法又可称为趋势曲线分析、曲线拟合或曲线回归。它是根据已有的历史数据资料来拟合一条曲线，使得这条曲线能够反映出研究对象本身的

3、增长趋势，然后按增长趋势曲线，对要求的未来的某一点进行估计，预测出该点该时刻的研究对象的预测值。能够正确并掌握认识销售量的变化规律，才能建立正确的销售量预测模型，做出准确的预报，为企业未来发展方向，企业发展定位，商品生产作出有利预测。下表为某地区1981年到1993年间高压锅的销售量列表（单位：万台）。表1-1： 高压锅的销售量（单位：万台）年份ty年份ty1981043.65198871238.7519821109.86198981560.0019832187.21199091824.2919843312.671991102199.0019854496.581992112438.891986

4、5707.651993122737.7119876960.25本文根据表1-1高压锅的销售量数据，分析销售量的变化规律，得到该地区高压锅的销售量的变化趋势的拟合曲线，建立销售量的模型，通过建立的销售量模型，得到模拟曲线，根据得到的拟合曲线，并对该地区的高压锅销售量进行预测。通过预测对该地区高压锅生产企业作出指导。当地企业可以根据预报的高压锅销售量，对高压锅生产量进行控制，避免了因为市场的盲从效应而造成的损失。2 基本假设1、产品的销售不受人为因素影响。2、高压锅的销售量随时间连续变化。3、任一单位时刻，高压锅的增长量与当时的高压锅总量成正比。4、销售量的增长律短时间内是不变的。5、在处理数据、

5、拟合曲线，得到模拟曲线的过程，都不考虑随机误差。6、在一段时间里，销售市场是平稳发展的，在高压锅的销售过程中，市场对高压锅的需求量稳定。7、假设说明：严格的说，讨论销售量所建立的模型属于离散型模型，但在销售量基数很大的情况下，突然地增加或减少的只是单一的个体或少数几个个体数，相对于全体数量而言，这种改变量是极其微小的，因此，销售量可以看作是可随时间的连续变化，这样，就可以采用微分方程的工具来研究这一问题。3 问题分析本文要求根据某地的1981年到1993年间高压锅销售量的数据，建立高压锅的销售量模型。根据 表1-1：高压锅的销售量（单位：万台）提供的数据，以时间t为衡轴，销售量y为纵轴，建立销

6、售时间t与销售量的关系图，如图表3-1：图表 3-1高压锅的销售量05001000150020002500300002468101214时间t销量y根据图表3-1，显然，高压锅的销售量随时间的变化呈指数增长。产品的销售平稳，经济发展的稳定，销售市场的平稳发展，高压锅在销售过程中，市场的外部环境总体稳定，在这一假设下，高压锅的销售量是随着时间的连续变化而变化的，也就是说，高压锅的销售数量是连续变化的，但总体分析可以得出增长量与当时的高压锅总量成正比，销售量的增长律是不变的。据此，建立高压锅的销售量指数增长模型。 符号说明 固有高压锅销售量增长率，即： 时段的高压锅销售量数 高压锅的最大销售量，显

7、然有4 模型建立4.1 建立高压锅指数增长模型假设商品是自然销售的，即不受人为因素影响，记时刻t的销售量为x(t)，在销售量基数很大的情况下，突然地增加或减少的只是单一的个体或少数几个个体数，相对于全体数量而言，这种改变量是极其微小的，可以忽略不计。即销售量可以看作是可随时间连续变化的，将x(t)视为连续、可微函数。记初始时刻(t=0)的销售量为，假设销售量增长率为常数r，即单位时间内的x(t)的增量等于r乘以x(t)。考虑到时间内高压锅销售量的增量，有令，得到满足微分方程 (1)由这个方程解出： （2）（2）式的参数r和可以用表1-1数据估计。为了利用最小二乘法，将（2）式取对数，可得 ，

8、（3）分别以1981年到1992年的数据和1981年到1993年的数据拟合（3）式，用matlab计算：t=0:11; x=43.65 109.86 187.21 312.67 496.58 707.65 960.25 1238.75 1560 1824.29 2199 2438.89;y=log(x); p=polyfit(t,y,1) r=p(1),x0=exp(p(2) Y=polyval(p,t); X=exp(Y); p = 0.3435 4.4914r = 0.3435x0 = 89.2424得到r0.3435，89.2424t=0:12; x=43.65 109.86 187.2

9、1 312.67 496.58 707.65 960.25 1238.75 1560 1824.29 2199 2438.89 2737.71;y=log(x); p=polyfit(t,y,1) r=p(1),x0=exp(p(2) Y=polyval(p,t); X=exp(Y);得到r0.3205，97.1060结果分析：用上面得到的参数r和代入（3）式，将结果与实际数据比较。X1是用1981年到1992年的数据拟合的结果，计算人口x2用的是全部数据的拟合的结果。表4-1年实际销售量计算销售量x1计算销售量想198143.6589.242497.10601982109.86125.820

10、5133.79421983187.21177.3910184.34391984312.67250.0989253.99211985496.58352.6078349.95451986707.65497.132348225700.8937664.346019881238.75988.1715915.346919891560.001393.21261.219901824.29196421737.719912199.00276932394.219922438.89390443298.819932737.71550474545根据表5-1的数据，用matlab制作指数增长型拟

11、合图形，图5-1.1表示用1981年到1992年的数据拟合的结果图，图4-1.2表示用全部数据（1981年到1993年）拟合的结果。图5-1.1、图5-1.2中曲线是计算结果，“*”表示实际数据。图表 4-1.2通过以上数据及图表，明显地，对于高压锅的销售量的分析指数模型并不适用，其原因是高压锅的销售是受市场的影响的，人的消费观念不会一成不变，购买欲是不随时间稳定变化的，由此可见，高压锅的销售量的增长律是不变的假设不成立。需重新建立模型分析。4.2 高压锅销售线性模型的假设由图表3-1，假设y和t满足线性关系，所以建立线性模型，设利用最小二乘法确定的具体值，并根据的值拟合高压锅的销售情况，与原

12、数据进行比较。MATLAB的程序实现如下：y=43.65 109.86 187.21 312.67 496.58 707.65 960.25 1238.75 1560.00 1824.29 2199.00 2438.89 2737.71;t=0:12;p=polyfit(t,y,1);yy=polyval(p,t,1) plot(t,y,'*',t,yy)画出原数据与拟合曲线，如图5-2图4-2线性模型的误差分析预测的标准误差为：.模型分析从图形及标准误差可以看出，线性模型虽然简单，但误差太大。并且当时，而高压锅销售量是有限的，也就是说高压锅的销售量是一个有限的数，不可能是一个

13、无限大的。所以，用线性模型不能完全反映高压锅的销售情况。必须寻找一个更好的模型去分析高压锅的销售情况。4.3 建立高压锅logistic模型由于高压锅在进入市场初期没有太多人家了解或清楚高压锅，致使人们对高压锅的购买数量不大，所以此时高压锅销售数量的增长率小；随着时间的推移，人们开始认识到使用高压锅的好处，高压锅的销售量增长率也逐渐提高；之后越来越多的人家都有高压锅了，而高压锅的经久耐用决定了高压锅销售量的增长率会逐渐减少；到最后该地区基本上所有人家都有高压锅了，高压锅的销售数量将趋于一个定值，即L，此时销售数量的增长率将趋于0.所以，综上所述，高压锅的销售情况满足Logistic模型。设高压

14、锅的销售数量的增长率为 r（t），高压锅的销售量的上限为L ， 销售量为 y（t）则有 ： . （4-1）. （4-2）建立模型. (4-3) 模型分析 当L与y(t)相比很大时，与ry相比可以忽略不计，Logistic模型可以转化为指数模型；而当L与y相比不是很大时，就不能忽略，其作用是使高压锅的销售量的增长速度减缓下来。作图结果，如图4-3图4-3对线性化Logistic增长模型并做非线性回归线性化Logistic增长曲线模型Logistic增长曲线模型为. (3-1)两边同时取倒数得 ，整理得. (3-2)两边同时取对数得 . (3-3)令，.即可得到一个线性关系式：. (3-4)所以L

15、ogistic增长曲线模型能线性化。用Matlab对Logistic模型做非线性回归Matlab程序如下：y=43.65 109.86 187.21 312.67 496.58 707.65 960.25 1238.75 1560.00 1824.29 2199.00 2438.89 2737.71;t=0:12;L=3000; y1=log(L./y-1); p=polyfit(t,y1,1);k=-p(1);a=exp(p(2);yy=L./(1+a*exp(-k*t);plot(t,y,'*',t,yy); 拟合Logistic模型，画出拟合图形，如图4-4 如图4-45

16、 拟合Gompertz模型线性化Gompertz模型Gompertz增长曲线模型为. (5-1)两边同时除以 L 得. (5-2)两边同时取对数得. (5-3)再两边同时取对数得. (5-4)令，. 得线性关系式. (5-5)拟合Gompertz模型用Matlab拟合Gompertz模型 画出Gompertz模型并与原数据比较，如图5-1 如图5-16 综合评价与结论两模型相同之处Logistic增长曲线模型，俗称“S曲线”，由Verhulst于1845年提出，当时主要目的是模拟人口的增长。其一般形式为. (5-1)也称狭义的逻辑增长曲线模型。增长曲线有两个重要特征。一是y随着t的增加直至+而

17、趋向于K，K即是Y的饱和值；反过来，当t-时，y0。二是增长曲线具有一个拐点，在拐点之前，y的增长速度越来越快；在拐点之后，y的增长速度越来越慢，逐渐趋近于0。在现实经济生活中，许多指标的增长过程具有这两个特征。例如，一种新产品、新技术的普及率，一种耐用品的存量，它们的增长过程都遵循逻辑增长曲线模型。所以，逻辑增长曲线模型在经济预测中有广泛的应用，是一种重要的预测模型。Gompertz曲线用于描述这样一类现象：初期增长缓慢，以后逐渐加快，当达到一定程度后，增长率又逐渐下降，最终接近一条水平线。Gompertz曲线通常用于描述事物的发展由萌芽、成长到饱和的周期过程。在现实生活中有许多现象符合Go

18、mpertz曲线形式，如工业生产的增长、产品的寿命周期和一定时期的人口增长等。可见Gompertz增长曲线与逻辑增长曲线相似，只是二者的拐点的位置不同。分析两种模型的拐点Logistic模型转化为微分形式是. (5-2)Gompertz模型转化为微分形式是. (5-3)在Logistic模型中dy/dt是关于y的二次函数，很容易看出，当y=L/2是增长率最大，即是模型的拐点。利用Matlab作图，去L=3000，很容易作出Logistic模型与Gompertz模型的拐点，估计处Gompertz模型的拐点在y=1015.2处，如图6-1 如图6-17 参考文献1姜启源，谢金星，叶俊，数学模型，北京：高等教育出版社，2005，13-216.2 中国人口增长预测模型3数学模型简介及其给予MATLAB的实现