牙膏销售量.docx
《牙膏销售量.docx》由会员分享，可在线阅读，更多相关《牙膏销售量.docx（10页珍藏版）》请在冰点文库上搜索。
牙膏销售量
实验报告
学院：
理学院系：
数学系专业：
应用数学实验名称：
计算实验
年级：
姓名：
学号：
组
号：
实验时间：
2010.12.1成绩：
指
导教师签字：
实验题目：
实验四回归分析与逐步回归
目的和要求：
1.学会对实际问题进行数学抽象。
2.熟悉用MINITAB软件中进行方差分析和回归分析的基本命令的操作；
3.学会用MINITAB进行多因素试验的回归分析的程序设计。
主要仪器设备：
多媒体计算机
实验内容：
牙膏的销售量
问题：
某大型牙膏制造企业为了更好的拓展产品市场，有效地管理库存，公司董事会要求销售部门根据市场调查，找出公司生产的牙膏销售量价格，广告投入等之间的关系，从而预测出在不同价格和广告费用下的销售量。
为此，销售部的研究人员收集了过去30个销售周期（每个销售周期为4周）公司生产的牙膏销量，销售价格，投入的广告费用，以及同期其他厂家生产的同类牙膏的平均销售价格，见表1。
试根据这些数据建立一个数学模型，分析牙膏销售量与其他因素的关系，为制定价格策略和广告投入策略提供数据依据。
销售
公司销售价
其他厂家平］广告费用
价格差
销售量
格（元）
均价格（元）
（百万兀）
（元）
（百万支）
1
3.85
3.80
5.50
-0.05
7.38
2
3.75
4.00
6.75
0.25
8.51
3
3.70
4.30
7.25
0.60
9.52
4
3.70
3.70
5.50
0
7.50
5
3.60
3.85
7.00
0.25
9.33
6
3.60
3.80
6.50
0.20
8.28
7
3.60
3.75
6.75
0.15
8.75
8
3.80
3.85
5.25
0.05
7.87
9
3.80
3.65
5.25
-0.15
7.10
10
3.85
4.00
6.00
0.15
8.00
11
3.90
4.10
6.50
0.20
7.89
12
3.90
4.00
6.25
0.10
8.15
13
3.70
4.10
7.00
0.40
9.10
14
3.75
4.20
6.90
0.45
8.86
15
3.75
4.10
6.80
0.35
8.90
16
3.80
4.10
6.80
0.30
8.87
17
3.70
4.20
7.10
0.50
9.26
18
3.80
4.30
7.00
0.50
9.00
19
3.70
4.10
6.80
0.40
8.75
20
3.80
3.75
6.50
-0.05
7.95
21
3.80
3.75
6.25
-0.05
7.65
22
3.75
3.65
6.00
-0.10
7.27
23
3.70
3.90
6.50
0.20
8.00
24
3.55
3.65
7.00
0.10
8.50
25
3.60
4.10
6.80
0.50
8.75
26
3.65
4.25
6.80
0.60
9.21
27
3.70
3.65
6.50
-0.05
8.27
28
3.75
3.75
5.75
0
7.67
29
3.80
3.85
5.80
0.05
7.93
30
3.70
4.25
6.80
0.55
9.26
表1牙膏销售量与销售价格，广告费用等数据
（其中价格差指其他厂家平均价格与公司销售价格之差）
1.实验设计方案
1）前期分析：
牙膏的销售量受到多种因素影响，例如：
产品销售价格、同类产
品销售价格、广告费用投入量、产品质量等因素。
但是我们只考虑两
个对结果有显著性影响的因素，广告费用投入量及同类价格产品。
在
考虑同类产品价格时不好处理，在这里我们仅考虑其他产品同本公司产品的价格差
2）模型假设：
1在一定时期内假设市场总需求量没有太大的变化。
2同类产品在一定时期内价格无明显变化。
3通过调节本公司的价格调整都能够达到理想的价格差
3）建立模型：
y=订「必「2X2「3X22■；
4）编写程序：
［b,bint,r,rint,stats］二regress（y,x,alpha）
5）对结果进行分析，讨论诸如：
结果的合理性、正确性，算法
的收敛性，模型的适用性和通用性，算法效率与误差等。
2.基本模型
为了大致分析y与xi和X2的关系，首先利用散点图观察销售量y与价格差xi及y与广告投入量X2之间的关系。
丫与x1的关系：
模型
y£讥1；
（1）
丫与x2的关系:
图2y对X2的散点图
拟合的（其中；是随机误差）在图2中，当X2增大时，y有向上弯曲增加的趋势，图中的曲线用二次函数模型：
y=札
（2）
拟合。
综上分析，结合模型
（1）和
（2）建立如下回归模型
y=\-ix<2X23x|（3）
其中，y是建立的模型，我们用讨二'-0■-^23x22对y进行估计，
其中-0,-1,-2,-3是我们待估计的参数。
3.模型求解
利用MATLA统计工具箱中的命令regress求解，使用格式为：
[b,bint,r,rint,stats]=regress（y,x,alpha）
具体实验操作过程如下：
x1=[-0.050.250.600.250.20.150.05-0.150.15
0.20.10.40.450.350.30.50.50.4-0.05-0.05-0.10.2
0.10.50.6-0.0500.050.55]';
x2=[5.56.757.255.576.56.755.255.2566.5
6.2576.96.86.87.176.86.56.2566.576.86.86.5
5.755.86.8]';
x3=[30.2545.562552.562530.254942.2545.562527.562527.5625
3642.2539.06254947.6146.2446.2450.414946.24
42.2539.06253642.254946.2446.2442.2533.062533.64
46.24]';
Y=[7.388.519.527.59.338.288.757.877.18
7.898.159.18.868.98.879.2698.757.95
7.657.2788.58.759.218.277.677.93
9.26]';
X=[ones（30,1）x1x2x3];
[b,bint,r,rint,stats]=regress（Y,X）
rcoplot（r,rint）
结果如下：
b=
17.3244
1.3070
-3.6956
0.3486
bint=
5.728228.9206
0.68291.9311
-7.49890.1077
0.03790.6594
r=
-0.0988
-0.0795
-0.1195
-0.0441
0.4660
-0.0133
0.2912
0.2735
-0.2351
0.1031
-0.4033
0.1747
0.0400
-0.1504
0.1284
0.1637
-0.0527
-0.1907
-0.0870
-0.0165
-0.1292
-0.3002
-0.2933
-0.1679
-0.2177
0.1116
0.3035
0.0693
0.2474
0.2270
rint=
-0.52700.3294
-0.53090.3718
-0.51060.2716
-0.47310.3848
0.08130.8507
-0.46090.4343
-0.13740.7197
-0.08700.6340
-0.59600.1258
-0.32800.5341
-0.81900.0125
-0.26180.6112
-0.40320.4832
-0.59330.2925
-0.32070.5775
-0.28410.6116
-0.48300.3776
-0.62480.2434
-0.53480.3609
-0.44230.4092
-0.56090.3024
-0.71810.1177
-0.72430.1377
-0.55480.2190
-0.64490.2095
-0.29940.5226
-0.10370.7106
-0.37140.5099
-0.18070.6755
-0.18900.6430
stats=
0.905482.94090.00000.0490
得到模型⑶的回归系数的估计值及其置信区间（置信水平:
=0.05）、检
验统计量r2,f,p的结果见下表1
参数
参数估计
参数置信区间
17.3244
[5.728228.9206]
1.3070
[0.68291.9311]
-3.6956
[-7.49890.1077]
0.3486
[0.03790.6594]
R2=0.9054F=82.9409p<0.0001
4.结果分析
表1
据矩阵［1XiX2X22］,alpha为置信水平：
•，缺省时〉=0.05;输出b为一：
的估计值，常记作1，bint为b的置信区间，r为残差向量，rint为r的
置信区间，stats为回归模型的检验统计量，有3个值回归方程决定系数R2、F统计量、F统计量对应的概率值P。
由表中的数据显示，R2=0.9054指因变量的y的90.54%可由模型确定，F值远远超过F检验的临界值，p远远小于：
•，因而模型（3）可用。
表2的回归系数给出了模型（3）中-0，'1，-2，飞的估计值
0=17.3244，「=1.3070，2=-3.6956，-3=0.3486。
检查它们的置信区
间发现，只有■-2的置信区间包含零点（但区间右端点距零点很近），表明回
归变量X2（对因变量y的影响）不是太显著的，但由于x；是显著的，我们
仍将变量x2保留在模型中。
5.销售量预测
经回归系数的估计值代入模型（3），即可预测公司未来某个销售周期牙膏
的销售量y，将预测值记为y，得到模型（3）的预测方程：
y=r」1花「2X2「3X；（4）
只需知道该销售周期的价格差X1和投入的广告费用X2，就可以计算预测值
A
y。
公司无法直接确定价格差X,，只能制定公司的牙膏销售价格人，但是其它厂家的平均价格一般可以通过根据市场情况及原材料的价格变化等估计。
模型中用
价格差做为回归变量的好处在于公司可以更灵活地来预测产品的销售量或市场需求量，因为其它厂家的平均价格不是公司所能控制的。
预测时只要调整公司的
牙膏销售价格达到设定的回归变量价格差X1的值。
回归模型的一个重要应用是，对于给定的回归变量的取值，可以以一定的
置信度预测因变量的取值范围，即预测区间。
6.模型改进
模型（3）中回归变量X1，X2对因变量y的影响是相互独立的，即牙膏销售
量y的均值和广告费用X2的二次关系由回归系数：
2，-3确定，而不依赖与价格
差Xi，同样，y的均值与Xi的线性关系由回归系数M确定，不依赖于X2。
根据经验可参想，Xi和X2之间的交互作用会对y有影响，简单的用Xi，X2的乘积代表他
（5）
们的交互作用，将模型（3）增加一项，得到:
y=0:
必亠』2x2•：
3X■"■;,4XiX2-；
在这个模型中，y的均值与x2的二次关系为2x^3x|4x1x2，由系数：
2，
-3，V确定，并依赖与价格差Xio
讨论：
对实验中存在的问题、进一步的想法等进行讨论。
升级会员 