佛山科学技术学院 上 机 报 告 课程名称 数学建模 上机项目 牙膏销售量模型 专业班级 姓 名 学 号 一、问题提出 根据牙膏销售量与价格、广告费等表格1中的数据，建立三个模型，要求： 1）画出散点图：y对x1的散点图1；y对x2的散点图2； 2）确定回归模型系数，求解出教程中模型（3）； 3）对模型进行改进，确定回归模型系数，求解出教程中模型（5）； 4）对模型进一步改进，求解出教程中模型（10）。 二、问题分析 由于牙膏是生活必需品，对大多数顾客来说，在购买同类产品的牙膏时更多地会在意不同品牌之间的价格差异，而不是它们的价格本身。因此，在研究各个因素对销售量的影响时，用价格差代替公司销售价格和其他厂家平均价格更为合适。 三、模型假设 记牙膏销售量为y，其他厂家平均价格和公司销售价格之差（价格差）为x1，公司投入的广告费用为x2，其他厂家平均价格和公司销售价格分别为x3和x4，x1?x3?x4。基于上面的分析，我们仅利用x1和x2来建立y的预测模型。 四、模型建立 （显示模型函数的构造过程） 1）、为了大致地分析y与x1和x2的关系，首先利用表1的数据分别作出y对x1和x2的散点图。 建立程序如下： y=[ ]; x1=[ 0 0 ]; x2=[ ]; a=polyfit(x1,y,1); y1=polyval(a,x1); b=polyfit(x2,y,2); x3=::; y2=polyval(b,x3); subplot(2,1,1);plot(x1,y,'*',x1,y1,'b');title('í?1 y??x1μ?é￠μ?í?'); subplot(2,1,2);plot(x2,y,'o',x3,y2,'b');title('í?2 y??x2μ?é￠μ?í?') 从图1可以发现，随着x1的增加，y的值有比较明显的线性增长趋势，图中的直线是用线性模型 y??0??1x1??(1) 拟合的（其中?是随机变量）。而在图2中，当x2增大时，y有向上弯曲增加的趋势，图中的曲线是用二次函数模型 2y??0??1x2??2x2??(2) 拟合的。 综合上面的分析，结合模型（1）和（2）建立如下的回归模型 2y??0??1x1??2x2??3x2??(3) 2（3）式右端x1和x2称为回归变量（自变量），?0??1x1??2x2??3x2是给定价差x1，广告费用x2时，牙膏销售量y的平均值，其中的参数?0,?1,?2,?3称为回归系数，由表1的数据估计，影响y的其他因素作用都包含在随机误差?中。如果模型选择合适，?应该大致服从均值为0的正态分布。 五、模型求解 （显示模型的求解方法、步骤及运算程序、结果） 2）、确定回归模型系数，求解出教程中模型（3）： 建立程序如下： x1=[ 0 0 ]'; x2=[ ]'; X=[ones(30,1) x1 x2 x2.^2]; Y=[ ]'; [b,bint,r,rint,stats]=regress(Y,X); b,bint,stats 结果如下： b = bint = stats = 表2模型（3）的计算结果 参数 参数估计值 参数置信区间 [，] [，] [,] [,] ?0 ?1 ?2 ?3 R2= F= p< s2= 结果分析：表2显示，R=指因变量y（销售量）的%可由模型确定，F值远远超过F检验的临界值，p远小于?，因而模型（3）从整体来看是可用的。 表?22的回归系数给出了模型（3）中?0,?1,?2,?3的估计值，即????0?17.3244,?1?1.3070,?2??3.6956,?3?0.3486。检查他们的置信区间发现，只有?2的置信区间包含零点（但区间右端点距零点很近），表明回归变量x2（对因变量y的影响）不是太显著，但由于x2是显著的，我们仍将变量x2保留在模型中。 2六、模型改进 3）对模型进行改进，确定回归模型系数，求解出教程中模型（5）： 模型（3）中回归变量x1和x2对因变量y的影响是相互独立的，即牙膏销售量y的均值与广告费用x2的二次关系由回归系数?2和?3确定，而不依赖于价格差x1，同样，y的均值与x1的线性关系由回归系数?1确定，而不依赖于x2。根据直觉和经验可以猜想，x1和x2之间的交互作用会对y有影响，不妨简单地用x1，x2的乘积代表它们的相互作用，于是将模型（3）增加一项，得到 2y??0??1x1??2x2??3x2??4x1x2??(5) 2在这个模型中，y的均值与x2的二次关系为(?2??4x1)x2??3x2，由系数?2,?3,?4确定，并依赖于价格x1。 建立程序如下： x1=[ 0 0 ]'; x2=[ ]'; X=[ones(30,1) x1 x2 x2.^2 x1.*x2]; Y=[ ]'; [b,bint,r,rint,stats]=regress(Y,X); b,bint,stats 结果如下： b = bint =
数学建模作业8牙膏销售量模型
