高教社杯全国大学生数学建模竞赛.docx
- 文档编号:27506796
- 上传时间:2023-07-02
- 格式:DOCX
- 页数:22
- 大小:227.84KB
高教社杯全国大学生数学建模竞赛.docx
《高教社杯全国大学生数学建模竞赛.docx》由会员分享,可在线阅读,更多相关《高教社杯全国大学生数学建模竞赛.docx(22页珍藏版)》请在冰豆网上搜索。
高教社杯全国大学生数学建模竞赛
2005高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们将受到严肃处理。
我们参赛的题目是:
B题:
DVD在线租赁
我们的参赛报名号为(如果赛区设置报名号的话):
03
所属学校(请填写完整的全名):
广西师范学院
参赛队员(打印并签名):
1.蒙智
2.陶国飞
3.张引俊
指导教师或指导教师组负责人(打印并签名):
陈建伟
日期:
2005年9月16日
赛区评阅编号(由赛区组委会评阅前进行编号):
2005高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
DVD在线租赁服务管理模型
摘要
本文建立对DVD在线租赁服务管理问题的研究模型,针对其中的问题分别进行研究:
针对问题一,我们运用概率统计中的条件概率公式,用期望去求表1中的每种DVD所需要准备的张数,以保证希望看到该DVD的会员中至少50%在一个月内能够看到该DVD,求得每种DVD需要准备的张数如表(1-1):
DVD名称
DVD1
DVD2
DVD3
DVD4
DVD5
需要准备的张数
6250
3125
1563
782
391
表(1-1)
保证在三个月内至少95%的会员能够看到该DVD,求得每种DVD需要准备的张数如表(1-2):
DVD名称
DVD1
DVD2
DVD3
DVD4
DVD5
需要准备的张数
3960
1980
990
495
391
表(1-2)
针对问题二,我们用excel进行表2的数据进行处理,根据会员的偏爱程度(用0…10的整数表示,数字越小表示偏爱程度越高,而数字“0”表示对应的DVD不在当前会员的在线订单中),为了研究方便我们把偏爱程度数“0”换成“11”,建立0——1线性规划模型,利用软件lingo8.0对DVD进行最优化分配,利用目标函数所求得的最小值来表达会员获得的最大的满意度,再根据lingo8.0所得的结果进行数据分析,列出前30位会员见表(2-1);
同样我们在问题二的基础上针对问题三,进一步改变模型的约束条件,在满足一个月内95%的会员得到他想看的DVD并且满意度最大的约束条件下,利用软件lingo8.0
求得目标函数的优化值为6000,并且得到DVD的分配的数据,对所得的数据进行统计得
到每种DVD购买量见表(3-1);
针对问题四,首先针对搜集到的数据运用软件Matlab6.5做出散点图,使用线性回归,研究DVD的需求预测、购买、分配以及一个月内的租赁次数等因素对DVD租赁的影响,对多个变量进行变换,使得非线性问题转换为线性问题并加以处理。
最后从网站经营管理人员出发,为会员提供更为周到的服务。
关键词:
最大满意度最优化解线性回归散点图
问题重述
顾客缴纳一定数量的月费成为会员,订购DVD租赁服务,会员对哪些DVD有兴趣,只要在线提交订单,网站就会通过快递的方式尽可能满足要求。
会员提交的订单包括多张DVD,这些DVD是基于其偏爱程度排序的。
网站会根据手头现有的DVD数量和会员的订单进行分发。
每个会员每个月租赁次数不得超过2次,每次获得3张DVD。
会员看完3张DVD之后,只需要将DVD放进网站提供的信封里寄回(邮费由网站承担),就可以继续下次租赁。
请考虑以下问题:
1)网站正准备购买一些新的DVD,通过问卷调查1000个会员,得到了愿意观看这些DVD的人数(表1给出了其中5种DVD的数据)。
此外,历史数据显示,60%的会员每月租赁DVD两次,而另外的40%只租一次。
假设网站现有10万个会员,对表1中的每种DVD来说,应该至少准备多少张,才能保证希望看到该DVD的会员中至少50%在一个月内能够看到该DVD?
如果要求保证在三个月内至少95%的会员能够看到该DVD呢?
2)表2中列出了网站手上100种DVD的现有张数和当前需要处理的1000位会员的在线订单,如何对这些DVD进行分配,才能使会员获得最大的满意度?
请具体列出前30位会员(即C0001~C0030)分别获得哪些DVD。
3)继续考虑表2,并假设表2中DVD的现有数量全部为0。
如果你是网站经营管理人员,你如何决定每种DVD的购买量,以及如何对这些DVD进行分配,才能使一个月内95%的会员得到他想看的DVD,并且满意度最大?
4)如果你是网站经营管理人员,你觉得在DVD的需求预测、购买和分配中还有哪些重要问题值得研究?
请明确提出你的问题,并尝试建立相应的数学模型。
问题背景
随着信息时代的到来,网络成为人们生活中越来越不可或缺的元素之一。
许多网站利用其强大的资源和知名度,面向其会员群提供日益专业化和便捷化的服务。
例如,音像制品的在线租赁就是一种可行的服务。
这项服务充分发挥了网络的诸多优势,包括传播范围广泛、直达核心消费群、强烈的互动性、感官性强、成本相对低廉等,为顾客提供更为周到的服务。
但是在线租赁必须要考虑到顾客的偏爱程度、现有数量、需求预测、购买和分配等因素,针对这些问题我们就有必要建立相应的数学模型进行研究,解决DVD的在线租赁问题!
模型假设:
1、网站允许用户租赁的最长时间为一个月;
2、同一张DVD租给同一个会员两次或两次以上,而且同种DVD不租给同一个会员2张或2张以上;
3、当天返回的DVD,只要订单需要,就在当天邮寄出;否则网站库储多余的DVD,从而造成网站的浪费;
4、网站订单一天处理一次,收到订单数目相对平稳,只有少许波动;
5、在每个会员密集区设立邮寄专区,负责DVD的回收,寄发,租赁中心和会员之间DVD传递的时间忽略不记;
6、设DVD在传送过程不会出现意外,即都能安全的到达会员和租赁中心。
符号说明:
(j=1,2,3,4,5);表示调查1000个会员愿意观看第j号DVD的人数;
(j=1,2,3,4,5);表示调查1000个会员愿意观看第j号DVD的比例;
(j=1,2,3,4,5);表示10万中至少50%在一个月内能够看到该DVD的人数;
(0——1)变量若第i个会员获得第j号DVD,则记
=1,否则
=0;
i(i=1、2……1000)表示会员数,
j(j=1、2……1000)表示DVD编号数
表示第i个会员对第j号DVD的偏爱程度;
表示第j号DVD的现有量;
模型的建立与求解
问题一:
根据表1给出的5种DVD数据:
C1=200,C2=100,C3=50,C4=25,C5=10;
根据
=
100%;
则
=20%,
=10%,
=5%,
=2.5%,
=1%;
=100000×
×50%;
则
=10000
=5000
设事件A=“为某人所获得到该DVD”
事件B1=“某人一个月内租第j号DVD(j=1,2,3,4,5)一次”,则P(B1)=0.4;
事件B2=“某人一个月内租第j号DVD(j=1,2,3,4,5)两次”,则P(B2)=0.6;
显然,B1,B2互不相容,且B1∪B2=
,可用全概率公式求解:
则
(A)=P(B1)P(A|B1)+P(B2)P(A|B2)=0.4×
+
(其中:
P(A|B1)=
P(A|B2)=
);
则P1=1.6x1/20000,P2=1.6x2/10000,P3=1.6x3/5000,P4=1.6x4/2500,P5=1.6x5/1000;
要保证至少50%的会员在一个月内能够看到该DVD
则愿意看DVDj的会员在一个月内能够看到DVDj的人数期望满足:
E(C)>=
*0.5;
E(C)=
代入数据得:
c1=6250;
c2=3125;
c3=1563;
c4=782;
c5=391,
(2)设需要DVD1为x1,DVD2为x2,DVD3为x3,DVD4为x4,DVD5为x5;
把借碟人分为A.B两类,
A类是60%的人,B类是40%的人
A类每月租碟两次,B类每月租碟一次,根据调查表的比例知:
10万人观看DVD1的人有10万*200/1000=20000;
观看DVD2的人有10万*100/1000=10000;
观看DVD3的人有10万*50/1000=5000;
观看DVD4的人有10万*25/1000=2500;
观看DVD5的人有10万*10/1000=1000,
要保证至少95%的会员三个月内看到该DVD,则需要该DVD碟数为:
3*(2x1*60%+x1*40%)>=20000*95%,
得出x1>=3960;
3*(2x1*60%+x1*40%)>=10000*95%,
得出x2>=1980;
3*(2x1*60%+x1*40%)>=5000*95%,
得出x3>=990;
3*(2x1*60%+x1*40%)>=2500*95%,
得出x4>=495;
3*(2x1*60%+x1*40%)>=1000*95%,
得出x5>=391;
因为保证至少,所以只能取等号:
x1=3960
x2=1980
x3=990
x4=495
x5=391
问题二:
从100种DVD的现有张数和当前需要处理的1000位会员的在线定单(表2),对这些DVD进行分配,根据每个会员每次租赁如果能够获得就只能获得3张DVD,如果不能获得就得到0张DVD,可以引进0——1变量表示会员是否获得DVD,从而建立这个问题的0——1线性规划模型,借助lingo8.0数学软件求解。
模型的建立与求解
引入0——1变量
,若第i个会员获得第j号DVD,记
=1;否则
=0。
根据DVD的服务管理要求,
应该满足两个约束条件:
第一、每次如果能够获得则只能获得3张DVD,即对于j=1,2……100,应有
;
第二、每种DVD都有一定的现有数量,对于1000会员要借该种DVD应该满足
;
用
表示第I个会员对第j号DVD的偏爱程度,当第I个会员获得第j号DVD时,
*
表示该会员的满意度,否则
*
=0。
从表2中会员的在线订单的情况来看,用数字1,2……10表示会员的偏爱程度,数字越小表示会员的偏爱程度越高,数字0表示对应的DVD当前不在会员的在线定单中,于是会员获得的满意度可表示为Z=
。
基于寻求使全体会员订单满意度最大,即是偏爱程度数最小的DVD分配方案,我们寻求的目标方案是使当前需要处理的所有订单的满意度数之和最小的处理方法。
综上,这个问题的0——1规划模型可写作
MinZ=
s.t
i=1,2……1000;
j=1,2……100;
xij={0,1};
但是在用lingo8.0对表2的处理中发现最小的数字“0”所示的意义和我们所寻求的目标方案有冲突,在综合数字0……10所表示的偏爱程度越来越大,而数字“0”所表示对应的DVD当前不在会员的在线订单中,我们按程度关系把数字“0”改换成数字“11”,即把表2处理成表2-1:
DVD编号
D001
D002
D003
D004
……
DVD现有量
10
40
15
20
……
C0001
6
11
11
11
……
C0002
11
11
11
11
……
C0003
11
11
11
3
……
C0004
11
11
11
11
……
……
……
……
……
……
……
运用lingo8.0进行最优化求解:
model:
sets:
w/1..1000/:
c;
v/1..100/:
d;
links(w,v):
a,x;
endsets
min=@sum(links:
a*x);
@for(v(j):
@sum(w(i):
x(i,j))<=d(j));
@for(w(i):
@sum(v(j):
x(i,j))=3);
@for(links(i,j):
@bin(x(i,j)));
data:
d=10,40,15,20,20,12,30,33,35,25,29,31,28,61,2,28,28,26,31,38,34,29,35,22,29,81,1,19,25,41,29,35,1,40,39,5,106,30,29,2,110,6,15,36,34,11,32,25,2,64,40,26,33,26,61,2,11,38,44,36,27,31,42,44,12,81,10,35,33,30,2,40,15,11,28,24,20,88,9,28,31,8,22,3,70,21,34,4,38,27,39,28,24,15,50,24,36,55,2,40;
a=6111111……
11111111……
1111113……
11111111……
…………;!
表2-1中的数据;
Enddata
End
将表(2-1)所转换来的数据代入这一模型,用lingo8.0软件求最优化的解,再对xij进行数据处理,得前30位会员所获得DVD如下表(2-2):
会员号
获得的DVD编号
C0001
D008
D041
D098
C0002
D006
D044
D062
C0003
D032
D050
D080
C0004
D007
D018
D041
C0005
D011
D066
D068
C0006
D019
D053
D066
C0007
D026
D066
D081
C0008
D031
D035
D071
C0009
D053
D078
D100
C0010
D041
D055
D085
C0011
D059
D063
D064
C0012
D002
D031
D041
C0013
D021
D078
D096
C0014
D023
D052
D089
C0015
D013
D052
D085
C0016
D010
D084
D097
C0017
D047
D051
D067
C0018
D041
D060
D078
C0019
D066
D084
D086
C0020
D045
D061
D089
C0021
D045
D050
D053
C0022
D038
D055
D057
C0023
D029
D081
D095
C0024
D037
D041
D076
C0025
D009
D069
D081
C0026
D022
D068
D095
C0027
D050
D058
D078
C0028
D008
D034
D082
C0029
D026
D030
D055
C0030
D037
D062
D098
问题三:
假设表2中DVD的现有数量全部为0,经营管理人员决定每种DVD的购买量为未知量
要使一个月内95%的会员得到他想看的DVD,对于j=1,2……100即要满足约束条件
;
在问题二的基础上,我们很快得到使一个月内95%的会员得到想看的DVD的最大满意度模型为:
MinZ=
;
s.t
;
i=1,2……1000;
xij={0,1};
运用lingo8.0建立模型求解:
model:
sets:
w/1..1000/:
c;
v/1..100/:
d;
links(w,v):
a,x;
endsets
min=@sum(links:
a*x);
@sum(links(i,j):
(x(i,j)))=0.95*3000;
@for(w(i):
@sum(v(j):
x(i,j))=3);
@for(links(i,j):
@bin(x(i,j)));
data:
a=a=6111111……
11111111……
1111113……
11111111……
…………;!
表2-1中的数据;
Enddata
End
利用lingo8.0求最优化解,再对所得到的
(j=1、2……100)进行统计得到DVD的分配情况如下:
现在将每种DVD的购买量绘制图表如下:
DVD编号
D001
D002
D003
D004
D005
D006
D007
D008
D009
D010
DVD购买量
21
36
27
38
21
28
30
33
35
25
DVD编号
D011
D012
D013
D014
D015
D016
D017
D018
D019
D020
DVD购买量
29
31
28
31
27
38
28
26
31
38
DVD编号
D021
D022
D023
D024
D025
D026
D027
D028
D029
D030
DVD购买量
34
29
35
22
29
31
26
19
25
41
DVD编号
D031
D032
D033
D034
D035
D036
D037
D038
D039
D040
DVD购买量
29
35
31
31
39
35
21
30
29
28
DVD编号
D041
D042
D043
D044
D045
D046
D047
D048
D049
D050
DVD购买量
53
35
26
35
34
25
32
25
32
34
DVD编号
D051
D052
D053
D054
D055
D056
D057
D058
D059
D060
DVD购买量
40
26
33
26
31
32
31
28
34
36
DVD编号
D061
D062
D063
D064
D065
D066
D067
D068
D069
D070
DVD购买量
27
31
32
34
32
31
30
35
33
30
DVD编号
D071
D072
D073
D074
D075
D076
D077
D078
D079
D080
DVD购买量
36
34
25
31
28
24
20
30
31
28
DVD编号
D081
D082
D083
D084
D085
D086
D087
D088
D089
D090
DVD购买量
29
18
22
19
33
21
34
24
24
27
DVD编号
D091
D092
D093
D094
D095
D096
D097
D098
D099
D100
DVD购买量
39
28
24
23
40
24
36
32
18
35
问题四:
1.一元线性回归预测:
(1)模型的建立:
作为网站的管理人员,需要对DVD的需求量和网站的购买量作出预测,这里我们提出了一个线性回归预测模型对未来的需求量和购买量作出预测。
根据第3问的表(3-1)所反映出的购买量与需求量的关系,可假设DVD的购买量为y与需求量为x之间有如下一元线性结构式:
y=a0+a1x+ε(4-1)
(其中a0,a1是两个未知参数,ε为其他随机因素对y的影响,x是非随机可精确观察的,ε是均值为零的随机变量是不可观察的,一般地称由(4-1)确定的模型为一元线性回归模型)
式
(1)两边同时取期望得:
Y=a0+a1X(4-2)
称式(4-2)为Y对X的回归直线方程,在该模型下,现假设该网站有某段时间内有n个历史观测值,该数据可以看作Yi=a0+a1xi+εi(这些样本相互独立但不同分布,i=1,2,……n)的实际抽样值,即样本值。
其中假定Y的方差为常数,a0和a1的回归系数,分别表示直线在Y轴的截距和直线的斜率。
这些可以根据历史观测值用最小二乘法求解,这使得实际数据与该直线的估计之间的误差最小。
给定s个样本或形如(x1,y1)(x2,y2)…....(xs,ys)的数据点,回归系数a0和a1可以用matlab6.5处理计算:
a1=
(4-3)
a0=
(4-4)
其中
是x1,x2,…xs的平均值,而
是y1,y2……ys的平均值。
系数a0和a1通常给出其他情况下复杂回归方程的很好近似。
(2)回归方程的显著性检验:
在实际工作中,事先我们并不能断定y与x之间的线性关系,式(4-1)只是一种假设。
当然,这个假设不是没有根据,我们可以通过专业知识和散点图作粗略判断。
但在求出回归方程之后,还须对这种线性回归方程同实际观测数据拟合的效果进行检验
由式(4-1)可知,当|a1|越大,y随x的变化的趋势就越明显;反之,当|a1|越小,y随x的变化是趋势就越不明显,特别当a1等于0时,则认为y与x之间不存在线性关系。
当a1不等于0时,则认为y与x之间有线性关系(4-1)。
因此,问题归结为对假设
H0:
a1=0;H1:
a1≠0
进行检验。
假设H0:
a1=0;被拒绝,则回归显著,认为y与x存在线性关系,所求的线性回归方程有意义;否则回归不显著,y与x的关系不能用一元线性回归模型来描述。
对问题3得出需求量和购买量,用matlab软件处理数据加以检验预测是模型的显著性。
在matlab输入:
x=[841841749978878710093909597851028494102911001169610110993891018783979710087918210997919487871191049390106949488919410791989297991087785103941031051089810590961051019510685829086889982989977728490789573949810794939010278951018086
]';
X=[ones(100,1)x];
Y=[213627382128303335252931283127382826313834293522293126192541293531313935213029285335
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高教 全国大学生 数学 建模 竞赛