牙膏销售量的数学模型课程设计

关于牙膏销售量的数学模型课程设计

————————————————————————————————作者：————————————————————————————————日期:

?

统计12-1 李本恩

一、 课题名称:牙膏销售量的影响因素 二、 课题条件

参考文献:：《ＭA ＴLAB 从入门到精通》 三、 设计任务

本文从收集有关牙膏的销售量开始，从牙膏销售量和价格、广告投入之间的关

系出 发，分别通过对这三个方面的深入研究从而制定出各自的最佳方案，最后再综合这 三个主要因素,进一步深入并细化，从而求得最优解。

四、论文内容

摘要内容 ：本文从收集有关牙膏的销售量开始，从牙膏销售量和价格、广告投入

之间的关系出发,分别通过对这三个方面的深入研究从而制定出各自的最佳方案， 最后再综合这三个主要因素,进一步深入并细化，从而求得最优解。

模块Ⅰ中，我们假设在x 1和x 2对y 的影响独立 ，从而得到了方程

εββββ++++=2

2322110x x x y

模块Ⅱ中,我们假设ｘ1和x ２对y 的影响有交互作用，进一步得到新的方程

2

0112232412y x x x x x βββββε=+++++

关键词：线性回归模型 相关系数

问题重述：某大型牙膏制造企业为了更好的拓展产品市场,有效地管理 库存，公司董事会要求销售部门根据市场调查,找出公司生产 的牙膏销售量价格，广告投入等之间的关系,从而预测出在不 同价格和广告费用下的销售量。为此,销售部的研究人员收集 了过去30个销售周期（每个销售周期为4周）公司生产的牙 膏销量，销售价格,投入的广告费用，以及同期其他厂家生产 的同类牙膏的平均销售价格,见表1。试根据这些数据建立一 个数学模型,分析牙膏销售量与其他因素的关系，为制定价格 策略和广告投入策略提供数据依据。

销售 周期 公司销售价格(元） 其他厂家平均价格(元） 广告费用（百万元）

价格差 (元) 销售量 （百万支) 1 2 3 4

3.８5 ３.75 3.70 3.70

3.８0 ４．00

4.30 3.70

５.50 6.75 7．25 5．50

-0.05 0．25 0.６０ 0

７.3８ 8.５1 9.52 7.50

5 ６7 8 ９１0

11

12

13

14 1５16 1７1８

19

20 ２1 22 ２3

24

25 ２6 27 ２８

29

30

3.60

3.60

3.６0

3.80

３.80

3.８5

3.90

3.９0

3．７0

3.75

3.75

3．８0

３．70

3.８0

３．７０

3．8０

3.80

3.7５

3.70

３.55

3.60

3．６5

３.70

3．75

３.8０

3.７０

３．85

3.８0

3.７5

3．8５

3．6５

4.0０

４.10

4.００

4．1０

4．20

4.1０

4.１0

４.20

4.３0

4.10

3.75

3．75

3.65

3.９0

3.65

４.10

４．25

3.65

３．75

3.85

４．２５

7.0０

６.5０

6.７5

5.25

5．2５

6.00

6.５0

6.25

7.０0

6.9０

６.８０

6.80

７．10

7.00

6.８０

６．50

6．25

6．0０

6.５0

7．00

6.80

6.80

6．50

5．7５

５.80

6.８0

0.2５

0.20

０．15

０.0５

－０．１5

0.15

0.２0

0．10

0.40

０.４5

0.35

0．30

0.５0

0．50

０.40

－0.０５

-0.0５

－０.10

0．20

0.１0

0．50

0．60

－0.０５

0.０5

０.55

9．３3

８.2８

８．75

7.８7

７.1０

8．０0

７.89

8．15

9.1０

8．86

８.90

8.87

9．26

9.00

8．７5

7.9５

7.６５

7.2７

8.０0

8.50

8．75

9．21

８.２7

7．67

7.93

９.２6 表1 牙膏销售量与销售价格,广告费用等数据

(其中价格差指其他厂家平均价格与公司销售价格之差）

1.实验设计方案

1) 前期分析:

由于牙膏是生活必需品,对大多数顾客来说,在购买同类产品的牙膏时更多的会在意不同品牌之间的价格差异,而不是它们的价格本身。因此,在研究各个因素对销售量的影响时，用价格差代替公司销售价格和其他厂家平均价格更为合适。

记牙膏销售量为y,其他厂家平均价格与公司销售价格之差（价格差）为x

1

公司投入的广告费用为ｘ2 ,其他厂家平均价格与公司销售价格分别为x 3和ｘ4 ，x 1=x 3-x ４。基于上面的分析,我们尽利用ｘ1和x 2来建立y 的预测模型。

２) 模型假设：

１ 在一定时期内假设市场总需求量没有太大的变化。 ２ 同类产品在一定时期内价格无明显变化。

3 通过调节本公司的价格调整都能够达到理想的价格差 3) 建立模型:εββββ++++=2

2322110x x x y

4) 编写程序:[b,bint,ｒ，r ｉnt,st ａt ｓ]=regres ｓ(y ，x,alpha ） 5) 对结果进行分析，讨论诸如:结果的合理性、正确性,算法的收敛性，模型的适用性和通用性，算法效率与误差等。

2.基本模型

为了大致分析ｙ与1x 和2x 的关系，首先利用散点图观察销售量y 与价格差1x 及ｙ与广告投入量2x 之间的关系。

ｙ与1x 的关系: ｙ 与2x 的关系:

图1 y 对1x 散点图（1) 图2 y 对2x 的散点图

从图（1)发现,随着1x 增加,y 的值有明显的线性增加趋势，图中直线用线性模型

011y x ββε=++ (1)

拟合的（其中ε是随机误差）在图2 中,当2x 增大时,y 有向上弯曲增加的趋势，图中的曲线用二次函数模型：

εβββ+++=2

23220x x y （2） 拟合。

综上分析,结合模型（1)和(2)建立如下回归模型

2

0112232

y x x x ββββε=++++ （３)

其中，y 是建立的模型,我们用=∧

y 2

2322110x x x ββββ+++对y 进行估计,其中

3210,,,ββββ是我们待估计的参数。

3．模型求解

利用ＭAT ＬA Ｂ统计工具箱中的命令regress 求解,使用格式为： [b,ｂｉnt ，ｒ，ｒin ｔ,st ａts]=reg ｒｅｓｓ(y,x,al ｐha)

得到模型(３)的回归系数的估计值及其置信区间（置信水平05.0=α）、检验统计量P F R ,,2，Ｓ2的结果见下表２ 参数

参数估计 参数置信区间 0β 17.３244 [5.72８2 28.9206] 1β

1．３070 [0.6８2９ 1.9311] 2β -3.６95６ [-7.４９8９ 0．107７] 3β

0.34８6

[0.03７9 0．6594]

2R =0.９０54 F =82．9４0９ p <0.０0０1 Ｓ2=0．0４9０

表 ２ ４.结果分析

由表中的数据显示,2R =0.90５4指因变量的y 的90.54%可由模型确定,F 值远远超过F 检验的临界值，p 远远小于α，因而模型（3)可用。

表2的回归系数给出了模型(3）中0β,1β,2β，3β的估计值

0β∧=1７.3244，1β∧=１.３０７０，2β∧=－３.6９５6，3β∧

＝0.3486。检

查它们的置信区间发现,只有2β的置信区间包含零点(但区间右端点距零点很近)，表明回归变量2x (对因变量y 的影响)不是太显著的,但由

于22x 是显著的，我们仍将变量2x 保留在模型中。

5.销售量预测

经回归系数的估计值代入模型(3),即可预测公司未来某个销售周期牙膏的销售量y ，将预测值记为y ∧

,得到模型(３）的预测方程:

y ∧

=2

0123122x x x ββββ∧∧∧∧

+++ (４） 只需知道该销售周期的价格差1x 和投入的广告费用2x ,就可以计算预测值

y ∧

。

公司无法直接确定价格差1x ，只能制定公司的牙膏销售价格4x ,但是其它厂家的平均价格一般可以通过根据市场情况及原材料的价格变化等估计。模型中用

价格差做为回归变量的好处在于公司可以更灵活地来预测产品的销售量或市场需求量，因为其它厂家的平均价格不是公司所能控制的。预测时只要调整公司的牙膏销售价格达到设定的回归变量价格差1x 的值。

回归模型的一个重要应用是,对于给定的回归变量的取值，可以以一定的置

信度预测因变量的取值范围，即预测区间。

6.模型改进

模型（3)中回归变量1x ,2x 对因变量y 的影响是相互独立的,即牙膏销售量y 的均值和广告费用2x 的二次关系由回归系数2β,3β确定,而不依赖与价格差1x ，同样，y 的均值与1x 的线性关系由回归系数1β确定，不依赖于2x 。根据经验可参想,1x 和2x 之间的交互作用会对y 有影响,简单的用1x ,2x 的乘积代表他们的交互作用,将模型(3）增加一项，得到:

20112232412y x x x x x βββββε=+++++ (５)

在这个模型中,y 的均值与2x 的二次关系为22232412x x x x βββ++，由系数

2β,3β,4β确定，并依赖与价格差1x 。

下面让我们用表１的数据估计模型(５）的系数。利用MATLAB 的统计工具箱得到的结果见表3.

验统计量P F R ,,2，Ｓ2的结果见下表１ 参数

参数估计 参数置信区间

0β ２９.1133 ［13．70１３ 44．５252]

1β

11.1342 [1.97７８ ２0.290６]

2β -7.６0８0 [－12.６932 -2.52２8] 3β

０．671２ [0.253８

1.0887]

β4

-1.4777

[-2.8518 -0．103７]

2R =０.９２0９ F ＝7２.7７７1 p <0．00０1 S 2=０.0.0426

表 ３

表3与表2的结果相比,Ｒ２有所提高，说明模型（5)比模型(3)有所改进，相信模型（5)更符合实际。

用模型（5）对公司的牙膏销售量做预测，仍设在某个销售周期中,维持产品的价格差X1=0.２元,并投入X2=6.5百万元的广告费用,则该周期牙膏销售量y 的估计值为=∧

y ∧

β

０

+ ∧β1ｘ1 ＋∧β２ｘ2+∧

β

３x 22 +∧

β4x 1x 2 =２９.１13３ + 11.134 ×

0．2 – 7.60８×6.5 + 0．6712 ×６.52 -１.4777 ×0．２ ×6.5 =8.3253百万支，置信度为95％的预测空间为[7.8953,8．7５9２],与模型（3)的结果相比,∧

y 略有增加，而预测区间长度短些。在保持广告费用x 2=6.5百万元不变的条件下,分别对模型（3)和（5）中牙膏销售量的均值∧

y 与价格差x 1的关系作图，见图３和图4

图３ 模型(3）∧

y 与ｘ1 的关系 图4 模型(５）∧

y 与ｘ1 的关系

在保持价格差ｘ1=０.2元不变的条件下，分别对模型(3）和(5）中牙膏销售量的均值∧

y 与广告费用x 2的关系作图,见图5和图6

图5 模型（3)∧

y 与2x 的关系 图6 模型(5）∧

y 与2x 的关系

可以看出,交互作用项1x 2x 加入模型,对∧

y 与1x 的关系稍有影响,而∧

y 与2x 的关系有较大变化，当2x < ６时∧y 出现下降,2x > ６以后∧

y 上升则快得多。

进一步讨论: 为了解1x 和2x 之间的相互作用,考察模型（5）的预测方程

∧

y ＝29.1133+11.134２1x －７.６0８02x +0.67122x ２ -1.４7771x 2x (6）

如果取价格差1x =0.1元，代入（6)可得

∧

y

1x =0.１ =30．2267

-7．７55８2x +0.67122x 2 (7） 再取1x =0．3元,代入（6）可得

∧

y

1x =0.3 =３2．453６ －8.05132x +0.67122

x 2

(8) 它们均为2x 的二次函数,其图形见图7，且

∧

y

1x =0.3 -

∧

y

1x =0．1

= 2.2269-0．２9552x (9)

由(9）式可得，当

2x ＜ ７.5360时，总有

∧

y

1x =０.3 ＞

∧

y

1x =０.1,即

若广告费用不超过大约7.5百万元，价格差定在０．3元时的销售量，比价格差定在0.1元的大 ，也就是说,这时的价格优势会使销售量增加。

图７

y与

2

x的关系（（7）与(８）的图形)

附录：源程序

ｘ１=[-0.０5；０．25;0.６0;0;0.25；０.20;0．15;０.05；－0.15;０.１5;０.20;0.１0；0.40;0.45；

0.３5;0.30；0.５0;０．50;0.40；-0.05；-0．0５;-0．10;0.20;0.10;０.5０;0.60;-0.05;0;0．０５;0.５５]；

y=［7.38;８.５１；9．52;7．5０;9.33;8.2８;8.75;7.8７;７．１0;8.00;7.89;8.１５；9.10;８.86;8．90;8.87;９.2６;９.０0;8．75；7.95；7.65；7.2７;8．00;8.50;8.75；9.２１;8.２7;7.６７;7．93；９．２6]；

aａ=ｐolｙfｉt(x1,y,１);

y1=polyval(aa,x1)；

plot（ｘ1,y1,x１,ｙ,＇ｒｏ')

%图1:Y对X1的散点图

x１=[-０.05;0．25;0.６0；0；0.25;0.２0；0．15；0．05;－０.15;０．１５;0.２0;0．10；0.４0;0.45;０．35；0．３0;0.50;0.50;0．40;-0．０5;－0.0５;－0.1０；0.2０;0.1０；0.50;0.60;-0．０5；0;0.05；0.55];

x2＝[5.50；６.７5;７.25;5．50;7.0０；6.50;6．７５;５．25;5．25;6.00；6．50;６.25;７.00;6.90;6．8０;6.80;７.10；7.00;6.８0;6．50;6．２5;６.00；6．５0；７．00;6.80；6.80;6.50;5．7５;5.80;6.８0];

ｙ＝[7.38;8.51;9.5２;7.５0;9.３3;8.28;８.75;7.8７;7．10;8.00；7.89;8.1５；９.1０；８.86;8．90；８．87;9．26;9.00；８.75;７.95；7.6５；７.2７;8.0０;8.50;8.75;9.２1;8.2７;7.67；７.93;9.２6］; ａa=polyｆit（x2,ｙ，２);

ｘ３=5.2５:0.05:７.２5;

y2=ｐolyval(ａa,x3);

pｌｏt（ｘ2，y,＇ｒo',x3，ｙ2)

%图２：Ｙ对X２的散点图d

x４=[onｅs（30，1),x1,x２，x2.^2];[b，biｎｔ,ｒ，rｉｎt,staｔs]=rｅgress(y,ｘ4）%表2

ｘ５=[oｎes(30,1),x1,x2,x2.^2,x1.*x２]；[ｂ,ｂint，r,rｉｎｔ,ｓtats］＝regrｅss（y，ｘ5) %表3

x1＝[-0.05;０.25;0.60；0;0．２5；0.20;0.１5;0.05;－0.１５;0.15;0.20;0.10；0.４0；0.４５;０．35;0.３0;0.５0;0.50；０.4０;-０．0５;-0.05;-０．1０；0.20；0.1０；0．50;0．６0；－0.05;0;0.05；

0.５5]；ytu３=17.３２44+1．307*x1+(-3.6９５6)*6．5＋０.３486*6.5*６.5；

plｏt(ｘ1，ｙtu3)

ｇrid oｎ%图3

x1＝[-0.05;0.25;0.６0；0;0.２5;0.20;0.１5；0.05；－0.15;0.15；0．20;0．10；０.40;0.45;０．35;0.30；０.５0；0.50;０.４0；－0.05;-0.０5;-0．1０;0.20;０.1０;0．５0;0．60;－0.05;0;０.０5；０.5５];

ytu4=29．1133+１1．13４2＊x1＋（-7.60８*６.5)+0.6７１２*６.5*6.5+(－１.477７）*6.５＊x１；

ｐｌoｔ（x１,ｙtu４)

ｇｒｉd on %图4

x2=[5.５0;６.7５;7.25;5.50;７.０0;6．50;6．7５；５.25;5.25；6.00；６.５0；6.25;7.０0;6.９0;６．80;6.80;7.10;７.00；6.８0；6.50；６.２5;6.００;6.５0；7.0０;6.8０;6.80;６.５0;５．75；5．80;6.８０]; ytu5=17.3２44+１.3０７*0．2+(-3.695６)*x2+0.３486*x2.*x2;

bｂ＝pｏｌyfit(x2，ytu5,2);

xtｕ5=5.25:0.05:７．2５;

ytu51=ｐolyｖａl(bb,ｘtu５);

pｌot(xtu5,ytu5１）

griｄoｎ%图5

ｘ２＝［5.50;6.７5;7．25;5.5０；7.０0;6.50;6.75;5.２5;5.25;6.00;6．5０;6.25;7．０0;6.９0;６．80;6.8０；7.10；7.00;６.８０;6.50；6.25；６.０0;6.50;7.0０;６．８0;6．80;６.50;5.７５;5．8０；6．8０];

ｙtu６＝2９.1１3３+１1.134２＊0.2+(－7.608*x2)+０.6712*x2.＊x2+(－1.4７7７)*ｘ2＊0．2; bｂ＝polyfit(x2,ytu６，2);

ｘｔｕ6=5.２5:0.05：7.25；

ytu6１=poｌyval(bb,xtｕ6)；

ploｔ（xｔu６,ytu６１)

gｒｉd on %图6

x２＝[５.５0;6.75;7．25;５．50;7.６5;6.50；6.7５;5．25;5．25;６．00；６.50；

6.2５;７.00;６.90；６.80;6.80;7．1０；7．00;6.80;６．50;6.２５;6．00;6.50;

7.00;6.8０;６.８０;６．5０；5.7５;5.８0；6．80];

xｔu7＝sort（x2);

yｔｕ7=3０.2267-７.75５8*xtu7+0.６712*xtu７.^2；

ｐlot(xtu７,ytｕ7)

grid on

hold on

ｙtｕ8=32.45３6-８.０５13*xｔu7＋0.6712*xtｕ７.^2;

pｌｏt（xtu7，ｙtu8)

hｏｌd ofｆ

%图7

经过讨论和模型的改进，可以预测出在不同价格和广告费用下的牙膏销售量。